[R-br] dados DTK.test

Ze Henrique jhguilhen em gmail.com
Domingo Março 10 12:27:24 BRT 2013


Olhe este segundo aqui.


library(geoR)
library(agricolae)
dados=read.geodata("c:/fsd.txt ",coords.col=13:14,data.col=4,sep =
"",header=T)
y=read.geodata("c:/fsd.txt ",coords.col=13:14,data.col=1,sep = "",header=T)
n=summary(dados)[[1]]
n
MSerror1<-aov(y$data ~ dados$data)
summary(MSerror1)
HSD.test(dados$data,y$data,n-1,0.540,alpha=0.05)


Parece que MSerror tem que aplicar a anova e ver o valor dentro dela e
colocar, não tenho certeza.

Ate..




Em 10 de março de 2013 12:05, Ze Henrique <jhguilhen em gmail.com> escreveu:

> Entendi o que faltou para vc conseguir rodar é o pacote agricolae.
>
> Use assim que vc conseguirá:
>
> library(geoR)
> library(agricolae)
> dados=read.geodata("c:/fsd.txt ",coords.col=13:14,data.col=4,sep =
> "",header=T)
> y=read.geodata("c:/fsd.txt ",coords.col=13:14,data.col=1,sep = "",header=T)
> n=summary(dados)[[1]]
> n
> MSerror=19 # não sei o valor disto e nem pra que serve.
> HSD.test(dados$data,y$data,n-1,MSerror,alpha=0.05)
>
> Se precisar de algo avise.
>
> Ate.
>
> Em 10 de março de 2013 11:48, Hélio Gallo Rocha <heliogallorocha em gmail.com
> > escreveu:
>
>> José Henrique,
>> Não sou expert em R, na realidade tenho usado apenas para geoestatistica,
>> não tenho feito a descritiva nele, uso o sisrvar.
>> Aprendi a carregar os dados no R na UFLA com o José Marcio no LEMAF.
>> Eu acho essa forma de "carregar" os atributos bastante fácil, pois
>> normalmente tenho vários atributos a serem analisados, dai só indico uqal a
>> coluna que vou analisar naquele momento.
>>
>> Fiz o seguinte, com om mesmo arquivo anexado no inicio "FSD.txt", que são
>> dados para previsão de safra em café, meço altura, diâmetro e numero de
>> frutos, que serão comparados em quatro modelos para previsão.
>>
>>
>> O script rodou, mas emperrei no MSerror, que chutei 19 script:
>>
>> dados=read.geodata("fsd.txt ",coords.col=13:14,data.col=4,sep =
>> "",header=T)
>> y=read.geodata("fsd.txt ",coords.col=13:14,data.col=1,sep = "",header=T)
>> n=summary(dados)[[1]]
>> n
>> MSerror=19
>> HSD.test(dados$data,y$data,n-1,MSerror,alpha=0.05)
>>
>> Abraço e bom domingo!
>>
>> Hélio
>>
>> Em 10 de março de 2013 08:44, José Henrique [via R-br] <
>> ml-node+s2285057n4658522h50 em n4.nabble.com> escreveu:
>>
>>> Um mais ou menos entendi o que o Sr. tentou fazer, mas o comando
>>> attach() tem apenas a função de vc não precisar ficar fixando os dados da
>>> tabela que vc necessita utilizar no exercício que está sendo realizado.
>>>
>>> Não respondi antes devido ontem o dia ser corrido.
>>>
>>> Porem vc consegui fazer o que precisava ou ainda não?
>>>
>>> Outra cousa nunca trabalhei com geoestatistica mas quando quero fazer
>>> teste de Tukey não utilizo "DTK.test", pois acho que talvez este comando
>>> faz mais que apenas só o teste que vc quer, caso vc queira para fazer o
>>> teste de Tukey utilize o pacote (agricolae) e o comando
>>> HSD.test(x,trat,Glres,Qmres,alpha=ns)
>>>
>>> Caso vc queira posso tentar te ajudar.
>>>
>>> Ate logo.
>>>
>>> Abraço.
>>>
>>>
>>>
>>> Em 9 de março de 2013 10:15, Hélio Gallo Rocha <[hidden email]<http://user/SendEmail.jtp?type=node&node=4658522&i=0>
>>> > escreveu:
>>>
>>>> Caro José Henrique
>>>> No script que enviei, na linha 4, faltou indicar qual a coluna que
>>>> estava analisando, no caso seria a 1, data.col=1
>>>> Então X seria a coluna 6 que é a relação diametro/altura
>>>> F seria a coluna 1 = Variedade ( que no caso são 3), cada uma com 14
>>>> repetições.
>>>>
>>>> Ma tentei  a sua forma de carregar os dados, não compreendi direito o
>>>> comando attach(fsd), ( depois eu estudo), mas retornou o mesmo erro.
>>>>
>>>> Abraço
>>>>
>>>>
>>>> --
>>>> Hélio Gallo Rocha
>>>> IFSULDEMINAS - Câmpus Muzambinho
>>>>
>>>> _______________________________________________
>>>> R-br mailing list
>>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658522&i=1>
>>>>
>>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>>> código mínimo reproduzível.
>>>>
>>>
>>>
>>>
>>> --
>>> *José Henrique Soler Guilhen*
>>> Graduando em Agronomia
>>> Universidade Federal do Espírito Santo - CCAUFES
>>>
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658522&i=2>
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>> código mínimo reproduzível.
>>>
>>> ------------------------------
>>>  If you reply to this email, your message will be added to the
>>> discussion below:
>>>
>>> http://r-br.2285057.n4.nabble.com/R-br-dados-DTK-test-tp4658513p4658522.html
>>>  To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4>
>>> .
>>> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>>>
>>
>>
>>
>> --
>> Hélio Gallo Rocha
>> IFSULDEMINAS - Câmpus Muzambinho
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
>
> --
> *José Henrique Soler Guilhen*
> Graduando em Agronomia
> Universidade Federal do Espírito Santo - CCAUFES
>
>


-- 
*José Henrique Soler Guilhen*
Graduando em Agronomia
Universidade Federal do Espírito Santo - CCAUFES
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130310/5703273e/attachment.html>


Mais detalhes sobre a lista de discussão R-br