- R-br - DINF UFPR

Classificação de distâncias a partir de um raio em uma matriz de distância
by ASANTOS 07 Jan '15

07 Jan '15

Caros listeiros, Tenho um objeto com coordenadas e um identificador (da) que representa um conjunto de pontos. Calculei a matriz de distância entre todos os pontos e gostaria de a partir de um raio arbitrário, criar uma nova coluna de classificação dos pontos (em da) que estão contidos nesta distância, ou seja, quais pontos estão próximos uns dos outros em uma distância escolhida, segue CRM: #Pacotes require(spatstat) require(sp) ## Gera um processo de pontos xypois=rpoispp(25,win=owin(c(0,1),c(0,1)))## ~25 pontos aleatórios xd=xypois$x yd=xypois$y cod<-seq(1,length(yd)) da<-cbind(xd, yd, cod) ## Crio um vetor com coordenada x e y e um identificador ##Transformando em um objeto da classe espacial pontos<- SpatialPoints(cbind(da[,1],da[,2])) pontos<- SpatialPoints(list(da[,1],da[,2])) pontos<- SpatialPoints(data.frame(da[,1],da[,2])) ##Cria a matriz de distâncias res<-spDists(pontos) res Aqui eu calculei a matriz de distância, mas agora preciso de alguma função que à partir de um raio, como R=0.1, faça a classificação que eu preciso, alguém teria alguma sugestão ou conhecimento de alguma função que faça isto? 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

2 2

Homogeneidade de variancias
by Andre Oliveira 06 Jan '15

06 Jan '15

Pessoal boa tarde, dado meu modelo em fatorial. modelo=aov(respsota~factor1*factor2) Como avaliar homogeneidade de variâncias? O que fiz aqui foi ..par(mfrow=c(2,2)) boxplot(resp~fat1) boxplot(resp~fat2) boxplot(resp~fat1*fat2) interaction.plot(fat1,fat2,resp,ylab="Médias",main="Interação Fat1 vs Fat2") Mas está muito subjetivo. André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

6 7

Suavização pela mediana (median polish)
by Felinto COSTA 06 Jan '15

06 Jan '15

Prezados. Por razões didáticas gostaria de decompor a variável regionalizada "v" no exemplo abaixo em duas partes: (i) uma parte atribuída a uma tendência associada às suas coordenadas (x,y) (ii) outra que seriam seus resíduos suavizados pela mediana (median polish) library(geoR) dados.mat=matrix(nrow=8, ncol=3,byrow=F, c(486259,483200,487021,478100,482646,481735,480789,484181, 7428046,7428217,7424717,7420978,7426862,7427070,7427206,7427091, 159,350,220,265,450,110,224,530)) colnames(dados.mat)=c("x","y","v") dados.mp=medpolish(dados.mat) dados.res=as.vector(dados.mp$residuals[,3]) dados.mat.mp[,3]=dados.res dados.mat.mp geo.dados.mat.mp=as.geodata(dados.mat.mp,coords.col = 1:2,data.col=3) Minhas dúvidas são: (i) a "extração" dos resíduos está correta? (ii) numa análise geoestatística, ao se modelar os resíduos (semivariograma experimental > ajuste dos parâmetros de uma função teórica) e usar esta função para estimar os valores numa grade de pontos (krigagem ordinária) ainda deve ser "adicionada" a parte atribuída à tendência removida dos dados experimentais no começo? (ou teoricamente não há razão para isso?) Antecipadamente agradeço a colaboração e desejo a todos um ótimo 2015. Felinto

1 0

Correlação
by Gustavo Faria 06 Jan '15

06 Jan '15

Boa tarde, estou com dificuldade (iniciantes ¬¬ hehe) em fazer uma correlação com dados qualitativos, meus dados são 3 metodologias de identificação de organismos, e eu gostaria de saber o quão parecido e/ou divergente eles foram nas identificações realizadas, como poderia montar essa matriz? ja que tenho somente as identificações realizadas pelas 3 metodologias.... -- --- Gustavo de Faria Silva Biólogo herpetólogo Crbio 87421/04-P Pós-Graduado em Consultoria e Licenciamento Ambiental Mestrando em Zoologia de Vertebrados Telefone: (31) 9579-8473 (Vivo)

3 3

Tratamento de colunas
by Angélica Jaconi 06 Jan '15

06 Jan '15

Prezados, Eu tenho uma planilha ocm numeros. Essa tem o cabecalho (head) indicando as repeticoes, por exemplo: A1, A2, A3...A10, média. Gostaria de uma ideia de como tratar essas colunas iniciadas em "A" como uma única coluna para aplicar posteriores tratamento de dados. Desde ja obrigada pela atencao AJ

