Dezembro 2018 - R-br

Re: [R-br] Digest R-br, volume 96, assunto 8
by Edmar Caldas 10 Dez '18

10 Dez '18

Não funcionou a função prediction. O engraçado é quem nem mensagem de erro mostra. lixo$previsao <-predict(simples, newdata=lixo,interval = "prediction")lixo$previsao <-predict(simples, newdata=lixo,interval = "prediction", level=0.95) segue meu codigo, attach(waste)names(waste)simples <-lm(waste$wastee ~ restrnts)simplesajuste<-fitted(simples)waste$ajuste <-ajusteerro <- (waste$ajuste - waste$wastee)waste$erro<-erroprova <- waste$ajuste-waste$errowaste$prova <-provahead(waste) y<- (0.14685 + 0.01014 * waste$restrnts)waste$y<-y plot(wastee ~ ajuste) head(waste) #Prevendo em novos dados. library(readxl)lixo <- read_excel("Dieretorio/lixo.xlsx")View(lixo) attach(lixo) lixo$previsao <-predict(simples, newdata=lixo,interval="prediction",level=0.95)lixo eu consigo fazer a previsão, mas preciso do intervalo de confiança para previsão. Edmar nsegue Em segunda-feira, 10 de dezembro de 2018 12:00:18 BRST, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Re: IC para previsao (Elias T. Krainski) 2. Problema para plotar IC para modelo glm Gamma usando stat_smooth() no ggplot2 (ASANTOS) ---------------------------------------------------------------------- Message: 1 Date: Sun, 9 Dec 2018 16:54:45 +0000 (UTC) From: "Elias T. Krainski" <eliaskrainski(a)yahoo.com.br> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] IC para previsao Message-ID: <815218723.1276327.1544374485235(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Coloca "prediction" On Fri, Dec 7, 2018 at 22:39, Edmar Caldas por (R-br)<r-br(a)listas.c3sl.ufpr.br> wrote: Preciso de uma ajuda para criar o IC para ajuste do modelo de regressão linear simples. não dá nenhuma mensagem de erro, cria somente a previsao do modelo em novos dados. lixo$previsao <-predict(simples, newdata=lixo,interval="confidence",level=0.95)lixo Edmar_______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181209/9931f2c9/atta…> -------------- Próxima Parte ---------- Um texto embutido e sem conjunto de caracteres especificado foi limpo... Nome: Untitled Url: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20181209/9931f2c9/atta…> ------------------------------ Message: 2 Date: Mon, 10 Dec 2018 01:47:12 -0300 From: ASANTOS <alexandre.santos(a)cas.ifmt.edu.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Problema para plotar IC para modelo glm Gamma usando stat_smooth() no ggplot2 Message-ID: <8d326872-10b9-06f7-cc1b-564caf54f3e7(a)cas.ifmt.edu.br> Content-Type: text/plain; charset=utf-8; format=flowed Prezados Membros do R-br, Estou tentando plotar um modelo glm Gamma com stat_smooth() no ggplot2 sem sucesso. Minha ideia é representar o intervalo de confiança para um modelo onde realizei a junção de dois níveis para o fator Feature. Segundo exemplo fictício do meu problema: #Simulação de um banco de dados, onde a variável resposta Production tem distribuição Gamma set.seed(123) d<-NULL N<-50 d$Production <- rgamma(N,10) d$Feature <- ifelse(d$Production >7 & d$Production<10, c("green"),ifelse(d$Production>11, c("red"), c("blue"))) d$Temp<-rnorm(N,20,5) d<-as.data.frame(d) # # Ajusto um modelo de Gamma completo mG<- glm(Production~ Feature + Temp, family= Gamma, data = d) summary(mG) anova(mG,test="Chi") # Comparo os níveis do fator Feature library(multcomp) PW<-summary(glht(mG, linfct = mcp(Feature= "Tukey"))) PW cld(PW) # Faço de conta que os níveis green = blue e realizo um novo ajuste Feature2<-d$Feature levels(Feature2) levels(Feature2)[1]<-"blue&green" levels(Feature2)[2]<-"blue&green" levels(Feature2) d$Feature2<-Feature2 mG2<- glm(Production~ Feature2 + Temp, family= Gamma, data = d) # Crio os valores de predição pred.data = data.frame( Feature2<-d$Feature2, Temp<-d$Temp ) pred.data$Production = predict(mG2, newdata=pred.data, type="response") #Realizo o plot library("ggplot2") ggplot(d, aes(Temp, Production, colour = Feature2)) + geom_point() + geom_line(data=pred.data) + stat_smooth(method = "glm", formula = y ~ x, family = Gamma) # E tenho uma representação errada do intervalo de confiança do meu modelo mG2 porque não consigo fazer a função stat_smooth() considerar o ajuste mG2. Como posso resolver isso? Obrigado -- ====================================================================== 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/ ====================================================================== ------------------------------ Subject: Legenda do Digest _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ------------------------------ Fim da Digest R-br, volume 96, assunto 8 ****************************************

