- R-br - DINF UFPR

Maximizar plot window no RStudio
by Elias Carvalho 13 Abr '17

13 Abr '17

Bom dia pessoal Alguém sabe como posso maximizar a janela plot do RSstudio via linha de comando ? Explicando: Estou fazendo uma rotina para verificação de outliers baseado no boxplot então ao examinar os dados eu gero os gráficos abaixo, no entanto, se a janela estiver pequena demais os gráficos não aparecem e ainda surge uma mensagem Error in plot.new() : figure margins too large in R <http://stackoverflow.com/questions/12766166/error-in-plot-new-figure-margin…> : -- *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

Combinar multiplos datasets - Inputação Múltipla com MICE
by Elias Carvalho 11 Abr '17

11 Abr '17

*Boa noite pessoal* *Eu preciso tratar dados missing do meu dataset, li bastante e o MICE me pareceu uma boa ferramente para isso, então p**rossegui da seguinte maneira:* mice(df.factor.variable, m=5, maxit = 50, method = 'polyreg', seed = 500) *e depois para juntar os 5 datasets usei:* complete.factor.data <- complete(factor.variable.imputed, 'long') *Como tenho 1999 registros no meu dataset e foram gerados 5 conjuntos de dados imputados, o resultado foi de 9995 observações (1999 * 5).* *Alguém sabe me dizer se o procedimento está correto ou há uma maneira de unir os 5 datasets imputados em um só gerando um dataset de 1999 registros sem NAs?* *ps: vou usar esse novo dataset sem missing para aprender uma rede Bayesiana* -- *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"*

1 0

código para construir tabela com dados pareados
by Mauricio Cardeal 09 Abr '17

09 Abr '17

Boa tarde! Por favor, como faço para construir uma tabela como a do modelo do final do script, onde o conteúdo das 4 caselas são contagens do número de pares que combinam melhora e não melhora de 2 tratamentos aplicados simultaneamente no mesmo individuo (perna direita e perna esquerda), por exemplo ? # CMR: par <- c(1,1,2,2,3,3,4,4,5,5,6,6,7,7) # refere-se ao número do individuo trata <- c(1,0,1,0,1,0,1,0,1,0,1,0,1,0) # tratamento 1 (A) ou 0 (B) melhora <- c(1,0,1,0,1,0,0,1,0,1,0,0,1,1) # resultado do tratamento 1 (melhora) 0 (não melhora) tratac <- factor(trata) levels(tratac)<-list("A"=1,"B"=0 ) melhorac <- factor(melhora) levels(melhorac)<-list("+"=1,"-"=0) data.frame(par,tratac,melhorac) par tratac melhorac 1 1 A + 2 1 B - 3 2 A + 4 2 B - 5 3 A + 6 3 B - 7 4 A - 8 4 B + 9 5 A - 10 5 B + 11 6 A - 12 6 B - 13 7 A + 14 7 B + # COMO FAÇO PARA GERAR ESSA TABELA A PARTIR DOS DADOS ACIMA? # tratamento B # melhora melhora # + - #tratamento A melhora + 1 3 # melhora - 2 1 Obrigado, Maurício Cardeal UFBA

2 3

Suavização do gráfico
by Mateus Dias Nunes 08 Abr '17

08 Abr '17

Bom dia, Eu estou com um problema para suavizar meu gráfico. Eu tenho dados horários e desejo plotar individualmente cada data (no período das 06h as 18h). Entretanto o gráfico tem ficado com muito ruído devido as observações serem horárias e existirem falhas (NA) nas observações. E por isso eu gostaria de suavizar o gráfico com as médias horárias, porém não estou conseguindo implementar isso. Espero que alguém possa me ajudar a resolver. dado utilizado: https://www.dropbox.com/s/b80zzzeaesnh8qp/dado.txt?dl=0 library(openair) mydata<-read.table("dado.txt", h = T, sep = ",") mydata$direta <- as.numeric(as.character(mydata$direta)) mydata$date <- as.POSIXct(strptime(mydata$date, format = "%d/%m/%Y %H:%M:%S",tz = "Etc/GMT-1")) means <- aggregate(mydata["direta"], format(mydata["date"],"%d-%m-%Y %H:%M:%S"), mean, na.rm = TRUE) means[ ,"date"] <- seq(min(mydata[ , "date"]), max(mydata[ , "date"]), length = nrow(means)) Figura1<-plot(mydata$date[2:721], mydata$direta[2:721], type = "l", col="red",xlab = "minutes",main = "01 de Janeiro de 2007", ylab = "Direct Solar Radiation (W/m²)") Figura1 dev.copy(device=jpeg,filename="Brasilia_01dejaneiro2007.jpg", width = 800, height = 500) dev.off () desde já, agradeço. _______________________________________________________________________________________ MATEUS DIAS NUNES Postgraduate Program in Meteorology - PGMET - CPTEC / INPE Center for Weather Forecasting and Climate Studies / National Institute for Space Research Cachoeira Paulista - SP - Brasil FONE: +55 (53) 98112-5154 SKYPE: omsndomateus(a)hotmail.com ____________________________________________ ___________________________________________ <https://mailtrack.io/> Sent with Mailtrack <https://mailtrack.io/install?source=signature&lang=en&referral=nunes.mateus…> <#>

