- R-br - DINF UFPR

Identidade de modelos de regresão via lm
by Marcelo Cardoso mello 18 Out '11

18 Out '11

Como posso comparar os dois modelo (ma e mb) na CMR abaixo? Tipo um teste de identidade de modelos em lm? Como testar se os tratamentos 'a' e 'b' tem (estatisticamente) a mesma curva de regressão? Obrigado por qualquer dica. a<- 'a'; b<- 'b' t<-c(a,a,a,a,a,a,a,a,b,b,b,b,b,b,b,b) x<-c(0,0,10,10,30,30,45,45,0,0,10,10,30,30,45,45) r<-c(2,3,4,5,7,9,11,13,5,4,7,8,8,9,10,11) d<-data.frame(t,x,r);d da<-subset(d, t=='a'); da db<-subset(d, t=='b'); db m<-lm(r~x, data=d) ma<-lm(r~x, data=da) mb<-lm(r~x, data=db) m;ma;mb deviance(m)/df.residual(m) deviance(ma)/df.residual(ma) deviance(mb)/df.residual(mb)

2 1

Erro em Script SEM
by Joao Ricardo Nickenig Vissoci 18 Out '11

18 Out '11

Pessoal, boa noite. Sou recente no fórum e também no estudo com o R, mas estou aprendendo. Bom, queria tirar uma dúvida, estou trabalhando em um script para analise de Modelos de Equações Estruturais, e sempre que tento rodar o script que fiz encontro o seguinte erro: Iteration limit exceeded. O Script é esse, desculpem-me pelo Inglês, mas a discussão do grupo é toda em inglês. # MODEL II # Install the package 'sem' # carry data library library (sem) # Change R's directory to that in which the data base .CSV file is # Load database teste<- read.csv("BaseSEM.csv",header=T) teste # sem package model specification mod.wh.1 <- specify.model()# Type these values that specify the model's relations (just use de Ctrl+R over each relation). Support->SP,NA,1 Support->SM,sm1,NA Rejection->RP,rp1,NA Rejection->RM,NA,1 Overprot->OP,op1,NA Overprot->OM,NA,1 MotI->MIAO,NA,1 MotI->MIEE,moti1,NA MotI->MIC,moti2,NA MotE->MERE,NA,1 MotE->MEIn,mote1,NA MotE->MEId,mote2,NA Support->MotI,NA,1 Rejection->MotI,rej1,NA Overprot->MotI,over1,NA Support->MotE,sup,NA Rejection->MotE,rej2,NA Overprot->MotE,NA,1 MotI->SCI,misci,NA MotE->SCI,mesci,NA Support<->Support,ss1,NA Rejection<->Rejection,rr1,NA Overprot<->Overprot,oo1,NA OP<->OP,opcov,NA OM<->OM,omcov,NA SP<->SP,spcov,NA SM<->SM,smcov,NA RP<->RP,rpcov,NA RM<->RM,rmcov,NA MotI<->MotI,micov,NA MotE<->MotE,mecov,NA MIAO<->MIAO,miaocov,NA MIEE<->MIEE,mieecov,NA MIC<->MIC,miccov,NA MEIn<->MEIn,meincov,NA MEId<->MEId,meidcov,NA MERE<->MERE,merecov,NA SCI<->SCI,scicov,NA Support<->Rejection,suprej,NA Support<->Overprot,supoover,NA Rejection<->Overprot,rejover,NA MotI<->MotE,mots1,NA # Insert de covariance matrix cov.teste <- cov(teste, y = NULL, use = "everything", method = c("pearson", "kendall", "spearman")) cov.teste # Estimate de model (Here is where we have been finding difficulties) sem.wh.1 <- sem(mod.wh.1, cov.teste, N=126) summary(sem.wh.1) Obrigado!!! João Ricardo -- *Joao Ricardo N. Vissoci* Psicólogo CRP 08/12469 Prof. Ms. Faculdade Ingá Doutorando em Psicologia Social - PUCsp Grupo Pro-Esporte UEM/CNPq Tel. 4499298078/4399368588 joaovissoci@gmail.com/jrvissoci@ig.com.br proesporteuem.blogspot.com.br

