- R-br - DINF UFPR

Dúvida summary com factor GLM
by Thiago Cesar Lima Silveira 20 Jan '15

20 Jan '15

Bom dia Pessoal, Tenho uma dúvida sobre como apresentar os dados de um GLM rodado no R. Estou modelando a distribuição e abundâncias de um peixe em relação a preditores ambientais. As coletas foram realizadas em diferentes épocas do ano, portanto usei a estação do ano como fator para controlar essa possível influência. O problema é: Segundo alguns posts no R-br, vi o comentário do Prof. Walmes citando que o glm no R tem por padrão zerar o primeiro nível do fator. Gostaria de apresentar os dados dos coeficientes do modelo ajustado, mas não faz sentido mostrar os coeficientes para três níveis somente. Alguém teria uma sugestão? Muito obrigado. Segue um exemplo código reproduzível: # Dados y<-c(1,33,500,700, 334,320, 703,303,3030,3002,200,0,100,100,169,12) v1<-c(12,33,544,600, 34,30, 03,3390,3030,302,20,108,170,101,2009, 22) v2<-c('s','s','s','s','a','a','a','a','w','w','w','w','f','f',’f, ‘f' ) df<-data.frame(y, v1, v2) #GLM glm1<-glm(y~v1 + factor(v2), data=df) summary(glm1) __________________________________ Thiago Cesar Lima Silveira Biólogo, Me. Biologia Animal Doutorando do PPG Zoologia - PUCRS e-mail: thiagoclsilveira(a)yahoo.com.br Skype: thiagocesarls

1 0

Barchart ordenado
by Odirley Campos 19 Jan '15

19 Jan '15

# Boa noite a todos!# Estou tentando criar graficos de barra utilizando a funções do pacote lattice KinA37.5<-KinA[[1]][5] #extraindo da lista de testes KinA37.5<-as.data.frame(KinA$`37.5`$groups) #de list para data.frame KinA37.5<-KinA37.5[sort(KinA37.5$trt,decreasing=T),] #ordenando por trt KinA37.5 #visualizando KinA50<-KinA[[2]][5] #extraindo da lista de testes KinA50<-as.data.frame(KinA$`50`$groups) #de list para data.frame KinA50<-KinA50[sort(KinA50$trt,decreasing=T),] #ordenando por trt KinA50 KinA62.5<-KinA[[3]][5] #extraindo da lista de testes KinA62.5<-as.data.frame(KinA$`62.5`$groups) #de list para data.frame KinA62.5<-KinA62.5[sort(KinA62.5$trt,decreasing=T),] #ordenando por trt KinA62.5 #---juntando em um único data----# resumoKinA<-rbind(KinA37.5,KinA50, KinA62.5) id<-as.factor(rep(c(37.5, 50, 62.5),each=5)) resumoKinA<-cbind(id,resumoKinA) resumoKinA ##--- graficando ---# require(lattice) barchart(means~trt|id, data=resumoKinA, horiz=FALSE, layout=c(3,1), ylim=c(0,max(resumoKinA$means*1.3)), panel=function(x, y, subscripts, ...){ panel.barchart(x, y, subscripts= subscripts,...) panel.text(x, y, label=resumoKinA[subscripts,"M"], pos=3) } ) # Mas as barras dentro de cada grafico não estão ordenadas como os dados em "means" no resumoKinA.# Existe alguma função para ordenar os dados dentro da função barchart?# Tentei sort e order, mas não consegui. # Antecipadamente agradeço! Odirley R. Campos Engenheiro Agrônomo UFV/MGDoutorando em Solos e Nutrição de Plantas UFV/MG Tell:(31)88452031

1 0

Obter a distância de um municipio a outro
by Daniel Marcelino 19 Jan '15

19 Jan '15

Eu estou com uma dúvida e gostaria muito de ouvir os expertos da área geografia/estatística espacial. Usualmente eu vejo estudos usando a distância euclidiana de centro a centro de municípios. Mas isso não faz tanto sentido já que estes são ligados por estradas sinuosas. Minha dúvida é: alguém saberia me dizer se alguém já fez, ou se eu consigo obter a distância de um município a outro usando a malha de estradas? Evidentemente o Google faz isso, mas seria um trabalho monstruoso verificar isso para 5656 municípios. att. Daniel

2 2

Prof. Assistente vs Prof. Adjunto
by andrebvs＠bol.com.br 19 Jan '15

19 Jan '15

6 5

Dificuldade em montar matriz de contraste com sapply()
by ASANTOS 18 Jan '15

18 Jan '15

