Não sei se a forma de inserir os dados está correta. Anexei o fsd.txt que tem objetivo de estudo de geoestatistica, mas estou tentando rodar o teste de tukey. script: library(geoR) library(DTK) dados=read.geodata("fsd.txt ",coords.col=13:14,data.col=6,sep = "",header=T) f=read.geodata("fsd.txt ",coords.col=13:14,data.col=,sep = "",header=T) DTK.test (x = "dados$data", f = "f$data", a = "0.05") grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
Caro colega verifique se isso te cerve: library(geoR) library(DTK) fsd<-read.table("fsd.txt", h=T) attach(fsd) dados=read.geodata("c:/fsd.txt ",coords.col=13:14,data.col=6,sep = "",header=T) f=read.geodata("c:/fsd.txt ",coords.col=13:14,data.col=,sep = "",header=T) DTK.test (x = "Alt.", f = "Variedade", a = 0.05) # porem tende que fazer umas por vez ou será necessário criar um lupo para vc fazer tude se uma vez. Bem se foi isso que eu entendi que vc queria fazer né. Caso precise mais alguma ajuda e eu possa avise. ate... Em 9 de março de 2013 09:11, Hélio Gallo Rocha <heliogallorocha@gmail.com>escreveu:
Não sei se a forma de inserir os dados está correta. Anexei o fsd.txt que tem objetivo de estudo de geoestatistica, mas estou tentando rodar o teste de tukey.
script: library(geoR) library(DTK) dados=read.geodata("fsd.txt ",coords.col=13:14,data.col=6,sep = "",header=T) f=read.geodata("fsd.txt ",coords.col=13:14,data.col=,sep = "",header=T) DTK.test (x = "dados$data", f = "f$data", a = "0.05")
grato
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list R-br@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
Não sei se a forma de inserir os dados está correta. Anexei o fsd.txt que tem objetivo de estudo de geoestatistica, mas estou tentando rodar o teste de tukey. script: library(geoR) library(DTK) dados=read.geodata("fsd.txt ",coords.col=13:14,data.col=6,sep = "",header=T) f=read.geodata("fsd.txt ",coords.col=13:14,data.col=,sep = "",header=T) DTK.test (x = "dados$data", f = "f$data", a = "0.05") grato Em 9 de março de 2013 09:41, José Henrique [via R-br] < ml-node+s2285057n4658514h2@n4.nabble.com> escreveu:
Caro colega verifique se isso te cerve:
library(geoR) library(DTK) fsd<-read.table("fsd.txt", h=T) attach(fsd) dados=read.geodata("c:/fsd.txt ",coords.col=13:14,data.col=6,sep = "",header=T) f=read.geodata("c:/fsd.txt ",coords.col=13:14,data.col=,sep = "",header=T) DTK.test (x = "Alt.", f = "Variedade", a = 0.05) # porem tende que fazer umas por vez ou será necessário criar um lupo para vc fazer tude se uma vez.
Bem se foi isso que eu entendi que vc queria fazer né.
Caso precise mais alguma ajuda e eu possa avise.
ate...
Em 9 de março de 2013 09:11, Hélio Gallo Rocha <[hidden email]<http://user/SendEmail.jtp?type=node&node=4658514&i=0>
escreveu:
Não sei se a forma de inserir os dados está correta. Anexei o fsd.txt que tem objetivo de estudo de geoestatistica, mas estou tentando rodar o teste de tukey.
script: library(geoR) library(DTK) dados=read.geodata("fsd.txt ",coords.col=13:14,data.col=6,sep = "",header=T) f=read.geodata("fsd.txt ",coords.col=13:14,data.col=,sep = "",header=T) DTK.test (x = "dados$data", f = "f$data", a = "0.05")
grato
-- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list [hidden email] <http://user/SendEmail.jtp?type=node&node=4658514&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=4658514&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-tp4658513p4658514.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
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
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 <heliogallorocha@gmail.com>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 R-br@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, 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@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
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@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@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@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
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@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@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@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@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
José Henrique, Funcionou Quanto a ideia do Tiago, vou testar até quarta.. agora uma boa macarronada Abraço e grato pela ajuda Hélio Em 10 de março de 2013 12:27, José Henrique [via R-br] < ml-node+s2285057n4658527h5@n4.nabble.com> escreveu:
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 <[hidden email]<http://user/SendEmail.jtp?type=node&node=4658527&i=0>
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 <[hidden email]<http://user/SendEmail.jtp?type=node&node=4658527&i=1>
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] <[hidden email]<http://user/SendEmail.jtp?type=node&node=4658527&i=2>
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. 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 [hidden email] <http://user/SendEmail.jtp?type=node&node=4658527&i=3> 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
_______________________________________________ R-br mailing list [hidden email] <http://user/SendEmail.jtp?type=node&node=4658527&i=4> 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-tp4658513p4658527.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
Hélio se o seu objetivo é fazer um teste de médias eu te indico a usar o pacode ExpDes disponivel no site https://sites.google.com/site/ericbferreira/unifal/downloads-1 com ele você usa os seguintes comandos. Ex:y #nome do seu arquivo de dados attach(y) require(ExpDes) trat #nome dos seus tratamentos bloco # nome da coluna com os blocos caso o seu experimento tenha sido em DBC caso contrario use o comando para DIC dbc(trat,bloco,variavel resposta,quali=TRUE,mcomp="tukey",sigT=0.05) #Exemplo para dbc. dic(trat,variavel resposta,quali=TRUE,mcomp="tukey",sigT=0.05) #Exemplo para dic Att. Tiago. Date: Sun, 10 Mar 2013 11:48:27 -0300 From: heliogallorocha@gmail.com To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] dados DTK.test 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]]nMSerror=19HSD.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@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]> 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/alturaF 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] 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] 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. NAML -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho _______________________________________________ R-br mailing list R-br@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.
participantes (3)
-
Hélio Gallo Rocha -
Tiago Souza Marçal -
Ze Henrique