Versión de HTML Básico
Josué M. Polanco-Martínez
Analíti a
k
8
Revista de Análisis Estadístico
Journal of Statistical Analysis
primer valor del tiempo en el archivo de entrada mientras
que
t
N
es el último valor) [41, 42]. La tercera pauta es que
aunque REDFIT fue programado para utilizarse en Win-
dows, es posible compilar el código fuente para trabajar
en otros sistemas operativos
5
. La cuarta, es que los auto-
res de REDFIT recomiendan utilizar el software Gnuplot
(
http://www.manfredmudelsee.com/soft/redfit/
) para
graficar la salida de REDFIT. Sin embargo, en este artícu-
lo presentamos un sencillo programa en lenguaje R [34]
(véase anexo B) para visualizar la salida de REDFIT.
ANEXO B
B Programa para graficar la salida de REDFIT
######################################################################
# Este script grafica la salida de REDFIT [Schulz & Mudelsee 2002]
######################################################################
#::: Josu\’ M. Polanco-Mart\’nez (josue.m.polanco@gmail.com)
#::: C\’digo liberado bajo la Licencia P\’blica General de GNU (GPL)
#::: http://www.gnu.org/licenses/gpl-3.0.html
#::: Bilbao, Agosto 2014.
######################################################################
#::: Esta function transforma a decibeles
source("./decibel.R")
######################################################################
######################################################################
#::: Se lee la salida n1000.red de REDFIT y se definen variables.
######################################################################
redfitdat <- read.table("./n1000.red")
len = length(redfitdat[,1])
vars <- 1:dim(redfitdat)[2]
names(vars) <- c("Freq", "Gxx", "Gxx_cor", "Gred_th",
"<Gred>", "CorrFac", "80%-Chi2",
"90%-Chi2", "95%-Chi2", "99%-Chi2",
"80%-MC", "90%-MC", "95%-MC",
"99%-MC")
######################################################################
#::: Se grafica y guarda
######################################################################
pdf_file <- ("grafica_salida_redfit.pdf")
ranche <- range(db(redfitdat[,3]))
Y1 <- ranche[1] - 40
Y2 <- ranche[2] + 20
pdf(file=pdf_file, height=4, width=5)
plot(redfitdat[,1], db(redfitdat[,3]), t="l",
xpd=FALSE, xlab="", ylab="", ylim=c(Y1, Y2))
title(main="", xlab="Frecuencia [ciclos por A\’o]",
ylab="Amplitud espectral [dB]")
points(redfitdat[,1], db(redfitdat[,10]), t="l", lty=2, col="red")
points(redfitdat[,1], db(redfitdat[,14]), t="l", lty=3, col="blue")
5
El autor de este artículo compiló el código fuente de REDFIT para trabajar en Linux; si tiene interés en cómo llevar a cabo esta tarea, póngase en
contacto vía correo-e.
22
Analítika,
Revista de análisis estadístico
, 4 (2014), Vol. 8(2): 7-23