- R-br - DINF UFPR

Aplicação de métodos de dados faltantes no ambiente R
by Renato Rocha Batista 31 Mai '19

31 Mai '19

Boa tarde Prezados. Estou a tratar dados experimentais em excel, sendo que alguns dados estão em falta em relação a algumas variáveis e objetos (entradas em Excel). Vou realizar a análise multivariada com dados de uma planilha (40 objetos versus 24 colunas de dados), mas antes, eu tenho que prever os dados em falta de algumas células. Por gentileza alguém teria recomendação de script no R com a implementação dos métodos? - Método do algoritmo E-M (Expectativa-Maximização) que pressupõe a normalidade dos dados referentes às variáveis de resposta e - Modelo de regressão logística. A equação gerada permitiria estimar os valores dos dados perdidos. Não assume uma distribuição normal das variáveis de resposta. Ao contrário da regressão linear cujos dados variam de - ∞ a + ∞, na regressão logística eles variam de 0 a 1 e, portanto, são apropriados para variáveis cujos valores variam de 0 a 1. Em vista do exposto, alguém poderia me ajudar ou recomendar alguém que saiba implementar esses métodos no R? Grato pela atenção. Obrigado

2 1

Problema na importação dos dados
by Ze Henrique 29 Mai '19

29 Mai '19

Bom dia senhores, estou com problema para importar dados. Estou utilizando o pacote R "RcppCNPy", pois o arqui está na extensão ".npy". Dos 3 arquivos que tenho, apenas 1 importar os outros dois quando vou importar o R fecha. Alguém sabe o motivo e como poderia resolver este problema? Obs: Arquivo que importa tem 103 KB e os que não consigo importar tem 106 KB e 465.000 KB. Desde já agradeço pela ajuda. -- *José Henrique Soler Guilhen* Doutorando em Genética e Melhoramento Universidade Federal do Espírito Santo - CCAUFES

2 2

Solução alternativa evitando uso de loop
by Paulo Eduardo de Mesquita 28 Mai '19

28 Mai '19

Prezado(a)s colegas da lista, Peço a vossa ajuda nessa questão. Embora eu já tenha uma solução, gostaria de explorar eventuais alternativas que evitem o uso do loop for. Eu tenho um vetor de números inteiros v1 com 10 elementos. Eu preciso de um código que produza um outro vetor numérico v2, também com 10 elementos, em que: - o valor do primeiro elemento seja igual ao valor do primeiro elemento de v1; - o valor do segundo elemento seja igual a soma do primeiro e do segundo elemento de v1; - o valor do terceiro elemento seja igual a soma do primeiro, do segundo e do terceiro elementos de v1; - todos os demais elementos (do 4o ao 10o) sejam iguais a soma dele e dos três elementos imediatamente anteriores a ele. Tentei soluções aplicando as funções da familia apply, mas acaba caindo no loop. Eu imagino que seja possível uma solução utilizando a função cumsum ou alguma variante que restrinja a soma aos últimos 4 elementos e defina uma regra para a exceção representada pelos 03 primeiros elementos. Qualquer ajuda nessa implementação é muito benvinda. Grato pela atenção, Código: v1 <- c(1,2,3,4,5,6,7,8,9,10) v2 <- numeric(length = length(v1)) for (i in seq_along(v1)) { v2[i] <- ifelse(i < 2, v1[i], ifelse(i < 3, v1[i - 1] + v1[i], ifelse(i < 4, v1[i - 2] + v1[i - 1] + v1[i], v1[i - 3] + v1[i - 2] + v1[i - 1] + v1[i]))) } v2 data.frame(v1, v2) -- Paulo Eduardo de Mesquita Disciplina de Infectologia - Faculdade de Medicina Universidade do Oeste Paulista Presidente Prudente - São Paulo - Brasil telefone: 5518 97718261

3 2

importar dados do Sistema IOS macbook para o R
by Jonatan Alexandre De Oliveira 28 Mai '19

