Re: [R-es] R encountered a fatal error

2018-02-18 Thread Carlos Ortega 
Hola,

¿Cómo es de grande el conjunto que estás modelizando?.
Aunque liquidSVM es menos intenso en necesidades de máquina que SVM, fíjate
que estás explorando un número muy grande de combinaciones y es muy
probable que te estés quedando sin recursos. Abre en la ejecución

*gammas = c(0.0001,0.001,0.04,0.1,50, 1,10, 100), c_values = c(0.1,
0.0001, 0.001,1,10,25,50)*

Te sugeriría que probaras varias cosas:

   - Primero reduce el número de gammas y c_values a tres valores cada uno
   o un par. Puedes utilizar un valor muy pequeño, otro intermedio y otro
   grande y así comenzar a ver cómo se comporta el modelo.
   - Y si lo anterior tampoco funciona, reduce el número del conjunto
   "trainning". Haz un muestreo quedándote con el 80% de filas y baja si no
   funciona.


Saludos,
Carlos Ortega
www.qualityexcellence.es

2018-02-18 11:22 GMT+01:00 Alberto :

> Hola a todos,
>
> cada vez que intento ejecutar el código adjuntado cuando llego al comando
> de predict() salta un pop-up de R diciendo 'R encountered a fatal error'.
> No se donde puede estar el fallo y no entiendo tampoco por qué pasa esto
> puesto que a mi parecer el código está bien. Alguna idea de cual puede ser
> el problema?
>
> Muchas gracias!
>
> set.seed(6)
>
> model <- init.liquidSVM(Casualty_Severity~Vehicle_Type+Vehicle_
> Manoeuvre+Junction_Location+Skidding_and_Overturning+Hit_
> Object_off_Carriageway+First_Point_of_Impact+Journey_
> Purpose_of_Driver+Sex_of_Driver+Age_Band_of_Driver+Propulsion_Code+Age_of_
> Vehicle+Driver_Home_Area_Type+Sex_of_Casualty+Age_Band_of_
> Casualty+Car_Passenger+Casualty_Type+Number_of_Vehicles+Hour_of_Day+First_
> Road_Class+Road_Type+Speed_limit+Junction_Detail+Light_
> Conditions+Weather_Conditions+Road_Surface_Conditions+Urban_
> or_Rural_Area+month+other_vehic, trainning, threads = -1, gammas =
> c(0.0001,0.001,0.04,0.1,50, 1,10, 100), c_values = c(0.1, 0.0001,
> 0.001,1,10,25,50))
>
> trainSVMs(model, threads = -1, gammas = c(0.0001,0.001,0.04,0.1,50, 1,10,
> 100), c_values = c(0.1, 0.0001, 0.001,1,10,25,50), solver = 'ls',
> command.args=list(L=2, T=-1, d=1))
>
> selectSVMs(model)
>
> svm.probs <- predict(model, type = 'response', newdata = tst)
>
>
>
> [[alternative HTML version deleted]]
>
>
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>


-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

[R-es] Estudio posthoc de la interacción

2018-02-18 Thread Rodrigo López Correa 
Estimados,

Estoy analizando el siguiente modelo para producción de leche y he
utilizado el paquete car para plantear el siguiente modelo de efectos fijos:


modelo<-Anova(lm(leche~tambo_ypart+mpart+lact+rcs_categ3+lact*:*
rcs_categ3,data=a),type="III")

El resultado dio siginificativa la interacción entre lact y  rcs_categ3.

Para analizar la interacción, utilicé el paquete phia:

interacc_medias<-interactionMeans(modelo)

Pero me produce el siguiente error:

*Error in terms.default(model) : no terms component nor attribute*


Me podrían ayudar a estudiar la interacción posthoc con este paquete u otro?


Muchas gracias!

Rodrigo.



Virus-free.
www.avast.com

<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>

