- R-br - DINF UFPR

Fatorial com 2 tratamentos adicionais
by Ze Henrique 21 Jun '19

21 Jun '19

Boa noite, gostaria de uma ajuda. To tentando fazer uma análise fatorial com tartamento adicional. Porém são 2 tratamentos adicionais, o positvo e outro o negativo. Alguem poderia me dar uma dica onde posso conseguir mais informaçções como fazer as analsie no r. Desde já agradeço a ajuda.

1 0

Problema para plotar intervalo de confiança no ggplot2 para modelo GLM de Poisson com junção de níveis
by ASANTOS 20 Jun '19

20 Jun '19

Prezados Membros, Estou tentando plotar sem sucesso o intervalo em um modelo GLM de Poisson onde fiz a junção de níveis para a variável categórica Gender. Alguém baseado no CRM abaixo poderia dar uma luz? #Pacotes library(ggplot2) library(dplyr) library(tidyverse) #Dados inventados com uma variável resposta (Consumption) e duas explicativas, sendo uma categórica e outra quali Consumption <- c(501, 502, 503, 504, 26, 27, 50, 56, 63, 60, 72, 93, 78, 43, 59, 70, 53, 80) Gender <- gl(n = 3, k = 6, length = 3*6, labels = c("Male", "Female","hermaphrodite"), ordered = FALSE) Income <- c(5010, 5020, 5030, 5040, 260, 270, 550, 560, 680, 690, 720, 550, 560, 680, 690, 720,500,512) df3 <- data.frame(Consumption, Gender, Income) df3 # GLM de Poisson fm1 <- glm(Consumption~Gender+Income, data=df3, family=poisson) summary(fm1) # ANOVA do modelo ajustado anova(fm1,test="Chi") #Comparo a variável Gender sort(tapply(df3$Consumption,df3$Gender,mean)) Gender2<-df3$Gender levels(Gender2) levels(Gender2)[2]<-"Fem_Her" levels(Gender2)[3]<-"Fem_Her" levels(Gender2) fm2<-glm(Consumption~Gender2+Income, data=df3, family=poisson) anova(fm1,fm2,test="Chi") # 0.7824 Female/Hermaphrodite são iguais então eu junto #Faço a predição sobre o modelo final e dos intervalos de confiança pred <- predict(fm2, type="response", se.fit = TRUE) df3 = cbind(df3, pred = pred$fit) df3 = cbind(df3, se = pred$se.fit) df3 = cbind(df3, ucl=df3$pred + 1.96*df3$se) df3 = cbind(df3, lcl=df3$pred - 1.96*df3$se) df3 = cbind(df3, Gender2) df<-df3 %>% dplyr::group_by(Income, Gender2) %>% dplyr::summarize(Consumption = mean(Consumption, na.rm = TRUE)) df<-as.data.frame(df) #Faço o plot usando o ggplot2 df3 %>% tidyr::gather(type, value, Consumption) %>% ggplot(mapping=aes(x=type, y=value, color = Gender2)) + geom_smooth(mapping=aes(ymin = lcl, ymax = ucl), stat = "identity") + geom_point(df,mapping=aes(x=Income, y=Consumption, color = Gender2)) + geom_line(mapping=aes(x=Income, y=pred)) # Obrigado, Alexandre -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 99686-6970 (VIVO) (+55) 65 3221-2674 (FIXO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 OrcID: orcid.org/0000-0001-8232-6722 - ResearcherID: A-5790-2016 Researchgate: www.researchgate.net/profile/Alexandre_Santos10 LinkedIn: br.linkedin.com/in/alexandre-dos-santos-87961635 Mendeley:www.mendeley.com/profiles/alexandre-dos-santos6/ ======================================================================

4 4

Experimento em faixas
by ana paula coelho madeira 15 Jun '19

15 Jun '19

Prezados colegas, Boa tarde! Estou analisando um experimento em DBC com dois fatores arranjados no esquema de faixas. Gostaria da ajuda de vocês para testar as hipóteses de normalidade e homocedasticidade. Desde já agradeço. Ana Paula

3 2

(sem assunto)
by Fabiana Goncalves Barbosa 12 Jun '19

12 Jun '19

Olá Pessoal Como faço um teste Kruskal Wallis para amostras não balanceadas?

2 1

Problema atualização.
by Diogo Jerônimo 08 Jun '19

08 Jun '19

Prezados, boa tarde. Acredito que seja ridícula, mas não achei no manual.. Atualizei meu R da versão 3.5.3 para a versão 3.6.0 (no "Ruindows"), e o programa está dando várias mensagens de erro em packages que antes funcionavam. Tentei retornar para a versão e o problema persiste. Alguém saberia o motivo, e a solução do problema? Obrigado!!!! Diogo JerônimoBacharel em Ciências Estatísticas - ENCE/IBGEMestre em Metrologia - PUC-Rio/PósMQICONRE: 8514 - SÉRIE Ahttp://lattes.cnpq.br/8996149312896520

2 1

Método de estimação do glm
by Elias Carvalho 04 Jun '19

04 Jun '19

Estou fazendo uma regressão logistica com o comando: modelo.rl1 <- glm(V1 ~ V2, data = data.to.work.training, family = binomial()) Como faço para saber o método de estimação utilizado? O resultado foi: summ(modelo.rl1) MODEL INFO: Observations: 24632 Dependent Variable: V1 Type: Generalized linear model Family: binomial Link function: logit MODEL FIT: χ²(1) = 478.25, p = 0.00 Pseudo-R² (Cragg-Uhler) = 0.08 Pseudo-R² (McFadden) = 0.07 AIC = 5954.13, BIC = 5970.35 Standard errors: MLE -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 1

Alternativas ao Rfacebook
by Jasmine Moreira 03 Jun '19

03 Jun '19

Olá pessoal, desenvolvi minha pesquisa de doutorado nos bons tempos do Rfacebook. Infelizmente, depois do escândalo da Cambridge Analytics, não há mais como coletar posts de páginas via API. Vi que o fb disponibiliza via parceiro um tal de CrowdTangle. Alguém já o utilizou? Poderia me dar umas dicas? Conhece alguma outra solução? Muito Obrigada, Jasmine

1 0

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