Julho 2018 - R-br - DINF UFPR

Re: [R-br] Pacote Chaid
by sznelwar＠uol.com.br 17 Jul '18

17 Jul '18

1 0

Pacote Chaid
by Edmar Caldas 16 Jul '18

16 Jul '18

Boa Tarde! Pessoal, não estou conseguindo executar esse código, podem me ajudar por favor? attach(credito1)str(credito1) install.packages("partykit")install.packages("CHAID", repos="http://R-Forge.R-project.org") library(partykit)library(CHAID) credit_rank<-as.factor(credit_rank)age <-as.factor(age)amex<-as.factor(amex)pay_week <-as.factor(pay_week)class <-as.factor(class) is.factor(credit_rank)is.factor(age)is.factor(class)is.factor(pay_week)is.factor(amex) > modelo <-chaid(credit_rank ~ ., data = credito1)Error: is.factor(x) is not TRUE

2 2

Eixo x
by Fabiana Goncalves Barbosa 16 Jul '18

16 Jul '18

Olá Pessoal Fiz um gráfico (plot) com os dados abaixo, porém no eixo x só apareceram alguns anos. Quero que apareçam todos os anos no eixo x. Como faço? qtidade=c(61,53,44,35,24,17,15,16,18,12,8,3,3,4,2,2,1,2,1,1,1,1,1) Anos=c(2018,2017,2016,2015,2014,2013,2012,2011,2010,2009,2008,2007,2006,2005,2003,2002,2001,1999,1998,1997,1992,1991,1981) plot(Anos,qtidade) Grata Fabiana

3 2

Modelos Mistos e Lattice (Grafico de predições)
by Fernando Souza 14 Jul '18

14 Jul '18

Caros colegas Desejo fazer um gráfico, semelhante ao retornado pela função augPred{nlme}, contendo os valores observados, os valores preditos pelo modelo fixo, e os valores individuais preditos. como esse exemplo (Só que para modelo não linear) https://rdrr.io/cran/nlme/man/plot.augPred.html por algum motivo a função augPred() não está funcionando com o meu modelo (Acredito por ser desbalanceado ) A função predict.nlme {nlme}, me retorna um dataframe com as informação que quero, porém como transformá-los em gráfico? Sei que é possível faze-lo através do xyplot{lattice} porém não estou conseguindo, devido a minha limitação no uso desse pacote. Alguém poderia me ajudar? segue CMR packages <- c("tidyr", "nlme""latticeExtra","lattice") #Instala pacotes necessários, se ja tiver instalado ativa-o for(p in packages){ if(!require(p,character.only = TRUE)) install.packages(p) library(p,character.only = TRUE) } ##==================================================================================================== ## **Preparação do banco de dados** ##==================================================================================================== DADOS.CURTO <- read.csv("https://www.dropbox.com/s/b4cckybgrcg86me/galinhas.csv?raw=1",head = TRUE) colnames(DADOS.CURTO) <- c("0","1","2","3","4","5","6","TRAT") DADOS.CURTO <-DADOS.CURTO[with(DADOS.CURTO,order(TRAT)),] DADOS.CURTO$ID <- c(1:93,94:280,281:370) DADOS<- DADOS.CURTO %>% gather('0','1','2','3','4','5','6',key ="SEMANA",value = "PESO", -TRAT,-ID) DADOS <- transform(DADOS, SEMANA = as.numeric(SEMANA),ID=factor(ID)) DADOS ##-------------------------------------------------------------------- GALINHA <- groupedData(PESO~SEMANA |ID,order.groups=FALSE,data = DADOS) ##------------------------------------------------------------------------------------------- modelo1.list <- nlsList(PESO ~ a * exp(-b * exp( -c * SEMANA)),start = c(a = 4.10,b = 4.53,c = 0.37 ),na.action=na.omit,data =GALINHA) ##--------------Definição do modelo de efeito Fixo-------------------------------- controle2<-nlmeControl(maxIter=300,msMaxIter = 300,niterEM= 100) modelo.nlme.0<-nlme(modelo1.list) modelo.nlme.2<-update(modelo.nlme.0,fixed=list(a~TRAT,b + c~1),start=c(4.10,0,0,4.77,0.34),control=controle2) ##------------------------------------------------------------------------------------------------------ controle3<-nlmeControl(maxIter=300,msMaxIter=1000) MM3.3<- update(modelo.nlme.2,weights = varExp(form = ~fitted(.)|SEMANA),control=controle3,correlation = corARMA(form = ~SEMANA |ID,q = 4)) summary(MM3.3) ##----------------------------------------------------------------- DADOS.ORIGINAL <- getData(MM3.3) predito<-predict(MM3.3) ##------------------------------------------------------------------ Sent from Mailspring (https://getmailspring.com/) the best free email app for work

2 1

Ajuda com Realocação de Inteiros
by Vanessa Felix 13 Jul '18

13 Jul '18

