Otra alternativa de resolución solo con la paquetería estándar:
rm(list = ls() )
## Tenemos un marco de 100 posibles participante
## que numeramos con un ID del 1 al 100 (su DNI en el estudio)
## Se asigna sexo aleatoriamente con una probablidad de 0.5 para
cada uno
marco <- data
Hola, aunque con un poco de retraso, hago una aportación a este hilo:
Si tienes un problema con dos escalas de tiempo sobre las que
evolucionan las lineas de seguimiento, creo que es mejor que hacer un
solo gráfico en forma de diagrama de Lexis, que dos gráficos
independientes. Muestro ejemp
Esta es otra forma:
CADENA_1<-c("caliente","frío","gordo","flaco")
CADENA_2<-c("bueno","malo","bueno","malo")
CADENA<-rbind(CADENA_1,CADENA_2)
CADENA<-t(CADENA)
LISTA<-c("caliente","perro","flaco","gato","coche")
# ---
CADENA <- CADENA_2
names(CADENA) <- CADENA_1
CADENA[LISTA] -> RES
A ver si esto te vale:
-
### Construyo objeto "m.sim" para pruebas
c( 0, 16, 10, 20, 16, 9, 14,
16, 0, 8, 4, 2, 13, 14,
10, 8, 0, 2, 13, 4, 6,
20, 4, 2, 0, 19, 3, 25,
16, 2, 13, 19, 0, 17, 31,
9, 13, 4, 3, 17, 0, 7,
14, 14, 6, 25, 31,
Si lo que quieres obtener es una valoraci�n global del efecto de cada
tratamiento, qu�zas lo mas adecuado es que use regresi�n de Cox (si
puedes asumir proporcionalidad de los efectos. La resoluci�n de modelo
creo que seria algo similar a esto
> s = Surv(c$tiempo, c$estado)
> summary(coxph(s ~
El error cuadrático medio es por la asunción de Homocedasticidad el
mismo para todos los estratos de predicción, la varianza del termino de
error de la regresión: sigma²
Por ejemplo un simulación "naive" con tus ejemplo, se podría hacer así:
##
y=c(8.35,12.
¿Será que en la simulación de nuevos valores no hayas incluido el
termino de error del modelo?. La predicción lineal de modelo es sobre la
media condicionada. Si lo que quieres es simular una posible observación
tendrás que añadir un termino aleatorio a la predicción condicionada de
la media. L
Una alternativa que creo resuelve la cuestión es mediante un
"precompilador", como "gcc -E" o cualquier otro, incluyendo irectiva
#include y #define para la parametrización del código. Si embargo me
parece una solución un tanto barroca, yo intentaría parametrizarlo
usando código R puro.
Fr
Creo que es mucho mas sencillo que te descargues el fichero completo con
toda la serie y luego selecciones los datos que te interesa en R. Por
ejemplo puedes descargar toda la serie en formato pc-Axis e importarlo
con el paquete "pxR". Adjunto ejemplo:
download.file('http://www.ine.es/jaxiT3/fi
Tal y como lo entiendo, después de algunos mensajes, se podría
considerar tu problema como un análisis de Supervivencia, y por lo tanto
emplear las funciones "Surv" y "coxph", si consideramos que "bmi"
medido en cada entrevista es un buen "proxi" de la edad cronológica del
individuo, tal como
Dentro de los Segundos Cursos de Estad�stica y Demograf�a de la ADEH
2015
<[http://www.adeh.org/?q=es/contenido/segundos-cursos-de-estadistica-y-demografia-de-la-adeh-2015-segundo-curso-de-estatistica-e>
que tendr�n lugar del 13 al 17 de Julio en Guimar�es, Portugal, se
incluyen este a�o dos
Para los que estéis por Sevilla o cerca, el próximo martes os recuerdo
que a las 20:00 en el "Bulebar Café" de la Alameda, tendremos la 2ª
reunión del grupo de usuario de R de Sevilla, con un participante de
excepción "Romain François" coautor de "dplyr" y "Rccp". Mas
información en este enlace:
o interaccion
y estratificación en la misma variable !!
El 05/02/15 a las 10:46, Francisco Viciana escribió:
> Lo que me resulta estraño del modelo que planteas es la combinación de
> la "interacion" con la "estratificación":
>
> S1.nuevo ~ mobilityPDurG2 * strata(Count
Lo que me resulta estraño del modelo que planteas es la combinación de
la "interacion" con la "estratificación":
S1.nuevo ~ mobilityPDurG2 * strata(Countryb) + strata(Q6),
Esta claro que "coxph" se lo ha tragado, pero ni idea de los que puede
significar.
En principio, yo creo que la estratific
La de Jorge, es la respuesta correcta, aunque le sobraba una "n" al
parámetro
encoding = 'latin1'
Los dos fichero que ajuntados por "Eric" son detectados como 'latin-1'
por mi emacs, luego la manera correcta de leerlos independiente del
operativo, GUI y la configuración del lenguaje de nuest
15 matches
Mail list logo