Versión de HTML Básico
Estimación espectral de datos ambientales no equiespaciados vía el periodograma...
Analíti a
k
8
Revista de Análisis Estadístico
Journal of Statistical Analysis
######################################################################
#::: Drawing the significant spectral peaks (above 99% CL)
the.points <- which(redfitdat[,3] >= redfitdat[,14])
cat("the points", the.points)
len.points <- length(the.points)
freq.points <- array(0, len.points)
psdDC.points <- array(0, len.points)
freqs
<- redfitdat[,1]
psdDC
<- redfitdat[,3]
for (i in 1:len.points) {
freq.points[i] <- freqs[the.points[i]]
psdDC.points[i] <- psdDC[the.points[i]]
}
points(freq.points, db(psdDC.points), pch=19, cex=0.5)
text(freq.points, db(psdDC.points)+.5, labels=c(1:len.points),
pos=4, cex=.5)
legend("bottomleft", paste(1:len.points, "->",
round(1/freq.points[1:length(freq.points)], digits=1),
"A\’os"), bty="n", cex=.75)
legend("topright", lty=c(2,3), c("99% NF v\’a AR1", "99% NF v\’a MC"),
col=c("red", "blue"), bty="n", cex=0.7)
dev.off()
Analítika,
Revista de análisis estadístico
, 4 (2014), Vol. 8(2): 7-23
23