Outubro 2016 - R-br

[OFF -Tópic] Cotas em concursos para negros e índios!
by Andre Oliveira 12 Jun '49

12 Jun '49

Pessoal, tem um colegar meu que vou fazer um concurso ele é descendente de negro ele é pardo. Ele veio com a pergunta como provar que sou negro, pois, tem vagas reservada para negros e índios? Eu sinceramente não sei e não tenho ideia pois nunca tive esta curiosidade. O que de fato comprova sua raça? Algum documento? O que são usados em concursos para provar? obrigado André Oliveira Souza

5 6

Legenda eixo X
by Luciano F. Sgarbi 10 Jun '49

10 Jun '49

Prezados, Estou tentando centralizar a legenda do eixo X em uma figura com vários gráficos. Estou conseguindo centralizar a legenda da linha inferior como no exemplo, mas as linhas do meio não. par(mfrow=c(2,2)) plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 1",outer=TRUE) ### o problema está aqui, pois a legenda desce para a parte inferior da figura. plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 2",outer=TRUE)## Esse dá certo Alguém podia dar uma luz de como eu posso fazer isso? Abraços -- *Luciano F. Sgarbi* Mestrando em Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás, campus II Goiânia-GO 74001-970 Brazil

2 2

Teste de independência dos resíduos
by Maurício Sangiogo 02 Nov '16

02 Nov '16

Bom dia listeiros! Estou iniciando no R, e gostaria de pedir ajuda a vocês sobre uma dúvida... Como eu faço no R para testar o pressuposto da anova de "independência dos resíduos"? Eu li algumas citações ao uso do Durbin-watson test para análise de resíduos em regressões, porém não sei se o mesmo se encaixaria (ou se ao menos está correto) para análise de experimentos com fatores qualitativos (em DBC ou DIC). Também vi algo em plotar um gráfico dos resíduos, mas não entendi o como. .. Se alguma boa alma estatística puder me ajudar, muito obrigado, se não puderem agradeço igualmente. Att .................. Eng. Agr. Maurício Sangiogo TEL:(55) 96822030

4 3

comparação de médias - problema com os pressupostos
by Maurício Lordêlo 02 Nov '16

02 Nov '16

Caros, Preciso fazer um teste de comparação de médias mas estou tendo problemas com os pressupostos. Tentei uma transformação Box Cox mas não consegui êxito. Alguma sugestão? Fenois = c(337.311, 344.874, 342.353, 325.546, 333.950, 330.588, 328.067, 328.067, 318.824, 331.429, 333.950, 334.790, 336.471, 338.151, 342.353, 259.160, 252.437, 268.403, 265.882, 266.723, 287.731, 88.571, 88.571, 90.252, 41.513, 52.437, 49.076, 88.571, 88.571, 90.252, 64.202, 60.000, 61.681) Cor = factor(c(rep("ambar",6),rep("ambar_claro",3),rep("ambar",6),rep("ambar_claro",6),rep("branco",6), rep("extra_ambar_claro",3),rep("branco",3))) Fenois; Cor boxplot(Fenois~Cor,ylab=expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro")) tapply(Fenois, Cor, mean) tapply(Fenois, Cor, sd) tapply(Fenois, Cor, var) #Teste de homogeneidade de variâncias bartlett.test(Fenois ~ Cor) #Gráficos para avaliar a variável resposta qqnorm(Fenois) qqline(Fenois) hist(Fenois) boxplot(Fenois) #Transformação Box-Cox require(MASS) boxcox(Fenois~Cor, lambda = seq(-2, 2, 1/10)) boxcox(Fenois~Cor, lambda = seq(0.5, 1.5, 1/100)) #Optando por lambda = 0.6 Fenois2 = ((Fenois^0.6) - 1)/0.6 par(mfrow=c(1,2)) boxplot(Fenois~Cor,ylab=expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro")) boxplot(Fenois2~Cor,ylab=expression(paste("Fenois"," ","(",mg," ",kg^-1,")")), xlab="cor",names=c("Âmbar", "Âmbar Claro", "Branco", "Extra Âmbar Claro")) bartlett.test(Fenois2 ~ Cor)

3 2

Dúvida para gerar gráfico
by Paulo Henrique Pimenta 31 Out '16

31 Out '16

Olá Listeiros, Tenho uma dúvida bem simples da qual gostaria de pedir ajuda a vocês. Criei um código que gera as evapotranspirações pelo método de Penman, Thronthwait e Blaney-Criddle. A partir desses resultados gostaria de gerar gráficos em dois modos: O primeiro seria pegando o resultado de Penman mostrando sua evolução com os meses do ano O segundo é pegando os métodos de Thornthwait e Blaney-Criddle juntos mostrando, também suas evoluções no ano. Sempre que eu tento plotar algo aparece no meu IDE (RStudio) uma série de caixas. Com este conjunto de dados, os .csv que possuo não conisgo criar os dados desejados. Os dados estão em: http://s000.tinyupload.com/?file_id=49893669914373744636 O código está em: http://s000.tinyupload.com/?file_id=96211791048541413747 Abraços a todos e obrigado desde já. Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435.

