Hola Isabel,
Mi primera idea ha sido pensar en un ajuste no-lineal (que creo que sigue
siendo válida la opción), pero en cualquier caso he mirado la lista del
R-Help (inglés) y hubo una entrada en el 2009 preguntando por lo mismo que
quieres.
http://tolstoy.newcastle.edu.au/R/e6/help/09/01/0626.html
En la respuesta verás que se referencia un artículo y aparece embebido el
código a utilizar en R.
Fuera de esto no sé si si hay alguna librería específica que lo contemple.
Voy a seguir mirando.
Y sobre lo que comentas que no sabes cómo obtener el valor de un parámetro.
Por favor, indícanos qué librería, función utilizas y vemos si es posible
obtenerlo.
Saludos,
Carlos Ortega
www.qualityexcellence.es
www.datanalytics.com/blog
2011/1/10 Isabel Martínez <ismartinez84@gmail.com>
> 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]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>
[[alternative HTML version deleted]]