Pessoal, Será que podem me ajudar a pensar em um problema. Tenho um grupo de pessoas que preciso redistribuir de acordo com uma porcentagem, mas não estou conseguindo distribui-las de forma mais justa e que não desequilibre muito a relevância. Segue um exemplo da minha base de dados: Relevância Total Grupo Nova Redistribuição Arredondamento Arredondamento para baixo Arredondamento para cima 0,90% 80 0,721 1,0 - 1,0 1,80% 80 1,441 1,0 1,0 2,0 25,23% 80 20,180 20,0 20,0 21,0 27,93% 80 22,342 22,0 22,0 23,0 44,14% 80 35,315 35,0 35,0 36,0 SOMA 100% 80 79 78 83 relev=c(0.009009009,0.018018018,0.252252252,0.279279279,0.441441441) n=80 relev*n sum(relev*n) ceiling(relev*n) sum(ceiling(relev*n)) floor(relev*n) sum(floor(relev*n)) round(relev*n) sum(round(relev*n)) -- ------------------------------------------------------------------------------------------ Vanessa Felix do Nascimento Sergio Graduada em Estatística - ICE/UFJF Mestre em Saúde Coletiva - IMS/UERJ http://lattes.cnpq.br/2516291371035474 (21) 97557-3956 -------------------------------------------------------------------------------------------

2 1

Uniao por coluna - chave
by Edmar Caldas 12 Jul '18

12 Jul '18

Bom Dia! Preciso de outra ajuda. Estou tentando unir 3 arquivo usando a variavel chave id e não estou conseguindo, ele só uni e arquivos. chave<-merge(Apend_A,Apend_B,Apend_C, by.x=c("id"), by.y=c("id"), by.z=c("id"), all.x=TRUE, all.y=FALSE, all.z=FALSE) chave Edmar

1 0

Soma duas colunas pela função apply
by Edmar Caldas 12 Jul '18

12 Jul '18

Boa Tarde! Estou tentando somar duas colunas e coloca-lá na variável soma_debitos e não estou conseguindo. Não pretendo usar função! > Banco$soma_debitos <-apply(Banco$cartao_credito,Banco$emprestimos,2,sum)Error in match.fun(FUN) : '2' não é uma função, caracter ou símbolo Edmar

3 6

Dados DIC
by fúlvia cristina oliveira 10 Jul '18

10 Jul '18

Boa noite, estou rodando uns dados em DIC, onde tenho dois tratamentos (dietas) e algumas variáveis a saída saiu assim:O que tenho que fazer? > crd(TRAT, Hemoglobina, quali = FALSE, mcomp = "tukey", sigT = 0.05, sigF = 0.05)------------------------------------------------------------------------Analysis of Variance Table------------------------------------------------------------------------ DF SS MS Fc Pr>FcTreatament 1 59.806 59.806 14.495 0.00081149Residuals 25 103.150 4.126 Total 26 162.956 ------------------------------------------------------------------------CV = 25.89 % ------------------------------------------------------------------------Shapiro-Wilk normality testp-value: 0.07082243 According to Shapiro-Wilk normality test at 5% of significance, residuals can be considered normal.------------------------------------------------------------------------ Adjustment of polynomial models of regression------------------------------------------------------------------------Error in treat^2 : non-numeric argument to binary operator> Fúlvia Cristina Oliveira Bióloga - UFMS Mestre em Zootecnia - UEMSDoutoranda em Ciência Animal - UFMS(67) 99610-2812 / 99129-6612

2 1

Construção de gráfico em ggplot2
by Carlos Peixoto Neto 09 Jul '18

09 Jul '18

Olá, busco um auxílio para aprimorar um script em R. Tenho alguns erros de execução que apresento logo abaixo. No estudo, avalio o efeito da interação entre girinos, var. preditora "INT", no uso de recurso espacial através de duas variáveis resposta: deslocamento em Posição Vertical (PV) e Posição Horizontal (PH). Erros: - Na execução da função summarySE: "/Warning message:In qt(conf.interval/2 + 0.5, datac$N - 1) : NaNs produced" - Na execução de legend.title=element_blank(): Error: unexpected symbol in:"legend.title=element_blank() legend.title" #Importar dados. Pacotes para o gráfico pmed<-read.table("https://drive.google.com/open?id=1cET5cROSb-_D-cZYud3yKBHGR-qamind",header=T, sep=',') require(Rmisc) require(ggplot2) require(EnvStats) require(sciplot) require(dplyr) summarypmed <- summarySE(pmed, measurevar="PV", groupvars=c("ID","INT"),na.rm=T)# gráfico 1 graf1 <- ggplot(summarypmed, aes(x=INT, y= PV, fill=INT)+ facet_wrap(~INT) + geom_dotplot(binwidth=0.05,binaxis="y",stackdir = "center") + geom_errorbar(aes (ymin=PV.y-se,ymax=PV.y+se),width = 0.25,size=0.25)+ xlab("Tratamentos") + ylab(" Posição Vertical (cm)") + geom_text(aes(label = paste("N", "==",N,sep = "")), parse = TRUE, y=-0.15) + geom_point(aes(y=PV.y), size=1,show.legend = F) + theme_bw () + scale_fill_manual(values=c("grey85","grey30")) + theme(axis.ticks.x=element_blank(), axis.title.x=element_blank(), axis.text=element_text(size=7), axis.text.y=element_text(size=10), axis.text.x=element_blank(), axis.title.y=element_text(size=11), plot.title=element_text(size=12, hjust = 0.5), panel.grid=element_blank(), plot.title=element_text(vjust=2), legend.title=element_blank()) *Carlos Alberto Algarves Peixoto Neto* Biólogo | Universidade Federal do Maranhão - UFMA Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.do? id=K4303028A6 ___ "O amor de Deus: Cristo morrer em nosso favor quando ainda éramos pecadores." Romanos cap. 5: vers. 8 (NVI)

1 0