- R-br - DINF UFPR

Re: [R-br] regressão multivariada
by Ivan Bezerra Allaman 28 Abr '11

28 Abr '11

Oi Ayuni, Nos últimos dois meses, estou estudando um pouco sobre análise multivariada, e até onde eu entendi, a análise de correlação canônica é uma extensão da análise de regressão múltipla. Neste tipo de análise você têm um grupo de variáveis dependentes e outro grupo de variáveis independentes. Dê uma olhada na teoria, acho que é o que você procura. O bom pacote que achei para fazer estes tipos de análises e a "candisc". Baixe o pdf do pacote e veja como utilizar as funções. Outro pacote que gera gráficos ótimos nestes casos é o "heplot". Só para lembrar, a candisc trabalha com variáveis canônicas. Veja em qual caso ela se aplique, pois pode ser que a utilização de componentes principais seja mais apropriado. Você poderia fazer também: modelo <- lm(cbind(var1, var2, etc..) ~ x1 + x2 + etc...) summary(modelo) Dê uma olhada na MASS para estes tipos de utilização. Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

4 4

regressão multivariada
by Ivan Bezerra Allaman 28 Abr '11

28 Abr '11

Opa!! Certinho Benilton!! Aquele caso se aplica quando o interesse é uma MANOVA. Ayuni, Desconsidere a parte da lm, ok? No entanto, os demais comentários se mantém. Apenas uma correção se vc me permite. A library "candisc", permite fazer um ajuste quando você tem um grupo de variáveis dependentes numéricas e um grupo de variáveis independentes como fatores ou do tipo dummy. Se vc quer fazer uma análise de regressão multivariada em que se tem dois grupos de variáveis numéricas propriamente ditas (matriz Y e X) aí você pode dar uma olhada no pacote "CCA" que ele tem as funções para isso. Compare as funcionalidades dos dois pacotes (candisc e CCA). Abraço. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Re: [R-br] R² ajustado de Modelos Não Lineares! - complementando mensagem
by Mauro Sznelwar 28 Abr '11

28 Abr '11

Eu copiei direto do corpo da mensagem Obrigado pelo retorno, mas originou uma só coluna em vez de duas, veja um trecho: Dias 1 Tambaqui 2 0 3 6.09 4 0 5 6.09 6 0 7 6.09 8 30 9 65.73 10 30 11 52.00 12 30 13 53.13 14 60 15 105.88 16 60 17 115.80 18 60 Olá Mauro. O comando "clipboard" extrai da área de trasnferência. Faça o seguinte. digite o comando abaixo e deixe o cursor no final do texto, após ) x=read.table('clipboard', h=T) Vá para a tabela em um arquivo de texto ou planilha e copie para a área de transferencia Volte para a área de trabalho do R e digite enter Pronto, o R recupera da área de transfêrência e vc cria o objeto x com a tabela que vc quiser. Valeu.

2 1

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 27 Abr '11

27 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