- R-br - DINF UFPR

Informação arquivo xls - Detalhes
by Eder David Borges da Silva 24 Abr '12

24 Abr '12

Boa noite pessoal, Venho mais uma vez recorrer as suas sabias ajudas, em: http://www.leg.ufpr.br/~eder/teste.xls deixei um arquivo xls, ###Caso queiram carregar ele: download.file('http://www.leg.ufpr.br/~eder/teste.xls','teste.xls') file.info('teste.xls') # Minha duvida é a seguinte, se clicarmos com o direito e ir nas abas: # detalhes -> descrição # Veremos que tem a opção de ter um titulo, assunto, marca etc. # minha missão é modificar este titulo (esta como MEU TITULO) para outro nome. # via comando... # Alguma sugestão? # condições de contorno: SO Windows, R 2.15. Obrigado pessoal

2 2

instalação de pacotes
by André Figueiras Rutz 24 Abr '12

24 Abr '12

Caros, Estou tendo dificuldades com o R na instalação de pacotes alguem poderia me ajudar? Parece que o pacote vai para o micro mas de alguma forma na hora de carregar nao funciona sugestões? Grato André Rutz > install.packages("rJava") --- Please select a CRAN mirror for use in this session --- tentando a URL ' http://cran-r.c3sl.ufpr.br/bin/windows/contrib/2.15/rJava_0.9-3.zip' Content type 'application/zip' length 746108 bytes (728 Kb) URL aberta downloaded 728 Kb package ‘rJava’ successfully unpacked and MD5 sums checked The downloaded binary packages are in C:\Users\user\AppData\Local\Temp\RtmpoBqDjS\downloaded_packages > require(rJava) Carregando pacotes exigidos: rJava Error : .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: fun(libname, pkgname) erro: JAVA_HOME cannot be determined from the Registry > require("rJava") Carregando pacotes exigidos: rJava Error : .onLoad falhou em loadNamespace() para 'rJava', detalhes: chamada: fun(libname, pkgname) erro: JAVA_HOME cannot be determined from the Registry

3 4

Como saber se é um 'bug' ou uma 'feature', no comportamento de 'plot'?
by Cleber N.Borges 24 Abr '12

24 Abr '12

olá a todos da lista, estava fazendo uns testes e me deparei com um comportamento do qual *acho* que pode ser incoerente. Vamos ao problema: > par(mfrow=c(2,2)) > plot(1:10) # grafico 01 > plot(1:10, cex.axis=0.5, las=1) # grafico 02 > plot(table(1:10)) # grafico 03 > plot(table(1:10), cex.axis=0.5, las=1) # grafico 04 a mudança dos parametros: 'cex.axis' e 'las' funcionam muito bem na mudança do gráfico 1 para 2, mas não há efeito na mudança de 3 para 4. Segunda pergunta: Como saber qual é a verdadeira função 'plot.ALGUMACOISAAQUI' que será chamada para fazer o grafico de um objeto da classe 'table'? dado que não há 'plot.table' como função... Grato pela atenção de todos. Cleber

2 3

Dúvida sobre ANOVA tipo III
by Olympio Neto 24 Abr '12

24 Abr '12

Pessoal, Tenho lido vários livros sobre R, muitos blogs (R-bloggers) e as excelentes páginas do Wiki-Leg da UFPR e as famosas Ridículas do Walmes e em muitas delas vejo que em muitos dos tutoriais, é necessária, para a obtenção da Anova Marginal (Tipo III), setar o contraste como "contr.sum" e outros que nada falam sobre isso, ficando setado o padrão do R, o "contr. treatment". A explicação para a alteração do contraste é que a soma deste deve ser nula, o que não acontece com o "contr.treatment". Caso não seja feita a alteração, os resultados podem estar errados. O artigo que tenta por luz neste assunto é aqui indicado: http://goanna.cs.rmit.edu.au/~fscholer/anova.php Normalmente eu uso o pacote "car" para cálculos que análises de variância do tipo III, mas nos exemplos que vejo, dados que simulo etc, não vejo diferenças. Acredito que nem poderiam haver, mesmo. A ANOVA trivial (tipo I) é invariante à parametrização, mas fico muito em dúvida se a tipo III também seria. O que vocês pensam? Já observaram diferenças importantes em suas análises? Abs e muito obrigado. Olympio

4 10

Paste
by Ricardo Fernandes Campos Junior 24 Abr '12

24 Abr '12

Olá, pessoal. Eu sei que essa pergunta é provavelmente muito simples, mas eu não consigo encontrar como fazê-lo. Eu tenho um objeto de caracteres de tamanho genérico, para fim de exemplo, vou usar 3. Eu quero criar um outro objeto com estes mesmos caracteres, mas concatenados. Porém não consigo. Tentei isso: exemplo = c("A", "B", "C") novoobjeto = paste(exemplo, " ") O que me é retornado é isso: "A" "B" "C" Mas o que eu quero é isso: "A B C" Alguém pode ajudar? -- Abraços, Ricardo Fernandes Campos Junior Departamento de Genética e Biologia Evolutiva Universidade de São Paulo

3 3

Notação
by Edson Lira 24 Abr '12

24 Abr '12

Tenho esse script abaixo: require(PASWR) trat<-rep(LETTERS[1:4],5) med<-c(25,31,22,33,26,25,26,29,20,28,28,31,23,27,25,34,21,24,29,28) exem1_4<-data.frame(med,trat) exem1_4 m1<-lm(exem1_4$med~exem1_4$trat) anova(m1) summary(m1) erro<-resid(m1) erro est<-predict(m1) med-est Quando gero o erro ele me dá essa saída. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 2.00e+00 4.00e+00 -4.00e+00 2.00e+00 3.00e+00 -2.00e+00 -3.33e-16 -2.00e+00 -3.00e+00 1.00e+00 2.00e+00 -1.11e-16 8.88e-16 1.11e-16 -1.00e+00 3.00e+00 -2.00e+00 18 19 20 -3.00e+00 3.00e+00 -3.00e+00 Gostaria de transformar estes números para a notação comum, tipo 2, 3, -4 .... Me lembro que certa vez alguém postou esse problema, mas faz muito tempo. Alguém tem alguma sugestão? Edson Lira Estatístico Manaus-Amazonas

3 2

[Dúvida] Números aleatórios.
by Pedro Rafael 23 Abr '12

23 Abr '12

