- R-br - DINF UFPR

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

Sorteios aleatórios de um cadastro (AAS ou amostra estratificada)
by Andre Oliveira 02 Nov '16

02 Nov '16

Bom dia, estou buscando informações de como realizar sorteios aleatórios de um cadastro (AAS ou amostra estratificada) com objetivo de inferir para média ou para a proporção. Alguma dica de pacote para tal trabalho? obg André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

1 0

Barplot
by Emerson Cotta Bodevan 01 Nov '16

01 Nov '16

Prezados, bom dia. Estou tentando construir um barplot horizontal tal que as barras representem a frequencia de uso de dez medicamentos. Gostaria de colocar os nomes de todos os medicamentos a esquerda das barras. Não estou conseguindo colocar todos os nomes. Somente alguns deles aparecem. Alguma sugestão? Agradeço qualquer ajuda. Att., *Emerson*

4 4

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 23 Out '16

23 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