- R-br - DINF UFPR

MANTEL
by Cézar Borges 05 Jun '13

05 Jun '13

Olá caros,Realizei um teste de mantel para 19 fragmentos para verificar a dependência espacial da comunidade arbórea, o resultado (r=-0.01; p=0.55) me mostrou que os fragmentos são independentes espacialmente. No entanto, estou com a seguinte dúvida: posso interpretar este resultado como se todos os fragmentos fossem diferentes entre si??? ou se não, como posso estar verificando quais os fragmentos que são iguais e diferentes entre si???Abraços! Att, Cézar Fonseca Borges Mestrando em Ciências Ambientais UFPA/MPEG/EMBRAPA - Amazônia Oriental Geógrafo Email: cezar-fonseca(a)hotmail.com/ caborges(a)museu-goeldi.br Lattes: http://lattes.cnpq.br/9805721240137485

1 0

teste de Jonckheere
by Cristina 04 Jun '13

04 Jun '13

Olá pessoal! Gostaria de saber se alguém já usou, no R, o teste de Jonckheere. tentei usar o teste mas consigo encontrar qual pacote está instalado. Obrigada Teresa Cristina

3 3

Coordenadas UTM para Lat Lon
by Tito Conte 04 Jun '13

04 Jun '13

Alguém poderia me explicar como faço a converção de utm para lat lon, estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão Tito Conte

4 4

erro na validação cruzada
by Hélio Gallo Rocha 04 Jun '13

04 Jun '13

Caros da Lista, Está dando este erro ao tentar rodar a validação cruzada na direção de 45 graus: xvalid: performing cross-validation at location ... 1, Erro em solve.default(Vcov, trend.d) : sistema é computacionalmente singular: condição recíproca número = 8.6266e-21 Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Remover compilação
by Sérgio Henrique almeida da silva ju 03 Jun '13

03 Jun '13

Amigos Estou usando a versão mais nova do R (3.0.1) no ubuntu. Tive que compilar a versão 2.14, para fazer um trabalho, agora quero remover essa compilação, como faço? Sempre quando chamo o R no terminal aparece a versão mais antiga, como faço pra voltar a chamar o R 3.0.1? 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) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

2 2

Separação de duas distribuições de dados
by ASANTOS 03 Jun '13

03 Jun '13

Boa tarde pessoal, Simulei duas distribuições normais e representei em um histograma que obviamente assumiu a forma bimodal, sendo: ### set.seed(10) dist1<-rnorm(n=400,mean=10,sd=5) dist2<-rnorm(n=400,mean=35,sd=8) dist.c<-rbind(dist1,dist2) # hist(dist.c) # Gostaria de saber se existe algum procedimento que me permita separar dist1 e dist2 com seus respectivos intervalos de confiança, imaginando que eu não conheça os parâmetros media e desvio padrão que os criaram? Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

2 1

Plotagem 3D
by Tiago Souza Marçal 03 Jun '13

03 Jun '13

Bom dia pessoal, estou querendo plotar um gráfico tridimensional com os escores canônicos que estão abaixo: dados<-structure(list(V1 = c(-18.7175, -20.1915, -21.1557, -17.5965, -18.142, -13.857, -18.0086, -14.0436, -12.571), V2 = c(38.3004, 36.4955, 34.4964, 39.9303, 40.4819, 38.3763, 34.9869, 35.0723, 35.9022), V3 = c(-20.9237, -20.12, -17.8586, -17.0183, -17.9961, -20.8457, -19.5551, -17.5703, -18.7393)), .Names = c("V1", "V2", "V3"), class = "data.frame", row.names = c(NA, -9L)) require(lattice) zlab <- "VC2" xlab <- "VC1" ylab <- "VC3" wireframe(V2~V1*V3, scales=list(arrows=FALSE), zlab=list(zlab, rot=90), xlab=list(xlab, rot=24), ylab=list(ylab, rot=-37),xlim=c(-12,-20),ylim=c(-18,-22),type=c("on","bottom")) No entanto, não estou conseguindo fazer com que os pontos apareçam no gráfico e gostaria de eliminar as 3 arestas frontais. Agradeço desde já pela ajuda. att. Tiago ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

