- R-br - DINF UFPR

distância entre municípios
by jairan dionizio 13 Mar '12

13 Mar '12

Bom dia a todos ! Estou querendo saber como faço para usar distancia entre municipíos na análise de agrupamento. Jairan Dionizio Graduando em Estatística Skype:jairan.dionizio Facebook: Jairan Dionizio Twitter: @jairandionizio Tel.:(84) 9423-5286 "A força não provém da capacidade física e sim de uma vontade indomável." (Mahatma Gandhi) Cuidados ao enviar e-mail Ao reenviar apague o meu e-mail e o meu nome, apague também os endereços dos amigos antes de reenviar qualquer e-mail, encaminhe como cópia oculta (Cco ou Bcc) aos SEUS destinatários.Agindo sempre assim dificultaremos a disseminação de vírus, spams, banners e outros.

2 2

Re: [R-br] Teste de Kruskal Wallis!
by Thiago de Paula Protásio 13 Mar '12

13 Mar '12

Prezados membros da R-Br; Como nas postagens anteriores estou tentando encontrar respostas sobre o teste t que a função kruskal() pacote agricolae executa após o teste não paramétrico em questão. Contudo, analisando as postagens antigas dos demais membros percebi que existe um outro pacote que faz o teste de Kruskal-Wallis: pgirmess função kruskalmc (). O problema é que essa função retornou um resultado diferente do obtido pela função kruskal() do agricolae. Alguém poderia me ajudar a esclarecer essa dúvida.? Muito obrigado!! Att., -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

2 1

Re: [R-br] Teste de Kruskal Wallis!
by Thiago de Paula Protásio 13 Mar '12

13 Mar '12

Obrigado Allaman pela dica! Olhei com atenção a documentação do pacote, mas não encontrei uma explicação sobre o teste t-Student que é aplicado após a função kruskal(). Não mencionam se é um teste t simples ou com algum tipo de proteção como Bonferroni. Muito obrigado! Att., -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

1 0

Re: [R-br] Ajuda com modelos mistos (medidas repetidas)
by Ivan Bezerra Allaman 13 Mar '12

13 Mar '12

continuando as perguntas, d) o intervalo de 48 h entre as medidas foi o suficiente para torná-las independentes? (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

Re: [R-br] Teste de Kruskal Wallis!
by Ivan Bezerra Allaman 13 Mar '12

13 Mar '12

Coisas interessantes em http://tarwi.lamolina.edu.pe/~fmendiburu/ talvez encontre resposta para sua pergunta! (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Ajuda com modelos mistos (medidas repetidas)
by Olympio Neto 13 Mar '12

13 Mar '12

Pessoal, gostaria de uma ajuda de vocês. Estou com um problema de estudo fatorial com medidas repetidas (assim considero) com as seguintes características: 1- Cada corpo de prova foi submetido a 1 tratamento com 3 níveis distintos, sendo que cada grupo continha 18 unidades amostrais (fator 1), 54 corpos de prova ao total; 2- Dentro de cada grupo definido pelo nível do tratamento no fator 1, metade das unidades amostrais do grupo (9) foram submetidas ao tratamento 2 com 2 níveis (fator 2), de modo que o i-ésimo indivíduo está submetido ao j-ésimo nível do fator 1 e k-ésimo nível do fator 2, ou seja, são independentes com relação ao fator 2. 3-A variável resposta foi medida em dois tempos distintos, tempo 0 e tempo 48h, por isso considerei medidas repetidas. Para cada unidade amostral, portanto, eu considerei um modelo assim: fixando-se a unidade amostra i,que foi somente submetida a 1 tratamento do fator 1 e 1 tratamento do fator 2. Logo, j e k também estão fixos e o que varia dentro de cada indivíduo é o tempo (l) em que a medida foi obtida. O modelo abaixo Y_l=b0+b1Trat1_l+b2*Trat2_l+b3*Tempo_l+todas as interações+a_i+e_il, l=1,2, a_i ~N(0, var_1) (efeito aleatório) e e_il~Normal(0, var_2) Desta maneira, as pseudorrepresentações são os tempos, com tamanho 2 para cada unidade e isso vem causando sérios problemas, porque acredito que o tratamento 2 deveria estar aninhado no tratamento 1, mas quando tento algo, os graus de liberdade são, obviamente insuficientes e um erro ocorre. Estimado o modelo acima, com a função lme, (algo do tipo lme(y~Trat1*Trat2*Tempo, random=1|Id)), os resíduos não são normais. Existe alguma maneira de utilizar a transformação de box-cox para tentar normalizar os resíduos. Eu fiz um aproximação usando o modelo linear comum com estes dados e a função boxcox do pacote MASS, e depois usei a transformação indicada no modelo linear neste misto, variando os valores de lambda em torno do valor estimado e vendo qual se aproximava mais da normalidade. Ficou muito ruim, ainda. Então, lanço as minhas perguntas: a) Com a experiência de vocês me digam, por favor; estou usando um modelo adequado para que quero encontrar? b) Existe algum pacote no R que calcule a transformação de box-cox em objetos lme? c) Será que o fator 2, dada natureza do experimento, deve ser aninhada ao fator 1 como efeito fixo? Algo como Y~Trat1*Trat1/Trat2*Tempo? Sei que estas perguntas fogem ao escopo deste fórum, mas logo imaginei que vocês pudessem me ajudar. Desde já muito obrigado. Olympio

