- R-br - DINF UFPR

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 11 Jul '18

11 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

Convite para a comunidade de R: submissões de propostas para o Mozilla Festival 2018
by Marcos Vital 29 Jun '18

29 Jun '18

Oi, pessoal da lista! A nona edição do Mozilla Festival (o MozFest) será realizada em Londres, sexta-feira, 26 de outubro a domingo, 28 de outubro. O festival é um evento super diverso, que tenta acolher pessoas do mundo todo que estejam de alguma forma interessadas em uma internet livre e acessível para todas as pessoas - isso passa por vários aspectos, que incluem iniciativas ligadas à Open Science e Open Source, de forma que a comunidade ligada ao R pode se interessar. No ano passado eu tive a chance de participar como facilitador de uma sessão sobre ensino de R, e neste ano sou um dos "wranglers", a equipe que organiza o festival. No momento a chamada para sessões está aberta, e vocês podem submeter uma proposta aqui: https://mozillafestival.org/proposals As inscrições de sessões ficam abertas até o dia 1o de Agosto - depois disso só será possível se inscrever como participante. Há recursos para que algumas propostas de sessões tenham ajuda de custo, no caso de candidatos que não tenham condições de viajar para Londres e se hospedar lá sem isso. Caso alguém queira saber mais, fiquem à vontade para entrar em contato comigo, ok? Abraços! Marcos -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

1 0

Eliminação de observações que não variam ao longo do período
by Alexandre Loures 28 Jun '18

28 Jun '18

Boa noite! Tenho uma base com a seguinte característica: Ano origem destino HS04 valor 1989 ARG BRA 0101 1 1989 ARG BRA 0102 0 1989 ARG BRA 0103 0 1989 ARG BRA 0104 1 2010 ARG BRA 0101 1 2010 ARG BRA 0102 1 2010 ARG BRA 0103 1 2010 ARG BRA 0104 1 Preciso eliminar as observações HS04 que não variaram ao longo do período. Ou seja, HS04 0101 e a HS04 0104 pois tanto no ano inicial quanto no final tinham valor 1. Desde já, muito obrigado pela atenção! -- *Alexandre Loures* Bolsista de Pós-doutorado Júnior do CNPq Centro de Ciências Sociais Aplicadas Universidade Federal da Paraíba www.ccsa.ufpb.br/ppge/ www.rodriguesloures.com

1 0

Faixa etária
by Elias Carvalho 27 Jun '18

27 Jun '18

Boa noite pessoal Qual a melhor forma de categorizar a variável "idade"? Existe uma pré-definição de faixa etária no Brasil ? Coisa do tipo: criança, jovem, adulto, idoso? -- *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"*

3 4

Equalizar duas colunas em dois data frames a fim de realizar posterior merge
by Jobenil Luiz Magalhães Júnior 26 Jun '18

26 Jun '18

Boa tarde pessoal, Estou tentando realizar o merge em dois dataframes, porém as respectivas colunas são ligeiramente diferentes: df1 <- (concatenado = c("1 _ 307", "22 _ 105", "12 _ 303", "6 _ 201"), var1 = 1:4, var2 = rnorm(4)) df2 <- (concatenar = c("500002_01 _ 307", "5001245_05 _ 303", "500033_22 _ 105", "501256_12 _ 303", "500345_06 _ 201", var3 = runif(5), var4 = rnorm(5)) Quero usar os 8 últimos caracteres da variável df2$concatenar e comparar com as da variável df1$concatenado. Noto que tenho dois problemas: 1) completar com 0 "zeros" à esquerda os valores à esquerda em df1$concatenado, quando forem menores do que 10. 2) Extrair os 8 caracteres à direita da variável df2$concatenar. Como posso equalizar as duas colunas? Agradeço desde logo qualquer ajuda que puderem dispensar. Atenciosamente Jobenil Júnior [image: Mailtrack] <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=sig…> Remetente notificado por Mailtrack <https://mailtrack.io?utm_source=gmail&utm_medium=signature&utm_campaign=sig…> 26/06/18 11:37:17

2 1