Fevereiro 2014 - R-br

API de Georreferenciamento do Google
by Helder Gramacho 15 Fev '14

15 Fev '14

Bom dia, Estou tentando conhecer um pouco mais do R por meio da publicação "Introdução à análise de dados espacialmente referenciados" do Elias Teixeira Krainski, nela existe uma etapa na qual é possível utilizar uma API do google.maps para fazer o georreferenciamento de endereços, ou seja, obter a latitude e longitude desses endereços o código é o que segue logo abaixo: require(sp) ltipo <- c("av", "rua", NA) lnome <- c("Brasil", "joubert", "prefeitura municipal") enum <- c(3500, 100, NA) fGetLatLonLog <- function(tipo, nome, num, mun, uf, pais) { end <- paste(tipo, nome, num, mun, uf, pais, sep="+") end <- gsub(" ", "+", end, fixed=TRUE) end <- gsub("NA", "", end, fixed=TRUE) end <- gsub("++", "+", end, fixed=TRUE) end <- paste('http://maps.google.com/maps/geo?q= ',end,'&output=csv&sensor=true_or_false&key=abcdefg',sep="") end <- sapply(end, readLines, warn=FALSE) end <- t(sapply(strsplit(end, ","), as.numeric)) colnames(end) <- c("Status", "Acuracia", "Latitude", "Longitude") rownames(end) <- 1:nrow(end) return(as.data.frame(end))} ll=fGetLatLonLog(ltipo, lnome, enum, "Maringa", "PR", "BR") ll Mas quando executo a função aparece o seguinte erro: Error in file(con, "r") : cannot open the connection In addition: Warning message: In file(con, "r") : cannot open: HTTP status was '403 Forbidden' Parece que o Google não permite realizar a conexão... Será que alguém sabe como me ajudar. Hélder Santos *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

3 4

Erro em lapply
by Ludmila Rattis 15 Fev '14

15 Fev '14

Olá, Rodando a função lapply como segue: uc_list<-lapply(sp_names, function(x) data$WDPAID[data$BINOMIAL==x]) Deu o seguinte erro (apesar de continuar rodando): Error in sprintf("^%s", makeRegexpSafe(text)) : promise already under evaluation: recursive default argument reference or earlier problems? Alguém sabe me dizer o que significa? -- *Me. Ludmila Rattis* Programa de Pós-graduação em Ecologia/UNICAMP Conservation Biogeography Lab http://www.wix.com/rdloyola/lab <http://sites.ffclrp.usp.br/ficus>

3 4

Dúvida no pacote ltm
by Felipe Buchbinder 14 Fev '14

14 Fev '14

Prezados professores, Estou utilizando o pacote ltm para análise de um modelo de TRI do tipo *Graded Response Model*. Todavia, eu tenho ocasionalmente encontrado a seguinte mensagem: Error in log.pr[xj, ] : subscript out of bounds Alguém pode me ajudar a saber por que isso ocorre e como consertar o problema? Envio anexo um CMR, onde eu gero um teste e analiso com base no pacote ltm. Desde já, obrigado, Felipe

1 0

bootstrap em correlação
by geovane barbosa 13 Fev '14

13 Fev '14

Olá pessoal boa tarde Tenho a seguinte situação. Possuo um conjunto de dados caracterizados por x1,x2,x3. e gostaria de realizar o bootstrap para calcular a correlação e efetuar intervalos de confiança bootsrap. A linha de pensamento é essa. Suponhamos que x1<-rnorm(100) x2<-rnorm(100) x3<-rnorm(100) dados<-data.frame(x1,x2,x3) rx1x2<-cor(x1,x2) rx1x3<-cor(x1,x3) rx2x3<-cor(x2,x3) theta1<- sqrt[(rx1x2+rx1x3)/rx2x3] theta2<- sqrt[(rx1x3+rx2x3)/rx2x1] Gostaria de ter como resultados final os intervalos bootstrap para theta1, theta2. obrigado

2 1

quebrar amostra através de sorteio aleatório sem reposição
by Jodavid Ferreira 13 Fev '14

13 Fev '14

Olá pessoal,estou com dúvida na hora de realizar este problema, O problema é: Gero uma amostra de tamanho 100 por exemplo com dist. poisson amostra<-rpois(100,lambda=2.5) amostra [1] 1 4 8 0 2 4 2 1 2 2 6 3 3 0 2 2 1 1 1 2 1 2 4 2 2 2 1 0 4 6 5 2 1 2 2 3 3 1 2 [40] 0 3 5 5 5 2 2 2 5 3 2 2 3 3 2 2 2 3 3 0 2 2 1 2 1 1 4 2 0 3 4 2 2 3 4 4 2 4 1 [79] 2 2 1 0 2 1 4 3 2 2 0 3 3 1 3 3 2 3 3 1 5 0 agora com esses valores quero podemos dizer que quebrar em três amostas, a primeira de tamanho 25, a segunda, 40, e a terceira 35. Só que o problema é que tenho que sortear valores para a primeira amostra e estes valores tem que sair da amostra original. ou seja, quando sortear para a primeira amostra, chamar amostra1, só deve restar 75 valores para que eu possa dividir e mais duas. Estava pensando em utilizar o sample, para sortear os valores e alocar a cada amostra, mas como faço para remover os valores entre uma amostra e outra? -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