1 0

Teste de Kruskal Wallis!
by Thiago de Paula Protásio 12 Mar '12

12 Mar '12

Prezados membros da R-Br; Novamente recorro ao conhecimento de vcs. Estou tendo algumas dúvidas sobre a função kruskal() do pacote agricolae. Essa função executa o teste de kruskal e posteriormente um teste t. Recorri ao help() para entender esse teste t, mas não obtive sucesso. Gostaria de saber se alguém poderia confirmar ou não se esse teste t tem algum tipo de proteção como Bonferroni ou se é apenas um teste t simples (nesse caso não seria o mais apropriado eu acho). Alguém poderia indicar uma função que executa o teste de kruskal e posteriormente o teste t com a proteção de Bonferroni? Grato pela valiosa ajuda. -- Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Departamento de Ciências Florestais Ciência e Tecnologia da Madeira (035) 9183 - 2246

1 0

problemas com função nls
by Fernando Antonio de souza 11 Mar '12

11 Mar '12

caros amigos estou tentando rodar um modelo não linear de Gormpertz da forma Log(Y) = a-b*exp(-c*X) utiliando a função nls. A função é como segue: gompertz<-nls(log(PESO.FETO)~a-b*exp(-c*Gest),start=list(a=3.4,b=22.99,c=0.018),data=gestsimp) dput(gestsimp) Acontece que toda vez que rodo, recebo o aviso de erro:" Erro em nls(log(PESO.FETO) ~ a - b * exp(-c * Gest), start = list(a = 3.4, : número de iterações excedeu o máximo de 50 Ja tentei elevar o máximo de interação utilizando o comando controle mas nunca converge. Alguem pode me esclarecer em que estou errando. os valore iniciais foram obtido sutilizando a função manipulate. um grande abraço

2 1

erro "initRefFields"
by D'Apolito 10 Mar '12

10 Mar '12

Olá, Alguém sabe como solucionar o erro: "Error in initFields(scales = scales) : could not find function "initRefFields" " ...sem ter que reinstalar o R 2.14.1? obrigado Carlos

2 7

Re: [R-br] problemas com função nls
by Ivan Bezerra Allaman 10 Mar '12

10 Mar '12

Provavelmente seus chutes estão errados! Acesse o site http://ridiculas.wordpress.com/ que provavelmente irás resolver o seu problema. Do contrário, nos forneça um CMR. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

3 8