Formula em LaTeX: F_{T}(t;\alpha,\beta) = \Phi(\frac{1}{\alpha}(\frac{t}{\beta}^(1/2) - \frac{\beta}{t}^(1/2))) Pessoal com gerar números aleatórios com distribuição Birbaum-Saunder em C? Tenho a distribuição acumulada da Birbaum-Saunder que é escrita como acumulada de uma normal e tenho a densidade da Birbaum. A formula da acumulada está escrita acima para t e beta maiores que zero. Pensei primeiramente em pegar números aleatórios normais pela biblioteca GSL e fazer ser igual à t na formula acumulada acima. Daria certo? Se for pela transformada inversa eu teria que gerar números uniformes[0,1] e fazer F(t) = U. Depois eu tinha que ver qual o t em F me daria U, sendo U o número uniforme gerado. Outra complicação é que tenho uma acumulada da normal que é uma integral de -inf a t. Existe outra forma mais fácil? -- Saudações, Pedro Rafael Diniz Marinho.

2 1

Ajuda para evitar um loop com elegância 2
by Gustavo Marcatti 23 Abr '12

23 Abr '12

Olá pessoal, tenho que executar uma função parecida com essa no R. Mas ela normalmente demora longos 500 segundos para executar. Fiz uma versão que chama um código C, e essa é bem eficiente (demora cerca de 0,5 segundos), porém gostaria de uma função totalmente em R. Caso alguém queira ajudar, desde já agradeço. set.seed(1) vetor1 <- rep(1:4, each = 5) var1 <- abs(rnorm(4,3)) * vetor1 cond1 <-abs(rnorm(4,10)) * vetor1 var2 <- matrix(0, nrow=length(vetor1)) cond2 <- matrix(0, nrow=length(vetor1)) for(i in 1:(length(vetor1) - 1)){ if(vetor1[i] == vetor1[i+1]){ var2[i] = var1[i+1] cond2[i] = cond1[i+1] }else { var2[i] = NA cond2[i] = NA } } var2[length(vetor1)] <- NA cond2[length(vetor1)] <- NA dados <- na.exclude(data.frame(var2, cond2)) Gustavo Marcatti Eng. Florestal

2 2

R servidor Windows?
by Vinicius Brito Rocha 23 Abr '12

23 Abr '12

Senhores, é possível instalar o R em um servidor Windows para a utilização de multiplus usuários? Att. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 2

bootstrap
by Maria Papa 23 Abr '12

23 Abr '12

Olá Pessoal, boa noite Estou tentando, sem sucesso, realizar uma reamostragem bootstrap para médias de 3 amostras estão organizados em uma matriz na qual, que sendo cada coluna representa uma amostra. Porém, a cada reamostragem, eu preciso que seja feita para as três amostras, ou seja, quando um dado for selecionado para formar a amostra bootstrap, deve ser retirada a linha toda. O resultado final deve ser 3 médias reamostradas, que foram obtidas em um único processo de reamostragem. Alguém sabe como fazer isso no R? Eu consegui somente quando eu realizo 3 procedimentos de reamostragem distintos. Desde já agradeço a todos célia

3 2

Dúvida de Geoestatística
by Mauro Sznelwar 21 Abr '12

21 Abr '12

Queria saber se a bibliotéca GeoR faz kriging external drift (KED), e se tem o recurso para análise multivariada, onde eu tenho mais de uma variável como variável resposta. Outra coisa, se ela é tridimensional, pois eu tenho as dita cujas em várias camadas, é uma exploração arqueológica. Caso não tenha estes recursos, quais bibliotécas as tem?

2 3

Problema com pacote parallel
by Vanessa Leotti Torman 20 Abr '12

20 Abr '12

Prezados membros da lista, Bom dia. Estou rodando uma simulação bastante complexa, pois para cada replicação tenho que rodar uma cadeia mcmc através do pacote brugs (que necessita do software OpenBUGS). Para que a execução seja mais rápida, estou usando o pacote paralell, mais especificamente a função parCapply deste pacote. Os testes com poucas replicações funcionam perfeitamente. Mas quando coloco para rodar todas (4000), mesmo após várias horas executando, o Windows dá uma mensagem que o R front end parou de funcionar. Eu percebo através do gerenciador de tarefas que apenas um dos núcleos parou de executar, então clico em fechar e acabo perdendo todos os resultados. Já tentei não usar todos os núcleos do processador, mas mesmo assim o problema acontece. Alguém teria alguma sugestão? Abraços, Vanessa

3 5

Dependencias no mclapply
by Junior Beleti 20 Abr '12

20 Abr '12

