Buenas,

Estoy intentando validar un modelo de regresión logística y varias fuentes consultadas indican que debe validarse la tanto la calibración como la discriminación.

Estoy en la fase de discriminación y quisiera representar la curva ROC y calcular el área bajo la misma (AUC). El caso es que veo que hay muchas funciones para hacerlo y que se puede hacer tanto desde el paquete ROCBC como RATTLE. Con cada una de ellas logro una cosa distinta y no estoy seguro de estar haciéndolo bien.

Lo que yo tengo como entrada son las probabilidades de pertenencia de cada individuo a la clase de referencia predichas por el modelo logístico. Esto es, las p una vez despejadas del modelo logit(p)=SUMA(beta_{i}*X_{i}).

¿Qué función me recomendáis?

Gracias

Buenas:

Estoy trabajando con la interfaz de Rattle puesto que me parece muy interesante para mi propósito. Me gustaría copiar el código que genera los outputs para poder incluirlo en mi archivo .Rnw de Sweave y así producir mis informes de manera automática.

¿Alguien sabe si es posible hacer esto?

Muchas gracias de antemano

http://jjgibaja.net/curso-basico-y-gratuito-de-introduccion-a-r/comment-page-1/

Curso básico y gratuito de introducción a R

Para ayudar a quienes quieran dar este paso, Carlos Gil Bellosta y yo hemos diseñado un curso de R básico, con las siguientes características:

• Es gratuito.
• No da derecho a diplomas o certificados de ningún tipo.
• No es presencial.
• Dura 5 semanas.
• Las plazas son ilimitadas.
• Está basado en el autoestudio: cada participante tendrá que leer y trabajar por su cuenta.
• Es colaborativo: hemos desarrollado una plataforma web para que quienes sigan el curso puedan plantear preguntas y, como parte fundamental del programa, traten de responder las de sus compañeros.
• Está supervisado por Carlos y yo: nos encargaremos del programa del curso así como de tratar de dar apoyo de última instancia a las preguntas planteadas por los estudiantes.
• Es genérico. Cada cual quiere aprender R por un motivo distinto: unos, para analizar encuestas; otros, por su interés en la minería de datos; algunos, para analizar series temporales,… El curso está pensado para llevar a cada cual hasta el mismo umbral de su tema de interés, de forma que pueda, después de él, avanzar en el tema por su cuenta. Pero sin hacer especial hincapié en ningún asunto concreto.
• El contenido estadístico será mínimo (se limitará a algo de estadística descriptiva y poco más).

Los interesados encontrarán más información en este documento.

Para registrarte, basta con que te des de alta en cursoRbasico.usaR.org.es. Y si conoces a alguien potencialmente interesado en él, ¡no dudes en invitarlo!

Esperamos que esta iniciativa resulte de tu interés.

Diseñar un algoritmo de inversion para generar la variable aleatoria max(X1;X2; ::;Xn),
con X1; :::;Xn variables aleatorias i.i.d. con funcion de distribucion F (supongamos
F1 conocida).



agradesco su ayuda

Hola, necesito hacer un gráfico tipo histograma pero que las barras se solapen una con otras parcialmente.

Buenas:

Buenas!

Estoy utilizando Sweave para generar documentos en Latex, todo esto a través de R-Studio.

Me interesa sobre todo automatizar las órdenes a R de tal forma que cuando cambie la base de datos a analizar Sweave me lo compile sin problemas y sin apenas tener que modificar el archivo .Rnw.

He logrado programar un resumen de estadísticos principales así como tablas de recuentos de tal forma que, introduciendo al principio del documento el conjunto de variables a analizar, Sweave me genera todo a través de un ciclo for, de la primera a la última variable.

En estos momentos estoy intentando automatizar la generación de gráficos para que me produzca un gráfico para cada una de las variables introduciento un ciclo for en un solo chunk.

La verdad es que he encontrado poca documentación al respecto, ahí va un ejemplo:

http://www.statistik.lmu.de/~leisch/Sweave/FAQ.html#x1-11000A.9

He hecho pruebas con el chunk que proponen pero no logro nada:

<<results=tex,echo=FALSE>>=
for(i in 1:4){
file=paste("myfile", i, ".eps", sep="")
postscript(file=file, paper="special", width=6, height=6)
plot(rnorm(100)+i)
dev.off()
cat("\\includegraphics{", file, "}\n\n", sep="")
}
@

Si alguien pudiera ayudarme con este tema le estaría tremendamente agradecido.

Eskerrik asko/Muchas gracias

Hola estoy empezando en esto de R. Tengo una tabla con 20-columnas y me pregunto como puedo automatizar para que cada columna la utilice para correlacionarla con cor.test().


¿Pongo las columnas concatenadas?
c(tabla[,5:24])->correl

¿No hace falta?
tabla[,5:24]->correl

¿Como se añade eso para automatizarlo mediante for (i) {}?

Como el cor.test() necesita dos vectores debería dar valor X e Y en cada análisis a cada columna y así hasta enfrentarlas todas entre sí.

¿Alguna pincelada para avanzar un pasito? Sé que me puede faltar teoría sobre la diferencia de data.frame(), apply()... Pero sigo buscando.

Un saludo.

Buenas de nuevo,