28 Mai '19

Caros, Estou tentando importar dados txt para o R, destaco que meu sistema é o IOS, o arquivo encontra-se numa pasta dentro do desktop. Tentei o seguinte comando. dados<read.table("desktop:/teste/arquivo.txt",h=T) O arquivo encontra-se no Desktop, pasta R, arquivo jonatanm.txt. ENTRETANTO, não obtive sucesso, alguém sabe como solucionar este caso? Agradeço desde já pela atenção jonatan -- *JONATAN OLIVEIRA | Phd Candidate Research Group - Food Security and Institutional ActionsFaculty of Agrarian and Veterinary Sciences of Jaboticabal - UNESP+5516996251552Institute of Geosciences and Exact Sciences of - UNESP/Rio Claro+551935269149*

3 2

Adequação modelo linear misto
by Lucas Oliveira 27 Mai '19

27 Mai '19

[OFF-TOPIC] Boa tarde a todos! Estou analisando um experimento onde quero avaliar a adaptabilidade e estabilidade de alguns clones em diferentes espaçamentos de plantio. O experimento foi montado em DBC (6 blocos) e com arranjo em parcelas sub-divididas. O modelo em questão é apresentado para delineamento em blocos completos em vários locais e várias observações por parcela – Método MHPRVG: Modelo 51 do software Selegen REML/BLUP. Modelo: Xb + Za + Wc + Tp + e. Onde, y é o vetor de dados, b é o vetor dos efeitos de repetição (assumidos como fixos) somados à média geral, a é o vetor dos efeitos genotípicos (assumidos como aleatórios), c é o vetor dos efeitos de parcela (aleatórios), p é vetor dos efeitos da interação genótipos x ambientes (aleatórios) e e é o vetor de erros ou resíduos (aleatórios). As letras maiúsculas representam as matrizes de incidência para os referidos efeitos. O experimento é em apenas um local, onde os diferentes espaçamentos testados foram considerados como locais/ambientes, blocos como repetições e as sub-parcelas foram consideradas como parcelas e desta forma não sendo aleatorizadas. Assim, ficaram os seguintes questionamentos: * Este modelo pode ser utilizado nesta situação? * Qual a implicação de utilizar este modelo para o referido experimento? * Qual a implicação da não aleatorização das parcelas? * É possível criar um modelo adequado para este experimento? Como processar os dados utilizando um modelo mais adequado no R? Agradeço a atenção! Lucas Guilherme Moura Oliveira Engenheiro Florestal - UFMG Mestrando Ciência Florestal - UFVJM Diamantina - MG (31) 96274633 / (31) 73521293 Lattes<http://lattes.cnpq.br/3332942776537748>

1 0

xyplot: incluir equação e r2 no gráfico
by Maurício Lordêlo 27 Mai '19

27 Mai '19

Caros, Preciso incluir as equações do modelo ajustado e o R2 no gráfico construído usando a função "xyplot" do pacote "lattice". Usei para o código mínimo reproduzível os dados do próprio pacote mas não deu certo. Agradeço quem puder auxiliar. Maurício library(lattice) str(iris) table(iris$Species) iris_setosa = subset(iris, Species =="setosa") mod1 = lm(Sepal.Length ~ Petal.Length, data=iris_setosa) iris_versicolor = subset(iris, Species =="versicolor") mod2 = lm(Sepal.Length ~ Petal.Length, data=iris_versicolor) iris_virginica = subset(iris, Species =="virginica") mod3 = lm(Sepal.Length ~ Petal.Length, data=iris_virginica) tp_mod1 = sprintf("%s = %.2f + %.2f%s", "Y", coef(mod1)[1], coef(mod1)[2], "X") r2_mod1 = cor(iris_setosa$Sepal.Length, fitted(mod1))^2 tp_mod2 = sprintf("%s = %.1f + %.2f %s", "Y", coef(mod2)[1], coef(mod2)[2], "X") r2_mod2 = cor(iris_versicolor$Sepal.Length, fitted(mod2))^2 tp_mod3 = sprintf("%s = %.1f + %.2f %s", "Y", coef(mod3)[1], coef(mod3)[2], "X") r2_mod3 = cor(iris_virginica$Sepal.Length, fitted(mod3))^2 tp_mod1;r2_mod1 tp_mod2;r2_mod2 tp_mod3;r2_mod3 ##Necessito em cada "tp_mod" incluir o r2 logo abaixo da respectiva equação### mypanel<-function(x,y,...) { panel.xyplot(x, y, ...) panel.abline(mod1) panel.abline(mod2) panel.abline(mod3) panel.text(4,7,labels=tp_mod1) panel.text(3,8,labels=tp_mod2) panel.text(2,7,labels=tp_mod3) } #O gráfico gerado: xyplot(Sepal.Length ~ Petal.Length | Species,panel = mypanel, data = iris, strip = strip.custom(par.strip.text = list(font = 3)))

