
Pessoal,este é um exemplo fictício apenas para exemplificação! require(survival)tempo <- c(1, 18, 21, 21, 22, 29, 29, 35, 37, 39, 40, 50, 52, 54, 60, 80, 80, 81, 83, 84, 8)status <- c(1,0,0,0,0,0,0,1,1,0,1,0,0,0,0,0,0,0,0,0,0)grupos=c("A","F","F","F","M","A","M","A","M","F","M","A","M","F","F","F","M","F","M","F","A")km <- survfit(Surv(tempo,status) ~ grupos) # Estratificação summary(km)plot(km) survdiff(Surv(tempo, status)~grupos, rho = 1) # Teste Peto A 10% as curvas diferem, mas como comparar as curvas 2 a 2 para ver tais diferenças? A survival faz isto? André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

Bom dia André, Tem como fazer sim, primeiro ajuste um modelo nulo km0 <- survfit(Surv(tempo,status) ~ 1) e depois compare anova(km, km0, test="Chi") Se for significativo, os dois grupos são diferentes, Espero ter ajudado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== On 05/06/2015 07:15, Andre Oliveira wrote:
Pessoal, este é um exemplo fictício apenas para exemplificação!
require(survival) tempo <- c(1, 18, 21, 21, 22, 29, 29, 35, 37, 39, 40, 50, 52, 54, 60, 80, 80, 81, 83, 84, 8) status <- c(1,0,0,0,0,0,0,1,1,0,1,0,0,0,0,0,0,0,0,0,0) grupos=c("A","F","F","F","M","A","M","A","M","F","M","A","M","F","F","F","M","F","M","F","A") km <- survfit(Surv(tempo,status) ~ grupos) # Estratificação summary(km) plot(km) survdiff(Surv(tempo, status)~grupos, rho = 1) # Teste Peto
A 10% as curvas diferem, mas como comparar as curvas 2 a 2 para ver tais diferenças? A survival faz isto?
André Oliveira Souza.
Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES
_______________________________________________ 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.

Comparar curvas não é exatamente o que se faz. Tem como testar hipóteses, de igualdade por exemplo, sobre os parâmetros desses curvas. Isso requer o uso de um modelo paramétrico, como a weibull, log normal, gamma. Havendo ajustado um desses, hipóteses podem ser facilmente testadas com a glht, embora, pelo fato da log-verossimilhança ser assimétrica nesses modelos, é preferível testar hipóteses por razão de verossmilhanças do que por procedimentos de Wald (baseados na aproximação quadrática da log-ver) À disposição. Walmes.

Walmes, Então quer dizer que ajustar as curvas de sobrevivência com algum modelo paramétrico e ir agregando os níveis e compara los por qui quadrado esta incorreto? Obrigado Alexandre Enviado do Yahoo Mail no Android De:"walmes ." <walmeszeviani@gmail.com> Data:6:42 Sáb, 6 de Jun de AM Assunto:Re: [R-br] Curvas sobrevivência - Comparação Comparar curvas não é exatamente o que se faz. Tem como testar hipóteses, de igualdade por exemplo, sobre os parâmetros desses curvas. Isso requer o uso de um modelo paramétrico, como a weibull, log normal, gamma. Havendo ajustado um desses, hipóteses podem ser facilmente testadas com a glht, embora, pelo fato da log-verossimilhança ser assimétrica nesses modelos, é preferível testar hipóteses por razão de verossmilhanças do que por procedimentos de Wald (baseados na aproximação quadrática da log-ver) À disposição. Walmes.

Certas as hipóteses são sobre parâmetros e não sobre curvas, embora em certas situações, as hipóteses sobre os parâmetros podem ser vistas por meio das curvas. Por exemplo, você ajusta equação para relação altura~dap para dois genótipos florestais. Em fórmula R seria alt~gen*dap (relação linear para simplificar). Se você ajustar o modelo anterior e compara lo ao alt~dap, a hipótese que tá sendo testada é H0: \beta_{0,a} = \beta_{0,b} interseção \beta_{1,a} = \beta_{1,b}, que é o mesmo que H0: as chuvas são iguais ou não existe efeito algum de genótipo. Duas curvas serão iguais se é somente se todos os parâmetros forem iguais. Combinar/aglutinar/fundir níveis e criar submodelos é uma forma viável de testar hipóteses sim. Na realidade, essa é a ideia fundamental por trás do procedimento Scott-Knott. À disposição. -- ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (4)
-
Alexandre Santos
-
Andre Oliveira
-
ASANTOS
-
walmes .