2 2

Associação conjunta
by Edson Lira 18 Out '11

18 Out '11

Boa noite a todos, tenho um banco de dados com algumas variávei: AR(sim, não), AIN(sim, não, talvez), T_HIV(menos de 250, entre de 250 e 300, acima de 300), HISTO_A(normal, nia, nia I). Gostaria de avaliar a associação conjunta da variável AR com o as demais variáveis conjuntamente. Seria um modelo log-linear? Alguém poderia me dar alguma sugestão ou material para eu dar uma olhada? []'s Edson Lira Estatístico Manaus-Amazonas

3 3

Como obter a estrutura de um data.frame (linux)
by Jose Claudio Faria 18 Out '11

18 Out '11

Pessoal, Alguém poderia me informar uma função (no Linux) para exportar (por exemplo para a área de transferência) a estrutura de um data.frame. No Windows existe uma função (trExport) que acompanha o pacote TinnR que faz isso e não estou encontrando (ainda) nada similar no Linux. O objetivo é poder compartilhar dados reais (data.frame) pela lista de discussão sem a necessidade de conexões, read.table, etc. O que preciso é algo assim como o resultado para o data.frame BOD que fica no pacote datasets: structure(list(Time = c(1, 2, 3, 4, 5), demand = c(8,3, 10,3, 19, 16, 15.6, 19.8), .Names = c("Time", "demand"), row.names = c(NA, -6L), class = "data.frame") Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 2

Média e desvio de muitas variáveis de um fatorial!
by Ivan Bezerra Allaman 18 Out '11

18 Out '11

Bom dia senhores! Tenho um experimento fatorial 14x4, com 12 variáveis, e gostaria de saber se tem como colocar de uma só vez em um único quadro a média e o desvio padrão dos tratamentos e de todas as variáveis. Fiz algo como: medias <- apply(da[,-c(1:2,15)],2,function(x) tapply(x,da$Trat,mean,na.rm=TRUE)) desvios <- apply(da[,-c(1:2,15)],2,function(x) tapply(x,da$Trat,sd,na.rm=TRUE)) descrit <- list() for(i in 1:12){ descrit[[i]] <- cbind(medias[,i],desvios[,i]) } Mais sai como lista não sendo o ideal que procuro. Desde já agradeço pela atenção! \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 3

Dois plots no mesmo frame
by Marcelo Claro de Souza 18 Out '11

18 Out '11

Gostaria de plotar duas regressões lineares no mesmo frame. Para tal plotei os pontos utilizando os comandos abaixo: x1=data$SLA1 y1=data$LLS 1 x2=data$SLA2 y2=data$LLS2 plot(x1, y1, type='p', xlim=range(x1,x2), ylim=range(y1, y2), xlab='x', ylab='y') points(x2, y2, type='p', col="red", xlab='x', ylab='y') reg1=lm(y1~x1) reg2=lm(y2~x2) Mas ao tentar aplicar as linhas de tendência utilizando abline(reg1) e abline(reg2) não obtive resultado algum. Como devo proceder para adicionar as linhas das regressões a essas figuras? Muito obrigado pela ajuda. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

2 1

Strings
by Rodrigo Sant'Ana 18 Out '11

18 Out '11

Caros Companheiros da R-BR, estou com um problema, que provavelmente, seja muito simples de resolver. Mas confesso, já quebrei a cabeça para resolver e num consegui chegar ao resultado que desejo. Abaixo, tento replicar meu problema para que vocês possam me ajudar! Muito obrigado pela atenção de todos! # Eis uma réplica simplificada da base de dados! dados <- data.frame(municipio = c("A", "B", "C", "D", "E"), pesc.emb = c("PESCADOR JOAO","EMBARCACAO ZE","PESCADOR FELIPE","EMBARCACAO AMARGURA II","EMBARCACAO STAR")) dados municipio pesc.emb 1 A PESCADOR JOAO 2 B EMBARCACAO ZE 3 C PESCADOR FELIPE 4 D EMBARCACAO AMARGURA II 5 E EMBACACAO STAR # Preciso criar uma coluna classificando a informação da linha em "PESCADOR" ou "EMBARCACAO", ou seja, preciso extrair a string "PESCADOR" ou "EMBARCACAO" da coluna pesc.emb. Já tentei utilizar todas as funções do pacote stringr, porém não obtive exito, pois preciso ambos na mesma coluna. library(stringr) dados$classify <- str_extract(dados$pesc.emb, "PESCADOR") dados$classify2 <- str_extract(dados$pesc.emb, "EMBARCACAO") dados municipio pesc.emb classify classify2 1 A PESCADOR JOAO PESCADOR <NA> 2 B EMBARCACAO ZE <NA> EMBARCACAO 3 C PESCADOR FELIPE PESCADOR <NA> 4 D EMBARCACAO AMARGURA II <NA> EMBARCACAO 5 E EMBARCACAO STAR <NA> EMBARCACAO sei que pode parecer simples, mas estou enfrentando dificuldades em solucionar este problema. Fico super agradecido pela atenção de todos, abs _______________ Rodrigo Sant'Ana Oceanógrafo

6 10

Diagrama de Vernn (Package eVenn)
by Íris Cardoso 17 Out '11

17 Out '11

Olá, pessoal Primeiramente, sou nova na lista e gostaria de dizer que é muito bom encontrar um grupo com uma troca tão grande e rica de informações. E, como toda mensagem tem um "interesse" por trás, gostaria de esclarecer uma dúvida: Apanhei do R um bom tempo para conseguir plotar um diagrama de Venn (package eVenn). Sempre acusava um erro do tipo ausência de "TRUE/FALSE" nos dados. Optei por alterar o único ponto diferente dos dados dos exemplos, os nomes das variáveis. Depois de alterar os nomes para ficarem iguais aos nomes da base do exemplo "Liste_1_.xxx." (onde xxx = soma de 1's na variável), consegui plotar o diagrama. Porém ele plotou a 1ª coluna (id das observações) como uma variável do diagrama. Excluí a referida coluna e obtive o Diagrama do jeito que queria. Gostaria de saber se alguém já conseguiu plotar usando o "nome verdadeiro" das variáveis ou se essa alteração é mesmo necessária. Achei estranho pois, se o nome da variável fosse mesmo necessário nesse formato, ele não deveria plotar a coluna de id. Desde já agradeço a atenção. P.S.: A coluna de id estava com o mesmo nome do exemplo "row.names" e o conteúdo estava configurado como character. Íris Cardoso Estatística e Metre em Logística

2 2

metamodelagem
by Edmar Alino 17 Out '11

17 Out '11

Srs, Estou tentando criar metamodelos para otimizacao e me recomendaram os seguintes pacotes: - Redes Neurais (rbf e mlp): AMORE - GaSP: Dicekriging. Como sou iniciante no R, qualquer sugestao/recomendacao/orientacao pertinente a esses pacotes ou outros pacotes mais apropriados serao muito bemvindos. Grato desde ja, Edmar

2 1

Ajuda
by Paloma Palmieri 17 Out '11

17 Out '11

Prezados, Boa tarde! Estou com o seguinte problema: Tenho um programa em MATLAB que gera um arquivo .prn que é utilizado por um arquivo compilado em C (que é chamado em um outro arquivo de extensão .bat). O que acontece é que eu tenho que rodar os programas individualmente e isso não é muito eficiente. Neste caso, gostaria de ver se alguém da lista sabe como eu chamo e rodo um executável (ou .bat) em C dentro do MATLAB. Desde já agradeço. Abraços, Paloma Palmieri Alves paloma.palmieri(a)gmail.com Cel.: (21) 8811-5230 (21) 7961-8718

2 1