Fevereiro 2014 - R-br

Função bitmap
by Roberto Guevara 10 Fev '14

10 Fev '14

Olá, usuários do R. Estou tentando utilizar a função bitmap para salvar figuras em formato TIFF com resolução própria para publicação. Quanto tento usá-la, me deparo com a seguinte mensagem de erro, Código: bitmap("Plot.tiff", height = 4, width = 4, units = 'in', type="tifflzw", res=300) Error in system(paste(gsexe, "-help"), intern = TRUE, invisible = TRUE) : 'gswin32c.exe' not found Procurei saber como solucionar esse problema através do google. Baixei um arquivo chamado Ghostscript e tentei utilizar esse código: Sys.setenv(GS_CMD = "C:/Program Files (x86)/gs/gs9.10/bin/gswin32c.exe") E a partir disso agora aparece outra mensagem de erro: Error in system(paste(gsexe, "-help"), intern = TRUE, invisible = TRUE) : 'c:\Program' not found Alguém tem ideia de como se resolve isso? Ou, pelo menos, sabe de alguma alternativa que eu possa utilizar para obter as figuras do jeito que preciso? Obrigado pela atenção -- Roberto Guevara Ferreira Lima Biólogo - CRBio 6ª Região: 73146/06-D Mestre em Zoologia Laboratório de Ecologia e Zoologia de Vertebrados Instituto de Ciências Biológicas Universidade Federal do Pará Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110 Currículo Lattes: http://lattes.cnpq.br/7363382159007600

3 5

R Session Aborted
by Fernando Gama 10 Fev '14

10 Fev '14

Olá galera, Instalei recentemente o pacote shiny no windows. Utilizo o RStudio como IDE. Meu problema é que quando executo o diretório contendo os scripts(cliente e servidor) um fatal error é lançado com a seguinte mensagem: R Session Aborted e aí ele me força a iniciar uma nova sessão. testei com esse código nativo da biblioteca, assim como outros também: library(shiny) runExample("01_hello") Alguém já passou por este problema? OBS: quem quiser visualizar tirei print(em anexo). Abs, Fernando

1 0

Intervalo de confiança para o valor estimado numa regressão múltipla.
by Helder Gramacho 08 Fev '14

08 Fev '14

Boa tarde, Estou trabalhando com dados de imóveis utilizando análise de regressão, consegui gerar o modelo que servirá para avaliar o preço do imóvel: > modelo2=lm(log_pu~area+Idade.+quartos+suites+cond+n_elev) > summary(modelo2) Call: lm(formula = log_pu ~ a_priv + Idade. + quartos + suites + cond + n_elev) Residuals: Min 1Q Median 3Q Max -0.072013 -0.041334 -0.008192 0.032892 0.091292 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 9.0404271 0.1336608 67.637 < 2e-16 *** area 0.0032140 0.0011938 2.692 0.01846 * Idade. -0.0110019 0.0036698 -2.998 0.01028 * quartos -0.1016121 0.0405304 -2.507 0.02624 * suites -0.0909447 0.0454335 -2.002 0.06664 . cond -0.0004982 0.0001582 -3.149 0.00769 ** n_elev 0.0555213 0.0368277 1.508 0.15557 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 0.0593 on 13 degrees of freedom Multiple R-squared: 0.8394, Adjusted R-squared: 0.7652 F-statistic: 11.32 on 6 and 13 DF, p-value: 0.0001628 Agora preciso a partir dos dados do imóvel que está sendo avaliado, substitui-los no modelo e encontrar o valor estimado para o mesmo, em seguida calcular o intervalo de confiança de 80% para este valor estimado. Abaixo tem a saída de um programa comercial, no qual encontrei os valores, mas gostaria de saber como calcula-los no R. [image: Imagem inline 2] Será que alguém poderia me ajudar a calcular estes intervalos de confiança para o valor estimado? Desde já agradeço, Hélder Santos *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 1

