- R-br - DINF UFPR

porque não recebo meus e-mails enviados para a lista
by Luciana Melo 28 Abr '11

28 Abr '11

2 1

R² ajustado de Modelos Não Lineares!
by Thiago De paula protásio 28 Abr '11

28 Abr '11

Pessoal; Gostaria de perguntar-lhe uma coisa: Como faço para extrair o R² ajustado dos modelos não lineares?Agradeço-lhes qualquer ajuda! Att., Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

4 5

Converter valores em numeros 1 correspondentes
by Alexandre dos Santos 27 Abr '11

27 Abr '11

Boa tarde amigos do R, Estou com uma conversão que imagino ser bem complicada de fazer, preciso pegar vários objetos como por exemplo um deles: > m.2 [1] 5 2 1 2 2 1 1 5 6 8 2 1 3 6 2 12 7 2 19 1 1 3 7 2 3 2 1 1 4 5 1 3 2 1 3 8 3 37 3 3 2 2 3 1 1 1 1 1 4 3 1 4 2 1 3 3 55 5 1 3 3 1 1 4 E converter para um novo vetor somente com números 1 correspondentes a cada um dos valores em m.2 , obviamente que os valores que já são 1 não seriam modicados, mas o primeiro valor o 5, seria substituido por 1 1 1 1 1, o 2 por 1 1 , assim successivamente para todos os valores, Alguém poderia me dar uma luz, se é que tem jeito, Obrigado, Alexandre dos Santos Ingenieur forestier, Msc. INRA- Biostatistique et Processus Spatiaux (BioSP) Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29

6 12

Comando para organizar os dados da planilhas
by Victor Eduardo 27 Abr '11

27 Abr '11

Como eu posso organizar uma planilha gigantesta de acordo com o dado apenas da primeira coluna? Por exemplo, tenho uma planilha com a primeira coluna com vários dados misturados e queria colocar em ordem, por exemplo, os dados são A801, A878, A602, A652, no caso queria que organizasse meus dados da planilha em A602, A652, A801, A878. Abraços

4 4

Problemas com leitura do arquivo
by Victor Eduardo 27 Abr '11

27 Abr '11

Galera, como faço para ler um arquivo que tem essa cara: A322 2011 04 24 00 12.5 23 20.6 20.7 20.6 93 93 92 19.4 19.5 19.3 922.0 922.1 921.5 2.1 103 3.9 -2.790 0.0 / //// ///// ///// = A811 2011 04 24 00 12.7 15 12.3 13.1 12.3 56 59 55 3.9 5.2 3.8 964.3 964.3 964.1 5.9 223 9.3 -3.540 0.0 / //// ///// ///// = A008 2011 04 24 00 12.4 23 20.5 21.1 20.1 64 66 62 13.4 13.6 13.4 888.8 888.8 888.4 2.0 98 3.3 -4 0.0 / //// ///// ///// = A025 2011 04 24 00 12.6 25 21.4 22.7 20.8 73 75 67 16.4 16.6 15.9 927.3 927.3 926.7 1.7 232 2.6 -3.536 0.0 / //// ///// ///// = Eu uso o read.table normal, mas fica dando o seguinte erro: Erro em scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : linha 53 não tinha 29 elementos Abraços

2 2

Dúvida para obtenção do poder do teste
by Olympio Neto 27 Abr '11

27 Abr '11

Pessoal, bom dia. Estou com uma dúvida. Pretendo calcular o poder do teste via simulação para dados de correlação intra-classe. Que eu saiba, não existe um método formal para isso. Assim, como pretendo fazê-lo via simulação gostaria de saber o seguinte: Devo replicar a amostra que tenho (como é feito no bootstrap) e tomar o número de vezes em que o teste foi menor que meu nível de significância, ou estes dados devem ser gerados aleatoriamente? O problema é que não há uma distribuição evidente para estes dados. Outra coisa, se for utilizar os dados que já disponho, posso alterar o número de elementos na reamostragem para verificar o efeito do tamanho da amostra no poder do teste e a fim de gerar um gráfico? Por exemplo, a amostra tem 100 avaliações, quero saber o poder do teste com 10, 30, 70 e 90, além das 100... é incorreto proceder dessa maneira? Muito obrigado pela atenção. Olympio

