similar to: problem with survreg and anova function

Hi, I make a weibull survival regression using suvreg function. Bu when I try to get the P values from anova, it give me NAs: I'm using R 2.4.0 and survival 2.29 Look: m <- survreg(Surv(tempo,censor)~grupo*peso) anova(m) Df Deviance Resid. Df -2*LL P(>|Chi|) NULL NA NA 148 966.6416 NA grupo -2 25.6334407 146 941.0081 NA

Hello fellows: I've een trying to set up a function that performs 100 loops producing the coresponding 100 series. I want to save all those datasets in a dataframe, so I wrote this... prep <- function() # Clase[1]/Categoria[2]/Phi[3]/Rf[4] peso <- c(.0,.03,.3,.6) # Extension del calculo for (i in 1:100) { # Calculos de todas las curvas # Variables (Valor Base) abase

Hola Fernando, Que pasa si usas subset(d, PESO == 76) ? Saludos, Jorge.- 2013/8/2 Fernando Macedo <fermace@gmail.com> > Estimados compañeros de lista, estábamos trabajando con unos datos que > adjunto y queríamos hacer una selección mediante índices pero nos deparamos > con algún problema. > > Por ejemplo si hacemos solamente: > > data[data$PESO==76] me da algo

Hi, I make an analysis and depending of the order of the variables, the significance change, look. m1 <- glm((infec/ntot)~idade+sexo+peso,family=binomial,weights=ntot) > anova(m1,test="F") Analysis of Deviance Table Model: binomial, link: logit Response: (infec/ntot) Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev F Pr(>F)

Hi, I have make an analise with presence and absence, y=(1 e 0). I have a area continuous data and a sp data with 25 levels. I have 300 points. When I make glm((presenca/peso)~area,weights=peso,family=binomial,maxit=1000) where presenca is 0 or 1. peso is the unit = 1. area is the continuous data. The analysis is OK. When I put the sp and interactions in analysis this warning appear.

Cuando haces... colIniAnalisis <- ncol(p.grupos)+ncol(p.pesos.ord)+2 ...sumas 2 al número de columnas de los dataframes pierdes el tipo "integer" y conviertes "colIniAnalisis" a coma flotante. Cuando construyes el rango ya no tienes enteros y no funciona. En teoría debería, pero no lo hace. Quizá alguien te pueda explicar por qué no funciona... Prueba con

Hi all, I try to understand the output of the summary(model). I make a simple model. > summary(m1anova) Call: glm(formula = peso ~ gen) Deviance Residuals: Min 1Q Median 3Q Max -4.3114 -2.3788 -0.9167 2.1581 5.7856 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 4.8660 0.9808 4.961 0.000101 *** gen2g 3.8504

XD... gajes del oficio. Te lo había dicho Carlos en el mensaje anterior... Un saludo Isidro Hidalgo Arellano Observatorio Regional de Empleo D.G. de Desarrollo de Estrategia Económica y Asuntos Europeos Avenida de Irlanda, 14 Tlf.: 925 28 80 98 ihidalgo en jccm.es Consejería de Empleo y Economía http://www.jccm.es > -----Mensaje original----- > De: Jorge Tornero - Listas

Pues este es el culpable... como soy nuevo en esto, disculpad si las estrategias de creación de datframes, etc. son poco ortodoxas. Y por cierto... esto del slicing con R es un poco... duro library(RPostgreSQL) library(reshape) #CARGA DE DATOS conn<-dbConnect("PostgreSQL",dbname="OFIDAT",user="antares") #consulta<-dbSendQuery(conn,"select

Hi! I`m trying to install xen 3.0.2 on my Debian system! But, the "make ARCH=xen menuconfig" don`t work! debian:/home/peso/xen-3.0.2/xen-3.0.2/linux-2.6.16-xen0# make ARCH=xen menuconfig Makefile:439: /home/peso/xen-3.0.2/xen-3.0.2/linux-2.6.16-xen0/arch/xen/Makefile: No such file or directory make: *** No rule to make target `/home/peso/xen-3.0.2/xen-3.0.2/linux-