formatação de gráfico
by Maurício Lordêlo 07 Fev '14

07 Fev '14

Saudações a todos. No CMR abaixo construo seis gráficos e estou tendo dificuldade para acrescentar uma legenda única e de posicionar um "título" no centro de cada dupla de gráficos. Fiz um desenho manualmente para verem como gostaria que ficasse o resultado final: https://www.dropbox.com/s/o7snb3i206bz4dr/graficos.jpg Consegui ter acesso ao livro "R Graphics - Paul Murrel" e encontrei esta página que talvez tenha a informação que preciso: https://www.dropbox.com/s/n96r3h1aw9az41q/pag91_92Rgraphics.jpg Porém, fiz várias tentativas (nas linhas com "comentário" estão as últimas mudanças) juntamente com as informações contidas no help da função "par" e não obtive sucesso. Grato se puderem me ajudar. Maurício par(mfrow=c(3,2) # ,plt = c(-10,-10,0,0) # ,pty="s" # ,mai = c(-5,1,-7,1) # ,mar = c(-5, 4, -4, 2) # ,oma = c(1,0,1,0) # ,pin = c(-3,-6) # ,xpd = T # ,din=c(0.5,0.5) # ,fin=c(0.5,0.5) ) #par(mfrow=c(3,2)) #layout(matrix(c(1, 2, 3, 4, 5, 6), byrow=TRUE, ncol=2)) tempo<-c(4,5,6,7,8,9,10,11) A1A10<-rnorm(8,90,4) A3A10<-rnorm(8,90,4) plot(A1A10 ~ tempo,axes = F,xlab = '',ylab ="medida",main = 'Intercepto',cex.main=1,type = 'b',xlim=c(4,11),ylim=c(80,100)) mtext('10 repetições',line = 2.9,cex = .8) lines(A3A10~tempo, type="o",pch="+",lty=2) axis(2) A1A20<-rnorm(8,90,4) A3A20<-rnorm(8,90,4) plot(A1A20 ~ tempo,axes = F,xlab = '',ylab ="medida",main = 'Intercepto',cex.main=1,type = 'b',xlim=c(4,11),ylim=c(80,100)) mtext('20 repetições',line = 2.8,cex = .8) lines(A3A20~tempo,type="o",pch="+",lty=2) axis(2) A1B10<-rnorm(8,90,4) A3B10<-rnorm(8,90,4) plot(A1B10 ~ tempo,axes = F,xlab = '',ylab ="medida",main = 'Tratamento B',type = 'b',xlim=c(4,11),ylim=c(80,100)) lines(A3B10~tempo,type="o",pch="+",lty=2) axis(2) A1B20<-rnorm(8,90,4) A3B20<-rnorm(8,90,4) plot(A1B20 ~ tempo,axes = F,xlab = '',ylab ="medida",main = 'Tratamento B',type = 'b',xlim=c(4,11),ylim=c(80,100)) lines(A3B20~tempo,type="o",pch="+",lty=2) axis(2) A1C10<-rnorm(8,90,4) A3C10<-rnorm(8,90,4) plot(A1C10 ~ tempo,axes = F,xlab = "tempo",ylab ="medida",main = 'Tratamento C',type = 'b',xlim=c(4,11),ylim=c(80,100)) lines(A3C10~tempo,type="o",pch="+",lty=2) axis(1, at = seq(4, 11, 1)) axis(2) #legend(9.5,60,c("Opção 1",""),lty = c(1,2),ncol = 1,horiz=TRUE) A1C20<-rnorm(8,90,4) A3C20<-rnorm(8,90,4) plot(A1C20 ~ tempo,axes = F,xlab = "tempo",ylab ="medida",main = 'Tratamento C',type = 'b',xlim=c(4,11),ylim=c(80,100)) lines(A3C20~tempo,type="o",pch="+",lty=2) axis(1, at = seq(4, 11, 1)) axis(2) #legend(0.5,60,c("","Opção 3"),lty = c(1,2),ncol = 1,horiz=TRUE) #legend(0,60,c("Opção 1","Opção 3"),lty = c(1,2),ncol = 1,horiz=TRUE)