3 3

Re: [R-br] Teste de independência dos resíduos
by Walmes Zeviani 31 Out '16

31 Out '16

Maurício, O primeiro gráfico até pode mostrar não independência quando ela é induzida por falta de ajuste. O código abaixo ilustra minha afirmação. library(lattice) library(latticeExtra) data(wtloss, package = "MASS") str(wtloss) xyplot(Weight ~ Days, data = wtloss) # Ajuste de um modelo linear. m0 <- lm(Weight ~ Days, data = wtloss) # Indica falta de ajuste e não independência. plot(m0, which = 1) # Mas eles exibem uma não independência devido à falta de ajuste. plot(acf(residuals(m0))) # Ajuste do modelo quadrático. m1 <- lm(Weight ~ poly(Days, degree = 2), data = wtloss) # Ausência de falta de ajuste. plot(m1, which = 1) # E sem correlação seria. plot(acf(residuals(m1))) start <- list(th0 = 185, th1 = -100, th2 = 140) with(start, { plot(Weight ~ Days, data = wtloss) curve(th0 + th1 * (1 - 2^(-Days/th2)), xname = "Days", add = TRUE, col = 2) }) # Ajuste do modelo não linear. m2 <- nls(Weight ~ th0 + th1 * (1 - 2^(-Days/th2)), data = wtloss, start = start) with(as.list(coef(m2)), { plot(Weight ~ Days, data = wtloss) curve(th0 + th1 * (1 - 2^(-Days/th2)), xname = "Days", add = TRUE, col = 2) }) # Valores preditos de cada modelo. pred <- with(wtloss, data.frame(Days = seq(min(Days), max(Days), length.out = 30))) pred$y0 <- predict(m0, newdata = pred) pred$y1 <- predict(m1, newdata = pred) pred$y2 <- predict(m2, newdata = pred) # Resíduos e valores ajustados de cada modelo. L <- lapply(list(m0, m1, m2), function(m) { data.frame(r = residuals(m), f = fitted(m)) }) names(L) <- 0:2 rf <- plyr::ldply(L, .id = "model") # Gráfico da curva predita com os valores observados. p1 <- xyplot(y0 + y1 + y2 ~ Days, outer = TRUE, type = "l", layout = c(NA, 1), data = pred) + layer(panel.points(x = wtloss$Days, y = wtloss$Weight)) # Gráfico dos resíduos vs valores ajustados. p2 <- xyplot(r ~ f | model, data = rf, layout = c(NA, 1), type = c("p", "smooth")) # Todos numa mesma "página". c(p1, p2, layout = c(NA, 2)) À disposição. Walmes.

1 0

Passar linhas de comandos do DOS para o script do R
by Adriele Giaretta Biase 28 Out '16

28 Out '16

Como executar um script do “R” com uma linha de comando do DOS passando parâmetros para o script. Tenho duas perguntas para essa dúvida: 1) Como executar uma linha de comando do DOS, passando parâmetros para o script do R? 2) Como reconhecer as linhas de comando com os parâmetros dentro do script do R? -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 98861-0619.

4 3

boxplot comparativos
by Elisa Henning 24 Out '16

24 Out '16

Boa noite Gostaria de construir boxplots múltiplos comparativos. *Exemplificando:* tenho dois arquivos com dados climáticos referentes a duas cidades. Gostaria de comparar, mês a mês, por meio de boxplots, as temperaturas destas duas cidades. Alguém poderia indicar algum material que auxilie a fazer estes gráficos? Desde já agradeço Elisa

4 4

não foi possível encontrar a função "%!in%"
by Elias Carvalho 23 Out '16

23 Out '16

Pessoal Executando uma função do pacote bn ocorreu o seguinte erro: check.bn.strength(strength) Error in check.bn.strength(strength) : não foi possível encontrar a função "%!in%" Alguém sabe porque dá esse erro no R ? -- *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

servidor cloud para rodar o R + shiny
by Marcio B 20 Out '16

20 Out '16

Olá pessoal boa tarde. Estou pesquisando um serviço de cloud para utilizar o R + Shiny. Minha duvida esta no custo e por isso gostaria de saber se alguém pode me ajudar para entender melhor: Quem é a empresa (cloud) que esta usando? AWS, Azure, Google ou IBM Qual o tipo de serviço para configurar o R + Shiny? Nestes servidores oferecem vários tipos, o que mais encontro na web é o EC2 - AWS, como vários tutoriais de instalação e configuração Se puder nos informar o uso e o custo deste serviço? Principalmente no uso de Memoria e Storage de dados, algum parâmetro para temos uma noção. Somos uma startup e queremos oferecer nos serviços utilizando o R + Shiny. Porem gostaríamos de um parâmetro para tomar a decisão de configurar este serviço e começar a utilizar. Espero que fui claro e toda a ajuda será bem vinda. Obrigado desde já. Marcio. -- Marcio B. Cel.: 19 99648.3949 tim 99798.0104 vivo skype: maborbaa http://www.linkedin.com/in/marciobar https://www.facebook.com/marciobar Empresa http://www.techagr.com/ ==== PAZ E BEM ====

6 8