3 2

como carregar somente algumas colunas de um arquivo csv?
by Cinara Santos 08 Dez '18

08 Dez '18

Olá! Continuo as voltas com arquivos realmente grandes... Desta vez, está em csv. Ocorre que não preciso trabalhar com todas as colunas. Tendo conhecimento das colunas que me interessam (sei seus nomes - "header"), gostaria de saber da possibilidade de carregar somente tais colunas ao invés de carregar o arquivo inteiro, pra depois manter somente as colunas que me interessam.... É possível? Ainda gostaria de saber se alguém tem alguma sugestão pra eu conseguir abrir aos pouquinhos o arqeuivo do stata que deu estouro de memória [Digest R-br, volume 95, assunto 21].... Desde já, grande abraço! ------------------------ Cinara Santos

4 3

análise de correlação entre variáveis qualitativas e quantitativas
by Flávia de Campos Martins 07 Dez '18

07 Dez '18

Prezad@s, O objetivo do meu trabalho é analisar a correlação entre o conhecimento que as pessoas tem sobre a fauna local (variável dependente) e dados sócio-econômicos (sexo, idade, renda, escolaridade, tempo de residência no local). Eu tenho uma amostra de 105 informantes e as variáveis independentes foram organizadas em frequências e para montar a matriz de dados usamos números que representam as frequências das variáveis. Ex.: Idade: 0 a 10 anos - 1; 10 a 20 anos - 2... Como estamos trabalhando com a variável dependente qualitativa ordinal, e as amostras não apresentaram distribuição normal, fizemos um teste de correlação de Spearman. Entretanto, minha dúvida é se os informantes não se distribuem igualmente nas classes (categorias) de idade, de renda, de escolaridade isso não vai afetar o resultado? Para esclarecer: em relação à renda mensal eu tenho quatro classes: não possui renda fixa (0); até um salário mínimo (1); até dois salários (2) e até 3 salários. Só que 58% dos informantes estão na classe 1. Ou seja eu tenho distribuição desigual entre as classes, isso pode afetar meu resultado? Se sim, como resolver? Obrigada! Professora Adjunta Universidade de Pernambuco (UPE), *Campus* Petrolina BR 203, km 2, S/N Vila Eduardo 56328903 - Petrolina, PE- Brasil.

2 1

COMPARAR TRATAMENTOS
by Edson Lira 02 Dez '18

02 Dez '18

Boa noite, os dados abaixo referem-se a 10 tratamentos e ao estagio de uma doença. O grau 2 é o estágio maior. doencaTrat 0 1 2 Total GM1 3 2 1 6GM2 0 2 4 6GM3 1 3 2 6GM4 0 3 2 5GM5 0 4 2 6GM6 0 1 5 6GM7 1 1 4 6GM8 0 0 6 6GM9 0 4 1 5GM10 0 0 6 6 Fiz um teste de associação (Exato de Fisher) e verifiquei que ha uma associação entre o tratamento e o estagio da doença. Agora gostaria de verificar qual dos tratamentos é o mais eficiente. Posso testar também: G1xG2G3xG4G5xG6G7xG8G9XG10 Alguém tem alguma sugestão? [ ]'s.Prof. Edson Lira, Me Estatístico Manaus-Amazonas

3 3