Janeiro 2014 - R-br

Re: [R-br] Licenciamento
by Juan C. R. Soto Sotelo 16 Jan '14

16 Jan '14

Leonardo, você pode informar que usa o R que é um software licenciado por meio da GPL-GNU, Licença Pública Geral detalhes sobre o uso da licença podem ser lidos na página http://www.gnu.org/licenses/gpl.html Existe um portal sobre software público brasileiro http://www.softwarepublico.gov.br/spb/ListaSoftwares geralmente os softwares listados nesse portal usam a licença GPL-GNU -- Juan C. R. Soto Sotelo

1 0

Analise variância fatorial com tratamento adicional 3x2+1
by Fernando Souza 14 Jan '14

14 Jan '14

Caros amigos Estou analisando uns dados de produção de gases onde foi avaliado o efeito da adição de 3 níveis de farelo de coco em uma ração base (feno+soja+Milho) sobre a produção de metano. Foi utilizado um DIC com arranjo fatorial 3 x 2 (nível de farelo adicionado x Industria de origem do farelo). Nesse experimento a produção de metano da dieta sem adição de farelo (nível controle tambem foi avaliada. Eu quero avaliar este experimento em relação a dieta controle? Como posso fazer isso, uma vez que não tenho um fatorial para o controle? Penso que este experimento deva ser avaliado como um fatorial com tratamento adicional 3 x 2 + 1, mas não sei como fazê-lo. Alguém pode me orientar? dados<-structure(list(Amostra = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), Industria = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), Metano = c(23.82046, 22.58261, 6.90982, 9.52078, 18.04296, 24.80488, 19.28557, 13.93736, 19.13685, 8.62221, 16.07414, 9.38481, 10.82267, 17.51822, 22.07641, 23.04412, 22.8202, 26.22151, 8.506, 19.75686, 2.87987, 14.44739, 13.23224, 2.10336, 9.61595, 1.94412, 8.66299, 4.49037, 14.61952, 15.54833, 17.3914, 1.89799, 15.31825, NA, 14.09188, 4.86393, NA, 10.15676, 4.80198, 14.74247, 18.07063, 8.45602, 15.9324, 8.86496, 13.61907, 4.17938, 16.34748, 14.99311, 12.23218)), .Names = c("Amostra", "Industria", "Metano"), class = "data.frame", row.names = c(NA, -49L)) Amostra= Nível de farelo adicionado (nível 0 é o controle) Industria= (Industria de origem do farelo), há duas industria avaliada dados<-transform(dados,Amostra=factor(Amostra),Industria=factor(Industria)) ANALISE<-aov(Metano~Amostra*Industria,data=dados[-c(32,3,21),]) drop1(ANALISE,~.,test="F") shapiro.test(resid(ANALISE) --- Este email está limpo de vírus e malwares porque a proteção do avast! Antivírus está ativa. http://www.avast.com

2 5

Tinn-R: eBook publicado
by Jose Claudio Faria 14 Jan '14

14 Jan '14

Caros usuários, Foi publicado hoje pela editora Editus um eBook do Tinn-R. O eBook pode ser baixado gratuitamente no link abaixo. http://www.uesc.br/editora/index.php?item=conteudo_livros_digitais.php Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

8 10

[OFF] Erro quando se examina toda poulação
by Andre Oliveira 14 Jan '14

14 Jan '14

bom dia, qual erro pode associar quando se examina toda uma população? André Oliveira Souza

2 1

Seleção de pixels em raster em função de coordenadas não inteiras
by ASANTOS 14 Jan '14

14 Jan '14

Caros Membros Gostaria de fazer a seleção e representar graficamente os pixels de um um raster (r) que contem pontos de interesse (pontos). No entanto, o raster possui coordenadas com valores inteiros (Ex. 6, 8, 9 ..) e as coordenadas dos pontos possuem subunidades(6.03, 8.22, 9.65), gostaria de saber como posso fazer para extrair os pixels que contem os pontos selecionados e sua representação gráfica, segue CRM que não funcionou, require(raster) require(sp) ## Criando um raster r <- raster(nc=10, nr=10) r <- setValues(r, round(runif(ncell(r))* 255)) ##Pontos selecionados x<-c(4.25,3.06,6.22,7.35) y<-c(7.00,8.90,1.18,7.03) pontos <- SpatialPoints(cbind(x,y))#cbind-junta duas colunas em uma, no caso as coordenadas ### Recorte dos pontos de interesse do raster r.crop <- crop(r, pontos, snap='out') image(r.crop, asp=1) ### Recorte dos pixels nas coordenadas de interesse r.masked <- mask(r.crop, pontos) image(r.masked, asp=1) ### 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 ======================================================================

1 0

Correlação canônica
by Marcos Toebe 10 Jan '14

10 Jan '14

Prezados! Alguém saberia me informar um pacote do r com o qual eu poderia fazer a análise de correlação canônica a partir de matrizes de correlação de dimensões diferentes? Por exemplo, um grupo X com matriz de correlação simétrica de dimensão 3x3 e um grupo Y com dimensão 4x4. Teria que ser um pacote ou comando que permitisse a análise a partir dessas matrizes, pois tenho problema de multicolinearidade num dos grupos e optei por adicionar uma constante *k *na diagonal dessa matriz para corrigir esse problema. Desde já agradeço pela ajuda *----------------------* *Marcos Toebe* *Engenheiro Agrônomo - UFSM* *Mestre em Agronomia - UFSM* *Doutorando em Agronomia - UFSMProfessor de Estatística e Experimentação da* *Universidade Federal do Pampa - Itaqui - RS* m.toebe(a)gmail.com 55-96617761

1 0

Vaga para doutorado em geoestatistica (Wagningen)
by Wenceslau Teixeira 09 Jan '14

09 Jan '14

Caro Moderador, veja a viabilidade de divulgar esta nota na Lista do R. Cordialmente, Wenceslau ___________________________________ Dear colleagues, We have a PhD vacancy on ‘Bayesian geostatistics for mapping soil and agronomic variables’.Perhaps you know suitable candidates, in which case we hope that you will forward them this email. Please also feel free to forward this message to anyone who might know of suitable candidates. The application deadline is January 20. More information and application instructions at https://www.academictransfer.com/employer/WUR/vacancy/21307/lang/en/. Best wishes, Gerard Gerard B.M. Heuvelink Wageningen University and ISRIC World Soil Information P.O. Box 47 6700 AA Wageningen The Netherlands -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

1 0

Re: [R-br] Análise de agrupamento
by Marcos Vital 08 Jan '14

08 Jan '14

Natália, é possível que o problema seja a métrica que você está utilizando. Aparece um erro quando você calcula a matriz de distância usando Bray Curtis? É isso? Ou é alguma outra coisa? Não sei se o coeficiente de Bray Curtis aceita linhas com zeros. Acredito que não aceite, pois acho que uma linha de zeros não deva fazer muito sentido nesta métrica. Normalmente usamos Bray Curtis para comparar comunidades biológicas de acordo com sua similaridade de composição e abundâncias de espécies (algo como a "estrutura" da comunidade), e neste contexto uma linha só de zeros não dá informações úteis para compreendermos as relações de similaridade. Se for isto mesmo, a questão aqui é encontrar uma métrica que possa responder a sua pergunta. Aliás, este é um ponto crucial: qual a sua pergunta? Por que os zeros importam? Dependendo da situação você pode até criar uma métrica própria usando a função designdist, também do pacote vegan. Se precisar de uma boa literatura sobre estas medidas, sugiro começar pelo Numerical Ecology (Legendre & Legendre). Abraços Marcos -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

2 1

Diferença no teste de Qui quadrado com e sem arredondamento
by ASANTOS 08 Jan '14

08 Jan '14

Boa tarde Pessoal, Gostaria de comparar duas variáveis var1 e var2 pelo teste de qui-quadrado e pergunto se é incorreto ou não fazer o arrendondamento com round() da var2, se nos dois casos os resultados são ligeiramente semelhantes, segue CRM abaixo: # Sem arredondamento var1<-c(11, 8, 8, 8, 7, 14) var2<-c(14.550682, 15.971429, 19.086132 , 3.403902 , 9.441030, 15.817173) teste1<-rbind(var1,var2) (Xsq <- chisq.test(teste1,correct=TRUE)) # # Com arredondamento var1<-c(11, 8, 8, 8, 7, 14) var2<-c(14.550682, 15.971429, 19.086132 , 3.403902 , 9.441030, 15.817173) teste2<-rbind(var1,round(var2)) (Xsq <- chisq.test(teste2,correct=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 ======================================================================

1 0

[OFF-TOPIC] Qual a diferença entre duas abordagens espaciais?
by ASANTOS 08 Jan '14

08 Jan '14

Boa tarde Pessoal, Estou estudando algumas abordagens utilizando GLM para a predição espacial e me deparei com uma dúvida, qual a diferença entre o ajuste de modelos GLM com o pacote geoRglm preconizado pelo Prof Paulo e o ajuste de modelos GLM com a função glm() e posterior uso dos resíduos no ajuste de semivariogramas como o pacote gstat preconizado pelo Prof. Tomislav Hengl, qundo o objetivo é a predição espacial da densidade de uma variável tendo variáveis explicativas a altitude, temperatura, etc, 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 ======================================================================

2 2