- R-br - DINF UFPR

Onde estão os autovalores do método "princomp"?
by Ivan Bezerra Allaman 30 Mar '11

30 Mar '11

Boa tarde senhores! Estou treinando uma análise multivariada com fins alométricos utilizando uma base de dados proveniente de um material (Multivariate allometry), e estou tendo dificuldades de encontrar alguns valores utilizando a função princomp. Antes que alguém questione ou recomende a utilização da função prcomp, estou usando a função princomp, porque ela calcula os componentes principais baseada na matrix de covariancia, que é obrigatório nos casos de estudos de crescimento alométricos. Pois bem, as dúvidas são as seguintes: 1) Não consegui achar nenhum comando para encontrar os valores absolutos dos autovalores. O summary no caso, só me retorna a proporção explicada por cada autovalor. Alguém sabe se é possível utilizando esta função ?(princomp). 2) Como é calculado os desvios padrões de cada componente ou autovetor? Já consultei o livro do Venables, e lá não tem nada a respeito. O livro do Mardia eu não tenho para fazer o estudo. Se alguém tiver este livro, talvez eu possa saber como foi feito estes cálculos. A referência é: Mardia, K. V., J. T. Kent and J. M. Bibby (1979). Multivariate Analysis, London: Academic Press. Desde já agradeço a atenção de todos. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Melhorias no Comando
by Victor Eduardo 30 Mar '11

30 Mar '11

Galera, estou tentando criar comandos que por eles, automaticamente abram pastas dentro de pastas para ir diretamente no arquivo (no caso txt), por exemplo, vou abrindo as seguintes pastas: Victor/teste/201005/20100520/2010052000.txt, sendo que 2010 é o ano, 05 o mes, 20 o dia e 00 a rodada (as rodadas vao de 00 até 23), no caso, esse é um exemplo, tenho que fazer para os demais meses, dias e rodadas. Queria saber uma estratégia de comando melhor, porque o que estou tentando montar não ficou muito legal, ficou assim: mes<-sprintf('%02d',seq(1,12))# sequencia de dois digitos de 1 ate 12 pasta<-paste(2010,mes,sep='') dia<-sprintf('%02d',seq(1,31)) # sequencia de 1 até 31 com dois dígitos dat.dir<-paste(dat.dir,pasta,'/',pasta,dia,'/',sep='') hh <-sprintf('%02d',seq(0,23)) # Cria uma sequência com 2 dígitos de 1 até 23 horario<-format(Sys.time(),'%Y%m') # formato ano/mes atual horario1<-format(Sys.time(),'%Y%m%d') #formato ano/mes/dia atual pa<-paste('2011',mes,sep='') #diretorio parcial pa1<-paste(dat.dir,pasta,'/',pasta,dia,'/',sep='') # diretorio completo do arquivo for(i in dia){ for(j in mes){ for(k in hh){ message(paste('C:/Users/max/Desktop/Victor/teste/',2010,j,'/',2010,j,i,'/',2010,j,i,k,sep='')) if(!file.exists(paste('C:/Users/max/Desktop/Victor/teste/',2010,j,'/',2010,j,i,'/',2010,j,i,k,sep=''))) next arquivo<-read.table(paste('C:/Users/max/Desktop/Victor/teste/',2010,j,'/',2010,j,i,'/',2010,j,i,k,sep='')) message('OK') } } } Abraços!

3 2

Res: [R_STAT] Teste de médias
by Walmes Zeviani 30 Mar '11

30 Mar '11

Douglas, A respeito do desbalanceamento, você terá que consultar materiais teóricos sobre o teste de Scott Knott e ver se ele se aplica à dados desbalanceados, se as taxas de erros são uniformes em relação ao desbalanceamento. Para o caso do teste de Tukey, se os efeitos são estimáveis (e são no seu caso), o contraste de médias pode ser aplicado sem problemas. A função agricolae::HSD.test() pode aplicar o teste usando a média harmônica do número de repetições. No entanto, essa opção é recomendada para desbalanceamentos amenos. Você pode também usar a função TukeyHSD(), que retorna o IC para a estimativa dos contrastes. Nessa última a representação não é com letras ao lado das médias porque o desbalanceamento implica que a diferença mínima significativa não seja constante. A [R-br] é a lista Brasileira oficial de usuários do R. Aconselho que faça inscrição nessa lista de e-mails conforme indicado pelo Ivan. À disposição. Walmes. ========================================================================== 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== ---------- Mensagem encaminhada ---------- De: walmes zeviani <walmeszeviani(a)yahoo.com.br> Data: 30 de março de 2011 14:37 Assunto: Enc: Res: [R_STAT] Teste de médias Para: walmeszeviani(a)gmail.com ----- Mensagem encaminhada ---- *De:* Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> *Para:* R_STAT(a)yahoogrupos.com.br *Enviadas:* Quarta-feira, 30 de Março de 2011 9:29:40 *Assunto:* Res: [R_STAT] Teste de médias Bom dia Fernando e Douglas! A maioria das pessoas que atuaram constantemente nesta lista, migraram para outro grupo do R, o "R-br". Caso tenham interesse em participar, segue o link com as instruções. Abraços. http://www.leg.ufpr.br/doku.php/software:rbr Allaman (S,f,P) *M.Sc Ivan Bezerra Allaman* Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA **email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924 ------------------------------ *De:* Fernando Henrique Toledo <fernandohtoledo(a)gmail.com> *Para:* R_STAT(a)yahoogrupos.com.br *Enviadas:* Quarta-feira, 30 de Março de 2011 9:09:09 *Assunto:* Re: [R_STAT] Teste de médias Douglas, Já que você está iniciando no programa vamos te mostrar o caminho das pedras! Veja a documentação do pacote "laercio", ou do "agricolae", em ambos há o teste de scott-knott implementado! Se você é novo no R, digite no prompt: >install.packages('laercio') e siga as instruções. Depois de intalado digite: >library('laercio') e depois procure a documentação da função através do comando: help(LScottKnott) Acho que você irá conseguir entender, de qualquer forma continuamos por aqui! att, FH 2011/3/30 Douglas Messias <dougmes(a)bol.com.br> > > > Pessoal, boa noite, esta é a primeira vez que me direciono ao grupo. > Estou com um problema aqui. Estou iniciando o uso do R, e preciso testar as > médias de uma anava que deu significativo, mas os dados são muito > desbalanceados e eu ainda estou fazendo tudo seguindo exemplos. > Gostaria de saber se alguem tem um exemplo de um teste de médias para eu > comparar meus dados, de preferencia o scott knott, que gosto muito por nao > encavalar valores. Meu experimento é em Dic com 5 tratamentos, mas um tem 32 > repetiçoes enquanto outro tem apenas 14. > Abraços... > > __._,_.___ <ivanalaman@yahoo.com.br?subject=Res%3A%20%5BR_STAT%5D%20Teste%20de%20m%C3%A9dias>| através de email<R_STAT@yahoogrupos.com.br?subject=Res%3A%20%5BR_STAT%5D%20Teste%20de%20m%C3%A9dias>| Responder através da web<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJya2hrYW01BF9TAzk3N…>| Adicionar um novo tópico<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmZTU2YmczBF9TAzk3N…> Mensagens neste tópico<http://br.groups.yahoo.com/group/R_STAT/message/15959;_ylc=X3oDMTM3NzNwcWxz…>( 3) Atividade nos últimos dias: - Novos usuários<http://br.groups.yahoo.com/group/R_STAT/members;_ylc=X3oDMTJnYXVpMTBqBF9TAz…> 8 Visite seu Grupo<http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmdnZ2MWQwBF9TAzk3NDkwND…> Tem muita gente querendo te conhecer! Que tal dar uma chance?<http://global.ard.yahoo.com/SIG=15m67uiub/M=758712.14532720.14422460.129601…> [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlNzlla2tvBF9TAzk3NDkwNDM1BGdycElkAz…> Trocar para: SÃ³ Texto<R_STAT-traditional@yahoogrupos.com.br?subject=Mudar%20Formato%20de%20Envio:%20Tradicional>, Resenha DiÃ¡ria<R_STAT-digest@yahoogrupos.com.br?subject=Envio%20de%20email:%20Resenha>• Sair do grupo <R_STAT-unsubscribe(a)yahoogrupos.com.br?subject=Sair%20do%20grupo> • Termos de uso <http://br.yahoo.com/info/utos.html> . __,_._,___

1 0

replica de data.frame
by Samuel Carvalho 30 Mar '11

30 Mar '11

Boa noite senhores(as) Por gentileza alguem poderia ajudar como faço para replicar um data.frame e "colar" estas observações replicadas abaixo das observações existentes? Por exemplo data <- data.frame(parcela=c(rep(1:5, each=4), rep(6, each=5))) data$resp <- rnorm(nrow(data)) data data2 <- rep(data, 2) dados <- cbind(data, data2) #aqui o data2 é inserido como colunas ao lado, porém gostaria que o data2 viesse colado abaixo do data dados Obrigado ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

