Hola Rubén (¿Lobo?),

Efectivamente la función demod estaba presente en S-Plus [1], y formaba
parte del paquete VR (Venables and Ripley) que luego pasó a ser MASS
(Modern Applied Statistics with S), pero en R esa función ya no se incluyó
aunque se creó su equivalente "spec" (Spectral Density Estimation) aunque
ya comentas que no te vale. Entiendo que no te es útil porque tu serie es
univariada y la función por tanto no devuelve la fase. ¿Es así?.

En cuanto al test de Thom, he visto en el paper que referencias que se
trata de un test de rachas, que tampoco he encontrado en R. Sí que está el
test de rachas: función runs.test (paquete lawstat).

Y sí, el "Run sequence Test" es una forma un tanto peculiar de indicar el
tipo de gráfico que indicas. Que por cierto, he sido incapaz de encontrar
como tal en el libro de Chambers...

[1]: Todavía se puede encontrar los volúmenes que formaban los manuales de
usuario de S-Plus. En este segundo se trata el uso de esta función (página
281 del pdf):
http://www.uni-koeln.de/themen/statistik/software/s/v70/statman2.pdf

Saludos,
Carlos Ortega
www.qualityexcellence.es

El 3 de febrero de 2012 00:19, Rubén Gómez Antolí <lobo@mucharuina.com>escribió:

> Hola a todos:
>
> Estoy intentado realizar un modelo senoidal de unos datos de radiación
> solar con el fin de afrontar el relleno de la serie y aplicar la prueba de
> Thom para verificar su homogeneidad [0].
>
> De momento me encuentro con los siguientes problemas:
>
> 1- ¿Existe la prueba de Thom en R? ¿O debo crearme mi propia función?
>
> 2- Para la realización del modelo senoidal estoy siguiendo los pasos que
> marcan en este texto:
>
> http://www.itl.nist.gov/**div898/handbook/eda/section4/**eda4253.htm<http://www.itl.nist.gov/div898/handbook/eda/section4/eda4253.htm>
>
> pero me encuentro con los siguientes problemas:
>
> 1a - La serie que estoy utilizando tiene valores perdidos (ND) y hay
> varias funciones, como spectrum o lag.plot, que no se pueden utilizar. Yo
> los he reconvertido a 0, lo que no se si sera la correcto o me puede
> desvirtuar el modelo.
>
> 1b - No soy capaz de encontrar como obtener la gráfica de demodulación
> compleja, tanto de fase [1] como de amplitud [2], en R. He buscado y lo
> único que he encontrado es una referencia a la función «demod» en S dentro
> de un extracto del libro de Venables y Ripley.
>
> 1c - ¿Estoy yo muy despistado o el «run-sequence plot» es una gráfica de
> la función, en este caso una línea temporal con los años en el eje X y el
> valor de radiación en el Y?
>
> Por ahora creo que es todo.
>
> Gracias por adelantado.
>
> Salud y Revolución.
>
> Lobo.
>
> [0] http://www.ub.edu/gc/**Documentos/Sevilla_SunDu.pdf<http://www.ub.edu/gc/Documentos/Sevilla_SunDu.pdf>
>
> [1] http://www.itl.nist.gov/**div898/handbook/eda/section3/**compdeph.htm<http://www.itl.nist.gov/div898/handbook/eda/section3/compdeph.htm>
>
> [2] http://www.itl.nist.gov/**div898/handbook/eda/section3/**compdeam.htm<http://www.itl.nist.gov/div898/handbook/eda/section3/compdeam.htm>
>
> [3] http://www.itl.nist.gov/**div898/handbook/eda/section3/**runseqpl.htm<http://www.itl.nist.gov/div898/handbook/eda/section3/runseqpl.htm>
> --
> Libertad es poder elegir en cualquier momento. Ahora yo elijo GNU/Linux,
> para no atar mis manos con las cadenas del soft propietario.
> ---------
> Desde El Ejido, en Almería, usuario registrado Linux #294013
> http://www.counter.li.org
>
> ______________________________**_________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/**listinfo/r-help-es<https://stat.ethz.ch/mailman/listinfo/r-help-es>
>



-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]