3 10

Teste de Dunn
by Michele Figueiró 05 Fev '14

05 Fev '14

Caros amigos Estou analisando um banco de dados e para a analise de comparação multipla foi utilizado o teste de kruskall wallis com teste post doc de Dunn utilizando o R. Eu gostaria de saber se alguém poderia me indicar alguma referencia sobre o teste de Dunn. Eu estou com dúvidas pois este teste fornece dois p valores, o usual e o ajustado. Eu já me deparei com uma situação na qual há uma diferença significativa entre as médias e no resultado do teste de Dunn, os p valores ajustados não são significativos mas os p valores usuais, como agir neste caso?. Também gostaria de saber se há outros testes não-paramétricos post-hoc que poderiam ser utilizados ao invés do teste de Dunn. Abs Michele *"To live is the rarest thing in the world. Most people just exist." Oscar Wilde*

1 0

Análise Espacial
by Paulo Vetz 05 Fev '14

05 Fev '14

Olá Pessoal, Eu quero avaliar o padrão espacial entre algumas espécies de animais. Eu tenho os dados (coordenadas) coletados durante deferentes épocas do ano. Vocês poderiam me indicar algumas referências de métodos estatísticos para poder realizar essa análise. Estou interessado em avaliar a interação no tempo, espaço e espaço-tempo. Desde já agradeço a ajuda. Saudações, Paulo

5 7

Pacote synbreed
by Giselle Davi 05 Fev '14

05 Fev '14

Prezados, alguém saberia me dizer o que aconteceu com as funções do pacote synbreed ??? Grande parte delas não podem ser encontradas quando se roda a rotina conforme o manual. Abaixo encontra-se o início da rotina descrita no manual do pacote synbreed, onde já existe o erro de função não encontrada. Desde já agradeço. library(synbreedData) data (maize) newDHpheno <- data.frame(Trait=1000,row.names="newDH") newDHpheno #simulating genotypic data newDHgeno <- matrix(sample(c(0,1), ncol(maize$geno), replace=TRUE), nrow=1) rownames(newDHgeno) <- "newDH" #new pedigree newDHpedigree <- data.frame(ID="newDH", Par1=0, Par2=0, gener=0) #new covar information newDHcovar <- data.frame(family=NA, DH=1, tbv=1000, row.names="newDH") # add individual maize2 <- add.individuals(maize,newDHpheno,newDHgeno,newDHpedigree,newDHcovar) summary(maize2) ##erro em add.individuals: função não pode ser encontrada

3 5

Dúvida com substituição
by Edson Lira 03 Fev '14

03 Fev '14

Tenho um banco com duas variáveis: bairro_resid bairro_ntabelado <NA> cidade nova <NA> SANTA ETELVINA Gostaria de subsituir o <NA> do bairro_resid pelo cidade nova na outra linha pelo SANTA ETELVINA O banco tem 58.000 linhas, em algumas delas acontece a situação acima mostrada, como substituir o NA do bairro residencial pelo conteúdo do bairro_ntabelado em todo o banco? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

3 4

Como colocar caracteres, letras, números em subscrito ou sobrescrito em gráfico no R
by Ederson Civardi 03 Fev '14

03 Fev '14

Prezado Fernando Toledo, vi que deixei de incluir um parêntesis no final e por isso dava erro. Porém corrigido este detalhe o problema não foi resolvido e ficou desta maneira: "Lesão ("*mm^2*")", no gráfico. -- *_________________________________________* *Ederson A. Civardi* *Doutorando Produção Vegetal - UFG - PPGA - Goiânia* Cel. (Claro): (62) 9237-7817 ________________________________________________

4 7