Estimado Luciano

Se podría y no se podría, separo en los datos y en el desarrollo, sobre el desarrollo hay algo de código que me llevo meses, donde adapto una librería que ya no existe, esa parte del código se podría compartir. Sobre los datos, unos son o eran públicos (cabras) y los otros que tengo pero son propiedad de un estudio y es confidencial.

Cuándo procese datos públicos (cabras) todo funciona perfecto, cuándo ejecuto el mismo código con los datos privados, aparece el error, como parte de las fórmulas tienen logaritmo, elimine todos los NA y 0, todos “rabautizados a 1”. Y el error continuaba pero solo en dos tipos de curvas. 

Como estaba utilizando lm, que era escrito en origen por el autor de la librería que nunca llego a ser pública (servidores r-project), incluso se abandonó el desarrollo, decidí utilizar gls, nlme, lme4, y ando probando, al “abandonar lm” el error decidió “abandonarme”. 

Se podría investigar la causa, es una mezcla de datos y lm, estadísticamente no es mi área pero hablando con mis palabras, las singularidades desaparecen al cambiar el método estadístico (si bien es por algebra ...). Entonces ¿para que hacerles perder tiempo en algo que tiene solución? 

Javier

From: Luciano Selzer 
Sent: Monday, April 16, 2012 10:06 AM
To: Marcuzzi, Javier Rubén 
Cc: r-help-es@r-project.org 
Subject: Re: [R-es] lm.fit(x, y, offset = offset, singular.ok = singular.ok

Podrías postear un ejemplo reproducible?
Luciano



El 15 de abril de 2012 00:44, Marcuzzi, Javier Rubén <javier.ruben.marcuzzi@gmail.com> escribió:

  Hola a todos Tengo un problema con lm, al procesar los datos informa el siguiente error: lm.fit(x, y, offset = offset, singular.ok = singular.ok Los datos tienen algunos NA, por eso uso na.action=..., el modelo es:lm(logY ~ Curva(Dias), data=C_L_J_1, na.action="na.exclude") Desconozco que puedo hacer, hay información, incluso un libro lo menciona en errores comunes en R, pero no pude solucionar mi problema. Aclaración: Curva(Dias) es una función, como seno + seno*2 para todos los valores de Dias, esa parte funciona perfecto con otros datos que no son los que uso en el modelo data=C_L_J_1 Javier
         [[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]]