Olá a todos. Na busca de melhorias no processamento paralelo, estou com uma dúvida quanto a função mclapply do pacote multicore. Benilton já me deu várias dicas de como proceder com o mesmo, mas ainda tenho a seguinte dúvida: - Quanto a questão de dependências dentro das operações do mclapply, por exemplo, nas seguintes instruções: mclapply(1:100, function(i) { g <- mvrnorm(n=1, mat.cokri[[1]],mat.cokri[[2]]) #1 seq1 <- seq(1,nlinhas,by=2) #2 seq2 <- seq(2,nlinhas,by=2) #3 y1 <- g[seq1] #4 y2 <- g[seq2] #5 gerado <- data.frame(y1,y2) #6 compos <- agl(gerado) #7 compos1 <- cbind(compos,compos1) #8 },mc.cores=4 ) - Em minha análise, cada um dos 4 núcleos vai executar todas as instruções dentro da função, ou seja, cada núcleo vai executar as 8 instruções dentro da função. Mas a dúvida surge quanto as dependências, visto que em uma execução sequencial, todas as instruções seriam executadas sequencialmente, e a próxima iteração só iniciaria após o resultado de "compos1" ter sido atribuído. E quanto ao processo paralelo, seria respeitado tal sequencia de instruções, ou seja, a segunda iteração inicia somente após "compos1" ter sido atribuído na primeira iteração? Como o mclapply possui tal controle? Estou ciente que esta questão pode ser mais da implementação do próprio pacote, mas se alguém puder me esclarecer ficaria muito grato. Att, Carlos Beleti.

6 13

Dúvida conectar R com Interbase/firebird
by Marco Pollo 20 Abr '12

20 Abr '12

Caros, Registro os agradecimentos do email anterior, mas, alguém saberia qual seria configuração do ODBC (na plataforma Windows 32bits) para essa base? E como ficaria o comando odbcConnect, considerando minhas tentativas (infrutíferas) registradas abaixo: Configuração do meu ODBC INTERBASE/FIREBIRD : TENTATIVAS NO R: 1ª tentativa: channel <- odbcConnect(SISATER, uid = "SYSDBA", pwd = "masterkey") Erro em paste("DSN=", dsn, sep = "") : objeto 'SISATER' não encontrado 2ª tentativa: channel <- odbcConnect(SISATER.GDB, uid = "SYSDBA", pwd = "masterkey") Erro em paste("DSN=", dsn, sep = "") : objeto 'SISATER.GDB' não encontrado 3ª tentativa: channel <- odbcConnect(C:\SISATER\SISATER.GDB, uid = "SYSDBA", pwd = "masterkey") Erro: entrada inesperado em "channel <- odbcConnect(C:\" Obrigado Marco Pollo

2 1

Re: [R-br] Dependencias no mclapply
by Junior Beleti 20 Abr '12

20 Abr '12

Okay Benilton, entendido.. Obrigado, Carlos Beleti.

1 0

STARMA
by Gilbert Queiroz 20 Abr '12

20 Abr '12

Alguém pode me ajudar a encontrar alguma coisa referente a isto???? Mesmo que ainda não tenha código p/ o R. Abs. Gilbert

3 3

(sem assunto)
by Paulo Ricardo Silva Moreira 19 Abr '12

19 Abr '12

Boa tarde, Estou precisando usar a função tapply. Alguém sabe como utilizá-la? Att, Paulo Ricardo.

3 2

Soma de quadrados de regressão..
by Julia Maria 19 Abr '12

19 Abr '12

Pessoal, como sempre estou com algumas dúvidas básicas. Eu ainda estou na graduação e estou me adaptando ao R. Será que alguem poderia me tirar uma dúvida? Quando faço uma analise de um modelo de regressão no R, por exemplo, do modelo: m1<-lm(Y~-1+x1+x2+x3+I(x1^2) +I(x2^2)+I(x3^2) (as variáveis estao definidas corretamente em um arquivo a parte) e peço anova ele solta as SQs de cada coeficiente. Gostaria de saber se há alguma forma de obter um quadro de análise de variância da regressão, com a soma de quadrados da regressao linear e da regressao quadratica, ao inves de cada coeficiente.

3 2

Fwd: [paleontropica] Macroevolution in R Short Course: Call for Applicants
by D'Apolito 19 Abr '12

19 Abr '12

We are pleased to announce an intensive short course on using R to perform comparative methods to be held in Santa Barbara on June 11th to June 15th. This course is funded by the National Science Foundation, and a number of stipends to cover or defray travel, room, and board are available to qualified students and post-docs. Topics covered will include an introduction to the R programming language, tree manipulation, independent contrasts and phylogenetic generalized least squares, ancestral state reconstruction, models of character evolution, diversification analyses, and community phylogenetic analysis. If you are interested please submit your CV along with a short (maximum 1 page) description of your research interests, background, and reasons for taking the course. Admission is competitive, and the best applications come from students with data sets to analyze. International applicants are welcome. To apply visit this URL: tinyurl.com/macro-in-R Application deadline: April 30th. Individuals from cultural, racial, linguistic, geographic and socioeconomic backgrounds that are currently underrepresented in science are especially encouraged to apply. Please contact the co-organizers, Michael Alfaro (michaelalfaro(a)ucla.edu) and Luke Harmon (lukeh(a)uidaho.edu) with any questions. -- David Bapst Dept of Geophysical Sciences University of Chicago 5734 S. Ellis Chicago, IL 60637 http://home.uchicago.edu/~dwbapst/ http://cran.r-project.org/web/packages/paleotree/index.html _______________________________________________ Paleonet mailing list Paleonet(a)nhm.ac.uk http://mailman.nhm.ac.uk/mailman/listinfo/paleonet __._,_.___ Responder a remitente<alc1006@hotmail.com?subject=RE%3A%20Macroevolution%20in%20R%20Short%20Course%3A%20Call%20for%20Applicants>| Responder a grupo<paleontropica@gruposyahoo.com?subject=RE%3A%20Macroevolution%20in%20R%20Short%20Course%3A%20Call%20for%20Applicants>| Responder mediante la Web<http://espanol.groups.yahoo.com/group/paleontropica/post;_ylc=X3oDMTJxbmVxa…>| Crear un tema nuevo<http://espanol.groups.yahoo.com/group/paleontropica/post;_ylc=X3oDMTJmNWo4a…> Mensajes con este tema<http://espanol.groups.yahoo.com/group/paleontropica/message/1633;_ylc=X3oDM…>( 1) Actividad reciente: Visita tu grupo<http://espanol.groups.yahoo.com/group/paleontropica;_ylc=X3oDMTJmMmg2N3E1BF…> [image: Yahoo! Grupos]<http://espanol.groups.yahoo.com/;_ylc=X3oDMTJldGQxdTA0BF9TAzk3NDkwNDU3BGdyc…> Cambiar a: Solo Texto<paleontropica-traditional@gruposyahoo.com?subject=Change+Delivery+Format:+Traditional>, ReseÃ±a Diaria<paleontropica-digest@gruposyahoo.com?subject=Email+Delivery:+Digest>• Cancelar suscripción<paleontropica-unsubscribe(a)gruposyahoo.xn--com?subject=cancelar+suscripcin-99c>• Condiciones de uso <http://e1.docs.yahoo.com/info/utos.html> . __,_._,___ -- _____________________________________ *Carlos D'Apolito Júnior - Intern* Center for Tropical Paleoecology and Archeology-CTPA Smithsonian Tropical Research Institute-STRI Panama City (507) 212-8019 CV: http://lattes.cnpq.br/4557754746424026

1 0

duvida comando read.fwf
by Alexandre Heinemann 18 Abr '12

18 Abr '12

Caros Colegas, Gostaria de saber se alguem poderia me ajudar a entender por que o comando abaixo não funciona? Abaixo segue exemplo com 1 arquivo, mas tenho mais de 100 arquivos para ler. http://www.datafilehost.com/download-7e562cad.html (link para download do arquivo utilizado como exemplo (somente 1 arquivo)) #codigo nasaAirTempFiles <- file.path("NASA", "Files", paste("temperature", 1:1, ".txt",sep="")) head(nasaAirTempFiles) [1] "NASA/Files/temperature1.txt" tempDF <- read.fwf(nasaAirTempFiles[1],skip=7,widths=c(-12, rep(7, 24))) #esse comando não funciona??? Error in file(file, "rt") : cannot open the connection In addition: Warning message: In file(file, "rt") : cannot open file 'NASA/Files/temperature1.txt': No such file or directory tempDF <- read.fwf("cloudhigh1.txt",skip=7,widths=c(-12, rep(7, 24))) # qdo substituo o nome do arquivo "nasaAirTempFiles[1]" para "cloudhigh1.txt" o nome original do arquivo (link para download) funciona head(tempDF) V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18 V19 V20 V21 V22 V23 V24 1 26.0 23 23.0 17.0 19.5 17.0 16.0 16.0 16.0 19.0 18.0 19.0 19.5 19.5 18.0 16.0 14.0 14.0 14.0 18.0 20.5 22.0 25.5 25.5 2 20.0 20 18.5 16.5 18.0 15.0 15.0 15.0 16.0 15.0 15.0 17.5 17.0 17.0 17.0 17.5 17.5 18.0 19.0 19.0 20.0 22.0 23.5 23.5 3 16.0 16 14.0 12.5 13.5 14.5 13.5 13.5 13.5 13.0 13.0 16.5 18.5 20.0 21.0 21.0 21.0 20.5 20.5 19.5 19.0 19.0 19.0 20.5 4 13.0 11 11.0 11.0 14.0 13.5 11.5 12.0 15.0 18.0 18.0 18.5 20.5 20.5 18.5 18.0 17.0 15.0 15.0 14.0 13.0 13.0 13.5 12.5 5 7.5 8 10.5 10.5 15.0 16.5 19.5 20.0 20.5 19.5 17.0 16.0 16.0 16.0 12.5 9.5 8.5 7.5 5.5 4.5 4.5 6.5 6.5 10.0 6 8.0 11 13.0 19.5 19.5 26.0 27.5 25.5 20.0 14.0 9.5 8.5 10.0 8.0 4.5 5.0 5.0 4.5 3.5 2.5 2.5 4.0 8.0 12.5 Att., Alexandre

4 4

Dúvida acessar base de dados
by Marco Pollo 18 Abr '12

18 Abr '12

Caros, Algum direcionamento para eu conseguir acessar um arquivo com extensão .GDB (FireBird Database) direto de uma base de dados. Marco Pollo Almeida

2 1

Re: [R-br] Duvida String
by Sérgio Henrique almeida da silva ju 18 Abr '12

18 Abr '12

---------- Mensagem encaminhada ---------- De: Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com> Data: 18 de abril de 2012 11:22 Assunto: Re: [R-br] Duvida String Para: r-br(a)listas.c3sl.ufpr.br Uma alternativa que encontrei foi f = function(v, dat, v1=sapply(gregexpr(dat,v),"[",i=1)){ n = ifelse(v1 > 1, nchar(v), NA) substr(v, v1,n) } dados1$causaa = f(dados1$linhaa, "A39") dados1$causab = f(dados1$linhab, "A39") dados1$causac = f(dados1$linhac, "A39") dados1$causad = f(dados1$linhad, "A39") dados1$causaii = f(dados1$linhaii, "A39") dados1$causa<- with(dados1, gsub('NA', '', paste(dados1$causaa,dados1$causab,dados1$causac,dados1$causad,dados1$causaii, sep=''))) Como faço para buscar mais de um código? Tem que fazer um monte de comando como esses? Existe uma forma mais simples de fazer isso? Segue o banco anexo Desculpe enviar o arquivo anexo, mas é uma tentativa para vocês me ajudarem, já que o dput não é uma boa saida! Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942 -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 3

Agregar colunas
by Sérgio Henrique almeida da silva ju 18 Abr '12

18 Abr '12

Olá Estou tentando agregar três colunas em uma, porém quando faço isso os NAs também se agregam, como faço para descosiderá-los Ex: # Dados tmp = data.frame(x=c('A1', NA, 'A2', NA, 'A3', NA, NA), y=c(NA, NA, NA, 'A4', NA, 'A5', NA), z=c(NA, 'A6', NA, NA, NA, NA, 'A7')) # Agregando com o PASTE tmp$W <-paste(tmp$x,tmp$y,tmp$z, sep = "") # Como os dados ficam dput(tmp) structure(list(x = structure(c(1L, NA, 2L, NA, 3L, NA, NA), .Label = c("A1", "A2", "A3"), class = "factor"), y = structure(c(NA, NA, NA, 1L, NA, 2L, NA), .Label = c("A4", "A5"), class = "factor"), z = structure(c(NA, 1L, NA, NA, NA, NA, 2L), .Label = c("A6", "A7"), class = "factor"), W = c("A1NANA", "NANAA6", "A2NANA", "NAA4NA", "A3NANA", "NAA5NA", "NANAA7")), .Names = c("x", "y", "z", "W"), row.names = c(NA, -7L), class = "data.frame") # Gostaria que minha coluna W ficasse assim e não como acima: W = c("A1", "A6", "A2", "A4", "A3", "A5", "A7") Como faço Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 3

Configuração de e-mails - Lista
by Danilo Cuzzuol Pedrini 17 Abr '12

17 Abr '12

Prezados, Gostaria de saber se é possível escolhermos de receber apenas um relatório diário sobre o que foi enviado à Lista? Pelo menos para mim, a caixa de e-mails fica cheia de assuntos os quais não posso ajudar e/ou tenho interesse. A coisa fica pior quando me ausento uns 4 dias do e-mail. Atenciosamente, Danilo Cuzzuol Pedrini www.danilocp.blogspot.com

2 1

Comparar duas regressões no mesmo plot
by Marcelo Claro de Souza 17 Abr '12

17 Abr '12

Plotei a seguinte regressão hipotética e gostaria de saber se existe alguma forma se verificar se as reg 1 e 2 são significativamente diferentes, talvez comparando ângulo de inclinação das retas? Tb gostaria de saber como faço p inserir na legenda o que significa cada reta, pois só consegui colocar o que significam os pontos. x1=c(1,2,3,4,5,6,7,8,9,10) x2=c(13,14,15,16,17,18,19) y1=c(30,40,50,60,70,80,90,100,110,120) y2=c(20,22,23,24,25,26,27) x=c(1,2,3,4,5,6,7,8,9,10,13,14,15,16,17,18,19) y=c(30,40,50,60,70,80,90,100,110,120,20,22,23,24,25,26,27) reg1=lm(x1~y1) reg2=lm(x2~y2) simbolos=c(rep(19,10),rep(2,7)) plot(x~y,pch=simbolos) abline(reg1,lty=2) abline(reg2,lty=1) legend("bottomright",c("D","W"),pch=c(19,2)) Muito obrigado pela ajuda. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

3 2

ajuda com sample
by DANIELA CUNHA DA SE 17 Abr '12

17 Abr '12

Olá Pessoal, Gostaria de saber se alguém pode me ajudar: Quero sortear uma amostra dentro de um data.frame com varias colunas.....porém quero resgatar da amostra selecionada as informações das outras colunas. Exemplo a1<-sample(dados$parc,10) [1] 21 51 59 3 30 65 10 2 13 47 exemplo a primeira unidade selecionada é a parcela 21 (da coluna dados$parc do data.frame), como relacionar a essas parcelas os valores das outras colunas? Obrigada pela atenção.

2 1

Ajuda com mapas
by Renato Braga 17 Abr '12

17 Abr '12

Prezados, Sou iniciante com mapas: vocês podem me encaminhar/indicar material de estudo? Como plotar as bordas de um mapa, definir limites (municípios, por exemplo)? Obrigado,

2 1

modelos mistos - funçao para nested ancova
by Catarina Jakovac 17 Abr '12

17 Abr '12

Olá, Alguém sabe qual a funçao para rodar uma ancova hierárquica? O modelo que quero testar é o crescimento de plantulas (RGR) em funçao do tipo de capoeira (Standtype), espécies de plantulas (spp), tratamento de abertura do dossel (treatment) e tendo como covariáveis a idade da capoeira (age) e os niveis de luz que atingem a plantula (light). Como o tratamento de abertura do dossel é pareado dentro de cada capoeira, preciso incluí-lo com um random effect. Desejo testar também todas as interaçoes possíveis. Pelo que eu estudei até agora, usando o pacote lme4, o modelo ficaria algo parecido com: ancovaRGR<-lme(RGR~Standtype + treat + spp + age + light + Standtype* treat *spp *age + Standtype* treat *spp *light, data = seedlings, random = ~ 1|num/treat/seedlingID) Mas esse modelo não está rodando corretamente. Alguém sabe o que está errado, e se essa é mesmo a função que eu devo usar? Obrigada, Catarina End -- MSc. Ana Catarina C. Jakovac PhD. Candidate Forest Ecology and Forest Management Group Centre for Ecosystem Studies Wageningen University and Research Centre P.O. Box 47 - Lumen Building room 1.212 6700 AA Wageningen, The Netherlands -- MSc. Ana Catarina C. Jakovac PhD. Candidate Forest Ecology and Forest Management Group Centre for Ecosystem Studies Wageningen University and Research Centre P.O. Box 47 - Lumen Building room 1.212 6700 AA Wageningen, The Netherlands -- MSc. Ana Catarina C. Jakovac PhD. Candidate Forest Ecology and Forest Management Group Centre for Ecosystem Studies Wageningen University and Research Centre P.O. Box 47 - Lumen Building room 1.212 6700 AA Wageningen, The Netherlands

4 9

(sem assunto)
by sznelwar 16 Abr '12

16 Abr '12

3 2

Problema com barplot múltiplo
by ASANTOS 15 Abr '12

15 Abr '12

Boa noite pessoal, Tenho 3 tratamentos hipotéticos avaliados em 3 tempos diferentes, porém no tempo dois não foi avaliado o tratamento 2. Para fazer a representação deles em um barplot múltiplo, vai ficar faltando uma barra, portanto vai ficar um espaço em branco entre o tratamento 1 e o tratamento 3 e gostaria de remover esse espaço (por motivo de uma força maior chamada reviewer do periódico) e deixar lado a lado as barras dos tratamentos 1 e 3 no tempo 2. Meu CMR seria: require(gplots) #tratmentos no tempo1 trat1<-rep(1:3,5) resp1<-rpois(15,25) med1<-tapply(resp1, trat1,mean)##Media dos tratamentos ep1<-tapply(resp1,trat1,sd)/sqrt(tapply(resp1,trat1,length))###Erros padrao da media #tratmentos no tempo2 trat2<-c(1,1,1,1,1,2,2,2,2,2,3,3,3,3,3) resp2<-c(7,7,11,12,8,NA,NA,NA,NA,NA,8,8,10,9,6)##Não foi avaliado o tratamento 2 med2<-tapply(resp2, trat2,mean, na.rm = TRUE)##Media dos tratamentos ep2<-tapply(resp2,trat2,sd, na.rm = TRUE)/sqrt(tapply(resp2,trat2,length))###Erros padrao da media #tratmentos no tempo3 trat3<-rep(1:3,5) resp3<-rpois(15,45) med3<-tapply(resp3, trat3,mean)##Media dos tratamentos ep3<-tapply(resp3,trat3,sd)/sqrt(tapply(resp3,trat3,length))###Erros padrao da media ### # Montando as tabelas para cada mêss: a<- med1 b<- med2 cc<-med3 a2<-ep1 b2<-ep2 c2<-ep3 # Unindo as tabelas TAB <- cbind(a,b,cc) TAB2<-cbind(a2,b2,c2) # Criando o barplot múltiplo mp <- barplot2(TAB, beside = TRUE, axisnames = FALSE, plot.ci=T, ci.u=TAB+TAB2,ci.l=TAB-TAB2, ylab="Number of insects", legend.text=T, ylim=c(0, 60),col =c("grey25","grey50","white")) # Adicionando os meses mtext(1, at = colMeans(mp), text = rep(c("April", "May", "June"), 7),line = 1, cex = 0.8) # Então eu pergunto se é possível fazer isso, pois a função barplot2() não me permite alterar isso uma vez que tenho no tempo 2 o tratamento 2 representado por NA, por outro lado se retiro o tratamento 2 na hora de unir as tabelas não posso fazer um cbind(), pois tenho vetores de média e erro padrão com comprimentos diferentes. Se algum puder me dar uma luz agradeço, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

4 9

[OFF] Linguagem Julia.
by Pedro Rafael 15 Abr '12

15 Abr '12

Linguagem Julia, que já ouviu falar e quem já a usa. Ouvi rumores de pessoas da estatística falando que ela pretende ser um "R" mais otimizado, isso é verdade? Vocês acreditam que ela venha a ser muito popular no meio estatístico? Link: http://julialang.org/ -- Saudações, Pedro Rafael Diniz Marinho.

2 3

Webinar: R for SAS and SPSS Users
by Mauro Sznelwar 15 Abr '12

15 Abr '12

Repassaram para mim, vejam se tem alguma utilidade. Webinar: R for SAS and SPSS Users This free webinar will provide a general introduction to R and also give insight into the differences between R, SPSS and SAS. The session will then offer a brief tutorial on how to create simple Bar Plots in R. The webinar will be sponsored by The Orange County R User Group and Lion Data Systems, LLC. If this webinar's start time does not work for you, please email info(a)liondatasystems.com to let us know a start time that would be better for your country. We can choose different start times for our future R webinars based on audience demand. You can also tweet about the event using a message similar to the following: #Rstats for #SAS and #SPSS users webinar with @LionDataSystems liondatasystems.com/events #in You and your colleagues can register for this webinar at: liondatasystems.com/events Webinar Date: May 9, 2012 Webinar Time: 9 am Pacific Time Webinar Duration: About 1 hour Please email info(a)liondatasystems.com with questions.

1 0

Problema pra criar Função - duvidas com rowSums
by Augusto Ribas 14 Abr '12

14 Abr '12

Ola pessoal, estou com problemas pra deixar uma função generica o suficiente, acho que provavelmente a solução deve ser bem simples, mas não estou conseguindo solucionar alguem poderia me dar uma ajudinha? Minha intenção a partir de uma matriz de network entre hospedeiros e parasitos conseguir calcular quantos parasitos morrem ao um ou mais hospedeiros serem extintos. em baixo eu dou um exemplo de matriz e a função que num funciona direito para quando so 1 hospedeiro morre. #vamos supor uma matrix onde 1 indica que aquele parasita(linha) vive #dentro do hospedeiro tal(coluna) #são 6 parasitos e 6 hospedeiros dados<-matrix(c(1,1,1,0,0,0, 0,0,1,0,0,0, 0,0,1,0,0,0, 0,0,0,1,0,0, 0,0,0,1,1,0, 0,0,0,0,0,1),byrow=T,ncol=6,nrow=6, dimnames=list(paste("P",1:6,sep=""),paste("H",1:6,sep=""))) dados #Separando os hospedeiros hospedeiros<-dimnames(dados)[[2]] hospedeiros #Sorteando uma ordem de extinção para os hospedeiros extinto<-sample(hospedeiros,ncol(dados),replace=F) extinto #função pra ver qts parasitas extinguem dados a morte dos hospedeiros extintos extcas<-function(dados,extintos) { interações<-dimnames(dados)[[1]][which(rowSums(dados[,extintos]>0)>0)] spextintas<-sum(rowSums(dados[interações,]>0)<=length(extintos)) return(spextintas) } #ela funciona mais ou menos #se todos os hospedeiros morrem, todos os parasitos morrem extcas(dados,extinto) #se morre 3 hospedeiros, os 3 primeiros morre nesse caso... extcas(dados,extinto[1:3]) #o mesmo acima pra 2 extcas(dados,extinto[1:2]) #Agora qd tento pra ver so pra 1 hospedeiro não funciona extcas(dados,extinto[1:1]) #o problema ta no rowSums #eu to usando esse comando pra ver qts vezes cada parasito interage com #os hospedeiros rowSums(dados[,extinto]>0) #a questão é até 2 ele funciona rowSums(dados[,extinto[1:2]]>0) #mas no 1 ele nao funciona rowSums(dados[,extinto[1:1]]>0) #pq ai nao fica uma matriz fica um vetor #então tem como eu faze o R entender isso dados[,extinto[1:1]]>0 #como uma matriz de 1 coluna e 6 linhas? Ai acho que o rowSums funcionaria #ou não sei se existe alguma solução mais interessante pra chegar no #mesmo resultado #mas tem que ser uma matrix que nao perca a informação do nome da coluna nem da linha pelo que entendo #minha meta final é fazer um grafico onde no eixo x é o numero de #especies extintas acumulado de hospedeiros e no Y qts parasitas #extinguiram juntos #então se essa função funcionar eu faço um loop pra cada número de #hospedeiros mortos, qts parasitas morrem juntos. #vetor pra receber o resultado resultado<-rep(NA,6) resultado #loop for(i in 2:6) { resultado[i]<-extcas(dados,extinto[1:i]) } resultado #grafico plot(resultado~c(1:6),type="b",xlab="Número de hospedeiros extintos", ylab="Número de parasitos extintos") #eu quero chegar nesse grafico, que fica faltando o resultado pro primeiro #hospedeiro morto. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Dados como fator para numérico
by Sérgio Henrique almeida da silva ju 14 Abr '12

14 Abr '12

Olá Estou trabalhando com alguns dados numéricos e que o R esta considerando como fator. > dput(dados) structure(list(NU_LESOES = c(12, 2, 1, 2, 2, 4, 2, 12, 12, 9, 1, 5, 6, 7, 1, 8, 10, 3, 4, 11)), .Names = "NU_LESOES", row.names = c(NA, -20L), data_types = "C", class = "data.frame") >str(dados2) 'data.frame': 20 obs. of 1 variable: $ NU_LESOES: Factor w/ 12 levels "1","10","15",..: 12 2 1 2 2 4 2 12 12 9 ... - attr(*, "data_types")= chr "C" >table(dados2$NU_LESOES) 1 10 15 2 20 3 30 5 7 8 80 NA 3 4 1 2 1 1 1 1 1 1 1 3 Porém quando uso o comando: dados2$NU_LESOES=as.numeric(dados2$NU_LESOES) os dados mudam > dados2$NU_LESOES=as.numeric(dados2$NU_LESOES) > table(dados2$NU_LESOES) 1 2 3 4 5 6 7 8 9 10 11 12 3 4 1 2 1 1 1 1 1 1 1 3 O que pode está acontecendo? -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

3 10

Cálculo de vértices à partir do centróide
by ASANTOS 14 Abr '12

14 Abr '12

Boa tarde pessoal, Bom descobri onde estou errando, quando faço: for(i in 1:(length(xcc))){ for(j in 1:(length(ycc))){ O loop vai pegar a primeira coordenada de xcc e fazer todas as combinações com ycc. Na verdade, eu queria uma ajuda para o loop e fazer com que ele pegue apenas o primeiro valor de xcc e comparar apenas com o primeiro valor de ycc, seguido de xcc[2] com ycc[2] ate xcc[11] com ycc[11], meu script completo é: xcc<-c(672036,672332,672674,672079,671386,671364,670753,671854,669986,670860,671366) ycc<-c(6783650,6784105,6783470,6782354,6782416,6782418,6782199,6783073,6781853,6781791,6781911) res=NULL for(i in 1:(length(xcc))){ for(j in 1:(length(ycc))){ xa<-xcc[i]-20 ya<-ycc[j]+20 v1<-cbind(xa,ya) xb<-xcc[i]+20 yb<-ycc[j]+20 v2<-cbind(xb,yb) xc<-xcc[i]+20 yc<-ycc[j]-20 v3<-cbind(xc,yc) xd<-xcc[i]-20 yd<-ycc[j]-20 v4<-cbind(xd,yd) lim.p<-rbind(v1,v2,v3,v4,v1) write.table(lim.p, file ="parc.mod", col.names=T, row.names=F,append=T) result=read.table("parc.mod",header=F) }} Obrigado -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 2

Para quem eh fa~ de "bleeding-edge tech" (em termos de R)
by Benilton Carvalho 14 Abr '12

14 Abr '12

As promessas sobre as quais comentei em algum ponto no passado (de que o R-core iria incorporar suporte a objetos mais longos que 2^31-1) comecam a ser cumpridas.... Eu imagino que no R-2.16.x eh possivel que jah tenhamos alguma dessa funcionalidade (lembrem-se que o release schedule do R mudou e agora eh apenas 1x por ano)... Enfim, o texto abaixo informa as mudancas que ja' foram feitas no codigo-fonte. b \subsection{LONG VECTORS}{ \itemize{ \item There are the beginnings of support for vectors longer than \eqn{2^{31}-1}{2^31 - 1} elements on 64-bit platforms. This applies to raw, logical, integer, double, complex and character vectors, as well as lists. (Elements of character vectors remain limited to \eqn{2^{31}-1}{2^31 - 1} bytes.) \item What can be done with such vectors is currently very limited, and most operations will return the error \sQuote{long vectors not supported yet}. They can be serialized and unserialized, \code{identical()} and \code{object.size()} work and means can be computed. \item \code{dist()} can produce dissimilarity objects for more than 65536 rows (but for example \code{hclust()} cannot process such objects). \item \code{serialize()} to a raw vector is no longer limited in size on 64-bit platforms. \item The C-level function \code{R_alloc} can now allocate \eqn{2^{34}}{2^34} or more bytes on 64-bit platforms. } }

2 1

Tinn-R com Windows 7 64 bit, R 2.15.0
by Jose Claudio Faria 14 Abr '12

14 Abr '12

Usuários do Tinn-R, Testemunho que praticamente não usei Windows nos últimos 2/3 anos. Nesse tempo estava ajudando o Jakson Aquino no projeto do Vim-R-plugin (que considero hoje uma das melhores alternativas de GUI para usuários intermediários e avançados de computador e R). Pois bem, esses dias decidi portar alguns recursos gráficos de um pacote que roda apenas sob Windows e com o R em 32 bit para todas as plataformas suportadas pelo R. Precisava então estudar o código fonte do pacote no Windows com o R rodando 32 bit. Meu velho Windows Vista 32 bit estava bem pesado, ele veio num notebook 64 bit que havia comprado a alguns anos. Instalei o Windows 7 64 bit, instalei o R (32/64) e o Tinn-R (confesso que estava com saudades...). Estou trabalhando a 2 semanas no projeto do pacote e gostaria de informar que nunca tinha usado esses programas nessa plataforma. Testemunho que o Tinn-R me surpeendeu... muito... ele vem atravessando bravamente a evolução do Windows e do R. Para quem usa Windows, acredito que ele ainda tem muito contribuir, e por muitos anos, pois sua plataforma de desenvolvimento (Object Pascal) foi/é ainda muito consistente. Bom, sugestões de melhorias são sempre bem vindas: Em 2 ou 3 meses deverá sair uma nova versão corrigindo umas coisinhas pequenas que não atrapalham seu uso efetivo. Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Ajuda para evitar um loop com elegância (agora correto)
by Jose Claudio Faria 14 Abr '12

14 Abr '12

Pessoal, Gostaria de evitar um loop no procedimento abaixo da forma mais elegante possível: # R help m <- matrix(runif(24), ncol=3) res <- matrix(0, nrow(m), ncol(m)) mg <- mean(m) # média geral mc <- colMeans(m) # média das colunas mr <- rowMeans (m) # médias das variáveis for (i in 1:nrow(m)) for (j in 1:ncol(m)) res[i,j] <- m[i,j] + mg - mc[j] - mr[i] res Toda ajuda é bem vinda... -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 7

SciELO Books
by Benilton Carvalho 14 Abr '12

14 Abr '12

Uma sugestao de um modo legal de distribuicao de livros (agradecam tbm `a licenca Creative Commons): http://books.scielo.org/ e algo um pouco mais "especifico": http://books.scielo.org/id/p2qh6 benilton

4 3

Ajuda para evitar loop com elegância
by Jose Claudio Faria 14 Abr '12

14 Abr '12

Pessoal, Gostaria de evitar um loop no procedimento abaixo da forma mais elegante possível: m <- matrix(1:24, ncol=3) res <- matrix(0, nrow(m), ncol(m)) mg <- mean(m) # média geral mc <- colMeans(matrixdata) # média das colunas mr <- rowMeans (matrixdata) # médias das variáveis for (i in 1:nrow(m)) for (j in 1:ncol(m)) res[i,j] <- m[i,j] + mg - mc[j] - mr[i] res Toda ajuda é bem vinda... -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 1

Juntando os dados de dois objetos um abaixo do outro
by ASANTOS 14 Abr '12

14 Abr '12

Bom dia Pessoal, Criei um objeto com duas variáveis: v1<-1:11 v2<-rnorm(11) tab1<-cbind(v1,v2) e outro: v1<-1:11 v2<-rnorm(11) tab2<-cbind(v1,v2) e gostaria de juntar o tab 1 e tab 2 um abaixo do outro em um novo objeto. Fiz algumas pesquisas e encontrei o pacote plyr que permite especificar o local onde se deseja a fusão, mas não funcionou não sendo: require(plyr) tab12<- join(tab1,tab2,type="inner") > tab12 v1 v2 [1,] 1 2.10533290 [2,] 2 -1.85441266 [3,] 3 0.95619644 [4,] 4 1.15970056 [5,] 5 -0.60433624 [6,] 6 0.73799947 [7,] 7 -1.28049285 [8,] 8 -2.63621583 [9,] 9 1.36940938 [10,] 10 -0.02570176 [11,] 11 0.99026465 Que não me dá o resultado esperado, pois esta faltando o tab2. Alguém teria alguma sugestão, Obrigado, Alexandre -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

5 8

Problemas r-br com yahoo mail
by ASANTOS 13 Abr '12

13 Abr '12

Caros colegas do r-br, Alguém saberia me dizer se existe algum problema do local onde esta hospedado o r-br@listas e o yahoo mail, não sei se alguém esta passando por isso, mas estou recebendo as 22h mensagens postadas as 11h e algumas mensagens minhas estão chegando a lista duplicadas. Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 1

Converter arquivos RData em csv dentro de um loop
by Marcos Silva 13 Abr '12

13 Abr '12

Olá Pessoal, Eu tenho, em um diretório, um conjunto de mais ou menos 30 arquivos RData que armazenam data frames. O que eu gostaria de fazer é escrever um scrip que lesse os nomes dos arquivos RData no referido diretorio e os exportasse no formato .csv. A abordagem que estou tentanto desenvolver é a seguinte: arquivos <- list.files() for(nomes in arquivos){ dframe <- load(nomes) load(nomes) write.csv2(????, file=paste(dframe, ".csv", sep=""), row.names=FALSE) rm(list=ls()) } O problema para o qual ainda não consegui vislumbrar uma solução é como usar o write.csv2() nesse contexto. Os arquivos RData tem nomes no seguinte padrão: "Emp_DETRAN_2009.2012fev.RData" e o data frame contido no arquivo tem nome "empDETRAN". Basicamente o que muda de um arquivo para o outro é o nome do do órgão. Possivelmente exista outra abordagem menos macarronica e mais direta ao ponto... Abs. -- Marcos F. Silva http://sites.google.com/site/marcosfs2006

4 3

arredondamento de números no R
by Vitor Aguiar 13 Abr '12

13 Abr '12

Caros, quero fazer o seguinte cálculo: pvalue = 1 - pchisq(sum(X2, na.rm = TRUE), dfreedom(X2)) #dfreedom( ) é uma função escrita por mim mas usando os meus dados isso dá sempre 0 (zero), porque a parte do pchisq( ) é um valor muito próximo de 1, então o R arredonda para 1, ficando 1 - 1 = 0. Mesmo que o verdadeiro valor desse cálculo seja algo muitíssimo pequeno, como zero seguido de vários zeros e eventualmente algum outro número, eu preciso saber qual é. Como fazer isso? Desde já agradeço a quem puder ajudar. Vitor Aguiar

5 10

Rscript
by Leandro Marino 13 Abr '12

13 Abr '12

Caros, alguem saberia me dizer aonde fica localizada a pasta Rscript em uma instalação padrão em Linux? Como posso localizá-la? Atenciosamente, Leandro Marino http://www.leandromarino.com.br (Fotógrafo) http://est.leandromarino.com.br/Blog (Estatístico) Cel.: + 55 21 9845-7707 Cel.: + 55 21 8777-7907

6 8

ajuda com integral
by DANIELA CUNHA DA SE 13 Abr '12

13 Abr '12

Olá a todos, Estou com o seguinte problema, quero calcular a integral de 0 até HT( valor da altura total dde cada arvore presente na base de dados em anexo). Rotininha: rm(list=ls(all=T)) setwd("C:/Documents and Settings/Daniela/Desktop/Ajustes no R") mdados<-read.table("INTEGRAL.txt", h=T, dec=",");mdados # funçao de afilamento fx<-function(hi) 0.91-(1.28*(hi/HT)^1)+(2.98*(hi/HT)^2)-(4.50*(hi/HT)^3)+(2.58*(hi/HT)^4)-(0.69*(hi/HT)^5) for (i in 1:nrow(mdados))##numero total de arvores do arquivo { integral<-integrate(fx,lower=0,upper=mdados$HT) volume<-with(mdados, ((pi*(mdados$DAP^2))/40000)*integral) mdados<-data.frame(mdados, volume) } Mas não está dando certo. A equação (polinomio do 5 grau) gerada é: 0.91-(1.28*(hi/HT)^1)+(2.98*(hi/HT)^2)-(4.50*(hi/HT)^3)+(2.58*(hi/HT)^4)-(0.69*(hi/HT)^5) Objetivo, calcular somente a integral para cada árvore (ex: 33 observações), para nao precisar faze-lo manualmente ao introduzir árvore a árvore (HT). O arquivo está no excel, a integral é de 0 até HT. O DAP servirá somente para o calculo do Volume. Alguém pode me ajudar? Meu computador é Windows. Agradeço a todos, desde já! Daniela

2 1

Problema para armazenar vértices de parcelas
by ASANTOS 13 Abr '12

13 Abr '12

Boa tarde pessoal, Tenho os centróides em UTM de 11 parcelas quadradas de 40x40m e estou tentando calcular os vértices, inclusive com a repetição do primeiro vértice para fechar o polígono. Ao final, quero armazenar os resultados. Para tanto fiz: xcc<-c(672036,672332,672674,672079,671386,671364,670753,671854,669986,670860,671366)##latitude dos centróides ycc<-c(6783650,6784105,6783470,6782354,6782416,6782418,6782199,6783073,6781853,6781791,6781911)##longitude dos centróides Fazendo o cálculo dos quatro vértices para cada par de coordenada: res=NULL for(i in 1:(length(xcc))){ for(j in 1:(length(ycc))){ xa<-xcc[i]-20 ya<-ycc[j]+20 v1<-cbind(xa,ya) xb<-xcc[i]+20 yb<-ycc[j]+20 v2<-cbind(xb,yb) xc<-xcc[i]+20 yc<-ycc[j]-20 v3<-cbind(xc,yc) xd<-xcc[i]-20 yd<-ycc[j]-20 v4<-cbind(xd,yd) lim.p<-rbind(v1,v2,v3,v4,v1) write.table(lim.p, file ="parc.mod", col.names=T, row.names=F,append=T) result=read.table("parc.mod",header=F) }} O arquivo result esta saindo todo bagunçado e com mais de 60 linhas que deveriam corresponder a cabeçalho e cinco pares de coordenadas para todos os centróides, porém quando inspeciono o objeto lim.p, os cálculos estão corretos para o último par de coordenadas do centróide. Algum poderia me dizer onde posso estar errando. Obrigado, Alexandre -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

1 0

Pergunta boba
by Sérgio Henrique almeida da silva ju 13 Abr '12

13 Abr '12

No meu conjunto de dados têm muitos dados em branco e o R não reconhece como NA, como posso fazer para que ele transforme o "em branco" para NA? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

3 12

Regressão não-linear
by Marcelo Claro de Souza 13 Abr '12

13 Abr '12

Existe algum pacote onde o R teste automaticamente vários modelos não-lineares sugerindo o mais adequado ao meiu conjunto de dados ou preciso testar manualmente diversos modelos como Gaussiano, Bolztman, Exponencial, etc... até encontral o ideal? Muito obrigado. Marcelo Claro de Souza Biologist, PhD student in Plant Biology Institute of Bioscience - UNESP, Brazil

3 2

cut
by Vinicius Brito Rocha 13 Abr '12

13 Abr '12

Pessoal, estou usando o comando cut, para criar a legenda de uma mapa do brasil que estou plotando. O problema é que minha variável tem uma escala muito alta e por isso ao executar o comando cut estou ficando com algo parecido com isso: cut(mapa$dados),10) [1] (397,1.32e+05] (397,1.32e+05] (397,1.32e+05] [4] (397,1.32e+05] (3.95e+05,5.27e+05] (3.95e+05,5.27e+05] [7] (1.32e+05,2.64e+05] (397,1.32e+05] (1.32e+05,2.64e+05] [10] (1.32e+05,2.64e+05] (397,1.32e+05] (397,1.32e+05] O que preciso é que de uma notação mais amigável para o cut. Alguém pode me ajudar? -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

3 2