Caros Listeiros, Estou tentando montar um matriz de contraste entre 3 níveis (colony1, colony2 e colony3) dentro de um fator chamado P1 com a função sapply, mas não consigo ver onde estou errando, pois os coeficientes entre as combinações de níveis não estão ficando corretos, alguém poderia dar uma ajuda, segue CRM: ## Cria uma matriz de contraste # Níveis dentro de um fator P1 P1 <- c("colony1", "colony2", "colony3") dados <- data.frame(P1) comb <- cbind(combn(unique(c(P1,P1)), 2)[2:1,])## Combinação entre os níveis comb # Comparação entre níveis com sapply() sapply(1:ncol(comb), function(x) with(dados, 2*( (P1==comb[1,x] & P1 == comb[2,x]) | (P1==comb[1,x] & P1 == comb[2,x])) - (P1 == P1 & P1 %in% comb[,x]) ) ) dados[ , 2:4] <- sapply(1:ncol(comb), function(x) with(dados, 2*( (P1==comb[1,x] & P1 == comb[2,x]) | (P1==comb[1,x] & P1 == comb[2,x])) - (P1 == P1 & P1 %in% comb[,x]) ) ) dados## Os coeficientes não estão corretos Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

1 0

Vaga para trabalhar com R em Porto Alegre, RS
by Fernando Henrique Ferraz Pereira da Rosa 17 Jan '15

17 Jan '15

Pessoal, Estamos com uma vaga para trabalhar com Estatística e R no RS. Se tiver interesse vejam mais detalhes em: https://www.linkedin.com/pulse/last-call-internship-position-data-science-f… , https://plus.google.com/+FernandoHenriqueFerrazPereiradaRosa/posts/p/pub ou https://www.facebook.com/photo.php?fbid=10152750369209601&set=a.10151426344… . Abcs, -- "Though this be randomness, yet there is structure in't." Fernando H Rosa - Statistician http://www.fernandohrosa.com.br / http://www.feferraz.net - Estatística, Matemática e Computação AprendaAlemao.net <http://aprendaalemao.net/> - Seu ponto de partida para melhorar seu Alemão! @fhrosa

1 0

artigo - geoestatistica
by Fernanda De Bastiani 16 Jan '15

16 Jan '15

Estimados, alguém tem acesso ao seguinte artigo: Diggle, P.J. \& Ribeiro Jr, P.J. (2002) Bayesian inference in Gaussian model-based geostatistics. Geographical and Environmental Modelling, Vol. 6, No. 2, 129-146. e poderia me enviá-lo em formato .pdf por gentileza.Não tenho acesso gratuito aos artigos publicados nesta revista. Desde já agradeço,Fernanda

2 1

Re: [R-br] Problema com download de arquivos no dropbox (ASANTOS)
by Rafael Tieppo 15 Jan '15

15 Jan '15

Tem outro meio: #solucao"http://stackoverflow.com/questions/14580850/reading-raw-data-in-r-to-be-sav…" #camnho do arquivo (nao sei se vc consegue acessar) #se precisar troque pelo caminho do seu arquivo https://www.dropbox.com/s/wo1whtogsz3wmlz/EXAMPLE_DONE.csv?dl=0 #instalar Pacote install.packages("httr") library(httr) FILE <- GET(url="https://www.dropbox.com/s/wo1whtogsz3wmlz/EXAMPLE_DONE.csv?dl=0") writeBin(FILE$content, "DADOS.csv") (read.csv("DADOS.csv", head=T, sep=";")) Testei aqui e deu certo.Saudações Rafael Tieppo UNEMAT - site: http://www2.unemat.br/rafaeltieppo https://sistemasagricolas.wordpress.com/ "Evite o desperdício: antes de imprimir pense na sua responsabilidade com o ambiente". On Thursday, January 15, 2015 12:00 PM, "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> wrote:

2 1

Material Regressão logística múltipla
by ELIAS MENDES COSTA MENDES 15 Jan '15

15 Jan '15

Boa tarde pessoal, Algum de vocês tem material sobre regressão logística múltipla? (De preferencia usando o R) Desde de já agradeço Att,Elias Mendes Costa

4 3

Distâncias no R
by Roberto Guevara 15 Jan '15

15 Jan '15

Olá, usuários de R. Minha dúvida é em relação ao cálculo de distâncias entre unidades amostrais para dados multivariados utilizando as funções dist e vegdist. Sabendo que o R calcula sempre as dissimilaridades, ao fazer uma análise de agrupamento baseado em algum índice de similaridade utilizando a função hclust eu devo diminuir os valores de 1 previamente ou a função já converte as dissimilaridades em similaridades por si só? 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

1 0