3 2

Movimentar caixa de legenda em graficos
by Victor Eduardo 29 Mar '11

29 Mar '11

Galera, como faço para movimentar a caixa de legendas de um barplot? Estou colocando varios graficos numa imagem pelo par(mfrow), so que a caixa de legendas de alguns saem por cima da barra do grafico. Abraços

4 4

Pacote diagnosis
by Edson Lira 29 Mar '11

29 Mar '11

Boa tarde, com esta rotina, require(DiagnosisMed) ele<-matrix(c(55,7,49,84),nrow=2,ncol=2,byrow=T, dimnames=list(Eletro=c("Presente","Ausente"), Etsenose=c("Presente","Ausente"))) ele diagnosis(ele,plot=F) Tenho a saída --------------------------------------------------------------- Etsenose Eletro Presente Ausente Sum Presente 55 7 62 Ausente 49 84 133 Sum 104 91 195 The test has the following parameters [95% confidence interval] --------------------------------------------------------------- Sample size: 195 Prevalence considered(%): 46.67 Sensitivity(%): 92.31 [ 84.96 - 96.22 ] Specificity(%): 52.88 [ 43.36 - 62.20 ] Positive predictive value(%): 63.16 [ 54.70 - 70.88 ] Negative predictive value:(%): 88.71 [ 78.48 - 94.42 ] Positive likelihood ratio: 1.96 [ 1.58 - 2.44 ] Negative likelihood ratio: 0.15 [ 0.07 - 0.30 ] Diagnostic odds ratio: 13.28 [ 5.47 - 37.30 ] Error trade off (FN : FP) 0.14 : 1 Error rate(%): 28.72 [ 22.83 - 35.43 ] Accuracy(%): 71.28 [ 64.57 - 77.17 ] Youden index: 0.4519 [ 0.4555 - 0.4484 ] Area under ROC curve: 0.726 ----------------------------------------------------------- Quando chequei os resultados, vi que sensibilidade e especificidade estao invertidos. Vejam: Sensibiidade = a/(a+c)=55/104=0,5288. Na sáida está como especificidade, ou seja, invertido. O que estou fazendo errado. Valores corretos: sensib=0,5288 especif=0,9231 Obrigado a todos. Edson Lira Estatístico Manaus-Amazonas

2 1

sugestao de melhorias em codigo?
by Benilton Carvalho 29 Mar '11

29 Mar '11

