Junho 2011 - R-br - DINF UFPR

Re: [R-br] Instalação do pacote rparallel
by Junior Beleti 02 Jun '11

02 Jun '11

Olá Benilton, nada especifico, mas pelo que eu vi sobre ele, me parece uma boa solução para aplicações paralelas pelo fato de "automatizar" o paralelismo. Porém não encontrei nenhum exemplo de sua utilização. O snow voce utiliza em maquina multicore mesmo ou cluster? Vou dar uma olhada nos outros pacotes. Obrigado -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

RES: RES: Download de dados pelo R
by Humberto Hazin (hotmail) 02 Jun '11

02 Jun '11

Olá Victor, Obrigado pela ajuda! Mais o problema não é esse! Não sei porque ele não esta reconhecendo meu diretório, porem nesse exemplo ele vai, só que tenho que colocar arquivo por arquivo! cho<-"http://aspera.jpl.nasa.gov/download/pub/sea_surface_temperature/avhrr/ pathfinder/data_v5/monthly/night/04km/1986/bsst/198601.s04m1pfv50-bsst-16b.h df" download.file(cho, destfile = "~/Dropbox/Humberto/Analises em R/ambientais/.hdf", method="auto", quiet = FALSE, mode = "w",cacheOK = FALSE) se alguem puder souber como resolver isso agradeço abs, Humberto De: r-br-bounces(a)listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Victor Eduardo Enviada em: Wednesday, June 01, 2011 9:47 AM Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] RES: Download de dados pelo R Outra coisa, estava olhando o seu erro e repare o seguinte: C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf' No finalzinho, tem um espaço entre 2001 e .hdf. De repente pode ser isso, confere direitinho o seu paste todo e ve se está encaminhando certo para o diretório. Em 1 de junho de 2011 09:44, Victor Eduardo <victorduca08(a)gmail.com> escreveu: Humberto, eu não consegui fazer o download aqui, de repente pode ser problema pelo site mesmo. Em 1 de junho de 2011 09:04, Humberto Hazin (hotmail) <hghazin(a)hotmail.com> escreveu: Olá Eduardo Ele deve estar fora de ar nesse instante ! Tenta com esse aqui http://aspera.jpl.nasa.gov/download/pub/sea_surface_temperature/avhrr/pathfi nder/data_v5/monthly/night/04km/1986/bsst/ ou http://oceandata.sci.gsfc.nasa.gov/SeaWiFS/Mapped/Monthly/9km/chlor/ Abs, Humberto De: r-br-bounces(a)listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Victor Eduardo Enviada em: Wednesday, June 01, 2011 8:58 AM Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Download de dados pelo R Humbero, não está dando para acessar o link enviado. Em 1 de junho de 2011 08:54, Humberto Hazin (hotmail) <hghazin(a)hotmail.com> escreveu: Bom dia pessoal! Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou conseguindo, alguém poderia me ajudar? ano<-2001 cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", method="auto") Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", : cannot open destfile 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such file or directory' Abs, Humberto _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br

4 8

Re: [R-br] Instalação do pacote rparallel
by Junior Beleti 02 Jun '11

02 Jun '11

Olá Ivan. Este pacote não é um pacote oficial, mas encontrei artigos sobre programação paralela com ele. Estou desenvolvendo minha pesquisa sobre pacotes paralelos no R. Conhece algum pacote que possa ser utilizado? Tentei utilizar o pnmath mas não tive sucesso. Até mais. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

graficos com dois eixos y
by Gilson Sanchez 01 Jun '11

01 Jun '11

Pessoal mais uma vez procuro ajuda de vc's, estou querendo montar um gráfico com dois eixos y. Então o seguinte como eu posso fazer o seguinte: 1. Colocar nomes dos anos no eixo X (jan, fev, mar, abr...dez). 2. oC quero colocar "o" sobreescrita. 3. Como coloco o nome do segundo eixo y. plot(dados$temp, type="b", pch=15, lwd=1.5, xlab="Meses (Anos 2010)", ylab="Temperatura (oC)", ylim=c(24,30), xlim=c(1,12), main= "A", font.main=2, font.lab=2) par(new=T) plot(dados$umidade, type="b", pch=14,axes=F, frame=T, ann=F) axis(4, mtext(text="Umidade (%))", 4, line = 3)) dados meses temp umidade rad chuva 1 Jan 26.49 86.58 795.88 0.36 2 Fev 26.65 88.49 710.24 0.34 3 Mar 27.19 86.16 772.99 0.21 4 Abr 26.28 89.75 574.88 0.67 5 Mai 26.62 89.22 614.02 0.31 6 Jun 26.13 87.83 680.08 0.26 7 Jul 25.83 86.57 675.97 0.15 8 Ago 27.05 83.14 756.44 0.07 9 Set 27.60 83.02 925.57 0.14 10 Out 27.44 85.16 927.71 0.17 11 Nov 26.56 88.18 788.87 0.19 12 Dez 25.87 90.63 703.94 0.33 Obrigado se puderem ajudar -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

2 4

Download de dados pelo R
by Humberto Hazin (hotmail) 01 Jun '11

01 Jun '11

Bom dia pessoal! Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou conseguindo, alguém poderia me ajudar? ano<-2001 cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", method="auto") Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", : cannot open destfile 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such file or directory' Abs, Humberto

3 5

modifica dataframe
by Samuel Carvalho 01 Jun '11

01 Jun '11