2 2

gráfico
by alanarocha＠sapo.pt 03 Jun '13

03 Jun '13

Boa tarde, porque será que o gráfico não é desenhado? pode ser por as medianas serem iguais? código: setwd("G:/SIN/Users/Ana Rua/dados não Panakeia e Panakeia") dados<-read.table("Dadoshgbnaophgbp2sem2011.csv",header=TRUE,sep=";", stringsAsFactors=TRUE, dec=",") dados #objectivo é: #verificar se há diferenças estatisticamente significativas entre os custos e os valores de hemoglobina para clinicas não panakeia no 2semestre de 2011. # Imprime as hipoteses do teste cat("\n\t Teste Wilcoxon\n") cat("Teste bicaudal:\n\tH0: há diferenças significativas entre os valores de hemoglobina para clínicas não Panakeia e Panakeia, mediana=0", "\n\tH1: Não há diferenças significativas entre os valores de hemoglobina para clínicas não Panakeia e Panakeia, mediana diferente de 0 \n") wilcox.test(dados$hgbnaopanakeia, dados$hgbpanakeia, paired=TRUE) boxplot(dados, axes=FALSE) # O parametro axes=FALSE faz com que os eixos não sejam desenhados title(col.main='blue',main= 'Valores de hemoglobina das clínicas não Panakeia e Panakeia no 2semestre de 2011') # Acrescenta o título do gráfico title(col.lab='blue',xlab= 'meses do 2semestre do ano 2011') # Acrescenta o titulo do eixo X title(col.lab="blue", ylab= 'Valores') # Acrescenta o titulo do eixo Y box() # Desenha a caixa em torno do gr?fico axis(2) # Coloca o eixo Y (o 2 ali representa o eixo Y) no gráfico axis(1, 1:6, c('Jul', 'Ago', 'Set', 'Out', 'Nov', 'Dez')) # Coloca no eixo X (representado pelo 1), nos pontos 1 a 6, as etiquetas Jul a Dez. dados: X2sem2011 hgbnaopanakeia hgbpanakeia 1 JUL 75.9 75.4 2 AGO 77.5 77.1 3 SET 79.4 77.0 4 OUT 76.8 77.2 5 NOV 76.1 78.0 6 DEZ 78.3 76.7 Ana

2 1

MObilidade urbana
by Edson Lira 03 Jun '13

03 Jun '13

Caros amigos, estou tentando estabalecer uma amostra para o seguinte problema: Tenho uma origem e quatro destinos. Na origem consigo coletar 1326 placas de carros em 60 minutos, São 4 destinos e consigo identificar 620 placas desta origem, como estabelecer uma amostra nos 4 destinos para identificar o origem? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 1

Seleção em objeto tipo list
by ASANTOS 02 Jun '13

02 Jun '13

Boa noite pessoal, Uso muito pouco objetos do tipo list e tenho uma pergunta trivial, quando faço: df <- data.frame(sample=rnorm(1,0,1),params=I(list(list(mean=0,sd=1,dist="Normal")))) df <- rbind(df,data.frame(sample=rgamma(1,5,5),params=I(list(list(shape=5,rate=5,dist="Gamma"))))) df <- rbind(df,data.frame(sample=rbinom(1,7,0.7),params=I(list(list(size=7,prob=0.7,dist="Binomial"))))) df <- rbind(df,data.frame(sample=rnorm(1,2,3),params=I(list(list(mean=2,sd=3,dist="Normal"))))) df <- rbind(df,data.frame(sample=rt(1,3),params=I(list(list(df=3,dist="Student-T"))))) # df[[2]] ## Objeto 2 novo<-df[[2]]=="Binomial" ## Não funciona novo<-df[[2=="Binomial",]] ## Tambem não funciona Quero selecionar somente "Binomial" dentro de [[2]] em novo objeto chamado novo e não estou conseguindo, como faço isto? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia 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 ======================================================================

2 1