Marcuzzi, Javier Rubén
2012-Apr-15 03:44 UTC
[R-es] lm.fit(x, y, offset = offset, singular.ok = singular.ok
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]]
Luciano Selzer
2012-Apr-16 13:06 UTC
[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]]
Marcuzzi, Javier Rubén
2012-Apr-16 13:45 UTC
[R-es] lm.fit(x, y, offset = offset, singular.ok = singular.ok
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]]