2 2

Teste de Aderência
by Cristiano Melo 27 Abr '11

27 Abr '11

É o seguinte: tenho em um arquivo txt um vetor que representa tempos até a falha de equipamentos. Gostaria de fazer alguns teste de aderência para verificar se estes dados se aproximam de algumas distribuições de probabilidade. Usei o lillie.test(dados) para verificar se dos dados aderem a uma distribuição nomal. No entanto, gostaria de verificar se estes mesmos dados (e algumas variações) se aderem a uma exponencial, gamma e weibull. Sei que a função é a ks.test(x, y,..., alternative=c("two.sided" "less" or "greater")) para o teste de Kolmogorov-Smirnov. Para testar normalidade com a função ks, fiz o seguinte: ks.test(vetor, "pnorm", sd=sd(vetor), mean=mean(vetor),alternative=c("two.sided")). Curioso que o resultado foi bem diferente da lillie.test, a seguinte mensagem foi apresentada: Warning message: In ks.test(vetor, "pnorm", sd=sd(vetor), mean=mean(vetor),alternative=c("two.sided")), : não é possível calcular os níveis descritivos corretos com empates. O que isso quer dizer??????? Quando tentei usar outra distribuição, pweibull por exemplo, o p-value foi menor que 2.2e-16, ou seja, nada a ver, e repetindo a mesma frase anterior. O mesmo resultado foi com as outras. Como não sei a sintaxe para weibull fiz o seguinte: ks.test(vetor, "pweibull", 1.129, 2,alternative=c("two.sided")) Onde estou errando? Vi que para montar uma pweibull são necessários o vetor de quantis e os parâmetros shape e scale. É necessário fazer separado e quardar em uma variável e depois jogar na ks? Como consigo esse vetor de quantis? Achei que seria automático. Estou correto se fizer assim: ks.test(vetor, "pweibull",10,2) E tem alguma forma de estimar os parâmetros shape e scale desta função?

5 16

Análise de Cluster em R
by karla Sabrina Magalhaes 27 Abr '11

27 Abr '11

> Prezados, > > Estou precisando rodar uma análise de Cluster em R, será que alguem sabe > como posso fazer isto, pois não consegui achar o tutorial na internet. > Alguem já fez este tipo de análise, pois não estou conseguindo, fazer os comandos. > > Aguardo respostas > Grata, > > -- > -- > Atenciosamente, > Karla Sabrina Magalhães > Bióloga - Mestranda em Bioengenharia Celular e Tecidual - UFSJ > MASP - 1240384 - 6 > > > > > -- -- Atenciosamente, Karla Sabrina Magalhães Bióloga - Mestranda em Bioengenharia Celular e Tecidual - UFSJ MASP - 1240384 - 6

2 1

Re: [R-br] R² ajustado de Modelos Não Lineares!
by Thiago De paula protásio 26 Abr '11

26 Abr '11

Ótima sugestão Walmes.....Já estou utilizando essa maneira simplificada de cálculo!!! Obrigado!Att., Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Teste de Autocorrelação Residual! Durbin-Watson!
by Thiago De paula protásio 26 Abr '11

26 Abr '11

Pessoal; Novamente conto com a experiência de vcs! Gostaria de fazer um teste para a comprovação ou não da presença de autocorrelação residual de modelos não lineares......uma vez que isso pode deixar as estimativas de y tendenciosas...... Para equações lineares eu sempre usei a função dwtest() do pacote lmtest, mas para modelos não lineares essa função não é válida..... Gostaria de perguntar-lhes se existe alguma função pronta no R para que eu possa realizar esse teste! Qualquer ajuda ou opinião é sempre muito bem vinda!Att., Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1