2 1

Letras gregas dentro do caption - xtable
by Andre Oliveira 04 Jan '15

04 Jan '15

Pessoas bom dia,estou tentando colocar, em grego, o lambda dentro do caption da tabela gerada pelo xtable, mas sem sucesso. Já olhei no ?plotmath tentei $\lambda$ ~lambda e não consegui. Alguma sugestão? <<keep.source=TRUE,results=tex>>=x=seq(0:12)p=dpois(0:12, 3)dados=data.frame(x,p)require(xtable)x=xtable(dados,caption="Distribuição Poisson ~lambda = 3",align="lll")print(x,caption.placement='top',table.placement="H",size="scriptsize",include.rownames=T, floating=TRUE)@ André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

2 2

ANOVA medidas repetidas no R
by Roberto Guevara 03 Jan '15

03 Jan '15

Caros usuários de R, minha dúvida na verdade é mais sobre estatística do que o R em si. Estou utilizando uma ANOVA de medidas repetidas com três fatores sendo que há repetição somente para dois destes fatores. O código e o resultado que estou obtendo segue o modelo abaixo: > summary(aov(dependente ~ fator1 * fator2 * fator3 + + Error(individuo/(fator2 * fator3)))) Error: individuo Df Sum Sq Mean Sq F value Pr(>F) fator1 1 80616 80616 0.260 0.645 fator2 1 809787 809787 2.611 0.205 Residuals 3 930281 310094 Error: individuo:fator2 Df Sum Sq Mean Sq F value Pr(>F) fator2 2 181649 90824 1.491 0.282 fator1:fator2 2 186811 93405 1.533 0.273 Residuals 8 487309 60914 Error: individuo:fator3 Df Sum Sq Mean Sq F value Pr(>F) fator3 2 474691 237345 14.526 0.00502 ** fator1:fator3 2 138023 69012 4.224 0.07163 . fator2:fator3 2 184239 92120 5.638 0.04189 * Residuals 6 98035 16339 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Error: individuo:fator2:fator3 Df Sum Sq Mean Sq F value Pr(>F) fator2:fator3 4 58378 14594 0.374 0.824 fator1:fator2:fator3 4 198072 49518 1.270 0.323 Residuals 16 624020 39001 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Residuals 51 2676429 52479 Todas as cinco tabelas são imprescindíveis para fins de avaliar o efeito dos fatores sobre a variável dependente? Como diferencio o que cada uma expõe? Por exemplo, o resultado da primeira tabela indica o efeito dos fatores 1 e 2 isoladamente ao retirar o efeito dos indivíduos? Qual a diferença entre os resultados expostos para o fator2 na primeira e na segunda tabela? Desde já, muito 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 1

Problema com Dates no Macbook pro
by Edson Lira 03 Jan '15

03 Jan '15

Estou com problemas na conversão de datas. Quero separar dia, mês e ano. Abaixo o link de uma amostra dos dados com as variáveis de datas Download TESTE DATA.csv | | | | | | | | | | | Download TESTE DATA.csvDownload TESTE DATA.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Link do Script Download CB.R | | | | | | | | | | | Download CB.RDownload CB.R free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Feliz 2015 para os amigos!!! Edson Lira Estatístico Manaus-Amazonas

2 2

Espaço entre texto e imagens - Sweave
by Andre Oliveira 30 Dez '14

30 Dez '14

Pessoas, boa noite. Depois de muita pesquisas eu recorro a vocês. O que posso colocar no preâmbulo de meu .rnw para reduzir os espaços entre imagens e texto? obrigado André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 7

Selecionando níveis dentro de fator
by ASANTOS 29 Dez '14

29 Dez '14

Caros listeiros, Tenho um tratamento constituído de 26 letras e gostaria de criar um vetor arbitrário para separar alguns níveis dentro deste tratamento, porém não quero utilizar a função & para cada nível (pois no meu banco de dados original ficaria muito longo), mas especificar diretamente um lista com eles sendo: samp.niv<-c("a","k","o","p","r","s","t")##Níveis escolhidos dentro de trat Porém, segundo o meu CRM se faço uma relação de igualdade dentro do banco de dados, no meu caso da, não funciona, alguém poderia me ajudar, segue CRM: ##Criando um banco artificial rm(list=ls(all=TRUE)) da <- expand.grid(trat=rep(paste(letters[1:26]),100), tempo=1:10) da$y1 <- rnorm(da$trat) samp.niv<-c("a","k","o","p","r","s","t")##Níveis escolhidos dentro de trat da.1<-da[da[,1]==samp.niv,]## Selecionando os níveis de interesse 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

2 2