1 0

Como extrair ano de data
by Emerson Cotta Bodevan 24 Mai '19

24 Mai '19

Prezados, boa tarde. Outra dúvida. Como extrair o ano de uma data? Por exemplo: de 03/12/2013 quero ficar apenas com 2013 Novamente agradeço qualquer ajuda. Atenciosamente, *Emerson*

2 2

[Loop em Análise de Tendência sobre série histórica]
by Renan Gon Ferreira 24 Mai '19

24 Mai '19

Prezados, boa tarde! Tenho numa planilha a série histórica de 56 estações pluviométricas. A primeira coluna refere-se à data (mês/ano) do evento e as demais colunas são as estações (cada coluna para uma estação da bacia hidrográfica). Quero analisar se existe tendência ou não na série de cada estação. Pra isso, estou isolando cada coluna (pra trabalhar separadamente com cada estação), excluindo os NA e aplicando o comando Pettitt.test (que pertence ao pacote "trend"). Mas apesar das colunas isoladas estarem com valores consistentes, sem problemas, o resultado final que estou obtendo não está fazendo sentido. Alguém conhece este pacote, ou outro de tendência com teste de Pettitt que poderia me indicar? Além disso, alguém teria um script com loop (utilização do for?) pra que eu possa obter o* p-value* e o *ano de mudança de tendência *referente a cada estação, caso contrário eu teria que realizar separadamente esta análise para 56 estações. A ideia seria fazer para isso pra primeira estação e automatizar a análise e exportação dos resultados para as demais estações. Muito obrigado, pessoal. *_______________* *Renan Gon Ferreira*

2 2

Converter dados
by Emerson Cotta Bodevan 24 Mai '19

24 Mai '19

Preados, boa tarde. Tenho um data frame em que uma das variáveis representa tempo de tartamento. Mas ela se encontra assim: *Tempo* 22m15d 10m10d ... Gostaria de converter e trabalhar apenas com meses (por exemplo): *Tempo* 22,50 10,33 ... Como fazer? Agradeço qualquer ajuda. *Emerson*

2 1

Duvida sobre tipos de categoria e análise estatística aplicável 2
by Tropicalnet 24 Mai '19

24 Mai '19

Prezado Cesar Obrigado pelo esclarecimento. Abordar a questão como você fez é a melhor resposta que podia ter tido. Sou da área de humanas, tenho trabalhado com o texto: PADUA, Jorge. Tecnicas de investigacion aplicadas a las ciencias sociales. 6a reimpressão - Cidade do Mexico: FCE, 1979. É um livro bem exaustivo para as tecnicas aplicáveis no âmbito das Ciências Humanas. Em seu capítulo IX o autor correlaciona categorias e testes estatisticos. Entretanto, dada a idade do texto, eu procurava algum embasamento mais atualizado. Agora acredito que continua um bom guia para orientar nossas atividades. Att. Alexandre

2 1