- R-br - DINF UFPR

Ajuda para elaboração de boxplot em painel
by Gilenio Borges Fernandes 28 Out '17

28 Out '17

Prezad@s Estou precisando de auxilio para elaborar boxplots em painel pelo menos aproximadamente iguais aos dos arquivos pdf anexos. O banco de dados é do pacote nlme do R; pode ser acessado com o CMR library(nlme) data(RatPupWeight) head(RatPupWeight) names(RatPupWeight) Grato Gilênio Borges Fernandes Professor Associado IV (Aposentado) Professor Adjunto A (Substituto, 03/2015 a 03/2017) Universidade Federal da Bahia Instituto de Matemática e Estatística Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 Skype: gilenio.fernandes Lattes: http://lattes.cnpq.br/6764860618464860

1 0

Identificação de variáveis
by Elias Carvalho 25 Out '17

25 Out '17

Como posso identificar se uma variável é 'ordered factor nas seguintes condições: - o read.csv leu ela como int: - o read.csv leu ela como factor e não ordered factor: Obrigado

2 1

ajuda
by Edmar Caldas 25 Out '17

25 Out '17

Turma tem hora que funciona e hora não. alguém pode me ajudar? a variável datanasc no Rstudio está como <dttm> o que significa isso? Edmar

4 3

ANOVA in R using the function lapply ans terms
by Euriann Yamamoto 24 Out '17

24 Out '17

Hi everyone. In my dataset I´ve got 36 genotypes in 5 locations in partially balanced incomplete design 6x6, simple. We are analyzing productivity. I need to use adjusted means to do the orders analysis. TO obtain ANOVA individual to the 5 locations I run command: mi <- lapply(split(da, f=da$local), aov, formula=producao ~ repeticao/bloco + cultivar) lapply(mi, summary)#not adjusted for blocks But to transform cultivar adjusted for blocks (to obtain adjusted means) I need use command terms(formula,keep.order=TRUE) I DO NOT KNOW HOW TO USE THIS IN THE CORRECT ORDER. Something like this? m0 <- lapply(split(da, f=da$local), aov(terms(producao ~ repeticao/bloco + cultivar,keep.order = TRUE),da)) lapply(m0, summary)#adjusted to blocks Thanks Euriann

2 1

ENC: Látice parcialmente balanceado em vários locais
by Euriann Yamamoto 24 Out '17

24 Out '17

Walmes, eu consegui fazer pelo comando lapply a anova individual dos 5 locais. Só que agora não estou conseguindo colocasr os termos da fonte de variação em ordem: # ajuste do modelo de efeitos fixos para bloco dentro de repeticao #Anova Individual mi <- lapply(split(da, f=da$local), aov, formula=producao ~ repeticao/bloco + cultivar) lapply(mi, summary)#nao ajustado pra bloco #Arrumar ordem dos termos e tornar cultivar ajustado pra bloco mi <- lapply(split(da, f=da$local), aov, terms(formula=producao ~ repeticao/bloco + cultivar, keep.order = TRUE)) Ocorre um erro: Error in terms.default(formula, "Error", data = data) : no terms component nor attribute Creio que seja na fórmula. Qual seria a ordenação correta para poder utilizar esse comando terms? Att., Euriann L. M. Yamamoto Doutoranda em Agronomia Universidade Federal da Grande Dourados ________________________________ De: Euriann Yamamoto Enviado: quinta-feira, 19 de outubro de 2017 17:05 Para: Walmes Zeviani Assunto: RE: [R-br] Látice parcialmente balanceado em vários locais Walmes, segue abaixo o script que alterei para o látice parc. balanc.. Consta nele a ANOVA com tratamentos não-ajustados e bloco/repetição não-ajustados. Acontece que eu tenho 5 experimentos em látice 6x6, cada um em 5 municípios na safra agrícola 2008-09. Como automatizar a análise para que eu utilize um mesmo .txt contendo os locais, tratamentos, blocos, repetições e a variável? Script: rm(list=ls()) ls() library(lattice) library(agricolae) library(contrast) library(multcomp) library(doBy) library(latticeExtra) library(nlme) #----------------------------------------------------------------------- # lendo arquivos de dados da <- read.table("C:/ADSR/Tese/teste.txt", header = TRUE, colClasses = rep(c("factor", "numeric"), c(3, 1))) str(da) #----------------------------------------------------------------------- # ver os dados xtabs(~repeticao + bloco, da) xtabs(~repeticao + cultivar, da) xyplot(producao ~ cultivar, groups = bloco, data = da) xyplot(producao ~ cultivar | repeticao, groups = bloco, data = da) #----------------------------------------------------------------------- # ajuste do modelo de efeitos fixos para bloco dentro de repeticao m0 <- lm(producao ~ repeticao/bloco + cultivar, data = da) par(mfrow = c(2, 2)) plot(m0) layout(1) # #Teste de normalidade # shapiro.test(residuals(m0)) # #Teste de Bartlett # bartlett.test(residuals(m0)~cultivar, data=da) anova(m0) # cultivar nao ajustado para blocos #----------------------------------------------------------------------- #Só arrumar ordem dos termos m0 <- lm(terms(producao ~ repeticao/bloco + cultivar, keep.order = TRUE), data = da) anova(m0) Euriann L. M. Yamamoto Doutoranda em Agronomia Universidade Federal da Grande Dourados ________________________________ De: Walmes Zeviani <walmeszeviani(a)gmail.com> Enviado: terça-feira, 17 de outubro de 2017 19:53 Para: Euriann Yamamoto Assunto: Re: [R-br] Látice parcialmente balanceado em vários locais Não sei te responder com segurança. Testa e vê no que dá, rs. É bem provável que precise de adaptações, ainda que mínimas. Walmes.

