Outubro 2017 - R-br

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

Tratamento não ajustado x tratamentos ajustados no látice e DBC
by Euriann Yamamoto 19 Out '17

19 Out '17

Boa tarde. Alguém possui algum material que fale um pouco a respeito de Tratamento não ajustado x tratamentos ajustados e quais as maneiras de utilização? Estou confusão sobre qual tipo utilizar na minha análise.

1 0

Pacotes
by Edmar Caldas 18 Out '17

18 Out '17

Boa Tarde! é possível saber o nome de todos os pacotes da versão 3.4.1?

3 3

Alocar elementos em lista...
by andrebvs＠bol.com.br 18 Out '17

18 Out '17

2 1

Merge tabelas não padronizados.
by Ronaldo Câmara de Araújo 18 Out '17

18 Out '17

Estou com um bando de dados que possui um texto ao qual, em algum local do campo tem o texto que necessito para fazer um merge com outra tabela que tem essa informação. Exemplo: SERVICO DE LOCACAO DE 02 GALPOES ONDE ABRIGARA AS INSTALACOES DO CDD ATALAIA, CONFORME DESCRICOES CONSTANTES DO DOCUMENTO REQUISITANTE. SERVICO DE LOCACAO DE IMOVEL TIPO GALPAO COM NO MINIMO 500 M², PARA FUNCIONAMENTO DO ARQUIVO GERAL DA DR/SE TRATA-SE DE PROCESSO DE RENOVACAO DO CONTRATO DE LOCACAO DE IMOVEL NAO RESIDENCIAL, LOCALIZADO NA AV. COLETORA A, S/N - CJ. MARCOS FREIRE I - N. SRª DO SOCORRO/SE, CEP 49160-000, DESTINADO AO FUNCIONAMENTO DA AC SHOPPING PREMIO. TRATA-SE DE PROCESSO DE RENOVACAO DO CONTRATO DE LOCACAO DE IMOVEL NAO RESIDENCIAL, LOCALIZADO NA AV. RUA F, Nº 70 E 70A, QD. 04, LOT. CHALES, BAIRRO AEROPORTO, ARACAJU/SE, CEP 49037-130, DESTINADO AO FUNCIONAMENTO DA CDD ATALAIA. Na outra tabela tenho o campo CDD ATALAIA 1015 DR/SE 0003 N. SRª DO SOCORRO 2105 CDD ATALAIA 1015 A minha chave vai passar a ser o nome da unidade. Alguém tem alguma sugestão de como fazer isso no R? Grato. Ronaldo Câmara de Araújo [https://mail.google.com/mail/?ui=2&ik=a38f3a957b&view=att&th=1335ff14b71079…]

2 1