2 1

Duvida em estatística no R
by Tiago G. Santos 08 Abr '17

08 Abr '17

Boa tarde, Alguém sabe como faço para apresentar os resultados dos testes tukey ou dunnet como as letras que indicam significância estatística ..... -- Tiago Gonçalves Santos http://lattes.cnpq.br/7645251131511453 Doutorando em Fìsica pelo Instituto Tecnológico de Aeronáutica - ITA Laboratório de Caracterização de Dispositivos Semicondutores Instituto de Estudos Avançados - IEAv Trevo Coronel Aviador José Alberto Albano do Amarante, no. 1 12.228-001 São José dos Campos – SP – Brasil

2 1

Retornando de NaN para raiz cúbica de números negativos alocados em vetores
by marcos paulo 07 Abr '17

07 Abr '17

Olá prezados boa tarde. Alguém saberia me explicar porque ao concatenar números negativos e positivos em um vetor não é possível obter a raiz cúbica dos números negativos? Como procederia para contornar esse problema, mantendo os dados no vetor ? Exemplo: > vec<-c(27,-27,27,27,-27) > vec^(1/3) [1] 3 NaN 3 3 NaN > -27^(1/3) [1] -3 Situação desejada: [1] 3 -3 3 3 -3 À disposição Marcos Paulo Técnico Agrícola em Zootecnia Bacharel em Agronomia Mestrando em Produção Vegetal - Universidade Federal de Goiás Estagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio de Goiás-GO Contato: (62) 85075783 http://lattes.cnpq.br/4322347592884852

4 3

Ajuda: If else
by Edmar Caldas 07 Abr '17

07 Abr '17

Pessoal estou com outra dúvida, preciso criar uma coluna de inconsistência que funciona assim. Zero quando a condição está ok e 1 quando tem erro de inconsistência.minha logica é simples. se estcivil == solteiro e felicidade_casamento == NRA coloca Zero caso contrario coloque 1. erro <-1if (estcivil == "solteiro"){ + else if (Felicidade_casamento != "NRA")else erro <- 0 Edmar

3 2

APSIM.R
by Simone toni ruiz corrêa 07 Abr '17

07 Abr '17

Bom dia a todos, Gostaria de saber se alguém já trabalhou com o packege APSIM.R Trabalho com o modelo APSIM e vi que existe essa ferramenta no R, para auxilio de calibração e análise de sensibilidade. Gostaria de tirar dúvidas... Obrigada, Simone

1 0

- Pacote bipartite
by Elvira De Bastiani 06 Abr '17

06 Abr '17

Boa noite pessoal! Sou iniciante no programa R e estou utilizando o pacote Bipartite para minhas análises. Eu estou tentando rodar um modelo com 1.000 amostras de redes de interação que são criadas aleatoriamente a partir de um banco de dados. E agora durante a análise, começou a aparecer o seguinte erro: Error in file(file, ifelse(append, "a", "w")) : cannot open the connection In addition: Warning message: In file(file, ifelse(append, "a", "w")) : cannot open file 'web-1-1.pairs': Permission denied Aparentemente, durante a análise é gerado um arquivo que neste caso é o 'web-1-1.pairs', e o mesmo tem acesso negado, não finalizando a análise para as 1.000 redes de interação. Pelo que andei pesquisando, este erro pode acontecer quando: Tem um arquivo aberto em outro programa, como o Excel, enquanto tenta gravá-lo..ou se o diretório de trabalho não tenha acesso para gravação a esse diretório.. mas acredito que isso não seja o meu caso. Alguém conhece ou já teve esse erro, e poderia compartilhar como solucioná-lo ? Muito obrigada Elvira D'Bastiani Mestranda em Ecologia e Conservação Laboratório de Ecologia e Evolução de Interações (Leei) Universidade Federal do Paraná - UFPR Curitiba, Paraná, Brasil (41) 9 9653-2860

2 2

Ajuda: Append
by Edmar Caldas 04 Abr '17

04 Abr '17

Pessoal estou com outra duvida. Estou tentando unir dois banco com a função append. tenho arquivo A que contém a variável id, nome, lazer e arquivo B que contém a variável id, nome, trabalho. quero unir os 02 arquivos e ficar com id, nome, lazer, trabalho mesmo que alguns casos fiquem NULL. append <- rbind(LAZER_1,LAZER_2) dá erro de numero de colunas ou argumentos não coincidem. estou tentando fazer isso: LAZER_1$lazer <-NULLLAZER_1 LAZER_2$trabalho <-NULLLAZER_2 depois append <- rbind(LAZER_1,LAZER_2) continua com o mesmo erro? Edmar

4 4