Prezados, existe algo obvio (em termos de R, visto que gostaria de evitar o uso de codigo compilado) que eu esteja esquecendo em termos de otimizacao para a seguinte operacao? Ela se resume a: R[i,j] = mean( X[i,j] < Y[,j]) O CMR segue: set.seed(1) X <- matrix(rnorm(100), nc=5) Y <- matrix(rnorm(50), nc=5) f = function(i, x, y){ n = ncol(x) sapply(1:n, function(j, v, ref) mean(ref[,j] > v[j]), x[i,], y) } R = t(sapply(1:nrow(X), f, X, Y) O problema que surge para mim refere-se ao tempo de execucao quando tenho q lidar com X com dimensoes 5e6 x 1e3.... pois, naturalmente, a solucao acima levara' um bom tempo para finalizar... obrigado, benilton

3 7

[OFF] [Dúvida] Atualizar o Ubuntu
by Pedro Rafael 28 Mar '11

28 Mar '11

Pessoal como sabemos a cada 6 meses é lançada uma nova versão do Ubuntu. Como faço para atualizar o sistema da versão 10.10 para a 11.04 a partir do momento que esta for lançada? Essa atualização pode me gerar problemas ou é completamente segura? Acho muito chato ter que sempre a cada versão formatar o computador para ter uma instalação limpa. Sou do tipo do usuário que sempre quero o SO mais recente contudo não acho vantajoso ter sempre que formatar o sistema. Tenho um certo medo quando o sistema vai atualizar o kernel possa dar alguns problemas de compatibilidade com alguns hardwares.... -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

3 4

[Dúvida] Ler dados direto no HD
by ivanalaman 28 Mar '11

28 Mar '11

Pessoal como faço para ler uma banco de dados com extensão .dbf ou qualquer outra base de dados sem que o R a jogue na RAM? Tenho um banco de dados com aproximadamente 700 MB e quando leio o arquivo com a função read.bdf()percebo que o R joga a base na memória RAM o que prejudica o desempenho do sistema. Softwares simples como o Tabwin por exemplo não demora à carregar essa base nem trava o computador. Será que o R não tem como carregar essa base sem tantos problemas? Saudações à todos, Pedro Rafael Diniz Marinho -- View this message in context: http://r-br.2285057.n4.nabble.com/Duvida-Ler-dados-direto-no-HD-tp3400886p3… Sent from the R-br mailing list archive at Nabble.com.

7 16

R2wd escrevendo resultados em .doc
by Pedro Emmanuel Alvarenga Americano do Brasil 28 Mar '11

28 Mar '11

Caros amigos, Gostaria de poder utilizar o R2wd para poder fazer um pequeno relatorio de inclusão de voluntarios de uma pesquisa a cada 15 dias. No entando me deparei com o seguinte problema... > library(R2wd) Carregando pacotes exigidos: rcom Carregando pacotes exigidos: rscproxy rcom requires a current version of statconnDCOM installed. To install statconnDCOM type installstatconnDCOM() This will download and install the current version of statconnDCOM You will need a working Internet connection because installation needs to download a file. > wdGet() Erro em wdGet() : tentativa de aplicar uma não-função sendo que essa é a inicialização após a instalação do installstatconnDCOM(), ou seja, toda vez que eu inicio o R, ele pede pra instala o statconnDCOM. O problema é que esse erro nao dá a menor pista pra mim do como eu poderia resolver a questão. Me parece que o aplicativo instalado com essa função nao está no local certo ou configurado de forma adequada para que o R2wd consiga chamar o Word. Ja tentei mandar um email para o desenvolvedor, mas nao recebi resposta até agora. Alguma luz? Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365 Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.openoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

3 4

Re: [R-br] Blocos incompletos no R
by Julia Maria 28 Mar '11

28 Mar '11

Walmes, O experimento da tese que tem os dados, é um pouco diferente do experimento que eu citei na mensagem acima. O autor utilizou aquele experimento acima como idéia pra obter o modelo, as equações normais, etc. O que ele analisou e tem resposta na tese é o seguinte: V=6 tratamentos principais (Dietas) k=4 tratamentos por bloco r=10 repetições por tratamento a=15 blocos lambda=6 ocorrencias de uma mesmo par de tratamentos u=3 tratamentos secundários, que são idades. Entao ele constrou o bloco incompleto utilizando os 6 tratamentos principais (utilizando somente 4 deles em cada bloco) e os tratamentos secundários aleatorizados dentro dos principais normalmente. O autor não apresenta as médias ajustadas, somente os efeitos ajustados dos tratamentos, que eu coloquei na sheet 2 conjuntamente com a analise de variancia. Pensei em calcular eu mesma, apenas somando os efeitos à média geral, mas para esse caso não sei se seria correto somar tanto para os efeitos principais quanto os secundários. De toda a forma, segue anexo os dados. Muito obrigada pela ajuda, eu gostaria muito de entender e aprender a fazer uma análise assim no R. --- Em sex, 25/3/11, Walmes Zeviani <walmeszeviani(a)gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani(a)gmail.com> Assunto: Re: [R-br] Blocos incompletos no R Para: r-br(a)listas.c3sl.ufpr.br Data: Sexta-feira, 25 de Março de 2011, 12:25 Júlia, Envie os dados em *.xls com colunas identificando as fontes de variação e reposta (sheet1). Vou tentar implementar e discutimos. Vou precisar da anova e médias ajustadas dos tratamentos para conferir. Pegue na tese e envie junto com os dados (sheet2). Procedimento para análise intra/interbloco de experimento com blocos incompletos tipo I, II e III então disponíveis nas Rídiculas do LEG. Confira os meus procedimentos com outros dados e outros programas. À diposição. Walmes. =========================================================== 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 =========================================================== -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

1 1

[Dúvida] Bibliotecas gráficasno R
by Pedro Rafael 27 Mar '11

27 Mar '11

Pessoal qual biblioteca para fazer gráficos descritivos simples no R? Muitas vezes preciso de gráficos simples descritivos como gráficos de setores, barras, linhas e para essas tarefas e análises simplórias contudo muitas vezes necessárias acabo por usar o Excel por gostar da qualidade dos gráficos e texturas. No R acho os gráficos por exemplo o de setores muito simples sem muita graça que as vezes é necessária para fazer uma apresentação ou coisas parecidas. No trabalho preciso muitas vezes mostrar graficos mais apresentáveis mais atraentes apesar de não gostar muito dessas frescuras. Outro motivo é que não queria usar um software pago como o Excel para isso. Qual ou quais bibliotecas no R eu posso utilizar para fazer gráficos mais sofisticados ou me indiquem outro programa que eu possa fazer isso...Fazer o que o Excel já faz na parte gráfica é o suficiente, ou seja, muito pouco, contudo suficiente para mim. -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

6 6

Calcular a área em Km2 a partir de um shapefile
by Carolina Pena 27 Mar '11

27 Mar '11

Alguém conhece algum pacote/função no R capaz de calcular a área em km2 de um mapa (shapefile). Eu li o mapa no R utilizando o comando 'readShapePoly'; consigo acessar a "area" (slots) do objeto através dos seguintes comandos: area = sapply(slot(map2, "polygons"), function(i) slot(i, "area")) No entanto, não sei como transformar essa informação em km2 Obrigada, Carolina

4 4

Re: [R-br] Pedindo arquivos de dados livro Johnson e Wichern
by Jose Claudio Faria 26 Mar '11

26 Mar '11

Olá, "Ao pessoal das listas muito obrigado! Minha intenção não era pedir coisa ilegal, acreditava que o livro deveria ter uma página onde permite-se que baixe dos data sets e os códigos fontes, como costuma acontecer, ou eventualmente alguém que já tivesse resolvido alguns exercícios de lá que os fornecesse!" Durante meu pós-doc sabático na ESALQ (2008/2009), um dos objetivos era preparar material didático em R sobre multivariada. O livro adotado pelo prof. Tadeu foi Johnson, R.A e Wichern, 6 ed. Uma parte considerável dos exercícios propostos dos capítulos 1-4 foram resolvidos em R (Sweave). Creio que este material seria bastante interessante e útil para quem está migrando ou se familiarizando com o R, no que diz respeito a mutlivariada, pois abrange os fundamentos. Seria interessante disponibilizar este material em um local onde fosse possível (aos interessados) terem acesso, assim como, poderem ir expandindo (acrescentando soluções faltantes, progressivamente). Alguém teria uma boa idéia a respeito no que diz respeito a um servidor? Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

3 2

Pedindo arquivos de dados livro Johnson e Wichern
by Mauro Sznelwar 26 Mar '11

26 Mar '11

Alguém teria os arquivos dos dados do livro: Johnson, R.A e Wichern, D.W., Applied Multivariate Statistical Analysis da 5ª ou 6ª edição? Tem algum lugar daonde baixá-los?

3 5

Como mudar o angulo do texto no eixo X do comando stripchart?
by Augusto Ribas 25 Mar '11

25 Mar '11

Opa, alguem pode dar uma mãozinha aqui. Eu to tentando fazer um grafico com o comando stripchart, mas o nome dos levels do fator são longos. Dai eu quero deixa eles inclinados, mudar o angulo com o qual o texto é escrito no eixo X, mas não estou conseguindo. No comando Text é com o argumento str, mas eu nao to conseguindo fazer dentro do stripchart, e com o comando text so da pra escrever dentro da area de plotagem. Tem algum argumento dentro do stripchart que faz isso, ou jeito facil de fazer isso? #Exemplo de Dados Tratamento<-rep(c("Controle","Ácido a concentração 10","Ácido a concentração 20","Estresse"),each=30) Velocidade<-rpois(120,18) #Grafico stripchart(Velocidade~Tratamento,vertical=T,pch=19,method="stack") -- Grato Augusto C. A. Ribas

3 4

continuando... existe algo como o argumento datevars do get.spss...
by Paula Mendes Luz 25 Mar '11

25 Mar '11

eu de novo agora com as vars de data... existe algo como o argumento datevars do get.spss para o read.csv eu sei que tem o as.Date para editar a var depois mas sempre levo um certo tempo para fazer isso funcionar... nesse momento o banco é lido e as vars de data sao importadas como: Factor w/ 911 levels "","1/1/1997 ",..: 127 697 756 421 694 402 111 430 150 611 ... se eu aplico o as.Date depois nao dá mto certo... (b1$dob <- as.Date(b1$dob, "%b/%d/%Y")) acho mais facil dizer que é data na importacao... tem como? obrigada pela ajuda.

2 1

[[OFF-TOPIC]] curso do Edzer no LEG e pré-curso
by Paulo Justiniano 25 Mar '11

25 Mar '11

Caros O curso do Edzer está confirmado com 1hr na terca dia 29 e 4 hrs na quarta dia 30. A fim de proporcionar um melhor rendimento e ter certeza que todos tenha condição de acompanhar com os recursos basicos e se familiarizar e testar o ambiente computacional faremos um pre-curso na segunda, 28/01 as 16:30 - 17:30 onde abordaremos 1. recursos computacionais apra o curso (voces podem trazer notebooks ou usar os terminais (linux do LABEST) 2. representacao de dados espaciais no R (alternativas, pacotes etc) aproveitando para falar de classes S3 e S4 3. Basico de analise espacial/geoestatistica de dados com sp, geoR e gstat 4. Duvidas/discussao Paulo Justiniano Ribeiro Jr LEG (Laboratorio de Estatistica e Geoinformacao) Universidade Federal do Parana Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR - Brasil Tel: (+55) 41 3361 3573 VOIP: (+55) (41) (3361 3600) 1053 1066 Fax: (+55) 41 3361 3141 e-mail: paulojus AT ufpr br http://www.leg.ufpr.br/~paulojus

2 1

usar letra minúscula ao ler um banco
by Paula Mendes Luz 25 Mar '11

25 Mar '11

oi pessoal alguém sabe dizer se há algo parecido com o argumento "lowernames = T" da função spss.get (que transforma o nome das variáveis em letra minúscula) para o read.csv ou read.xls obrigada! paula.

4 4

formatar datas e horas
by Guilherme Moraes Ferraudo 25 Mar '11

25 Mar '11

Boa noite a todos, por favor, peço a ajuda de vcs. estou tentando plotar a latitude e longitude filtrando por dia e período (horas), mas não consigo trabalhar com as horas. Por exemplo, gostaria de plotar 2 gráfico para o dia 28/04. Um contendo os pontos menor que 6horas e o outro contendo os pontos superiores a 6 horas. Será que fui claro o suficiente??? Obrigado. Abs -- Guilherme Moraes Ferraudo Jaboticabal/Campinas - SP http://lattes.cnpq.br/2096118558794430

4 5

Fwd: Res: [R_STAT] [Dúvida] Como responder e postar mensagens na R-br
by Pedro Rafael 24 Mar '11

24 Mar '11

Obrigado pelas informações. ---------- Mensagem encaminhada ---------- De: Pedro Rafael <pedro.rafael.marinho(a)gmail.com> Data: 24 de março de 2011 17:36 Assunto: Re: [R-br] Res: [R_STAT] [Dúvida] Como responder e postar mensagens na R-br Para: "Walmes Zeviani 2 [via R-br]" < ml-node+3403695-178076951-223914(a)n4.nabble.com> Obrigado pela dica Walmes. Em 24 de março de 2011 17:27, Walmes Zeviani 2 [via R-br] < ml-node+3403695-178076951-223914(a)n4.nabble.com> escreveu: Pedro, > > A postagem na lista é por e-mail. > O nabble é um portal que amazenamento busca de mensagem. Permite que > usuário faça conta, upload de avatar e inclusive que poste por lá. No > entanto, eu alterei as opções de membro para que ninguém poste via nabble. > Fiz isso porque tava dando confusão. As postagens via nabble só são aceitas > pela lista se o e-mail usado para criar a conta no nabble for o mesmo e-mail > de inscrição na lista. Caso contrário o administrador recebe e-mail de > notificação e o post via nabble não é distribuido para o e-mail dos > incritos, ficando presente só no nabble. Ou seja, implica que nabble, > mail-archive e lista tenham conteúdo diferentes, e não queremos isso. Então, > poste sempre via e-mail. > > Essas instruções estão na pagina www.leg.ufpr.br/rbr. No momento, estamos > editando um guia de postagem que para a lista. > > À disposição. > Walmes. > > > -- > =========================================================== > 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 > VoIP: (3361 3600) 1053 1173 > e-mail: [hidden email]<http://user/SendEmail.jtp?type=node&node=3403695&i=0&by-user=t> > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > =========================================================== > > *De:* pedrorafaelmarinho <[hidden email]<http://user/SendEmail.jtp?type=node&node=3403695&i=1&by-user=t> > > > *Para:* [hidden email]<http://user/SendEmail.jtp?type=node&node=3403695&i=2&by-user=t> > *Enviadas:* Quinta-feira, 24 de Março de 2011 15:46:36 > *Assunto:* [R_STAT] [Dúvida] Como responder e postar mensagens na R-br > > > > Pessoal como funciona a postagem de mensagens nessa nova comunidade? É via > e-mail ou por esta página (http://r-br.2285057.n4.nabble.com/<http://r-br.2285057.n4.nabble.com/?by-user=t> > )? > Estou atrapalhado com a lista :D. Pela página não estou conseguindo postar > as mensagens...como devo fazer? > > Abraço, > Pedro Rafael > > > _______________________________________________ > R-br mailing list > [hidden email]<http://user/SendEmail.jtp?type=node&node=3403695&i=3&by-user=t> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > If you reply to this email, your message will be added to the discussion > below: > > http://r-br.2285057.n4.nabble.com/R-br-Res-R-STAT-Duvida-Como-responder-e-p… > To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. > > -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6 -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

2 2

Res: [R_STAT] [Dúvida] Como responder e postar mensagens na R-br
by Walmes Zeviani 24 Mar '11

24 Mar '11

Pedro, A postagem na lista é por e-mail. O nabble é um portal que amazenamento busca de mensagem. Permite que usuário faça conta, upload de avatar e inclusive que poste por lá. No entanto, eu alterei as opções de membro para que ninguém poste via nabble. Fiz isso porque tava dando confusão. As postagens via nabble só são aceitas pela lista se o e-mail usado para criar a conta no nabble for o mesmo e-mail de inscrição na lista. Caso contrário o administrador recebe e-mail de notificação e o post via nabble não é distribuido para o e-mail dos incritos, ficando presente só no nabble. Ou seja, implica que nabble, mail-archive e lista tenham conteúdo diferentes, e não queremos isso. Então, poste sempre via e-mail. Essas instruções estão na pagina www.leg.ufpr.br/rbr. No momento, estamos editando um guia de postagem que para a lista. À disposição. Walmes. -- =========================================================== 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 =========================================================== *De:* pedrorafaelmarinho <pedrorafaelmarinho(a)yahoo.com.br> *Para:* R_STAT(a)yahoogrupos.com.br *Enviadas:* Quinta-feira, 24 de Março de 2011 15:46:36 *Assunto:* [R_STAT] [Dúvida] Como responder e postar mensagens na R-br Pessoal como funciona a postagem de mensagens nessa nova comunidade? É via e-mail ou por esta página (http://r-br.2285057.n4.nabble.com/)? Estou atrapalhado com a lista :D. Pela página não estou conseguindo postar as mensagens...como devo fazer? Abraço, Pedro Rafael

1 0

escala gráfica e Norte
by Mario Azevedo 24 Mar '11

24 Mar '11

Eu estou usando o R, com o pacote PBSmapping, para produzir uns Mapas de Moran e LISA. Eu gostaria de saber se é possível e como adicionar uma escala gráfica e um símbolo indicando o norte nesses mapas. Alguém pode me ajudar? Existe outro pacote melhor para essa tarefa? Obrigado. Mário

1 0

Testando a criação de um novo tópico.
by Pedro Rafael 24 Mar '11

24 Mar '11

Esse tópico é um simples teste. -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

2 1

Dúvida
by Pedro Rafael 24 Mar '11

24 Mar '11

Testando envio -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 0

Ler linhas
by Victor Eduardo 24 Mar '11

24 Mar '11

Como faço para ler linhas de arquivos que começam somente pelos numeros 82 ou 83? Por exemplo, vou colocar alguns dados aqui: V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 8345 2010 5 20 15 0 -38.20 -12.70 11 -9999 1016 1017 -9999 -9999.0 309.2 301.4 5122 2010 5 20 15 0 -43.80 -21.20 1000 -9999 887 998 284 2.5 292.1 289.1 8222 2010 5 20 15 0 -43.30 -22.60 33 -9999 1012 1016 131 0.4 294.8 292.6 8343 2010 5 20 15 0 -37.00 -10.90 5 -9999 1012 1013 -9999 -9999.0 301.6 296.9 4107 2010 5 20 15 0 -38.20 -9.40 253 -9999 984 1012 109 3.8 300.6 293.9 No caso, o meu interesse está na V1, quero pegar os numeros que comecam por 82 ou 83 e ao pegar, queria que os dados viessem dessa maneira, por exemplo: 8345 2010 5 20 15 0 -38.20 -12.70 11 -9999 1016 1017 -9999 -9999.0 309.2 301.4 8222 2010 5 20 15 0 -43.30 -22.60 33 -9999 1012 1016 131 0.4 294.8 292.6 8343 2010 5 20 15 0 -37.00 -10.90 5 -9999 1012 1013 -9999 -9999.0 301.6 296.9

5 7

Email teste (Como beber um litro de jurubeba numa golada só? )
by Ivan Bezerra Allaman 24 Mar '11

24 Mar '11

Senhores, este é um email teste visto que alguns membros não estão conseguindo postar suas dúvidas via email. Por favor, o mantenedor do grupo apague-a após recebe-la. Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Dúvida no rotina rcorr no pacote Hmisc
by Edson Lira 24 Mar '11

24 Mar '11

Como faço para colocar os nomes das variáveis na rotina abaixo. Faço a rotina para qa correlação:require(Hmisc) rcorr(cbind(corre$prevalência,corre$precipitação,corre$temp_max, + corre$temp_min) + ,type=c("pearson")) Esta saída é gerada. Como faço para colocar os nomes da variável. [,1] [,2] [,3] [,4] [1,] 1.00 -0.04 -0.15 -0.28 [2,] -0.04 1.00 -0.80 -0.41 [3,] -0.15 -0.80 1.00 0.67 [4,] -0.28 -0.41 0.67 1.00 n= 63 P [,1] [,2] [,3] [,4] [1,] 0.7841 0.2516 0.0249 [2,] 0.7841 0.0000 0.0010 [3,] 0.2516 0.0000 0.0000 [4,] 0.0249 0.0010 0.0000 Edson Lira Estatístico Manaus-Amazonas

2 2

Ajuda com modelo Log-linear (agora com dados)
by Diogo Borges Provete 24 Mar '11

24 Mar '11

Olá a todos novamente e desculpem-me pela submissão anterior, O problema que tenho é o seguinte: gostaria de uma ajuda para ajustar um modelo log-linear para dados com 3 categorias de pigmentação na superfície de órgãos da cavidade abdominal de anfíbios. Os dados em anexo referem-se ao coração. Os dados brutos consistem das categorias de pigmentação (0, 1 e 2) de 10 indivíduos divididos nos grupos experimentais (nas colunas). Eu gostaria de saber se há diferença nas categorias de pigmentação de acordo com os tratamentos experimentais. A tabela de contingência tem a informação da frequência das categorias de pigm distribuídas nos grupos experimentais. Eu usei a função loglm do pacote MASSque permite que se insira a fórmula parecida com o lm e glm. Eu inclui somente os dados pq a função diz que se a fórmula não for especificada pode-se colocar somente a origem dos dados. Mas não soube interpretar os resultados. Agradeço desde já Atenciosamente, Diogo Borges Provete ============================== Biólogo Mestre em Biologia Animal (UNESP) Doutorando PPG Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás - UFG Goiânia-GO CP: 131 74001-970 Brazil Tel. +55 62 3521-1732 Cel. +55 62 8231-5775 Skype: diogoprovete MSN: diogoprov(a)yahoo.com.br Personal web page Traduza conosco: <<<American Journal Experts>>> <<<D-Lang Soluções linguisticas>>> <<<Perfil no ProZ.com como tradutor freelancer>>> ==============================

1 0

Estimativa do intervalo de confianca em um modelo ZINB
by Humberto Hazin (hotmail) 24 Mar '11

24 Mar '11

Olá a todos, Tenho 2 problemas se alguém puder me ajudar (Um exemplo dos dados encontra-se em anexo): 1- como faço para estimar os intervalos de confiança dos coeficientes de um modelo ZINB. Estou usando o package (pscl). Já tentei por bootstrap mais ele não estima os intervalos. Segue abaixo o script que usei data(bum2) F1a<-formula(BUM ~ Year + Quarter + area + Strategy+offset(log(effort1))) options(contrasts=c("contr.treatment","contr.poly")) ZINB<-zeroinfl(F1a,dist="negbin",link = "logit",data=bum2) boot.zinb <- function(data, indices){ require(pscl) data <- data[indices,] try(mod <- zeroinfl(F1a, data = data, link = "logit", dist = "negbin", trace = FALSE EM = FALSE)) if (exists("mod")) { coef(mod) } else { NA } } system.time(zinb.boot.out <- boot(data = bum2, statistic = boot.zinb, R = 99,sim="permutation")) zinb.boot.out boot.ci(zinb.boot.out,type="bca") 2- como faço para estimar os predicts e os intervalos de confiança apenas da variável YEAR de um modelo ZINB. Abaixo o código usado: x<-seq(1980,2008) newdata=cbind.data.frame(Year=as.factor(x),Quarter=as.factor(rep(1,length(x) )),area=as.factor(rep(1,length(x))),Strategy=as.factor(rep(4,length(x))),eff ort1=rep(1000,length(x))) newdata fm_zip <- zeroinfl(F1a,dist="negbin",link = "logit",data=bum2) fit <- predict(fm_zip) Pearson <- resid(fm_zip, type = "pearson") VarComp <- resid(fm_zip, type = "response") / Pearson year <- bum2$Year bootstrap <- replicate(999, { yStar <- pmax(round(fit + sample(Pearson) * VarComp, 0), 0) predict(zeroinfl(yStar ~ year | 1), newdata = newdata) }) newdata0 <- newdata newdata0$fit <- predict(fm_zip, newdata = newdata, type = "response") newdata0<- t(apply(bootstrap, 1, quantile, c(0.025, 0.975))) Desde Ja agradeço a vossa atenção Humberto

1 0

Re: [R-br] [Dúvida] Ler dados direto no HD
by beniltoncarvalho 23 Mar '11

23 Mar '11

Eu nao mencionei converter p Excel :) . Ha' vezes que conversao para txt e' um mal necessario, mas se vc estiver em Windows pode ser que haja problemas tambem (depende do filesystem que vc estiver usando). De qualquer forma, voce pode usar o pacote RODBC com os comandos odbcConnectDbase e sqlFetch. Verifique o arquivo de ajuda para os detalhes. -- View this message in context: http://r-br.2285057.n4.nabble.com/Duvida-Ler-dados-direto-no-HD-tp3400886p3… Sent from the R-br mailing list archive at Nabble.com.

1 0

Ajuda com modelo Log-linear
by Diogo Borges Provete 23 Mar '11

23 Mar '11

Olá a todos, gostaria de uma ajuda para ajustar um modelo log-linear para dados com 3 categorias de pigmentação na superfície de órgãos da cavidade abdominal de anfíbios. Os dados seguem em anexo, a planilha é referente ao coração, a parte da esquerda são os dados originais, que consistem das categorias de pigmentação (0, 1 e 2) de 10 indivíduos divididos nos grupos experimentais (nas colunas). Eu gostaria de saber se há diferença nas categorias de pigmentação de acordo com os tratamentos experimentais. O pacote MASS tem uma função loglm que permite que se insira a fórmula parecida com o lm e glm, mas não sei como usar isso neste caso. Agradeço desde já Atenciosamente, Diogo Borges Provete ============================== Biólogo Mestre em Biologia Animal (UNESP) Doutorando PPG Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás - UFG Goiânia-GO CP: 131 74001-970 Brazil Tel. +55 62 3521-1732 Cel. +55 62 8231-5775 Skype: diogoprovete MSN: diogoprov(a)yahoo.com.br Personal web page Traduza conosco: <<<American Journal Experts>>> <<<D-Lang Soluções linguisticas>>> <<<Perfil no ProZ.com como tradutor freelancer>>> ==============================

2 1

visita Edzer pebesma ao LEG e mini-curso
by Paulo Justiniano 23 Mar '11

23 Mar '11

Caros Na proxima semana, entre 29-31/03/2011 estaremos recebendo no LEG o Edzer Pebesma Edzer Pebesma Institute for Geoinformatics (ifgi) University of Münster Germany. Phone: +49 251 O Edzer é o autor do pacote gstat, co-autor do pacote pacote sp e co-autor do livro: Analysis of Spatial Data with R (ASDAR) (com Roger Bivand e Virgilio G. Rubio) Durante a visita serah ministrado um curso de "Spatial temporal analysis" 1 hr na terca a tarde (29) 4 hrs na quata de manha (30) e se necessario 1 hr extra na quinta (31) (na segunda 28/03) rodaramos uma intro sobre recursos antes do inicio do curso Paginas relevantes: Instituto: http://ifgi.uni-muenster.de Pagina do Edzer no Instituto: http://ifgi.uni-muenster.de/en/staff/7-mitarbeiterinnen/3-pebesma-edzer-pro… Lista de Publicações: http://ifgi.uni-muenster.de/~epebe_01/publ.html Projeto/Comunidade http://www.52north.org/geostatistics Paulo Justiniano Ribeiro Jr LEG (Laboratorio de Estatistica e Geoinformacao) Universidade Federal do Parana Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR - Brasil Tel: (+55) 41 3361 3573 VOIP: (+55) (41) (3361 3600) 1053 1066 Fax: (+55) 41 3361 3141 e-mail: paulojus AT ufpr br http://www.leg.ufpr.br/~paulojus

1 0

Problemas ao instalar o pacote rgl
by ivanalaman 23 Mar '11

23 Mar '11

Olá senhores de novo. Estou com problemas há alguns dias ao instalar pacotes que dependem do "rgl". Este pacote ("rgl") não está instalando de jeito nenhum. Já tentei até via terminal (R CMD blbaalb) e nada. As mensagens são: * installing to library ‘/home/allaman/R/i486-pc-linux-gnu-library/2.10’ * installing *source* package ‘rgl’ ... checking for gcc... gcc -std=gnu99 checking for C compiler default output file name... a.out checking whether the C compiler works... yes checking whether we are cross compiling... no checking for suffix of executables... checking for suffix of object files... o checking whether we are using the GNU C compiler... yes checking whether gcc -std=gnu99 accepts -g... yes checking for gcc -std=gnu99 option to accept ISO C89... none needed checking how to run the C preprocessor... gcc -std=gnu99 -E checking for gcc... (cached) gcc -std=gnu99 checking whether we are using the GNU C compiler... (cached) yes checking whether gcc -std=gnu99 accepts -g... (cached) yes checking for gcc -std=gnu99 option to accept ISO C89... (cached) none needed checking for libpng-config... yes configure: using libpng-config configure: using libpng dynamic linkage checking for X... libraries , headers checking GL/gl.h usability... no checking GL/gl.h presence... no checking for GL/gl.h... no checking GL/glu.h usability... no checking GL/glu.h presence... no checking for GL/glu.h... no configure: error: missing required header GL/gl.h ERROR: configuration failed for package ‘rgl’ * removing ‘/home/allaman/R/i486-pc-linux-gnu-library/2.10/rgl’ Agradeço mais uma vez se puderem ajudar. Allaman (S,f,P) -- View this message in context: http://r-br.2285057.n4.nabble.com/Problemas-ao-instalar-o-pacote-rgl-tp3399… Sent from the R-br mailing list archive at Nabble.com.

3 4

graficos de regressão linear
by gilsonsch 23 Mar '11

23 Mar '11

Pessoal, to precissando fazer um grafico de analises de regressão de tres espécies dende, caiaue, e hi, eu já conseguir fazer os grafico, mas eu quero que as três linhas estejam num grafico só. o mesmo quando eu quero grafico de dispersão, quero colocar mas de uma variave no mesmo gráfico, como eu faço isso. exemplo: dende 2 dende 2 dende 3 caiuae 4 caiaue 3 caiaue 3 hi 2.5 hi 3.4 hi 5.4 Gilson -- View this message in context: http://r-br.2285057.n4.nabble.com/graficos-de-regressao-linear-tp3400025p34… Sent from the R-br mailing list archive at Nabble.com.

3 4

Ajuda com a função biplot.
by ivanalaman 23 Mar '11

23 Mar '11

Boa tarde senhores!! Estou com dois problemas com a função biplot. O primeiro é o seguinte: Quando tento colorir os pontos referentes a cada grupo, só sai uma cor em todos os grupos. Vejam. data(iris)#Já está na base do R pca <- prcomp(~Petal.Length + Sepal.Length + Petal.Width + Sepal.Width, data=iris, scale=TRUE) biplot(pca, xlabs=rep(c("1","2","3"),c(50,50,50)),col=c("red","blue","green"))#só é plotada a cor "red". Um outro problema, é que quando tento fazer a mesma coisa na minha base de dados, a função solta o seguinte erro. head(dmof) LIN IDADE BL REP TEMPO ESTA CP CC AC A1 A2 A3 D1 D2 L1 1 Verm 85 1 1 1 outono 8.3 25.9 26.4 32.1 28.1 11.8 38.2 11.9 16.5 2 Verm 85 1 2 1 outono 10.5 31.1 31.3 42.2 37.9 14.3 45.2 18.4 19.0 3 Verm 85 1 3 1 outono 9.2 28.1 28.4 36.7 30.1 11.8 42.5 16.2 17.1 4 Verm 85 1 4 1 outono 9.8 28.2 27.4 38.1 30.7 12.4 39.7 16.3 16.8 5 Verm 85 1 5 1 outono 9.1 28.4 29.2 36.2 31.2 12.4 42.7 18.2 18.4 6 Verm 85 1 6 1 outono 8.5 25.3 26.1 34.1 29.8 11.7 35.5 20.4 17.2 L2 L3 1 9.2 4.6 2 10.8 6.0 3 9.3 4.6 4 9.2 4.5 5 9.6 5.0 6 10.2 4.5 dmof$LIN Levels: Com Tai UFLA Verm pca1 <- prcomp(~CP+CC+AC+A1+A2+A3+D1+D2+L1+L2+L3, data=dmof, scale=TRUE) biplot(pca1, xlabs=rep(c("1","2","3","4"),c(200,200,200,200)), xlim=c(-0.1,0.1)) Erro em dimnames(x) <- list(xlabs, dimnames(x)[[2L]]) : comprimento de 'dimnames' [1] não é igual ao tamanho do array Calls: biplot -> biplot.prcomp -> biplot.default Não entendo o erro, pois estou obedecendo as dimensões da minha base de dados. dim(dmof) 800 17 Se alguém puder me ajudar eu agradeço. Allaman (S,f,P) -- View this message in context: http://r-br.2285057.n4.nabble.com/Ajuda-com-a-funcao-biplot-tp3399760p33997… Sent from the R-br mailing list archive at Nabble.com.

2 1

[OFF_TOPIC] - Renomear arquivos linux
by Eder David Borges da Silva 22 Mar '11

22 Mar '11

Pessoal, sei que não é duvida R, mas creio que vcs podem me ajudar. em uma pastas tenho uma sequencia de arquivos (a.nc ,b.nc,...,z.nc) não necessariamente em ordem alfabetica. se fosse para mudar *.nc para *.nc1 isso funcionaria. rename 's/\.nc/\.nc1/' *.nc só que querro algo mais, querro que fique aB.nc, bB.nc, ...xB.nc, ou seja antes do ponto da extensão do arquivo querro por um "B" por exemplo. Alguma dica? Atenciosamente ----- ------------------------------ Éder David Borges da Silva Eng° Agrônomo ederdbs(a)gmail.com eder(a)leg.ufpr.br "Tudo é relativo" Albert Einstein -- View this message in context: http://r-br.2285057.n4.nabble.com/OFF-TOPIC-Renomear-arquivos-linux-tp33973… Sent from the R-br mailing list archive at Nabble.com.

3 3

Teste média Subdividida
by Rafael Tieppo 22 Mar '11

22 Mar '11

Bom dia, Estou com dificuldade para realizar o teste de média para um experimento em parcelas subdivididas, conforme segue abaixo: summary(aov(Ds ~ Bloco + Amb*Prof + Error(Bloco/Amb), data = ATRIB)) Error: Bloco Df Sum Sq Mean Sq Bloco 5 0.0091 0.00182 Error: Bloco:Amb Df Sum Sq Mean Sq F value Pr(>F) Amb 2 1.20123 0.60062 34.396 3.293e-05 *** Residuals 10 0.17462 0.01746 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Prof 3 0.077294 0.0257648 8.8629 9.972e-05 *** Amb:Prof 6 0.045089 0.0075148 2.5850 0.03082 * Residuals 45 0.130817 0.0029070 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Até essas análises os resultados estão ok, conforme confronto de resultados da literatura. Porém quando insiro o comando de para análise de médias, os seguintes erros são apresentados: comando: TukeyHSD - TukeyHSD(aov(Ds ~ Bloco + Amb*Prof + Error(Bloco/Amb), data = ATRIB)) Erro: Erro em UseMethod("TukeyHSD") : no applicable method for 'TukeyHSD' applied to an object of class "c('aovlist', 'listof')" comando: LTukey (pacote Laercio) Erro em UseMethod("LTukey") : no applicable method for 'LTukey' applied to an object of class "c('aovlist', 'listof')" Como devo proceder para fazer a análise de médias? Obrigado Atenciosamente, Eng. Agrícola Rafael Cesar Tieppo UNEMAT - Campus Tangará da Serra - MT Departamento de Agronomia Fone:+55-65 3311-4922 site:http://www2.unemat.br/rafaeltieppo "Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente".

3 2

Re: [R-br] Teste média Subdividida
by Rafael Tieppo 22 Mar '11

22 Mar '11

Muito obrigado. -- View this message in context: http://r-br.2285057.n4.nabble.com/R-br-Teste-media-Subdividida-tp3396775p33… Sent from the R-br mailing list archive at Nabble.com.

1 0

Fwd: Popularity of R, SAS, SPSS, Stata, Statistica, S-PLUS updated
by Henrique Dallazuanna 22 Mar '11

22 Mar '11

Caros, Esta página mostra uma análise da popularidade de várias ferramentas estatísticas. Att ---------- Forwarded message ---------- From: Muenchen, Robert A (Bob) <muenchen(a)utk.edu> Date: Tue, Mar 22, 2011 at 7:37 AM Subject: Popularity of R, SAS, SPSS, Stata, Statistica, S-PLUS updated Greetings, I've just put out the latest version of "The Popularity of Data Analysis Software" at http://r4stats.com/popularity. This update includes complete data for 2010, the addition of number of blogs for each software, more coverage of Statistica, and, where possible, measures regarding the implementations of the SAS Language: Carolina and the World Programming System (WPS). Cheers, Bob ========================================================= Bob Muenchen (pronounced Min'-chen), Manager Research Computing Support Voice: (865) 974-5230 Email: muenchen(a)utk.edu Web: http://oit.utk.edu/research, News: http://oit.utk.edu/research/news.php ======================================================== ===================== To manage your subscription to SPSSX-L, send a message to LISTSERV(a)LISTSERV.UGA.EDU (not to SPSSX-L), with no body text except the command. To leave the list, send the command SIGNOFF SPSSX-L For a list of commands to manage subscriptions, send the command INFO REFCARD -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O

1 0

Re: [R-br] Como estimar os standard deviations de estimativas de herdabilidade em um teste
by Paulo Ricardo Gherardi Hein 22 Mar '11

22 Mar '11

Olá pessoal, Obrigado pelas dicas (função lmer do "lme4" e "SEI.GE"), mas o meu problema ainda não foi solucionado. Calcular o sd de herbabilidades estimadas a partir de modelos com efeito aleatório parece não ser uma tarefa simples no R. Os editores das revistas exigem o valor H² acompanhado pelo standard deviation/intervalo de confiança. Como é que os colegas que trabalham com melhoramento apresentam as suas estimativas de herdabilidade? Meus orientadores me ajudaram a calcular isso no SAS e no ASREML para outro design experimental (factorial mating design), mas eu não tenho licença para esses softwares e pretendo migrar 100% para o R até o ano que vem. Me desculpem por insistir no assunto, mas alguém da lista já calculou essas estatísticas no R? Abraços, Paulo Ricardo Gherardi Hein PhD candidate at University of Montpellier 2 CIRAD - PERSYST Department Research unit: Production and Processing of Tropical Woods TA B-40/16 bur.145 Adress: 73 rue Jean-François Breton 34398 Montpellier Cedex 5 - France phone: +33 4 67 61 44 51 / mobile: +33 6 47 90 57 60 email: paulo.hein(a)cirad.fr / skype: paulo_hein ------------------------------ *De:* paulo_hein <paulo_hein(a)yahoo.com.br> *Para:* R_STAT(a)yahoogrupos.com.br *Enviadas:* Terça-feira, 15 de Março de 2011 13:52:47 *Assunto:* [R_STAT] Re: Como estimar os standard deviations de estimativas de herdabilidade em um teste Prezado Fábio, obrigado pela sugestão. Prezado Fernando, você tem razão. Esse script retorna o valor do desvio padrão de cada componente da variância. No entanto, ao elevarmos os valores ao quadrado, a estimativa da herdabilidade é calculada apartir das variâncias, não é? O scrip retorna os valores de StdDev do efeito clone, interação clone x site e erro. Aqui vai um exmplo: 2581 > Resume1 2582 [[1]] 2583 Linear mixed-effects model fit by REML 2584 Data: Data 2585 Log-restricted-likelihood: -521.7632 2586 Fixed: Y ~ 1 + Site 2587 (Intercept) Site302 Site303 2588 60.79 -3.38 6.52 2590 Random effects: 2591 Formula: ~1 | Clone 2592 (Intercept) 2593 StdDev: *2.348988* 2595 Formula: ~1 | Site %in% Clone 2596 (Intercept) Residual 2597 StdDev: *5.427933* *7.052537* 2599 Number of Observations: 150 2600 Number of Groups: 2601 Clone Site %in% Clone 2602 10 30 Assim, estou calculando a broad-sense heritability da propriedade [[1]] como: H² = 2.348988^2 / (2.348988^2 + 5.427933^2 + 7.052537^2) O que eu preciso é calcular o standard deviation de cada estimativa de herdabilidade para apresentar os resultados como exigem as principais revistas da área. Algo do tipo: H² [[1]] = 0.065 (+ ou - SE) Meu orientador já calculou os SE usando o software ASREML, mas eu quero fazer isso no R. Abraços, Paulo Ricardo Gherardi Hein PhD candidate at University of Montpellier 2 CIRAD - PERSYST Department Research unit: Production and Processing of Tropical Woods TA B-40/16 bur.145 Adress: 73 rue Jean-François Breton 34398 Montpellier Cedex 5 - France phone: +33 4 67 61 44 51 / mobile: +33 6 47 90 57 60 email: paulo.hein@... / skype: paulo_hein > > > > > > ------------------------------ > > *De:* paulo_hein paulo_hein@... > > *Para:* R_STAT(a)yahoogrupos.com.br > > *Enviadas:* Segunda-feira, 14 de Março de 2011 5:24:20 > > *Assunto:* [R_STAT] Como estimar os standard deviations de estimativas de > > herdabilidade em um teste > > > > > > > > Olá, > > Trabalho com melhoramento de árvores (Eucalyptus) para producão de madeira. > > Atualmente estamos estimando os componentes de variância e herdabilidade de > > algumas propriedades da madeira proveniente de um teste clonal. O design > > experimental é composto por 10 clones, plantados em 3 sites. Cada clone é > > representado por 5 individuos, totalizando uma amostragem de 150 árvores (10 > > clones x 5 repetições x 3 sites). > > > > A maioria das estimativas do standard deviation de herdabilidades são > > feitas com base nas séries de Taylor, mas não sei como calculá-las. Alguém > > trabalha com esse assunto e conhece alguma solução? > > > > Para estimar os componentes de variância das 120 variáveis (estimados por > > espectroscopia de NIR) estamos usando dois modelos (com e sem interação). > > Para selecionar os modelos, estamos nos baseado no valor de BIC > > (selecionamos os menor BIC). Os cálculos são feitos a partir dos seguintes > > scripts: > > > > Data <- read.table("CNB150predNIR.csv", dec=".", sep=",", header=T) > > Data$Site <- as.factor(Data$Site) > > Data$Clone <- as.factor(Data$Clone) > > > > ## Estimating variance components > > Resume <- NULL > > Resume1 <- NULL > > for (i in 4:124){ > > Y <- Data[,i] > > lme <- lme(Y~ 1 + Site, random=~1|Clone, data=Data, method="REML") > > lme1 <- lme(Y~ 1 + Site, random=~1|Clone/Site, data=Data, method="REML") > > Resume <- c(Resume,list(lme)) > > Resume1 <- c(Resume1,list(lme1))} > > Resume > > Resume1 > > > > ## For calculating BIC > > Selection <- NULL > > Selection1 <- NULL > > for (i in 4:124){ > > Y <- Data[,i] > > lme <- lme(Y~ 1 + Site, random=~1|Clone, data=Data, method="ML") > > lme1 <- lme(Y~ 1 + Site, random=~1|Clone/Site, data=Data,method="ML") > > Selection <- c(Selection,list(lme)) > > Selection1 <- c(Selection1,list(lme1))} > > > > ## To display BIC results > > Resultat <- NULL > > for (i in 1:121){ lme0 <- Selection[[i]] > > lme1 <- Selection1[[i]] > > Resultat <- c(Resultat,list(anova(lme0,lme1)))} > > Resultat > > > > e as herdabilidades são estimadas da seguinte forma: > > > > H² (mle) = sdev(Clone)^2 / [sdev(Clone)^2 + sdev(Error)^2] ou > > H² (lme1) = sdev(Clone)^2 / [sdev(Clone)^2 + sdev(interaction)^2 + > > sdev(Error)^2] > > > > A maioria das estimativas do standard deviations de herdabilidades são > > feitas com base nas séries de Taylor, mas não sei como calculá-las. Alguém > > trabalha com esse assunto e conhece alguma solução? > > > > Atenciosamente, > > > > Paulo Ricardo Gherardi Hein > > PhD candidate at University of Montpellier 2 > > CIRAD - PERSYST Department > > Research unit: Production and Processing of Tropical Woods TA B-40/16 > > bur.145 > > Adress: 73 rue Jean-François Breton 34398 Montpellier Cedex 5 - France > > phone: +33 4 67 61 44 51 / mobile: +33 6 47 90 57 60 > > email: paulo.hein@... / skype: paulo_hein > > > > > > > > > > > __._,_.___ <ivanalaman@yahoo.com.br?subject=Res%3A%20%5BR_STAT%5D%20Re%3A%20Como%20estimar%20os%20standard%20deviations%20de%20estimativas%20de%20herdabilidade%20em%20um%20teste>| através de email<R_STAT@yahoogrupos.com.br?subject=Res%3A%20%5BR_STAT%5D%20Re%3A%20Como%20estimar%20os%20standard%20deviations%20de%20estimativas%20de%20herdabilidade%20em%20um%20teste>| Responder através da web<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJycDAydGVhBF9TAzk3N…>| Adicionar um novo tópico<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmZGxhb2hzBF9TAzk3N…> Mensagens neste tópico<http://br.groups.yahoo.com/group/R_STAT/message/15816;_ylc=X3oDMTM3ajNyaGtk…>( 5) Atividade nos últimos dias: - Novos usuários<http://br.groups.yahoo.com/group/R_STAT/members;_ylc=X3oDMTJnZ3A1MjJ1BF9TAz…> 9 Visite seu Grupo<http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmZGEzb2RnBF9TAzk3NDkwND…> Você pode ser quem está faltando para alguém.<http://global.ard.yahoo.com/SIG=15ma38fhn/M=758712.14532720.14422460.129601…> [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlbTY0dXIzBF9TAzk3NDkwNDM1BGdycElkAz…> Trocar para: SÃ³ Texto<R_STAT-traditional@yahoogrupos.com.br?subject=Mudar+Formato+de+Envio:+Tradicional>, Resenha DiÃ¡ria<R_STAT-digest@yahoogrupos.com.br?subject=Envio+de+email:+Resenha>• Sair do grupo <R_STAT-unsubscribe(a)yahoogrupos.com.br?subject=Sair+do+grupo> • Termos de uso <http://br.yahoo.com/info/utos.html> . __,_._,___

3 2

Criar um ambiente LaTeX que possui lista no sumário
by Walmes Zeviani 22 Mar '11

22 Mar '11

Saudações, Fui questionado por um aluno se é possível criar um novo ambiente no LaTeX para inserir fotos e ter a lista de fotografias assim como temos a lista de figuras e lista de tabelas. Alguém pode citar algum pacote? Desde já grato e desculpem pela dúvida fora do escopo da [R-br]. Walmes. -- ========================================================================== 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

4 3

Enc: [r-br] Função para quebra de linha em título
by Fabio Mathias Corrêa 21 Mar '11

21 Mar '11

----- Mensagem encaminhada ---- De: Henrique Dallazuanna <wwwhsd(a)gmail.com> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Sábado, 25 de Outubro de 2008 20:54:23 Assunto: Re: [R_STAT] Função para quebra de linha em título Marcio, uma alternativa talvez seja: nm <- c("Aaaaaa aaaaaaa aa", "Bbbbbb bbbbbb", "cccc") y <- c(10,4,9) barplot(y, names.arg = nm) strP <- .5 nm[strwidth( nm) > strP] <- gsub("( \\w*$)", "\n\\1", nm[strwidth( nm) > strP]) 2008/10/25 Marcio Nicolau <marcio.nicolau@ gmail.com> Prezados, > >A algum tempo, questionei na lista sobre a existência de uma função >para quebra de linha >automática para títulos de gráficos. > >Após a sugestão do Benilton (função gsub), segue abaixo o código... >caso necessitem. > > >## Pre-Processa titulos > >pretitulo <- function(orig, cut=55) { > nparts = ceiling(nchar( orig) / cut) > cuts = c(1,1:nparts * cut) > spt = strsplit(orig,"\\s")[[1]] > cs = cumsum(nchar( spt)) > > ids <- matrix(FALSE, nparts,length( nchar(spt) )) > for(i in 2:length(cuts) ) { > ids[i-1,] <- (cs >= cuts[i-1] & cs < cuts[i]) > } > > ch <- vector("character", nparts) > for(i in 1:nparts) { > ch[[i]] <- paste(spt[ids[ i,]],collapse = " ") > } > pretitulo <- paste(ch,collapse="\n") >} > >Não é o melhor, computacionalmente falando, mas funciona :) > >Abraços > >Marcio Nicolau >marcio.nicolau@ gmail.com > > > > >------------ --------- --------- ------ > >Links do Yahoo! Grupos > ><*> Para visitar o site do seu grupo na web, acesse: > http://br.groups. yahoo.com/ group/R_STAT/ > ><*> Para sair deste grupo, envie um e-mail para: > R_STAT-unsubscribe@ yahoogrupos. com.br > ><*> O uso que você faz do Yahoo! Grupos está sujeito aos: > http://br.yahoo. com/info/ utos.html > > > -- Henrique Dallazuanna Curitiba-Paraná -Brasil 25° 25' 40" S 49° 16' 22" O __._,_.___ Mensagens neste tópico (2) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 8 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

4 3

Função para quebra de linha em título
by Fabio Mathias Corrêa 21 Mar '11

21 Mar '11

Fábio Mathias Corrêa Departamento de Estatística Universidade Estadual de Santa Cruz Tel.: 73-3680-5076 Cel.: 73-9991-8155 ----- Mensagem encaminhada ---- De: Marcio Nicolau <marcio.nicolau(a)gmail.com> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Sábado, 25 de Outubro de 2008 19:47:08 Assunto: [R_STAT] Função para quebra de linha em título Prezados, A algum tempo, questionei na lista sobre a existência de uma função para quebra de linha automática para títulos de gráficos. Após a sugestão do Benilton (função gsub), segue abaixo o código... caso necessitem. ## Pre-Processa titulos pretitulo <- function(orig, cut=55) { nparts = ceiling(nchar( orig) / cut) cuts = c(1,1:nparts * cut) spt = strsplit(orig, "\\s")[[1] ] cs = cumsum(nchar( spt)) ids <- matrix(FALSE, nparts,length( nchar(spt) )) for(i in 2:length(cuts) ) { ids[i-1,] <- (cs >= cuts[i-1] & cs < cuts[i]) } ch <- vector("character" , nparts) for(i in 1:nparts) { ch[[i]] <- paste(spt[ids[ i,]],collapse = " ") } pretitulo <- paste(ch,collapse= "\n") } Não é o melhor, computacionalmente falando, mas funciona :) Abraços Marcio Nicolau marcio.nicolau@ gmail.com __._,_.___ Mensagens neste tópico (1) Responder (através da web) | Adicionar um novo tópico Mensagens | Arquivos | Fotos | Links | Banco de dados | Enquetes Alterar configurações via web (Requer Yahoo! ID) Alterar configurações via e-mail: Alterar recebimento para lista diária de mensagens | Alterar formato para o tradicional Visite seu Grupo | Termos de uso do Yahoo! Grupos | Sair do grupo Atividade nos últimos dias * 8 Novos usuáriosVisite seu Grupo Yahoo! Mail Conecte-se ao mundo Proteção anti-spam Muito mais espaço Yahoo! Barra Instale grátis Buscar sites na web Checar seus e-mails . Yahoo! Grupos Crie seu próprio grupo A melhor forma de comunicação . __,_._,___

6 6

Transformando dados com 2 fatores em uma tabela
by Augusto Ribas 21 Mar '11

21 Mar '11

Oi galera. Minha primeira postagem na nova lista que emoção :) Gostaria de ajuda com um problema simples que estou tendo aqui. #Eu queria fazer um plot assim: barplot(VADeaths, beside = TRUE, col = c("lightblue", "mistyrose", "lightcyan", "lavender", "cornsilk"), legend = rownames(VADeaths), ylim = c(0, 100)) title(main = "Death Rates in Virginia", font.main = 4) #Mas meu dados estão na forma de um data frame dessa forma... exemplo<-dados<-data.frame(cbind( Tratamento=rep(c("Controle","Ac10","Ac20","Estresse"),each=3) ,Placa=rep(1:3,4) ,Germinacao=round(rnorm(12,5,3),2) )) #Tem alguma maneira simples de fazer meu dataframe de exemplo virar uma tabela como a do "VADeaths" que é usado no help do barplot #pra fazer um plot daquele jeito, separando grupos de barrinhas tanto por tratamento como placa, no #meu caso, ou então fazer um barplot do jeito que quero(o do "VADeaths") direto com um dataframe assim do meu jeito? Obrigado -- Grato Augusto C. A. Ribas