Boa noite caros(a) senhores(as) membros Por gentileza alguém poderia me auxliar em um código que a princípio parece ser simples mas nao consegui uma solução Abaixo mostro o codigo e a seguir venho comentando ### Codigo R dados<-data.frame(idfuste=seq(1,5), dap=c(8.5, 11.0, 10.5, 17.5, 12.5), ht=c(10, 13, 12.5, 18, 14)) dados dados$hi<-seq(0.1, dados$ht, by=0.5) ### O que preciso é gerar uma sequencia para cada "idfuste" de 0.1 até a "ht" que corresponde aquele "idfuste", replicando também cada observação para cada "idfuste". A sequencia é de 0.5 em 0.5 até a altura total da árvore - 1 pois esta pode nao ser um número inteiro dividido por 0.5 sendo que a ultima observação de cada idfuste deve ser igual ao valor de ht Em resumo queria um data.frame de saída a partir dos dados de entrada, nesta estrutura idfuste dap ht hi 1 8,5 10 0 1 8,5 10 0,5 1 8,5 10 1 1 8,5 10 1,5 1 8,5 10 2 1 8,5 10 2,5 1 8,5 10 3 1 8,5 10 3,5 1 8,5 10 4 1 8,5 10 4,5 1 8,5 10 5 1 8,5 10 5,5 1 8,5 10 6 1 8,5 10 6,5 1 8,5 10 7 1 8,5 10 7,5 1 8,5 10 8 1 8,5 10 8,5 1 8,5 10 9 1 8,5 10 9,5 1 8,5 10 10 2 11 13,2 0 2 11 13,2 0,5 2 11 13,2 1 2 11 13,2 1,5 2 11 13,2 2 2 11 13,2 2,5 2 11 13,2 3 2 11 13,2 3,5 2 11 13,2 4 2 11 13,2 4,5 2 11 13,2 5 2 11 13,2 5,5 2 11 13,2 6 2 11 13,2 6,5 2 11 13,2 7 2 11 13,2 7,5 2 11 13,2 8 2 11 13,2 8,5 2 11 13,2 9 2 11 13,2 9,5 2 11 13,2 10 2 11 13,2 10,5 2 11 13,2 11 2 11 13,2 11,5 2 11 13,2 12 2 11 13,2 12,5 2 11 13,2 13 2 11 13,2 13,2 e assim por diante para todos idfustes Desde já agradeço Para quem é da area florestal estou tentando estruturar uma base de cubagem a partir de dados de parcelas Abraço e boa noite a todos Att, Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

6 8

Dúvida em argumentos do histograma
by Victor Eduardo 01 Jun '11

01 Jun '11

