Janeiro 2015 - R-br

Pca
by Fernando Zuchello 07 Jan '15

07 Jan '15

Boa tarde, estou com uma pequena dificuldade em realizar uma pca, tenho 49 variaves nas colunas e 4 areas nas linhas, consigo realizar a pca para as variaveis, mas como separar as areas (linhas). Agradeco a ajuda.

3 2

Discretizando uma variável contínua no chisq.test
by ASANTOS 07 Jan '15

07 Jan '15

Caros Listeiros, Gostaria de comparar valores entre parcelas observadas e estimadas para saber a similaridade ou não entre dados estimados vs observados. No entanto, as variáveis observadas são dados de contagem e os dados estimados são variáveis contínuas, então minha dúvida é se violo algum pressuposto estatístico ao discretizar os dados estimados com a função round() para realizar o teste de Chi-quadrado? Segue CRM: ## Comparação entre parcelas estimadas e observadas # #Dados das parcelas observadas - Dados de contagem data.obs<-c(11, 8, 8, 8, 7, 14) # #Medias estimadas para as parcelas observadas - Dados contínuos data.est<-c(14.6,16.0, 19.1, 3.4, 9.4, 15.8) # #Distretizando as médias estimadas tr.data.est<-round(data.est) # #Comparando parcelas observadas e estimadas chisq.test(x=tr.data.est, y=data.obs, simulate.p.value=TRUE) 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 3

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 07 Jan '15

07 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