2 2

Res: [R_STAT] Regressão de proporções com erro normal ou binomial
by Walmes Zeviani 21 Mar '11

21 Mar '11

Guilherme, Na regressão logística você observa significância porque o parâmetro de dispersão não é estimado mas sim conhecido/fixado como 1 (note que é aplicado o teste z). No modelo normal, além da estimação dos efeitos/coeficientes/parâmetros, é estimado a variância do erro aleatório e os erros padrões das estimativas estão em função disso (teste t). Se você usar uma family=quasibinomial você notará diferença. Esperamos você na [R-br]. À disposição. Walmes. -- ========================================================================== 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

Estudo de função!
by Elton Gean Araújo Araújo 21 Mar '11

21 Mar '11

Boa tarde! Gostaria de aplicar os conceitos de crescimento, decrescimento, máximos e mínimos de funções no R, poderiam me informar o pacote para isso? Desde já agradeço pela informação e pela nova lista! Abraço. Elton.

1 0

Especificação do modelo em lm e glm
by Diego Pujoni 21 Mar '11

21 Mar '11

#Olá pessoal, #Existem formas diferentes de se especificar um modelo na função #"lm" e "glm". Mas obtive os seguintes resultados: #Criando o data frame# X1=c(rep("A",4),rep("B",4)) X2=c(rep("C",2),rep("D",2),rep("C",2),rep("D",2)) Y=c(3,5,5,7,11,13,13,15) Table=as.data.frame(cbind(X1,X2,as.numeric(Y))) Table #Os quatro modos de se especificar o mesmo modelo fornecem a mesma deviance,# #ou seja, são iguais# glm(Y ~ X1+X2+X1:X2, data=Table, family=poisson())$deviance glm(Y ~ (X1+X2)^2, data=Table, family=poisson())$deviance glm(Y ~ X1*X2, data=Table, family=poisson())$deviance glm(Y ~ X1:X2, data=Table, family=poisson())$deviance #Mas há um deslocamento de -1.25 nos coeficientes no último modelo round(glm(Y ~ X1+X2+X1:X2, data=Table,family=poisson())$coefficients,digits=2) round(glm(Y ~ (X1+X2)^2, data=Table,family=poisson())$coefficients,digits=2) round(glm(Y ~ X1*X2, data=Table,family=poisson())$coefficients,digits=2) round(glm(Y ~ X1:X2, data=Table,family=poisson())$coefficients,digits=2) #O coeficiente de X1B:X2D sai como NA# #O quarto modo de especificar o modelo está errado? #Por que isso acontece? #Nos modelos abaixo: #Incluindo mas uma variável preditora em Table# X3=c("E","F","E","F","E","F","E","F") Table.X3=as.data.frame(cbind(X1,X2,X3,as.numeric(Y))) Table.X3 #As deviances e os coeficientes X1A:X3F e X1A:X3E também são iguais #além disso alguns coeficientes são deslocados de -1.4 (X1A:X2C,X1B:X2C,X1A:X1D), #mas os outros coeficientes não são iguais# glm(Y~X1+X2+X3+X1:X2+X1:X3,data=Table.X3,family=poisson())$deviance glm(Y~X1:X2+X1:X3,data=Table.X3,family=poisson())$deviance round(glm(Y~X1+X2+X3+X1:X2+X1:X3,data=Table.X3,family=poisson())$coefficients,digits=2) round(glm(Y~X1:X2+X1:X3,data=Table.X3,family=poisson())$coefficients,digits=2) #Especificando os coeficientes do primeiro modelo # Intercept = 1.16 # X1/X2| C D | X1/X3| E F | # A | 0 0.41| A | 0 0.41| # B |1.24 -0.25| B |1.24 -0.25| #Especificando os coeficientes do segundo modelo # Intercept = 2.56 # X1/X2| C D | X1/X3| E F | # A |-1.40 -0.99| A | 0 0.41| # B |-0.15 NA | B | 0 0.15| #Os coeficientes X1B:X3E e X1B:X3F mudam. #Alguém sabe aonde se encontra o meu erro? #Será que sou obrigado a incluir o termo de interação X2:X3? #Mas e se essas duas variáveis forem independentes? #Muito Obrigado # Diego PJ