3 2

Erro em script de apostila Embrapa Arroz e Feijão
by Euriann Yamamoto 24 Out '17

24 Out '17

Bom dia, Walmes. Estou rodando a análise 20 Análise de experimentos em vários locais (p.114) da apostila Embrapa Arroz e Feijão - Estatística experimental Aplicada 2011 e ocorreu o seguinte erro em anexo. Quando executo o comando mmqs, surge o seguinte erro: Error in qt((1 + conf.int)/2, dfdm) : Argumento não-numérico para função matemática Como resolver este problema visto que apenas copiei e colei todo o script. Att., Euriann

3 2

[OFF TOPIC] Dúvida diferente número de plantas por UA, mas, mesma área da UA
by Marcelo Laia 23 Out '17

23 Out '17

Caros colegas, Nesse final de semana me surgiu uma dúvida. Quando da análise de um experimento cuja unidade amostral possui o mesmo número de plantas, essas entram na análise normalmente. Por outro lado, tenho uma situação que, devido aos diferentes espaçamentos, a quantidade de plantas por UA é diferente. Para se fazer um teste de média eu posso utilizar as UA com números de plantas diferentes ou devo usar a média das parcelas? O experimento visa identificar qual é o melhor espaçamento de plantio para cada genótipo. Assim, a área amostral é sempre a mesma, mas, o N difere. O experimento conta com 6 blocos e os tratamentos foram sorteados dentro de cada bloco. Não sei me fiz entender, mas, lá vai minha dúvida que nada tem a ver com R. Se puder ajudar, agradeço muito. -- Marcelo

4 4

Perda de parcela - Lista R
by ana paula coelho madeira 20 Out '17

20 Out '17

Bom dia Walmes! A um tempo atrás você postou na lista de discussão do R uma sugestão acerca de experimentos em blocos com parcela perdida. Sugeriu, na ocasião, que a pessoa desse uma olhada no script http://www.leg.ufpr.br/~walmes/cursoR/cnpaf/cap08dbc-iso.R a partir do chunck number 6. Estou com o mesmo problema... tenho um DBC com perda de duas parcelas. No entanto, tentei seguir o script citado acima e aparece a seguinte msg de erro: > sim <- expand.grid(blc=gl(4,1, labels="bl"), trt=gl(5,1, labels="tr")) Error in as.character.factor(1:4) : malformed fator Você poderia me ajudar? Desde já Agradeço. PS: Peço desculpas por direcionar o email a você, porém como foi uma sugestão (script) dada por você, achei mais viável.

2 1

Preencher área entre duas linhas
by Delcio Rudinei Bortolanza 19 Out '17

19 Out '17

Olá! Segue um exemplo hipotético: #=========================================== x<-1:5 set.seed(4) y0<-sample(x) set.seed(2) y1<-sample(x)+5 matplot(x,data.frame(y0,y1),pch = c(0,1),col = c(1,1),las=1) xspline(x,y0,-0.5) xspline(x,y1,-0.5) #=============================================== Como preencher com cor entre as linhas suavizadas? Obs: a função polygon() não gera o resultado pretendido, pois une os pontos com linhas retas. Att., Delcio Bortolanza.

2 3

Efeitos diretos e indiretos no lavaan
by Elias Carvalho 19 Out '17

19 Out '17

Alguem sabe me dizer se é possível extrair os efeitos diretos e indiretos de uma variável específica de um grafo de Path Analysis gerado pelo Lavaan ? Explicando melhor, tennho um estudo avaliando o relacionamento de várias variáveis com um ou mais desfechos, então gostaria de criar uma tabela com esses valores e depois um grafico que permita analisar melhor os efeitos diretos e indiretos sobre a variável de desfecho. -- *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