[R-es] Gráfico tiempos de supervivencia

2019-07-18 Thread Griera-yandex 
Buenos días a todos:

Alguien me puede ayudar a hacer (si se puede) con unos datos similares a:

set.seed(20)
DATOS <- data.frame (
  ID  = c (1:10)
, TIEMPO  = sample(1:40, 10, replace=F)
, DEF = sample(0:1, 10, replace=T)
);DATOS

un gráfico que muestre los tiempos de supervivencia similar a:
https://miro.medium.com/max/700/1*yHhG4TVcAAi29Ln88t0_5Q.png

Lo he intentado con la función followup.plot del paquete
"epiDisplay" (https://cran.r-project.org/web/packages/epiDisplay/index.html),
pero no encuentro la forma.

Muchas gracias y saludos. 

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

Re: [R-es] Gráfico tiempos de supervivencia

2019-07-18 Thread PEDRO CONCEJERO CEREZO 
Hola, te vale esto? Es forma estandar de representar graficos supervivencia

Basado en esto:

https://rviews.rstudio.com/2017/09/25/survival-analysis-with-r/

set.seed(20)

DATOS <- data.frame (
  ID  = c (1:10)
  , TIEMPO  = sample(1:40, 10, replace=F)
  , DEF = sample(0:1, 10, replace=T)
)

DATOS

library(survival)

DATOS$DEF <- as.numeric(DATOS$DEF)
DATOS$TIEMPO <- as.numeric(DATOS$TIEMPO)

s <- Surv(DATOS$TIEMPO, DATOS$DEF)
head(s)

## Kaplan-Meier estimator.
km <- survfit(s ~ 1,
  data = DATOS,
  conf.type = "log-log")

## Show object
km

summary(km)

plot(km)

# Instala las librerías necesarias:
library(ranger)
library(ggplot2)
library(dplyr)
library(ggfortify)

autoplot(km)


El 18/07/2019 a las 12:00, 
r-help-es-requ...@r-project.org 
escribió:

Envíe los mensajes para la lista R-help-es a
r-help-es@r-project.org

Para subscribirse o anular su subscripción a través de la WEB
https://stat.ethz.ch/mailman/listinfo/r-help-es

O por correo electrónico, enviando un mensaje con el texto "help" en
el asunto (subject) o en el cuerpo a:
r-help-es-requ...@r-project.org

Puede contactar con el responsable de la lista escribiendo a:
r-help-es-ow...@r-project.org

Si responde a algún contenido de este mensaje, por favor, edite la
linea del asunto (subject) para que el texto sea mas especifico que:
"Re: Contents of R-help-es digest...". Además, por favor, incluya en
la respuesta sólo aquellas partes del mensaje a las que está
respondiendo.


Asuntos del día:

   1. Gráfico tiempos de supervivencia (Griera-yandex)

--

Message: 1
Date: Thu, 18 Jul 2019 10:18:42 +0200
From: Griera-yandex 
To: Lista R 
Subject: [R-es] Gráfico tiempos de supervivencia
Message-ID: <20190718101842.280414d5@debian-dde>
Content-Type: text/plain; charset="utf-8"

Buenos días a todos:

Alguien me puede ayudar a hacer (si se puede) con unos datos similares a:

set.seed(20)
DATOS <- data.frame (
  ID  = c (1:10)
, TIEMPO  = sample(1:40, 10, replace=F)
, DEF = sample(0:1, 10, replace=T)
);DATOS

un gráfico que muestre los tiempos de supervivencia similar a:
https://miro.medium.com/max/700/1*yHhG4TVcAAi29Ln88t0_5Q.png

Lo he intentado con la función followup.plot del paquete
"epiDisplay" (https://cran.r-project.org/web/packages/epiDisplay/index.html),
pero no encuentro la forma.

Muchas gracias y saludos.




--

Subject: Pié de página del digest

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


--

Fin de Resumen de R-help-es, Vol 125, Envío 7
*


--
Pedro Concejero
Telefónica CDO - 4th Platform - Internal Use Cases
E-mail: 
pedro.concejerocer...@telefonica.com
skype: pedro.concejero
twitter @ConcejeroPedro
linkedin pedroconcejero
eRReRo feliz, me puedes encontrar en gRupo R madRid 




Este mensaje y sus adjuntos se dirigen exclusivamente a su destinatario, puede 
contener información privilegiada o confidencial y es para uso exclusivo de la 
persona o entidad de destino. Si no es usted. el destinatario indicado, queda 
notificado de que la lectura, utilización, divulgación y/o copia sin 
autorización puede estar prohibida en virtud de la legislación vigente. Si ha 
recibido este mensaje por error, le rogamos que nos lo comunique inmediatamente 
por esta misma vía y proceda a su destrucción.

The information contained in this transmission is privileged and confidential 
information intended only for the use of the individual or entity named above. 
If the reader of this message is not the intended recipient, you are hereby 
notified that any dissemination, distribution or copying of this communication 
is strictly prohibited. If you have received this transmission in error, do not 
read it. Please immediately reply to the sender that you have received this 
communication in error and then delete it.

Esta mensagem e seus anexos se dirigem exclusivamente ao seu destinatário, pode 
conter informação privilegiada ou confidencial e é para uso exclusivo da pessoa 
ou entidade de destino. Se não é vossa senhoria o destinatário indicado, fica 
notificado de que a leitura, utilização, divulgação e/ou cópia sem autorização 
pode estar proibida em virtude da legislação vigente. Se recebeu esta mensagem 
por erro, rogamos-lhe q

Re: [R-es] Gráfico tiempos de supervivencia

2019-07-18 Thread Griera-yandex 
Hola Pedro:

Gracias por la ayuda. No conocía esta manera más elegante de mostrar las curvas 
de Kaplan-Meier. Te la compro.

En realidad quería mostrar un gráfico con la longitud de les tiempos de 
seguimiento y al final un símbolo para indicar el estado. Seria un gráfico 
similar a:

https://miro.medium.com/max/700/1*yHhG4TVcAAi29Ln88t0_5Q.png

Pero no encuentro la manera de hacerlo. Igual no se puede.

Muchas gracias y saludos.

On Thu, 18 Jul 2019 10:28:37 +
PEDRO CONCEJERO CEREZO  wrote:

> Hola, te vale esto? Es forma estandar de representar graficos supervivencia
> 
> Basado en esto:
> 
> https://rviews.rstudio.com/2017/09/25/survival-analysis-with-r/
> 
> set.seed(20)
> 
> DATOS <- data.frame (
>   ID  = c (1:10)
>   , TIEMPO  = sample(1:40, 10, replace=F)
>   , DEF = sample(0:1, 10, replace=T)
> )
> 
> DATOS
> 
> library(survival)
> 
> DATOS$DEF <- as.numeric(DATOS$DEF)
> DATOS$TIEMPO <- as.numeric(DATOS$TIEMPO)
> 
> s <- Surv(DATOS$TIEMPO, DATOS$DEF)
> head(s)
> 
> ## Kaplan-Meier estimator.
> km <- survfit(s ~ 1,
>   data = DATOS,
>   conf.type = "log-log")
> 
> ## Show object
> km
> 
> summary(km)
> 
> plot(km)
> 
> # Instala las librerías necesarias:
> library(ranger)
> library(ggplot2)
> library(dplyr)
> library(ggfortify)
> 
> autoplot(km)
> 
> 
> El 18/07/2019 a las 12:00, 
> r-help-es-requ...@r-project.org 
> escribió:
> 
> Envíe los mensajes para la lista R-help-es a
> r-help-es@r-project.org
> 
> Para subscribirse o anular su subscripción a través de la WEB
> https://stat.ethz.ch/mailman/listinfo/r-help-es
> 
> O por correo electrónico, enviando un mensaje con el texto "help" en
> el asunto (subject) o en el cuerpo a:
> 
> r-help-es-requ...@r-project.org
> 
> Puede contactar con el responsable de la lista escribiendo a:
> r-help-es-ow...@r-project.org
> 
> Si responde a algún contenido de este mensaje, por favor, edite la
> linea del asunto (subject) para que el texto sea mas especifico que:
> "Re: Contents of R-help-es digest...". Además, por favor, incluya en
> la respuesta sólo aquellas partes del mensaje a las que está
> respondiendo.
> 
> 
> Asuntos del día:
> 
>1. Gráfico tiempos de supervivencia (Griera-yandex)
> 
> --
> 
> Message: 1
> Date: Thu, 18 Jul 2019 10:18:42 +0200
> From: Griera-yandex 
> To: Lista R 
> Subject: [R-es] Gráfico tiempos de supervivencia
> Message-ID: <20190718101842.280414d5@debian-dde>
> Content-Type: text/plain; charset="utf-8"
> 
> Buenos días a todos:
> 
> Alguien me puede ayudar a hacer (si se puede) con unos datos similares a:
> 
> set.seed(20)
> DATOS <- data.frame (
>   ID  = c (1:10)
> , TIEMPO  = sample(1:40, 10, replace=F)
> , DEF = sample(0:1, 10, replace=T)
> );DATOS
> 
> un gráfico que muestre los tiempos de supervivencia similar a:
> https://miro.medium.com/max/700/1*yHhG4TVcAAi29Ln88t0_5Q.png
> 
> Lo he intentado con la función followup.plot del paquete
> "epiDisplay" (https://cran.r-project.org/web/packages/epiDisplay/index.html),
> pero no encuentro la forma.
> 
> Muchas gracias y saludos.
> 
> 
> 
> 
> --
> 
> Subject: Pié de página del digest
> 
> ___
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
> 
> 
> --
> 
> Fin de Resumen de R-help-es, Vol 125, Envío 7
> *
> 
> 
> --
> Pedro Concejero
> Telefónica CDO - 4th Platform - Internal Use Cases
> E-mail: 
> pedro.concejerocer...@telefonica.com
> skype: pedro.concejero
> twitter @ConcejeroPedro
> linkedin pedroconcejero
> eRReRo feliz, me puedes encontrar en gRupo R madRid 
> 
> 
> 
> 
> Este mensaje y sus adjuntos se dirigen exclusivamente a su destinatario, 
> puede contener información privilegiada o confidencial y es para uso 
> exclusivo de la persona o entidad de destino. Si no es usted. el destinatario 
> indicado, queda notificado de que la lectura, utilización, divulgación y/o 
> copia sin autorización puede estar prohibida en virtud de la legislación 
> vigente. Si ha recibido este mensaje por error, le rogamos que nos lo 
> comunique inmediatamente por esta misma vía y proceda a su destrucción.
> 
> The information contained in this transmission is privileged and confidential 
> information intended only for the use of the individual or entity n

[R-es] EpiLinux en DistroTest

2019-07-18 Thread miguel.angel.rodriguez.muinos 
Hola.

Ya se puede utilizar EpiLinux desde DistroTest.

https://distrotest.net/EpiLinux/5.0

Para los que no lo sepáis, EpiLinux es una distribución de linux especializada 
en software de bioestadística y epidemiología y viene con R + RCommander + 
RStudio instalados.
Es una oportunidad de poder "trastear" con R (y linux) sin instalar ni "romper" 
nada, directamente desde un navegador.

Más info: https://www.sergas.es/Saude-publica/Epilinux?idioma=es


Un Saludo,
--
Miguel Ángel Rodríguez Muíños
Coordinador del Proyecto EpiLinux
Dirección Xeral de Saúde Pública
Consellería de Sanidade
Xunta de Galicia
http://dxsp.sergas.es

















Nota: A información contida nesta mensaxe e os seus posibles documentos 
adxuntos é privada e confidencial e está dirixida únicamente ó seu 
destinatario/a. Se vostede non é o/a destinatario/a orixinal desta mensaxe, por 
favor elimínea. A distribución ou copia desta mensaxe non está autorizada.

Nota: La información contenida en este mensaje y sus posibles documentos 
adjuntos es privada y confidencial y está dirigida únicamente a su 
destinatario/a. Si usted no es el/la destinatario/a original de este mensaje, 
por favor elimínelo. La distribución o copia de este mensaje no está autorizada.

See more languages: http://www.sergas.es/aviso-confidencialidad

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

Re: [R-es] Gráfico tiempos de supervivencia

2019-07-18 Thread Carlos Ortega 
Hola,

Sí, lo puedes hacer de esta forma...

#-
set.seed(20)
DATOS <- data.frame (
ID  = c (1:10)
  , TIEMPO  = sample(1:40, 10, replace=F)
  , DEF = as.factor(sample(c(0,1), 10, replace=T))
)

library(ggplot2)

  ggplot( data = DATOS ) +
   geom_point( aes(x = TIEMPO, y = ID , shape = DEF, color = DEF), size = 5
) +
   geom_segment( aes( x = 0, y = ID,  xend = TIEMPO, yend = ID ) ) +
guides(colour = FALSE) +
labs(shape = 'LEGEND') +
scale_y_discrete() +
   theme_minimal()
#-

E incluso puedes reproducirlo usando fuentes parecidas a la de los comics
con el paquete "xkcd".

Saludos,
Carlos Ortega
www.qualityexcellence.es



El jue., 18 jul. 2019 a las 13:05, Griera-yandex ()
escribió:

> Hola Pedro:
>
> Gracias por la ayuda. No conocía esta manera más elegante de mostrar las
> curvas de Kaplan-Meier. Te la compro.
>
> En realidad quería mostrar un gráfico con la longitud de les tiempos de
> seguimiento y al final un símbolo para indicar el estado. Seria un gráfico
> similar a:
>
> https://miro.medium.com/max/700/1*yHhG4TVcAAi29Ln88t0_5Q.png
>
> Pero no encuentro la manera de hacerlo. Igual no se puede.
>
> Muchas gracias y saludos.
>
> On Thu, 18 Jul 2019 10:28:37 +
> PEDRO CONCEJERO CEREZO  wrote:
>
> > Hola, te vale esto? Es forma estandar de representar graficos
> supervivencia
> >
> > Basado en esto:
> >
> > https://rviews.rstudio.com/2017/09/25/survival-analysis-with-r/
> >
> > set.seed(20)
> >
> > DATOS <- data.frame (
> >   ID  = c (1:10)
> >   , TIEMPO  = sample(1:40, 10, replace=F)
> >   , DEF = sample(0:1, 10, replace=T)
> > )
> >
> > DATOS
> >
> > library(survival)
> >
> > DATOS$DEF <- as.numeric(DATOS$DEF)
> > DATOS$TIEMPO <- as.numeric(DATOS$TIEMPO)
> >
> > s <- Surv(DATOS$TIEMPO, DATOS$DEF)
> > head(s)
> >
> > ## Kaplan-Meier estimator.
> > km <- survfit(s ~ 1,
> >   data = DATOS,
> >   conf.type = "log-log")
> >
> > ## Show object
> > km
> >
> > summary(km)
> >
> > plot(km)
> >
> > # Instala las librerías necesarias:
> > library(ranger)
> > library(ggplot2)
> > library(dplyr)
> > library(ggfortify)
> >
> > autoplot(km)
> >
> >
> > El 18/07/2019 a las 12:00, r-help-es-requ...@r-project.org r-help-es-requ...@r-project.org> escribió:
> >
> > Envíe los mensajes para la lista R-help-es a
> > r-help-es@r-project.org
> >
> > Para subscribirse o anular su subscripción a través de la WEB
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
> >
> > O por correo electrónico, enviando un mensaje con el texto "help" en
> > el asunto (subject) o en el cuerpo a:
> > r-help-es-requ...@r-project.org r-help-es-requ...@r-project.org>
> >
> > Puede contactar con el responsable de la lista escribiendo a:
> > r-help-es-ow...@r-project.org r-help-es-ow...@r-project.org>
> >
> > Si responde a algún contenido de este mensaje, por favor, edite la
> > linea del asunto (subject) para que el texto sea mas especifico que:
> > "Re: Contents of R-help-es digest...". Además, por favor, incluya en
> > la respuesta sólo aquellas partes del mensaje a las que está
> > respondiendo.
> >
> >
> > Asuntos del día:
> >
> >1. Gráfico tiempos de supervivencia (Griera-yandex)
> >
> > --
> >
> > Message: 1
> > Date: Thu, 18 Jul 2019 10:18:42 +0200
> > From: Griera-yandex 
> > To: Lista R 
> > Subject: [R-es] Gráfico tiempos de supervivencia
> > Message-ID: <20190718101842.280414d5@debian-dde>
> > Content-Type: text/plain; charset="utf-8"
> >
> > Buenos días a todos:
> >
> > Alguien me puede ayudar a hacer (si se puede) con unos datos similares a:
> >
> > set.seed(20)
> > DATOS <- data.frame (
> >   ID  = c (1:10)
> > , TIEMPO  = sample(1:40, 10, replace=F)
> > , DEF = sample(0:1, 10, replace=T)
> > );DATOS
> >
> > un gráfico que muestre los tiempos de supervivencia similar a:
> > https://miro.medium.com/max/700/1*yHhG4TVcAAi29Ln88t0_5Q.png
> >
> > Lo he intentado con la función followup.plot del paquete
> > "epiDisplay" (
> https://cran.r-project.org/web/packages/epiDisplay/index.html),
> > pero no encuentro la forma.
> >
> > Muchas gracias y saludos.
> >
> >
> >
> >
> > --
> >
> > Subject: Pié de página del digest
> >
> > ___
> > R-help-es mailing list
> > R-help-es@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
> >
> >
> > --
> >
> > Fin de Resumen de R-help-es, Vol 125, Envío 7
> > *
> >
> >
> > --
> > Pedro Concejero
> > Telefónica CDO - 4th Platform - Internal Use Cases
> > E-mail: pedro.concejerocer...@telefonica.com pedro.concejerocer...@telefonica.com>
> > skype: pedro.concejero
> > twitter @ConcejeroPedro

[R-es] predict multinomial model con nnet

2019-07-18 Thread patricio fuenmayor 
Hola todos
Cuando realizo las predicciones del modelo multinomial con el paquete nnet,
estas cambian cada vez que lo ejecuto ... saben por qué pasa esto ??

Gracias por la ayuda.

[[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] predict multinomial model con nnet

2019-07-18 Thread Carlos Ortega 
Por la semilla.

Cada vez que inicias la red, los pesos comienzan con unos valores
aleatorios.
Si fijas la semilla, de ejecución en ejecución no debieras de ver variación.

Saludos,
Carlos Ortega
www.qualityexcellence.es

El jue., 18 jul. 2019 a las 14:58, patricio fuenmayor (<
patricio.fuenma...@gmail.com>) escribió:

> Hola todos
> Cuando realizo las predicciones del modelo multinomial con el paquete nnet,
> estas cambian cada vez que lo ejecuto ... saben por qué pasa esto ??
>
> Gracias por la ayuda.
>
> [[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

Re: [R-es] Gráfico tiempos de supervivencia

2019-07-18 Thread Griera-yandex 
Hola Carlos:

Como siempre, perfecta esta solución! Es lo que necesitaba.

Muchas gracias por esta ayuda y por el tiempo dedicado.

Saludos

On Thu, 18 Jul 2019 14:21:24 +0200
Carlos Ortega  wrote:

> Hola,
> 
> Sí, lo puedes hacer de esta forma...
> 
> #-
> set.seed(20)
> DATOS <- data.frame (
> ID  = c (1:10)
>   , TIEMPO  = sample(1:40, 10, replace=F)
>   , DEF = as.factor(sample(c(0,1), 10, replace=T))
> )
> 
> library(ggplot2)
> 
>   ggplot( data = DATOS ) +
>geom_point( aes(x = TIEMPO, y = ID , shape = DEF, color = DEF), size = 5
> ) +
>geom_segment( aes( x = 0, y = ID,  xend = TIEMPO, yend = ID ) ) +
> guides(colour = FALSE) +
> labs(shape = 'LEGEND') +
> scale_y_discrete() +
>theme_minimal()
> #-
> 
> E incluso puedes reproducirlo usando fuentes parecidas a la de los comics
> con el paquete "xkcd".
> 
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
> 
> 
> 
> El jue., 18 jul. 2019 a las 13:05, Griera-yandex ()
> escribió:
> 
> > Hola Pedro:
> >
> > Gracias por la ayuda. No conocía esta manera más elegante de mostrar las
> > curvas de Kaplan-Meier. Te la compro.
> >
> > En realidad quería mostrar un gráfico con la longitud de les tiempos de
> > seguimiento y al final un símbolo para indicar el estado. Seria un gráfico
> > similar a:
> >
> > https://miro.medium.com/max/700/1*yHhG4TVcAAi29Ln88t0_5Q.png
> >
> > Pero no encuentro la manera de hacerlo. Igual no se puede.
> >
> > Muchas gracias y saludos.
> >
> > On Thu, 18 Jul 2019 10:28:37 +
> > PEDRO CONCEJERO CEREZO  wrote:
> >
> > > Hola, te vale esto? Es forma estandar de representar graficos
> > supervivencia
> > >
> > > Basado en esto:
> > >
> > > https://rviews.rstudio.com/2017/09/25/survival-analysis-with-r/
> > >
> > > set.seed(20)
> > >
> > > DATOS <- data.frame (
> > >   ID  = c (1:10)
> > >   , TIEMPO  = sample(1:40, 10, replace=F)
> > >   , DEF = sample(0:1, 10, replace=T)
> > > )
> > >
> > > DATOS
> > >
> > > library(survival)
> > >
> > > DATOS$DEF <- as.numeric(DATOS$DEF)
> > > DATOS$TIEMPO <- as.numeric(DATOS$TIEMPO)
> > >
> > > s <- Surv(DATOS$TIEMPO, DATOS$DEF)
> > > head(s)
> > >
> > > ## Kaplan-Meier estimator.
> > > km <- survfit(s ~ 1,
> > >   data = DATOS,
> > >   conf.type = "log-log")
> > >
> > > ## Show object
> > > km
> > >
> > > summary(km)
> > >
> > > plot(km)
> > >
> > > # Instala las librerías necesarias:
> > > library(ranger)
> > > library(ggplot2)
> > > library(dplyr)
> > > library(ggfortify)
> > >
> > > autoplot(km)
> > >
> > >
> > > El 18/07/2019 a las 12:00, r-help-es-requ...@r-project.org > r-help-es-requ...@r-project.org> escribió:
> > >
> > > Envíe los mensajes para la lista R-help-es a
> > > r-help-es@r-project.org
> > >
> > > Para subscribirse o anular su subscripción a través de la WEB
> > > https://stat.ethz.ch/mailman/listinfo/r-help-es
> > >
> > > O por correo electrónico, enviando un mensaje con el texto "help" en
> > > el asunto (subject) o en el cuerpo a:
> > > r-help-es-requ...@r-project.org > r-help-es-requ...@r-project.org>
> > >
> > > Puede contactar con el responsable de la lista escribiendo a:
> > > r-help-es-ow...@r-project.org > r-help-es-ow...@r-project.org>
> > >
> > > Si responde a algún contenido de este mensaje, por favor, edite la
> > > linea del asunto (subject) para que el texto sea mas especifico que:
> > > "Re: Contents of R-help-es digest...". Además, por favor, incluya en
> > > la respuesta sólo aquellas partes del mensaje a las que está
> > > respondiendo.
> > >
> > >
> > > Asuntos del día:
> > >
> > >1. Gráfico tiempos de supervivencia (Griera-yandex)
> > >
> > > --
> > >
> > > Message: 1
> > > Date: Thu, 18 Jul 2019 10:18:42 +0200
> > > From: Griera-yandex 
> > > To: Lista R 
> > > Subject: [R-es] Gráfico tiempos de supervivencia
> > > Message-ID: <20190718101842.280414d5@debian-dde>
> > > Content-Type: text/plain; charset="utf-8"
> > >
> > > Buenos días a todos:
> > >
> > > Alguien me puede ayudar a hacer (si se puede) con unos datos similares a:
> > >
> > > set.seed(20)
> > > DATOS <- data.frame (
> > >   ID  = c (1:10)
> > > , TIEMPO  = sample(1:40, 10, replace=F)
> > > , DEF = sample(0:1, 10, replace=T)
> > > );DATOS
> > >
> > > un gráfico que muestre los tiempos de supervivencia similar a:
> > > https://miro.medium.com/max/700/1*yHhG4TVcAAi29Ln88t0_5Q.png
> > >
> > > Lo he intentado con la función followup.plot del paquete
> > > "epiDisplay" (
> > https://cran.r-project.org/web/packages/epiDisplay/index.html),
> > > pero no encuentro la forma.
> > >
> > > Muchas gracias y saludos.
> > >
> > >
> > >
> > >
> > > --
> > >
> > > Subject: Pié de página del digest
> > >
> > >

Re: [R-es] predict multinomial model con nnet

2019-07-18 Thread patricio fuenmayor 
Gracias Carlos por responder.
Te comento que si tengo fija la semilla.
Entendería yo que eso variaría en la estimación, pero en la predicción
debería tomar los resultados del modelo y aplicar los coeficientes. Lo
extraño es que si ejecuto varias veces solo las predicciones ... estas
tienen variaciones.

p1 <- data.table(predict(multinom.fit,newdata=train,"class"),1)
p2 <- data.table(predict(multinom.fit,newdata=train,"class"),2)
p3 <- data.table(predict(multinom.fit,newdata=train,"class"),3)
p4 <- data.table(predict(multinom.fit,newdata=train,"class"),4)
p5 <- data.table(predict(multinom.fit,newdata=train,"class"),5)
p6 <- data.table(predict(multinom.fit,newdata=train,"class"),6)
p7 <- data.table(predict(multinom.fit,newdata=train,"class"),7)

p <- rbindlist(list(p1,p2,p3,p4,p5,p6,p7))
dcast(p,V1~V2)
Using 'V2' as value column. Use 'value.var' to override
Aggregate function missing, defaulting to 'length'
  V1   1   2   3   4   5   6   7
1:  ALTO 773 785 772 785 786 783 772
2:  BAJO 551 533 537 539 532 555 551
3: MEDIO 512 518 527 512 518 498 513
>

Es extraño, voy a leer un poco mas como funciona la función predict, para
entender que pasa
Gracias por la ayuda.






El jue., 18 jul. 2019 a las 8:09, Carlos Ortega ()
escribió:

> Por la semilla.
>
> Cada vez que inicias la red, los pesos comienzan con unos valores
> aleatorios.
> Si fijas la semilla, de ejecución en ejecución no debieras de ver
> variación.
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
> El jue., 18 jul. 2019 a las 14:58, patricio fuenmayor (<
> patricio.fuenma...@gmail.com>) escribió:
>
>> Hola todos
>> Cuando realizo las predicciones del modelo multinomial con el paquete
>> nnet,
>> estas cambian cada vez que lo ejecuto ... saben por qué pasa esto ??
>>
>> Gracias por la ayuda.
>>
>> [[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