[R-br] Extração de valor-p na função hlme do pacote lcmm

[Éder Comunello]

Felipe, boa tarde!

Pelo que pude avaliar, o problema é que o código da função não prevê o
armazenamento da saída em objetos. Acredito que a solução mais correta
seria editar o código e definir uma nova função. Você pode sugerir essas
alterações ao autor do pacote.

Caso seu uso seja eventual, dá pra se virar com uma solução paliativa,
apresentada abaixo. Pode ser que os colegas da lista tenham uma ideia
melhor...

### <code r>
require(lcmm)
time <- rep(seq(1,10),50)
A <- as.data.frame(cbind(ID=rep(1:50,
each=10),time1=time,time2=time^2,Y=rnorm(500)))
gmm <- hlme(Y~time+time2, ng=3, mixture=Y~time+time2, subject='ID', data=A)
summary(gmm)

tmp <- summary(gmm) ### imprime na tela, mas gera objeto vazio
str(tmp)            ### retorna NULL

dput(summary.hlme) ### lista código da função
### observe ausência de instruções de retorno (return)
### uma saída é editar o código da função e inserir as instruções...

### uma solução paliativa! ;D
res  <- capture.output(summary(gmm)); res ### identifique o número das
linhas que interessam
sel  <- res[39:47]; sel ### selecione as linhas que interessam
pval <- read.table(textConnection(sel), head=F)[6]; pval ### a sexta coluna
refere-se ao p-valor
### <code r>

Éder Comunello <c <comunello.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140311/3e98f110/attachment.html>