Hi, My name is Jos?? Mar??a G??mez, and I am pretty new in R. Thus, I apologize deeply if my questions are extremmely na??ve.I have checked several available books and URL's, without finding any answer. I'm trying to fit Generalized Linear Mixed Models via PQL. Below I provide the structure of my data set. Year and Plot are random variables. Fate is the binomial dependent. I have severe

I've run an ancova, edadysexo is a factor with 3 levels,and log(lcc) is the covariate (continous variable) I get this results > ancova<-aov(log(peso)~edadysexo*log(lcc)) > summary(ancova) Df Sum Sq Mean Sq F value Pr(>F) edadysexo 2 31.859 15.9294 803.9843 <2e-16 *** log(lcc) 1 11.389 11.3887 574.8081 <2e-16 ***

Hola a todos, Me gustaria saber si en un modelo mixto se puede usar el tamaño de muestra como peso o es incorrecto hacerlo. Por ejemplo en el comando 'lmer' del paquete lme4 hay una opcion 'weights' (igual que en 'lm' de stats). Si tengo datos de una medida (por ejemplo el peso) de 10 especies de aves (rango, de 50 a 3000 datos segun la especie) ¿puedo usar esta n en

Hola a todos, Tengo que construir un gráfico para una base de datos que tengo en el que representar para los individuos de esta base su altura vs peso en función de su edad. Para ello, hay que utilizar los gráficos condicionados. En primer lugar usé la función coplot pero cargando el paquete lattice se pueden utilizar funciones que generan gráficos más completos, como es el caso de la función

Hola Buen día, estoy ajustando una regresión multiple, pero sucede que tengo un factor de expanción, digamos que deseo ajustar el modelo ingreso ~ edad + ocupacion mi problema es que cada observación tienen un factor de expanción, es decir cada variable no representa una observación sino más, que no es un número fijo para cada observación. por lo que a mí se me ocurrio usarlas como peso

Les consulto por lo siguiente: Se me ocurrió usar datos de MCMCglmm en otros procedimientos de R, comparar los resultados y el tiempo como facilidad de uso, para lo cual use LME4, los resultados fueron parecidos, entendiendo que internamente los algoritmos tienen diferencias, es lógico que los resultados no coincidan exactamente, pero aproximadamente tienen que ser compatibles, en este caso

Hola, paso el mini programita q estoy viendo, lo q me llama la atencion es una parte donde se definen las funciones. Probe primero meter adentro del boots la estadistica a estimar usando directamente gini(varible, pesos) pero no me dejo. Vi q en el ej del manual de boots, siempre define antes la funcion, entonces probe definir antes una funcion haciendo grini<-function(x) {gini(variable,

Hi, I know the function LETTERS, but, now I have some letters to convert it in numbers, like A=1,B=2, etc... Is any function to do that? Atenciosamente, Leandro Lins Marino Centro de Avalia??o Funda??o CESGRANRIO Rua Santa Alexandrina, 1011 - 2? andar Rio de Janeiro, RJ - CEP: 20261-903 R (21) 2103-9600 R.:236 ( (21) 8777-7907 ( leandro at cesgranrio.org.br "Aquele que suporta o peso da

Buenas tardes, Estoy intentando conseguir la media de la variable "numerocaracteres" por "producto" de mi base de datos, datos. Lo estoy haciendo con la función aggregate de este modo: AggregatedData<- aggregate(numerocaracteres ~ producto, data=datos, FUN=mean) El problema me viene porque quiero utilizar una variable de ponderación de modo que para construir la media

Hola: Si aún hay alguien que no está de vacaciones, igual me pueden ayudar. Quiero ajustar unos modelos: REG_LOG <- glm (low ~ X, family = "binomial", data = DATOS) Ejemplo: library(MASS) data(birthwt, package="MASS") birthwt$low <- factor(birthwt$low) birthwt$race <- factor(birthwt$smoke) REG_LOG <- glm (low ~ smoke, family =