Eu fiz desse jeito e deu essa mensagem:
clogit(as.numeric(con$status)~con$tnf+strata(con$cond))
Call:
clogit(as.numeric(con$status) ~ con$tnf + strata(con$cond))
coef exp(coef) se(coef) z p
con$tnfG/A 1.3084 3.700 67935274 1.93e-08 1
con$tnfG/G -0.0248 0.976 69421626 -3.57e-10 1
Likelihood ratio test=0 on 2 df, p=1 n= 85, number of events= 85
(170 observations deleted due to missingness)
Mensagens de aviso perdidas:
1: In Surv(rep(1, 255L), as.numeric(con$status)) :
Invalid status value, converted to NA
2: In fitter(X, Y, strats, offset, init, control, weights = weights, :
Ran out of iterations and did not converge
Em 5 de novembro de 2012 21:45, Victor Eduardo
<victorduca08@gmail.com> escreveu:
Oi professor. Se eu simplesmente transformar a variável para numérica com o comando as.numeric eu vou estar fazendo o procedimento certo? Eu fiz e rodou, mas da uma mensagem de aviso.
Abraços e boa semana.
Em 5 de novembro de 2012 21:07, Paulo Justiniano
<paulojus@leg.ufpr.br> escreveu:
entao. a msg de erro é clara:
Invalid status value, must be logical or numeric
o que nao é obedecido pois:
> class(con$status)
[1] "factor"
On Mon, 5 Nov 2012, Victor Eduardo wrote:
Pessoal, estou tentando rodar o script abaixo, mas está apresentando o seguinte erro:
clogit(con1$status~con1$tnf+strata(con1$cond),data=con1)
Erro em Surv(rep(1, 143L), con1$status) :
Invalid status value, must be logical or numeric
Tentei rodar sem o argumento data, mas apresentou o mesmo erro.
OBS: con1 é o meu banco de dados
A classe de cada uma das variáveis são:
> class(con$status)
[1] "factor"
> class(con$tnf)
[1] "genotype" "factor" > class(con$cond)
[1] "integer"
Abraços!
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.