Pessoal, estou com um loop e dentro do loop estou querendo utilizar o seguinte para a construção do histograma: hist(paste(var1,'_',var,sep=''),main=paste('tit.',var1,sep=''),xlim=c(0,10),ylim=c(0,85)) sendo que está dando o seguinte erro: Erro em hist.default(paste(var1, "_", var, sep = ""), main = paste("tit.", : 'x' deve ser numérico O primeiro paste me daria, por exemplo, o seguinte: obs_clcl Com isso, construiria o histograma para essa variável. Como eu poderia resolver esse problema? Abraços, Victor Eduardo

2 2

Re: [R-br] graficos com dois eixos y
by Ivan Bezerra Allaman 01 Jun '11

01 Jun '11

Gilson, Repare que no script que te mandei, antes do plot eu utilizo o parâmetro "par(mar=c(0,0,0,0))". Este parâmetro serve para alterar as margens do gráfico. Vá alterando um por um para descobrir a funcionalidade. Teste 5,0,0,0; depois 0,5,0,0 e assim por diante. Com certeza achará um valor adequado para que o nome do eixo y apareça. Pelo o código que te mandei o nome estava aparecendo, mais isso deve depender de máquina pra máquina. As vezes a margem que vc deve utilizar dentro do "mar" deve ser maior do que aqueles que estão no exemplo que te mandei. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] graficos com dois eixos y
by Ivan Bezerra Allaman 01 Jun '11

01 Jun '11

Enquanto o Walmes não responde, vai quebrando o galho com essa solução abaixo. lines <- 'meses temp umidade rad chuva Jan 26.49 86.58 795.88 0.36 Fev 26.65 88.49 710.24 0.34 Mar 27.19 86.16 772.99 0.21 Abr 26.28 89.75 574.88 0.67 Mai 26.62 89.22 614.02 0.31 Jun 26.13 87.83 680.08 0.26 Jul 25.83 86.57 675.97 0.15 Ago 27.05 83.14 756.44 0.07 Set 27.60 83.02 925.57 0.14 Out 27.44 85.16 927.71 0.17 Nov 26.56 88.18 788.87 0.19 Dez 25.87 90.63 703.94 0.33' da <- read.table(textConnection(lines), header=TRUE) par(mar=c(5,4,4,5)) plot(da$temp, type="b", pch=15, lwd=1.5, xlab="Meses (Anos 2010)", ylab="Temperatura (ºC)", ylim=c(24,30), xlim=c(1,12), axes=F, main= "A", font.main=2, font.lab=2) axis(1, at=1:12,labels=c("Jan.","Fev.","Mar.","Abr.","Maio","Jun.", "Jul.", "Ago.", "Set.", "Out.", "Nov.", "Dez.")) axis(2) par(new=T) plot(da$umidade, type="b", pch=14,axes=F, frame=T, ann=F) axis(4) mtext(text="Umidade (%)", 4, line = 3) Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

RStudio
by Diogo Borges Provete 01 Jun '11

01 Jun '11

Olá, alguém que usa o RStudio está conseguindo exibir as mensagens de erro etc em Português? Abraço Atenciosamente, Diogo Borges Provete ============================== Biólogo Mestre em Biologia Animal (UNESP) Doutorando PPG Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás - UFG Goiânia-GO CP: 131 74001-970 Brazil Tel. +55 62 3521-1732 Cel. +55 62 8231-5775 Skype: diogoprovete MSN: diogoprov(a)yahoo.com.br Personal web page Traduza conosco: <<<American Journal Experts>>> <<<D-Lang Soluções linguisticas>>> <<<Perfil no ProZ.com como tradutor freelancer>>> ============================== ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quarta-feira, 1 de Junho de 2011 9:04 Assunto: Digest R-br, volume 4, assunto 1 Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Re: modifica dataframe (Gustavo Henrique de Carvalho) 2. Re: teste de Homogeneidade e Normalidade (Ivan Bezerra Allaman) 3. Re: teste de Homogeneidade e Normalidade (Eduardo Vargas) 4. Re: teste de Homogeneidade e Normalidade (Mauro Sznelwar) 5. Re: teste de Homogeneidade e Normalidade (Eduardo Vargas) 6. Re: teste de Homogeneidade e Normalidade (Mauro Sznelwar) 7. Re: teste de Homogeneidade e Normalidade (Walmes Zeviani) 8. Re: modifica dataframe (Henrique Dallazuanna) 9. Re: teste de Homogeneidade e Normalidade (Eduardo Vargas) 10. Re: modifica dataframe (Samuel Carvalho) 11. Re: teste de Homogeneidade e Normalidade (Walmes Zeviani) 12. Re: teste de Homogeneidade e Normalidade (Eder David Borges da Silva) 13. Re: teste de Homogeneidade e Normalidade (Eduardo Vargas) 14. Download de dados pelo R (Humberto Hazin (hotmail)) 15. Re: Download de dados pelo R (Victor Eduardo) 16. RES: Download de dados pelo R (Humberto Hazin (hotmail)) ---------------------------------------------------------------------- Message: 1 Date: Tue, 31 May 2011 21:24:23 -0300 From: Gustavo Henrique de Carvalho <gustavo.bio(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Cc: Samuel Carvalho <samukajm(a)yahoo.com.br> Subject: Re: [R-br] modifica dataframe Message-ID: <BANLkTikLCzKcVdHqgWWUSv8f+Z8UMayxKA(a)mail.gmail.com> Content-Type: text/plain; charset=UTF-8 Não sei pra que colar esse negócio grotesco aqui, mas a primeira forma que me veio à cabeça foi essa: bleh <- function(x) {unique(c(seq(0, x, by = 0.5), x))} data.frame(lapply(dados, rep, unlist(lapply(sapply(dados$ht, bleh), length))), hi = unlist(sapply(dados$ht, bleh))) 2011/5/31 Benilton Carvalho <beniltoncarvalho(a)gmail.com> > > enqto o henrique nao responde, algo como o seguinte pode lhe servir: > hi = mapply(seq, 0.1, dados[['ht']], .5) > idx = mapply(rep, 1:length(hi), sapply(hi, length)) > dados2 = cbind(dados[unlist(idx),], hi=unlist(hi)) > b > > 2011/5/31 Samuel Carvalho <samukajm(a)yahoo.com.br> >> >> Boa noite caros(a) senhores(as) membros >> Por gentileza alguém poderia me auxliar em um código que a princípio parece ser simples mas nao consegui uma solução >> Abaixo mostro o codigo e a seguir venho comentando >> ### Codigo R >> dados<-data.frame(idfuste=seq(1,5), dap=c(8.5, 11.0, 10.5, 17.5, 12.5), >> ht=c(10, 13, 12.5, 18, 14)) >> dados >> dados$hi<-seq(0.1, dados$ht, by=0.5) >> ### >> O que preciso é gerar uma sequencia para cada "idfuste" de 0.1 até a "ht" que corresponde aquele "idfuste", replicando também cada observação >> para cada "idfuste". A sequencia é de 0.5 em 0.5 até a altura total da árvore - 1 pois esta pode nao ser um número inteiro dividido por 0.5 sendo que a ultima observação de cada idfuste deve ser igual ao valor de ht >> Em resumo queria um data.frame de saída a partir dos dados de entrada, nesta estrutura >> idfuste dap ht hi >> 1 8,5 10 0 >> 1 8,5 10 0,5 >> 1 8,5 10 1 >> 1 8,5 10 1,5 >> 1 8,5 10 2 >> 1 8,5 10 2,5 >> 1 8,5 10 3 >> 1 8,5 10 3,5 >> 1 8,5 10 4 >> 1 8,5 10 4,5 >> 1 8,5 10 5 >> 1 8,5 10 5,5 >> 1 8,5 10 6 >> 1 8,5 10 6,5 >> 1 8,5 10 7 >> 1 8,5 10 7,5 >> 1 8,5 10 8 >> 1 8,5 10 8,5 >> 1 8,5 10 9 >> 1 8,5 10 9,5 >> 1 8,5 10 10 >> 2 11 13,2 0 >> 2 11 13,2 0,5 >> 2 11 13,2 1 >> 2 11 13,2 1,5 >> 2 11 13,2 2 >> 2 11 13,2 2,5 >> 2 11 13,2 3 >> 2 11 13,2 3,5 >> 2 11 13,2 4 >> 2 11 13,2 4,5 >> 2 11 13,2 5 >> 2 11 13,2 5,5 >> 2 11 13,2 6 >> 2 11 13,2 6,5 >> 2 11 13,2 7 >> 2 11 13,2 7,5 >> 2 11 13,2 8 >> 2 11 13,2 8,5 >> 2 11 13,2 9 >> 2 11 13,2 9,5 >> 2 11 13,2 10 >> 2 11 13,2 10,5 >> 2 11 13,2 11 >> 2 11 13,2 11,5 >> 2 11 13,2 12 >> 2 11 13,2 12,5 >> 2 11 13,2 13 >> 2 11 13,2 13,2 >> e assim por diante para todos idfustes >> Desde já agradeço >> Para quem é da area florestal estou tentando estruturar uma base de cubagem a partir de dados de parcelas >> Abraço e boa noite a todos >> Att, >> Samuel >> >> ==================================== >> Samuel P. C. Carvalho >> Mestre em Ciências Florestais [UFLA] >> Doutorando em Recursos Florestais [ESALQ/USP] >> ============================================= >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> > > > > -- > Successful people ask better questions, and as a result, they get better answers. (Tony Robbins) > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ------------------------------ Message: 2 Date: Tue, 31 May 2011 17:34:14 -0700 (PDT) From: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> To: R Brasil <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <947015.46346.qm(a)web161820.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Se você quer transformar as variáveis que já estão em porcentagem, deve utilizar arcsen(x%)^0,5. Agora, se as variáveis estiverem entre 0 e 1 (em decimal), aí sim vc deve utilizar arcsen(x/100)^0,5. No entanto, como a variável "germinação" é do tipo Bernouli, aconselho-o a utilizar uma análise mais rebuscada como modelos lineares generalizados por exemplo. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/ef67c401/atta…> ------------------------------ Message: 3 Date: Tue, 31 May 2011 21:46:46 -0300 From: Eduardo Vargas <duvargas(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br, Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTinXDFS=T7UD7GrEKSxy-0z1-i77FA(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Ok eu entendi as transformações com relação quando já estão em porcentagem.Quando as variáveis por exemplo,não estiverem entre ) e 1, vc pode me informar que transformação eu utilizo?Por que pelo que eu entendi,as transformações com arcsen(x/110)^0.5 só entram nesse intervalo,entre 0 e 1. Em 31 de maio de 2011 21:34, Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br>escreveu: > Se você quer transformar as variáveis que já estão em porcentagem, deve > utilizar arcsen(x%)^0,5. Agora, se as variáveis estiverem entre 0 e 1 (em > decimal), aí sim vc deve utilizar arcsen(x/100)^0,5. > > No entanto, como a variável "germinação" é do tipo Bernouli, aconselho-o a > utilizar uma análise mais rebuscada como modelos lineares generalizados por > exemplo. > > Allaman > (S,f,P) > > *M.Sc Ivan Bezerra Allaman* > Zootecnista > Doutorando em Produção Animal/Aquicultura - UFLA > **email e msn - ivanalaman(a)yahoo.com.br > Tel: (35)3826-6608/9900-2924 > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/7c9a93de/atta…> ------------------------------ Message: 4 Date: Tue, 31 May 2011 22:11:34 -0300 From: "Mauro Sznelwar" <sznelwar(a)uol.com.br> To: <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <00f101cc1ff8$ddcd1fe0$6700a8c0@mauro> Content-Type: text/plain; charset="iso-8859-1" Onde eu pego os arquivos para rodar? Mauro José, Sem o seu codigo fica meio complicado, porem neste link tem uma serie de analises de estatistica experimental com onde se testa homogenidade e normalidade. http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R att Em 31 de maio de 2011 20:13, Eduardo Vargas <duvargas(a)gmail.com> escreveu: Boa noite pessoal,preciso de uma ajuda.Estou fazendo uma comparação de médias de germinação de sementes, e as varianveis são tratamentos(trat) fator e Germinação (germ) numerica.Preciso fazer o teste de homogeneidade de variancias com o Bartlett e normalidade com o teste de shapiro-Wilk, porém dá um mensagem de erro(Erro em bartlett.test(germ, trat) : objeto 'germ' não encontrado) e no Shapiro -wilk essa mensagem Erro em `[.data.frame`(x, complete.cases(x)) : undefined columns selected Será que alguem pode me ajudar? -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------------------------------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/ae8f9462/atta…> ------------------------------ Message: 5 Date: Tue, 31 May 2011 22:16:25 -0300 From: Eduardo Vargas <duvargas(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTin=XzAn+6AVxUAnEyeNqDCZQ1Ykug(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" São meus resultados de trabalho,preciso fazer essas transformações para normalizar. Em 31 de maio de 2011 22:11, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: > Onde eu pego os arquivos para rodar? > Mauro > > José, > Sem o seu codigo fica meio complicado, porem neste link tem uma serie de > analises de estatistica experimental com onde se testa homogenidade e > normalidade. > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > att > > Em 31 de maio de 2011 20:13, Eduardo Vargas <duvargas(a)gmail.com>escreveu: > >> >> Boa noite pessoal,preciso de uma ajuda.Estou fazendo uma comparação de >> médias de germinação de sementes, e as varianveis são tratamentos(trat) >> fator e Germinação (germ) numerica.Preciso fazer o teste de homogeneidade de >> variancias com o Bartlett e normalidade com o teste de shapiro-Wilk, >> porém dá um mensagem de erro(Erro em bartlett.test(germ, trat) : objeto >> 'germ' não encontrado) e no Shapiro -wilk essa mensagem Erro em >> `[.data.frame`(x, complete.cases(x)) : >> undefined columns selected >> Será que alguem pode me ajudar? >> -- >> José Eduardo Vargas Lopes de Aráujo >> Engenheiro Agrônomo >> Contato:(38)8811-6985 >> duvargas(a)gmail.com >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> >> > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/74c6398a/atta…> ------------------------------ Message: 6 Date: Tue, 31 May 2011 22:20:04 -0300 From: "Mauro Sznelwar" <sznelwar(a)uol.com.br> To: <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <002201cc1ffa$0decce90$6700a8c0@mauro> Content-Type: text/plain; charset="iso-8859-1" Eu me referi aos arquivos da página dos Exercícios do R no leg, lá tem exemplos do código mas não diz onde estão os arquivos dos dados que foram usados! São meus resultados de trabalho,preciso fazer essas transformações para normalizar. Em 31 de maio de 2011 22:11, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: Onde eu pego os arquivos para rodar? Mauro -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/c7a2ed12/atta…> ------------------------------ Message: 7 Date: Tue, 31 May 2011 22:53:19 -0300 From: Walmes Zeviani <walmeszeviani(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTikem0jdfVDhKZ8N-n8b83QyNkgXvw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Eduardo, Salve engano da minha parte, essa transformação é estabilizadora da variância, ou seja, o foco é corrigir heterocedasticidade e não olha para normalidade dos dados. Uma vez disse em um curso que essas transformações eram coisas do passado e fui por alguns instantes odiado (por aqueles mais tradicionais obviamente). Depois passou. Acontece que essa transformação é a recomendada quando seu dado é do tipo binomial (ou seja, já se sabe que o dado não é normal), dai parte algumas álgebras, integrais, até que se obtém essa função estabilizadora da variância. Transformação estabilizadora da variãncia existe para distribuição binomial e Poisson. Não lembro das funções porque nunca usei. Considero isso metodologia do passo porque, se você sabe que seu dado é binomial (ou Poisson), hoje já existe implementado nos aplicativos estatísticos modelos capazes de considerar essa característica inerente do dado. Esse conjunto de métodos, enfim, se chama modelos lineares generalizados, do qual a distribuição normal é um caso particular (com propriedades bem interessantes). Com essa metodologia você pode obter todos os resultados experimentais que obteria com a distribuição normal, diga-se a dobradinha anova e teste de médias (com algumas adaptações). Então, não há complicações e motivos para não se usar modelos lineares generalizados. Ainda é possível, caso o número de sementes colocadas para germinar seja grande (n, normalmente é 25 ou 50), e caso a probalidade de germinação não esteja na borda (próximo de 0 ou 1), de assumir que o dado é normal e analisa-lo assim. Porém, caso p mude muito entre tratamentos pode as variâncias amostrais serem muito discrepântes. Você pode tentar uma transformação Box-Cox. Verificar os resíduos. Em termos de funções do R, você pode usar a glm(, family=binomial), fazer a análise de resíduos usuais, obter teste (sequencial) para os efeitos fixos via anova() (que não é anova mas quadro de análise de deviance), se for comparar "médias" você pode usar a multcomp::glht() e contrast::contrast(). Verifique a documentação dessas funções para saber como usa-las. Para um começo, você pode rodar o CMR disponível nas Ridículas do LEG. http://www.leg.ufpr.br/doku.php/ridiculas?&#analise_de_dados_de_proporcao_u… À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/849b27e7/atta…> ------------------------------ Message: 8 Date: Tue, 31 May 2011 23:10:54 -0300 From: Henrique Dallazuanna <wwwhsd(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br, Samuel Carvalho <samukajm(a)yahoo.com.br> Subject: Re: [R-br] modifica dataframe Message-ID: <BANLkTi=77BMKw2WfJopo7xjR--he5P03VQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Samuel, Tente assim tbm: x <- sapply(dados$ht, seq, from = 0, by = 0.5) transform(dados[rep(seq(nrow(dados)), times = sapply(x, length)),], hi = unlist(x)) 2011/5/31 Samuel Carvalho <samukajm(a)yahoo.com.br> > Boa noite caros(a) senhores(as) membros > Por gentileza alguém poderia me auxliar em um código que a princípio parece > ser simples mas nao consegui uma solução > Abaixo mostro o codigo e a seguir venho comentando > ### Codigo R > dados<-data.frame(idfuste=seq(1,5), dap=c(8.5, 11.0, 10.5, 17.5, 12.5), > ht=c(10, 13, 12.5, 18, 14)) > dados > dados$hi<-seq(0.1, dados$ht, by=0.5) > ### > O que preciso é gerar uma sequencia para cada "idfuste" de 0.1 até a "ht" > que corresponde aquele "idfuste", replicando também cada observação > para cada "idfuste". A sequencia é de 0.5 em 0.5 até a altura total da > árvore - 1 pois esta pode nao ser um número inteiro dividido por 0.5 sendo > que a ultima observação de cada idfuste deve ser igual ao valor de ht > Em resumo queria um data.frame de saída a partir dos dados de entrada, > nesta estrutura > > idfuste dap ht hi 1 8,5 10 0 1 8,5 10 0,5 1 8,5 10 1 1 8,5 10 1,5 1 > 8,5 10 2 1 8,5 10 2,5 1 8,5 10 3 1 8,5 10 3,5 1 8,5 10 4 1 8,5 10 4,5 > 1 8,5 10 5 1 8,5 10 5,5 1 8,5 10 6 1 8,5 10 6,5 1 8,5 10 7 1 8,5 10 > 7,5 1 8,5 10 8 1 8,5 10 8,5 1 8,5 10 9 1 8,5 10 9,5 1 8,5 10 10 2 11 > 13,2 0 2 11 13,2 0,5 2 11 13,2 1 2 11 13,2 1,5 2 11 13,2 2 2 11 13,2 > 2,5 2 11 13,2 3 2 11 13,2 3,5 2 11 13,2 4 2 11 13,2 4,5 2 11 13,2 5 > 2 11 13,2 5,5 2 11 13,2 6 2 11 13,2 6,5 2 11 13,2 7 2 11 13,2 7,5 2 > 11 13,2 8 2 11 13,2 8,5 2 11 13,2 9 2 11 13,2 9,5 2 11 13,2 10 2 11 > 13,2 10,5 2 11 13,2 11 2 11 13,2 11,5 2 11 13,2 12 2 11 13,2 12,5 2 > 11 13,2 13 2 11 13,2 13,2 > > e assim por diante para todos idfustes > Desde já agradeço > Para quem é da area florestal estou tentando estruturar uma base de cubagem > a partir de dados de parcelas > Abraço e boa noite a todos > Att, > Samuel > > *====================================* > *Samuel P. C. Carvalho > *Mestre em Ciências Florestais [UFLA] > Doutorando em Recursos Florestais [ESALQ/USP] > ============================================= > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/4937eeab/atta…> ------------------------------ Message: 9 Date: Tue, 31 May 2011 23:34:58 -0300 From: Eduardo Vargas <duvargas(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTiniir0GPKzWcMm_=PsYwrBfmeJTjA(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Walmes,obrigado pela sua resposta,foi bastante pertinente,e acho que para eu fazer esse modelo linear generalizado é bastante complexo.Acho que meu experimento é bastante simples,fiz um experimento com 11 tratamentos + 0 o controle,com 25 sementes cada com 4 repetições.Quero comparar as médias dos tratamentos com a média do controle,Acredito eu,usando Dunnet.Meu orientador pediu que eu fizesse esse teste de normalidade(Shapiro-wilk) e homogeneidade de variâncias(bartlett).As transformações,que acho serem necessarias,eu andei olhando nos artigos que usei como referencia e todos fizeram as transformações.Mais uma vez obrigado pela ajuda,estou meio perdido. Em 31 de maio de 2011 22:53, Walmes Zeviani <walmeszeviani(a)gmail.com>escreveu: > Eduardo, > > Salve engano da minha parte, essa transformação é estabilizadora da > variância, ou seja, o foco é corrigir heterocedasticidade e não olha para > normalidade dos dados. Uma vez disse em um curso que essas transformações > eram coisas do passado e fui por alguns instantes odiado (por aqueles mais > tradicionais obviamente). Depois passou. > > Acontece que essa transformação é a recomendada quando seu dado é do tipo > binomial (ou seja, já se sabe que o dado não é normal), dai parte algumas > álgebras, integrais, até que se obtém essa função estabilizadora da > variância. Transformação estabilizadora da variãncia existe para > distribuição binomial e Poisson. Não lembro das funções porque nunca usei. > > Considero isso metodologia do passo porque, se você sabe que seu dado é > binomial (ou Poisson), hoje já existe implementado nos aplicativos > estatísticos modelos capazes de considerar essa característica inerente do > dado. Esse conjunto de métodos, enfim, se chama modelos lineares > generalizados, do qual a distribuição normal é um caso particular (com > propriedades bem interessantes). > > Com essa metodologia você pode obter todos os resultados experimentais que > obteria com a distribuição normal, diga-se a dobradinha anova e teste de > médias (com algumas adaptações). Então, não há complicações e motivos para > não se usar modelos lineares generalizados. > > Ainda é possível, caso o número de sementes colocadas para germinar seja > grande (n, normalmente é 25 ou 50), e caso a probalidade de germinação não > esteja na borda (próximo de 0 ou 1), de assumir que o dado é normal e > analisa-lo assim. Porém, caso p mude muito entre tratamentos pode as > variâncias amostrais serem muito discrepântes. Você pode tentar uma > transformação Box-Cox. Verificar os resíduos. > > Em termos de funções do R, você pode usar a glm(, family=binomial), fazer a > análise de resíduos usuais, obter teste (sequencial) para os efeitos fixos > via anova() (que não é anova mas quadro de análise de deviance), se for > comparar "médias" você pode usar a multcomp::glht() e contrast::contrast(). > Verifique a documentação dessas funções para saber como usa-las. > > Para um começo, você pode rodar o CMR disponível nas Ridículas do LEG. > > http://www.leg.ufpr.br/doku.php/ridiculas?&#analise_de_dados_de_proporcao_u… > > À disposição. > Walmes. > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/43f2d612/atta…> ------------------------------ Message: 10 Date: Tue, 31 May 2011 20:02:08 -0700 (PDT) From: Samuel Carvalho <samukajm(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] modifica dataframe Message-ID: <369328.15725.qm(a)web161206.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Obrigado caros pela gentileza de sempre. Vou verificar as sugestões []'s Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] ============================================= ________________________________ De: Henrique Dallazuanna <wwwhsd(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br; Samuel Carvalho <samukajm(a)yahoo.com.br> Enviadas: Terça-feira, 31 de Maio de 2011 23:10 Assunto: Re: [R-br] modifica dataframe Samuel, Tente assim tbm: x <- sapply(dados$ht, seq, from = 0, by = 0.5) transform(dados[rep(seq(nrow(dados)), times = sapply(x, length)),], hi = unlist(x)) 2011/5/31 Samuel Carvalho <samukajm(a)yahoo.com.br> Boa noite caros(a) senhores(as) membros >Por gentileza alguém poderia me auxliar em um código que a princípio parece ser simples mas nao consegui uma solução >Abaixo mostro o codigo e a seguir venho comentando >### Codigo R >dados<-data.frame(idfuste=seq(1,5), dap=c(8.5, 11.0, 10.5, 17.5, 12.5), >ht=c(10, 13, 12.5, 18, 14)) >dados >dados$hi<-seq(0.1, dados$ht, by=0.5) >### > >O que preciso é gerar uma sequencia para cada "idfuste" de 0.1 até a "ht" que corresponde aquele "idfuste", replicando também cada observação >para cada "idfuste". A sequencia é de 0.5 em 0.5 até a altura total da árvore - 1 pois esta pode nao ser um número inteiro dividido por 0.5 sendo que a ultima observação de cada idfuste deve ser igual ao valor de ht > >Em resumo queria um data.frame de saída a partir dos dados de entrada, nesta estrutura > > >idfuste dap ht hi >1 8,5 10 0 >1 8,5 10 0,5 >1 8,5 10 1 >1 8,5 10 1,5 >1 8,5 10 2 >1 8,5 10 2,5 >1 8,5 10 3 >1 8,5 10 3,5 >1 8,5 10 4 >1 8,5 10 4,5 >1 8,5 10 5 >1 8,5 10 5,5 >1 8,5 10 6 >1 8,5 10 6,5 >1 8,5 10 7 >1 8,5 10 7,5 >1 8,5 10 8 >1 8,5 10 8,5 >1 8,5 10 9 >1 8,5 10 9,5 >1 8,5 10 10 >2 11 13,2 0 >2 11 13,2 0,5 >2 11 13,2 1 >2 11 13,2 1,5 >2 11 13,2 2 >2 11 13,2 2,5 >2 11 13,2 3 >2 11 13,2 3,5 >2 11 13,2 4 >2 11 13,2 4,5 >2 11 13,2 5 >2 11 13,2 5,5 >2 11 13,2 6 >2 11 13,2 6,5 >2 11 13,2 7 >2 11 13,2 7,5 >2 11 13,2 8 >2 11 13,2 8,5 >2 11 13,2 9 >2 11 13,2 9,5 >2 11 13,2 10 >2 11 13,2 10,5 >2 11 13,2 11 >2 11 13,2 11,5 >2 11 13,2 12 >2 11 13,2 12,5 >2 11 13,2 13 >2 11 13,2 13,2 > >e assim por diante para todos idfustes >Desde já agradeço >Para quem é da area florestal estou tentando estruturar uma base de cubagem a partir de dados de parcelas >Abraço e boa noite a todos >Att, >Samuel > >==================================== >Samuel P. C. Carvalho >Mestre em Ciências Florestais [UFLA] >Doutorando em Recursos Florestais [ESALQ/USP] >============================================= >_______________________________________________ >R-br mailing list >R-br(a)listas.c3sl.ufpr.br >https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40" S 49° 16' 22" O -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/7d0c4086/atta…> ------------------------------ Message: 11 Date: Wed, 1 Jun 2011 00:52:54 -0300 From: Walmes Zeviani <walmeszeviani(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTimC0Ux64PEz2DChLzsedsWigrmoTQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Eduardo, O código abaixo mostra como não é difícil analisar os seus dados pela abordagem que discuti #------------------------------------------------------------------------------------------ # experimento com 11 tratamentos e 1 controle, resposta binomial n=25, 4 repetições trat <- c("controle", paste("trat", 1:11, sep="")) da <- expand.grid(trat=trat, rept=1:4) da$germinadas <- rbinom(nrow(da), size=25, prob=0.8) levels(da$trat) # controle é o nível de referência # ajuste um modelo linear generalizado g0 <- glm(cbind(germinadas, 25-germinadas)~trat, data=da, family=binomial) anova(g0, test="Chisq") # testa hipótese sobre os efeitos principais summary(g0) # cada efeito é o contraste controle vs trat_i # usar correção de bonferroni nas comparações com controle #------------------------------------------------------------------------------------------ Se todo mundo continuar fazendo as coisas como todos já fizeram não vamos evoluir, certo? À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/5331e054/atta…> ------------------------------ Message: 12 Date: Wed, 1 Jun 2011 08:13:13 -0300 From: Eder David Borges da Silva <eder(a)leg.ufpr.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTin6WSj1N+aP3h=J+Z9caY15_+pBKg(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Mauro, a pagina completa com os dados estão em: http://www.leg.ufpr.br/~eder/cursoR/ <http://www.leg.ufpr.br/~eder/cursoR/>e o wiki com as explicações em: http://www.leg.ufpr.br/doku.php/cursos:ragronomia <http://www.leg.ufpr.br/doku.php/cursos:ragronomia>Att Em 31 de maio de 2011 22:11, Mauro Sznelwar <sznelwar(a)uol.com.br> escreveu: > Onde eu pego os arquivos para rodar? > Mauro > > José, > Sem o seu codigo fica meio complicado, porem neste link tem uma serie de > analises de estatistica experimental com onde se testa homogenidade e > normalidade. > http://www.leg.ufpr.br/~eder/cursoR/Exercicios.R > att > > Em 31 de maio de 2011 20:13, Eduardo Vargas <duvargas(a)gmail.com>escreveu: > >> >> Boa noite pessoal,preciso de uma ajuda.Estou fazendo uma comparação de >> médias de germinação de sementes, e as varianveis são tratamentos(trat) >> fator e Germinação (germ) numerica.Preciso fazer o teste de homogeneidade de >> variancias com o Bartlett e normalidade com o teste de shapiro-Wilk, >> porém dá um mensagem de erro(Erro em bartlett.test(germ, trat) : objeto >> 'germ' não encontrado) e no Shapiro -wilk essa mensagem Erro em >> `[.data.frame`(x, complete.cases(x)) : >> undefined columns selected >> Será que alguem pode me ajudar? >> -- >> José Eduardo Vargas Lopes de Aráujo >> Engenheiro Agrônomo >> Contato:(38)8811-6985 >> duvargas(a)gmail.com >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> >> > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/8b5cb003/atta…> ------------------------------ Message: 13 Date: Wed, 1 Jun 2011 08:37:36 -0300 From: Eduardo Vargas <duvargas(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] teste de Homogeneidade e Normalidade Message-ID: <BANLkTincK6wOSuWSgLyYEv_A4wmBXQ+q5A(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Obrigado mais uma vez pela informação. Em 1 de junho de 2011 00:52, Walmes Zeviani <walmeszeviani(a)gmail.com>escreveu: > Eduardo, > > O código abaixo mostra como não é difícil analisar os seus dados pela > abordagem que discuti > > > #------------------------------------------------------------------------------------------ > # experimento com 11 tratamentos e 1 controle, resposta binomial n=25, 4 > repetições > > trat <- c("controle", paste("trat", 1:11, sep="")) > da <- expand.grid(trat=trat, rept=1:4) > da$germinadas <- rbinom(nrow(da), size=25, prob=0.8) > levels(da$trat) # controle é o nível de referência > > # ajuste um modelo linear generalizado > g0 <- glm(cbind(germinadas, 25-germinadas)~trat, data=da, family=binomial) > anova(g0, test="Chisq") # testa hipótese sobre os efeitos principais > summary(g0) # cada efeito é o contraste controle vs trat_i > # usar correção de bonferroni nas comparações com controle > > > #------------------------------------------------------------------------------------------ > > Se todo mundo continuar fazendo as coisas como todos já fizeram não vamos > evoluir, certo? > > À disposição. > Walmes. > > > ========================================================================== > Walmes Marques Zeviani > LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) > Departamento de Estatística - Universidade Federal do Paraná > fone: (+55) 41 3361 3573 > VoIP: (3361 3600) 1053 1173 > e-mail: walmes(a)ufpr.br > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -- José Eduardo Vargas Lopes de Aráujo Engenheiro Agrônomo Contato:(38)8811-6985 duvargas(a)gmail.com -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/833c1853/atta…> ------------------------------ Message: 14 Date: Wed, 1 Jun 2011 08:54:31 -0300 From: "Humberto Hazin (hotmail)" <hghazin(a)hotmail.com> To: <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Download de dados pelo R Message-ID: <COL108-DS18964BBC6881C7A035B8CEB47D0(a)phx.gbl> Content-Type: text/plain; charset="iso-8859-1" Bom dia pessoal! Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou conseguindo, alguém poderia me ajudar? ano<-2001 cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", method="auto") Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", : cannot open destfile 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such file or directory' Abs, Humberto -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/7cafde66/atta…> ------------------------------ Message: 15 Date: Wed, 1 Jun 2011 08:58:14 -0300 From: Victor Eduardo <victorduca08(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Download de dados pelo R Message-ID: <BANLkTimuvALu_kvZ3iFTzho3P4WCd+bfpQ(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Humbero, não está dando para acessar o link enviado. Em 1 de junho de 2011 08:54, Humberto Hazin (hotmail) <hghazin(a)hotmail.com>escreveu: > Bom dia pessoal! > > > > > > Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou > conseguindo, alguém poderia me ajudar? > > > > ano<-2001 > > cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" > > > > download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em > R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), > > "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", > method="auto") > > > > Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em > R/ambientais_", : > > cannot open destfile > 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais > ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such > file or directory' > > > > > > Abs, > > > > Humberto > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/d4a14280/atta…> ------------------------------ Message: 16 Date: Wed, 1 Jun 2011 09:04:35 -0300 From: "Humberto Hazin (hotmail)" <hghazin(a)hotmail.com> To: <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] RES: Download de dados pelo R Message-ID: <COL108-DS15CEF73034BD68424BCDFFB47D0(a)phx.gbl> Content-Type: text/plain; charset="iso-8859-1" Olá Eduardo Ele deve estar fora de ar nesse instante ! Tenta com esse aqui http://aspera.jpl.nasa.gov/download/pub/sea_surface_temperature/avhrr/pathfi nder/data_v5/monthly/night/04km/1986/bsst/ ou http://oceandata.sci.gsfc.nasa.gov/SeaWiFS/Mapped/Monthly/9km/chlor/ Abs, Humberto De: r-br-bounces(a)listas.c3sl.ufpr.br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de Victor Eduardo Enviada em: Wednesday, June 01, 2011 8:58 AM Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Download de dados pelo R Humbero, não está dando para acessar o link enviado. Em 1 de junho de 2011 08:54, Humberto Hazin (hotmail) <hghazin(a)hotmail.com> escreveu: Bom dia pessoal! Estou tentando baixar uns arquivos do site abaixo pelo R porem não estou conseguindo, alguém poderia me ajudar? ano<-2001 cho<-"http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/" download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), destfile=paste(getwd(), "~/Dropbox/Humberto/Analises em R/ambientais_", ano, ".hdf"), mode="w", method="auto") Error in download.file(paste(cho, ano, "~/Dropbox/Humberto/Analises em R/ambientais_", : cannot open destfile 'C:/Users/Familia/Documents/Dropbox/Humberto/Analises em R/ambientais ~/Dropbox/Humberto/Analises em R/ambientais_ 2001 .hdf', reason 'No such file or directory' Abs, Humberto _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/18aad303/atta…> ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 4, assunto 1 ***************************************

1 0