####Rotina de preparação do conjunto de dados######Ide2009 dados=read.table('c:\\dadosepv1.txt',header=T) dados=na.omit(dados) attach(dados) data2=rbind(dados,dados,dados)[,1:5] sujeito=rep(seq(1:nrow(dados)),3) tempo=c(rep('t30',nrow(dados)),rep('t60',nrow(dados)),rep('t90',nrow(dados))) resposta=c(t30,t60,t90) novosdados=data.frame(sujeito,tempo,data2,resposta) # dicotomizando detach(dados) attach(novosdados) novoresp=resposta novoresp[which(resposta%in%c(1,2))]=1 novoresp[which(resposta%in%c(3,4))]=0 novosdados=data.frame(novosdados,novoresp) novosdados<-novosdados[order(sujeito),] # criando três conjuntos separados por diagnóstico datA<-split(novosdados,novosdados[,"diag"])$A datB<-split(novosdados,novosdados[,"diag"])$B datC<-split(novosdados,novosdados[,"diag"])$C ####Obs: ver função reshape (PJ)