[[alternative HTML version deleted]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Estudio posthoc de la interacción

2018-02-18 Thread Javier Marcuzzi 
Estimado Rodrigo López Correa

Yo realizo ese tipo de análisis en mi vida profesional.

¿Qué título tiene o que estudia?

Lo que usted escribe tiene sentido pero biológicamente está con errores, en
su modelo, la producción de leche depende del tambo, mes de parto, año de
parto, lactancia, recuento de células somáticas y relación entre lactancia
y células somáticas. Luego dice que tiene relación significativa entre la
lactancia y células somáticas.

Biológicamente las células somáticas afectan la producción de leche, la
producción de leche afecta las lactancias, entrando muy fino puede haber
diferencias en la histología, pero un cambio por ejemplo ambiental afecta
el recuento de células somáticas. Si mira algún artículo del Journal Dairy
Science su modelo no tiene algo típico como leche = hys + ..., una simple
pregunta como si al envejecer la vaca aumentan las células somáticas, es lo
único a lo que apunta claramente su modelo a responder desde mi punto de
vista, el resto de su modelo tiene problemas biológicos, le recomiendo
definir un poco más ese aspecto de acuerdo a lo que buscas, R le puede
procesar información y dar resultados, pero las vacas se manejan con otras
reglas.

Usted tiene un error anterior, el análisis de las variancias puede ser
bueno, pero usted no está reflejando la producción de leche, no le falta
mucho para escribir el modelo en una forma adecuada, luego lo puede
complicar todo lo que desea, pero le recomiendo ver ese aspecto antes de
continuar.

Javier Rubén Marcuzzi



El 18 de febrero de 2018, 12:13, Rodrigo López Correa 
escribió:

> Estimados,
>
> Estoy analizando el siguiente modelo para producción de leche y he
> utilizado el paquete car para plantear el siguiente modelo de efectos
> fijos:
>
>
> modelo<-Anova(lm(leche~tambo_ypart+mpart+lact+rcs_categ3+lact*:*
> rcs_categ3,data=a),type="III")
>
> El resultado dio siginificativa la interacción entre lact y  rcs_categ3.
>
> Para analizar la interacción, utilicé el paquete phia:
>
> interacc_medias<-interactionMeans(modelo)
>
> Pero me produce el siguiente error:
>
> *Error in terms.default(model) : no terms component nor attribute*
>
>
> Me podrían ayudar a estudiar la interacción posthoc con este paquete u
> otro?
>
>
> Muchas gracias!
>
> Rodrigo.
>
>
>  source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=icon>
> Virus-free.
> www.avast.com
>  source=link&utm_campaign=sig-email&utm_content=webmail&utm_term=link>
> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
> [[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]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Microsoft R Open

2018-02-18 Thread Lucas Fernandez Seivane 
Hola Jesús
Es la antigua Revolution R. Está muy acelerada, usa MKL por defecto Yo
la uso desde hace mucho tiempo.
Un saludo.
ᐧ

2018-02-18 23:50 GMT+01:00 Jesús Para Fernández <
j.para.fernan...@hotmail.com>:

> Buenas,
>
> Alguno ha probado a usar Microsfot R Open? Segun leo, se han centrado en
> hacer que R por defecto sea multicore (cosa que podemos hacer usando
> parallel, foreach...)
>
> https://mran.revolutionanalytics.com/documents/rro/multithread
>
> Home [mran.revolutionanalytics.com] documents/rro/multithread>
> This site uses cookies for analytics, personalized content and ads. By
> continuing to browse this site, you agree to this use. Learn more
> mran.revolutionanalytics.com
>
>
> De ser así, me llama la atención lo inadvertida que pasa esta versión..
>
> Un saludo
> Jesús
>
> [[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]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Microsoft R Open

2018-02-18 Thread Jesús Para Fernández 
Me la recomiendas por  encima de la version de CRAN?

De: Lucas Fernandez Seivane 
Enviado: domingo, 18 de febrero de 2018 23:52
Para: Jesús Para Fernández
Cc: r-help-es@r-project.org
Asunto: Re: [R-es] Microsoft R Open

Hola Jesús
Es la antigua Revolution R. Está muy acelerada, usa MKL por defecto Yo la 
uso desde hace mucho tiempo.
Un saludo.
[https://mailfoogae.appspot.com/t?sender=acXVldmVkaW5AZ21haWwuY29t&type=zerocontent&guid=fedf1dc9-793b-40ba-88d9-2fea19ed8af4]ᐧ

2018-02-18 23:50 GMT+01:00 Jesús Para Fernández 
mailto:j.para.fernan...@hotmail.com>>:
Buenas,

Alguno ha probado a usar Microsfot R Open? Segun leo, se han centrado en hacer 
que R por defecto sea multicore (cosa que podemos hacer usando parallel, 
foreach...)

https://mran.revolutionanalytics.com/documents/rro/multithread

Home 
[mran.revolutionanalytics.com]>
This site uses cookies for analytics, personalized content and ads. By 
continuing to browse this site, you agree to this use. Learn more
mran.revolutionanalytics.com


De ser así, me llama la atención lo inadvertida que pasa esta versión..

Un saludo
Jesús

[[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]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Microsoft R Open

2018-02-18 Thread Lucas Fernandez Seivane 
En general sí (es sencilla de instalar, aunque si eres hábil compilando y
usas las MKL de intel podrías obtener las mismas o incluso superiores
ganancias de velocidad).
https://mran.microsoft.com/documents/rro/multithread#mt-bench

ᐧ

2018-02-19 0:02 GMT+01:00 Jesús Para Fernández :

> Me la recomiendas por  encima de la version de CRAN?
> --
> *De:* Lucas Fernandez Seivane 
> *Enviado:* domingo, 18 de febrero de 2018 23:52
> *Para:* Jesús Para Fernández
> *Cc:* r-help-es@r-project.org
> *Asunto:* Re: [R-es] Microsoft R Open
>
> Hola Jesús
> Es la antigua Revolution R. Está muy acelerada, usa MKL por defecto Yo
> la uso desde hace mucho tiempo.
> Un saludo.
> ᐧ
>
> 2018-02-18 23:50 GMT+01:00 Jesús Para Fernández <
> j.para.fernan...@hotmail.com>:
>
> Buenas,
>
> Alguno ha probado a usar Microsfot R Open? Segun leo, se han centrado en
> hacer que R por defecto sea multicore (cosa que podemos hacer usando
> parallel, foreach...)
>
> https://mran.revolutionanalytics.com/documents/rro/multithread
> 
>
> Home [mran.revolutionanalytics.com
> 
> ] 
> >
> This site uses cookies for analytics, personalized content and ads. By
> continuing to browse this site, you agree to this use. Learn more
> mran.revolutionanalytics.com
> 
>
>
> De ser así, me llama la atención lo inadvertida que pasa esta versión..
>
> Un saludo
> Jesús
>
> [[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]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Microsoft R Open

2018-02-18 Thread Javier Marcuzzi 
Estimado Jesús Para Fernández

Usted puede tener una confusión que yo tuve en su momento, ahora utilizo
Microsoft R Open.

Algunos puntos, el compilador, desconozco si alguna vez compilo fortran con
gfortran y el compilador de Intel, el segundo es mucho más grande, demora
más, pero realiza algunas optimizaciones, la diferencia, es comercial.

El segundo punto es el código, este puede ser escrito por el autor o desde
una librería que tenga optimizado respecto a versiones libres válidas, pero
que la ciencia de la computación pudo mejorar. En el caso de Microsoft
compró tecnología, por lo cuál tiene código más eficiente y compiladores
que tienen mejor performance.

Lo nuevo es R, luego Microsoft tiene su versión que va con tecnología
probada, en otras palabras, si requiere una función nueva puede ser que la
única alternativa es el R GNU, de lo contrario Microsoft ofrece la
característica de poder fijar una versión para siempre tener reproducciones
exactas.

Lo que se agrega, por ejemplo, la librería R español por decir algo, si
está en los repositorios de Microsoft R Open tengo la duda si fue compilada
con los compiladores comerciales o los GNU, si son los repositorios GNU
seguro fueron compilados con los compiladores GNU.

Javier Rubén Marcuzzi

El 18 de febrero de 2018, 19:50, Jesús Para Fernández <
j.para.fernan...@hotmail.com> escribió:

> Buenas,
>
> Alguno ha probado a usar Microsfot R Open? Segun leo, se han centrado en
> hacer que R por defecto sea multicore (cosa que podemos hacer usando
> parallel, foreach...)
>
> https://mran.revolutionanalytics.com/documents/rro/multithread
>
> Home [mran.revolutionanalytics.com] documents/rro/multithread>
> This site uses cookies for analytics, personalized content and ads. By
> continuing to browse this site, you agree to this use. Learn more
> mran.revolutionanalytics.com
>
>
> De ser así, me llama la atención lo inadvertida que pasa esta versión..
>
> Un saludo
> Jesús
>
> [[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]]

___
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es