2 2

Re: [R-br] R Session Aborted
by Fernando Gama 12 Fev '14

12 Fev '14

Pronto, Agora imagem em um link externo. http://imageshack.com/a/img593/951/hekz.png Abs, Fernando. Em 11 de fevereiro de 2014 12:15, Fernando Gama <f.fabiogama88(a)gmail.com>escreveu: > Valeu Ivan, > > no linux é na hora da instalação do package. Veja na imagem em anexo. > > abs, > > Fernando. > > > Em 11 de fevereiro de 2014 11:27, Ivan Bezerra Allaman < > ivanalaman(a)yahoo.com.br> escreveu: > >> Bom dia Fernando! >> >> Eu uso a shiny no linux sem problemas. Me relate detalhadamente qual o >> problema estás tendo no linux para eu tentar te ajudar. >> >> Abraços! >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> > >

4 9

Dúvida coeficiente de regressão partial
by Fernando Souza 11 Fev '14

11 Fev '14

Caros colegas. Rodei um modelo de regressão misto multivariado por meio da função lme (pacote nlme) eu obter o coeficiente de regressão parcial das variáveis selecionadas no modelo (R2 parcial). Estou com algumas dúvidas sobre como se calcula esse coeficiente. Alguém pode me explicar ou me indicar algum bom texto (tipo for Dummies) para que eu possa me inteirar sobre o assunto? att

2 1

Re: [R-br] R Session Aborted
by Ivan Bezerra Allaman 11 Fev '14

11 Fev '14

Bom dia Fernando! Eu uso a shiny no linux sem problemas. Me relate detalhadamente qual o problema estás tendo no linux para eu tentar te ajudar. Abraços!

1 0

Re: [R-br] R Session Aborted
by Ivan Bezerra Allaman 11 Fev '14

11 Fev '14

Boa noite Fernando! Há dias estou tendo solucionar um problema parecido só que com a função "runApp" no R console. Já procurei em "quase tudo" na internet mais não encontrei nenhuma solução. Tente o seguinte, talvez funcione para você. Vá no menu iniciar e clique com o botão direito em "computador". Depois vá em "propriedades". Em seguida clique em "Configurações avançadas do sistema". Depois em "Avançado" - "Configurações". Depois na terceira aba "Prevenção de Execução de dados" e clique em "Ativar a DEP para todos os ....". Clique em "adicionar" e procure pelo o executável "Rgui.exe". Detalhe: Se o seu sistema for 64bits acredito que não dará certo, pois segundo os blbababa da microsoft, os programas 64bits não estão sob DEP. Logo, se algum santo tiver a solução para esse problema eu agradeço, senão, tanto faz, pois estou no windows temporariamente e logo voltarei para o linux onde não EXISTE problemas imbecis como este. Abraço Fernando!

2 1

Gráfico de Pareto
by Ari Clecius 10 Fev '14

10 Fev '14

Alguém poderia me dizer uma forma mais fácil de fazer um gráfico de pareto para los efeitos significativos de um experimento fatorial? Segue código abaixo. ATtt, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 #Código resp=scan(nmax=8) 82.8 55.1 74.3 72.2 64.7 51.3 69.9 53.3 #replicated rad=rep(c(+1,-1),time=4,each=1) rad pH=rep(c(+1,-1),time=2,each=2) pH area=rep(c(+1,-1),each=4) area order=(1:8) order fat=data.frame(rad,pH,area,resp,order) fat aov.2=aov(resp~rad+pH+area, data=fat) summary(aov.2) #Gráfico de pareto MSE=41.7 #Obtido a partir de coef1=c(aov.2$coef[2:4]) coef1 n=8 f1=t(abs(coef1))/(sqrt(MSE/n));f1 z=data.frame(f1) names(z) v1=qt(0.975,8) coef3=sort(abs(coef1)) coef3 g=barplot(sort(f1),horiz=TRUE,names.arg=c( names(coef3)),col=ifelse(sort(f1)>v1,"blue","gray"))##xlim, ylim box("plot", col="blue") abline(v=qt(0.975,8),col="red",lty=2) s=format(qt(0.975,8),digits=2) text(2.5,1,label=s,cex=1) mtext("Efeito padronizado", side=3, line=2, cex=1.0, col="blue", outer=FALSE) p=format(f1,scientific = TRUE) legend(x=6,y=4,cex=0.6,bg="yellow",box.col="blue", legend=c( "Rad=7.475","area=5.650","pH=-1.975" ), col="blue", lwd=1:2, lty=c(NA,NA,NA))

2 1