Hola a todos,

es la primera vez que escribo a la lista, espero explicarme correctamente.

Trabajo en cronobiología de mamíferos y quiero comenzar a analizar mis datos
con R (haciendo análisis de varianza ya llevo unos años en el entorno). Me
interesa saber si cada una de mis variables podría ajustarse a una función
sinusoidal; tengo medidas cada dos horas para cada una. Querría extraer el
mesor y la acrofase, y en definitiva obtener una p de bondad de ajuste.

He encontrado que las funciones 'cosinor()' de los paquetes 'season' (ritmos
de periocidad mayor a 24h) y 'psych' (de 24h) podrían valer, pero no sé cómo
extraer tal p. Con la receta que aparece en 'psych.cosinor()' puedo obtener
las correlaciones entre mis datos y unos en los que artificialmente he
introducido ruido sumándoles un 'rnorm()' de las dimensiones adecuadas, pero
¿de ahí cómo extraigo la significación y los parámetros restantes?

Quizá mi problema es que hasta ahora he empleado software cerrado que no sé
qué hace exactamente, así que no sé a dónde acudir: ¿es aplicable un glm? ¿o
quizá se trataría de regresiones circulares como circ.reg() del paquete
'CircStats', tomando las horas como ángulos?

En Matlab (no lo he probado) creo que existe algo que hace lo que busco: ?1.

Me he extendido demasiado así que resumo: ¿hay algún paquete para
cronobiología que haga un test de ritmicidad circadiana? Si no es así, ¿qué
debería hacer, glm, regresión circular...?

Cordialmente,

Isabel Martínez

?1
http://www.mathworks.com/matlabcentral/fileexchange/20329-cosinor-analysis

	[[alternative HTML version deleted]]