2 1

Solução de Equação não Linear
by Carolina Pena 21 Mar '11

21 Mar '11

Boa tarde! Alguém sabe como consigo resolver uma equação não linear no R ? Por exemplo suponha que eu tenho a seguinte equação: y = 2,5*x + 3,6*x^2 + 7,8*x^4+ 10,0*x^12 Para um dado valor de y (exemplo y = 5 ) qual é o valor de x ? Obrigada, Carolina

4 4

Re: [R-br] [R_STAT] [Dúvida] Carregar dados sem ser na memória RAM
by Fernando Mayer 20 Mar '11

20 Mar '11

O pacote SOAR pode resolver. De uma olhada no manual e na vignette em http://cran-r.c3sl.ufpr.br/web/packages/SOAR/index.html []s, --- Fernando Mayer URL: http://sites.google.com/site/fernandomayer e-mail: fernandomayer [@] gmail.com 2011/3/20 pedrorafaelmarinho <pedrorafaelmarinho(a)yahoo.com.br>: > Pessoal estou com um problema ou talvez seja do R. Tenho que ler uma base de dados de 600 mb em formato DBF (DBase). O R pelo que ando lendo ele ler a base de dados e joga na memória RAM. Toda vez que carrego essa base de dados o pc trava a ponto de o mouse não funcionar. Existe alguma forma de fazer com que o R não carregue toda a base de dados diretamente na memória? Por exemplo a base de dados ficar no HD e somente carregar algumas variáveis a medida que vou precisando? Softwares simples como o Tabwin eu consigo acessar a base de dados e o pc não trava...O que está acontecendo? > O comando que usei para carregar a base foi read.dbf() usando o pacote foreign. > > > > ------------------------------------ > > Links do Yahoo! Grupos > > <*> Para visitar o site do seu grupo na web, acesse: > http://br.groups.yahoo.com/group/R_STAT/ > > <*> Para sair deste grupo, envie um e-mail para: > R_STAT-unsubscribe(a)yahoogrupos.com.br > > <*> O uso que você faz do Yahoo! Grupos está sujeito aos: > http://br.yahoo.com/info/utos.html > > >

3 4