[R-br] RStudio
Diogo Borges Provete
diogoprov em yahoo.com.br
Quarta Junho 1 10:22:00 BRT 2011
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 em 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 em listas.c3sl.ufpr.br" <r-br-request em listas.c3sl.ufpr.br>
Para: r-br em 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 em 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 em listas.c3sl.ufpr.br
Você poderá entrar em contato com a pessoa que gerencia a lista pelo
endereço
r-br-owner em 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 em gmail.com>
To: r-br em listas.c3sl.ufpr.br
Cc: Samuel Carvalho <samukajm em yahoo.com.br>
Subject: Re: [R-br] modifica dataframe
Message-ID: <BANLkTikLCzKcVdHqgWWUSv8f+Z8UMayxKA em 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 em 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 em 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 em 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 em 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 em yahoo.com.br>
To: R Brasil <r-br em listas.c3sl.ufpr.br>
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <947015.46346.qm em 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 em 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/attachment-0001.html>
------------------------------
Message: 3
Date: Tue, 31 May 2011 21:46:46 -0300
From: Eduardo Vargas <duvargas em gmail.com>
To: r-br em listas.c3sl.ufpr.br, Ivan Bezerra Allaman
<ivanalaman em yahoo.com.br>
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <BANLkTinXDFS=T7UD7GrEKSxy-0z1-i77FA em 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 em 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 em yahoo.com.br
> Tel: (35)3826-6608/9900-2924
>
> _______________________________________________
> R-br mailing list
> R-br em 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 em gmail.com
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/7c9a93de/attachment-0001.html>
------------------------------
Message: 4
Date: Tue, 31 May 2011 22:11:34 -0300
From: "Mauro Sznelwar" <sznelwar em uol.com.br>
To: <r-br em listas.c3sl.ufpr.br>
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <00f101cc1ff8$ddcd1fe0$6700a8c0 em 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 em 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 em gmail.com
_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
------------------------------------------------------------------------------
_______________________________________________
R-br mailing list
R-br em 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/attachment-0001.html>
------------------------------
Message: 5
Date: Tue, 31 May 2011 22:16:25 -0300
From: Eduardo Vargas <duvargas em gmail.com>
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <BANLkTin=XzAn+6AVxUAnEyeNqDCZQ1Ykug em 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 em 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 em 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 em gmail.com
>>
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>
>>
> ------------------------------
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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 em gmail.com
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/74c6398a/attachment-0001.html>
------------------------------
Message: 6
Date: Tue, 31 May 2011 22:20:04 -0300
From: "Mauro Sznelwar" <sznelwar em uol.com.br>
To: <r-br em listas.c3sl.ufpr.br>
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <002201cc1ffa$0decce90$6700a8c0 em 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 em 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/attachment-0001.html>
------------------------------
Message: 7
Date: Tue, 31 May 2011 22:53:19 -0300
From: Walmes Zeviani <walmeszeviani em gmail.com>
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <BANLkTikem0jdfVDhKZ8N-n8b83QyNkgXvw em 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_usando_modelo_linear_generalizado
À 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 em 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/attachment-0001.html>
------------------------------
Message: 8
Date: Tue, 31 May 2011 23:10:54 -0300
From: Henrique Dallazuanna <wwwhsd em gmail.com>
To: r-br em listas.c3sl.ufpr.br, Samuel Carvalho <samukajm em yahoo.com.br>
Subject: Re: [R-br] modifica dataframe
Message-ID: <BANLkTi=77BMKw2WfJopo7xjR--he5P03VQ em 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 em 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 em 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/attachment-0001.html>
------------------------------
Message: 9
Date: Tue, 31 May 2011 23:34:58 -0300
From: Eduardo Vargas <duvargas em gmail.com>
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <BANLkTiniir0GPKzWcMm_=PsYwrBfmeJTjA em 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 em 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_usando_modelo_linear_generalizado
>
> À 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 em ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em 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 em gmail.com
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110531/43f2d612/attachment-0001.html>
------------------------------
Message: 10
Date: Tue, 31 May 2011 20:02:08 -0700 (PDT)
From: Samuel Carvalho <samukajm em yahoo.com.br>
To: "r-br em listas.c3sl.ufpr.br" <r-br em listas.c3sl.ufpr.br>
Subject: Re: [R-br] modifica dataframe
Message-ID: <369328.15725.qm em 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 em gmail.com>
Para: r-br em listas.c3sl.ufpr.br; Samuel Carvalho <samukajm em 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 em 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 em 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/attachment-0001.html>
------------------------------
Message: 11
Date: Wed, 1 Jun 2011 00:52:54 -0300
From: Walmes Zeviani <walmeszeviani em gmail.com>
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <BANLkTimC0Ux64PEz2DChLzsedsWigrmoTQ em 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 em 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/attachment-0001.html>
------------------------------
Message: 12
Date: Wed, 1 Jun 2011 08:13:13 -0300
From: Eder David Borges da Silva <eder em leg.ufpr.br>
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <BANLkTin6WSj1N+aP3h=J+Z9caY15_+pBKg em 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 em 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 em 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 em gmail.com
>>
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>
>>
> ------------------------------
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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/attachment-0001.html>
------------------------------
Message: 13
Date: Wed, 1 Jun 2011 08:37:36 -0300
From: Eduardo Vargas <duvargas em gmail.com>
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br] teste de Homogeneidade e Normalidade
Message-ID: <BANLkTincK6wOSuWSgLyYEv_A4wmBXQ+q5A em 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 em 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 em ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em 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 em gmail.com
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/833c1853/attachment-0001.html>
------------------------------
Message: 14
Date: Wed, 1 Jun 2011 08:54:31 -0300
From: "Humberto Hazin (hotmail)" <hghazin em hotmail.com>
To: <r-br em listas.c3sl.ufpr.br>
Subject: [R-br] Download de dados pelo R
Message-ID: <COL108-DS18964BBC6881C7A035B8CEB47D0 em 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/attachment-0001.html>
------------------------------
Message: 15
Date: Wed, 1 Jun 2011 08:58:14 -0300
From: Victor Eduardo <victorduca08 em gmail.com>
To: r-br em listas.c3sl.ufpr.br
Subject: Re: [R-br] Download de dados pelo R
Message-ID: <BANLkTimuvALu_kvZ3iFTzho3P4WCd+bfpQ em 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 em 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 em 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/attachment-0001.html>
------------------------------
Message: 16
Date: Wed, 1 Jun 2011 09:04:35 -0300
From: "Humberto Hazin (hotmail)" <hghazin em hotmail.com>
To: <r-br em listas.c3sl.ufpr.br>
Subject: [R-br] RES: Download de dados pelo R
Message-ID: <COL108-DS15CEF73034BD68424BCDFFB47D0 em 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 em listas.c3sl.ufpr.br
[mailto:r-br-bounces em listas.c3sl.ufpr.br] Em nome de Victor Eduardo
Enviada em: Wednesday, June 01, 2011 8:58 AM
Para: r-br em 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 em 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 em 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/attachment.html>
------------------------------
_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Fim da Digest R-br, volume 4, assunto 1
***************************************
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110601/cdce1f1e/attachment-0001.html>
Mais detalhes sobre a lista de discussão R-br