Pequeno detalhe na cálculo do coeficiente de variacao! 100*sqrt(qms)/m Gilenio Fernandes Em 17 de março de 2014 12:00, <r-br-request@listas.c3sl.ufpr.br> escreveu:
Enviar submissões para a lista de discussão R-br para r-br@listas.c3sl.ufpr.br
Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request@listas.c3sl.ufpr.br
Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner@listas.c3sl.ufpr.br
Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..."
Tópicos de Hoje:
1. Coeficiente de determinação (cv): experimentos com mais de uma estimativa do erro (Jose Claudio Faria) 2. Re: Coeficiente de determinação (cv): experimentos com mais de uma estimativa do erro (walmes .) 3. Re: Coeficiente de determinação (cv): experimentos com mais de uma estimativa do erro (Jose Claudio Faria)
----------------------------------------------------------------------
Message: 1 Date: Mon, 17 Mar 2014 07:41:15 -0300 From: Jose Claudio Faria <joseclaudio.faria@gmail.com> To: r-br@listas.c3sl.ufpr.br Subject: [R-br] Coeficiente de determinação (cv): experimentos com mais de uma estimativa do erro Message-ID: < CAN+Emd8zJtNHF1CKM_6oQcmVV2UX1Mmn0iEuXdQjzgZSCcOF9w@mail.gmail.com> Content-Type: text/plain; charset=ISO-8859-1
Pessoal,
Alguém teria um bom algoritmo em R para determinar a estimativa do coeficiente de variação (cv - medida estatística básica e clássica sobre a precisão) dos experimentos com mais de uma estimativa do erro (por exemplo: parcelas subdivididas, subsubdividads, etc).
Tenho um feito apenas para experimentos com apenas uma estimativa do erro:
#=============================================================================== # Name : cv # Author : José Cláudio Faria/UESC/DCET # Date : 26/7/2006 15:25:07 # Version : v3 # Aim : Calculate coefficient variation from aov objects
#===============================================================================
# Arguments: # av aov object
cv <- function(av) { if(is.null(av)) stop('Please, check the parameter!') qmee <- with(av, sum(residuals^2) / df.residual) cv <- 100 * sqrt(qmee) / mean(av$fitted.values) return(round(cv, 2)) }
# Exemplo de uso:
av <- aov(Sepal.Length ~ Species, data=iris)
cv(av) [1] 8.81
Ab, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9966.9100 - VIVO 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
------------------------------
Message: 2 Date: Mon, 17 Mar 2014 10:29:17 -0300 From: "walmes ." <walmeszeviani@gmail.com> To: "r-br@listas.c3sl.ufpr.br" <r-br@listas.c3sl.ufpr.br> Subject: Re: [R-br] Coeficiente de determinação (cv): experimentos com mais de uma estimativa do erro Message-ID: <CAFU=Eka7XYyAJZbGpXc7dP2aP7NKzXS=uyPeqcgyRj-Hg= i9KQ@mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1"
Prezado Faria,
Embora eu seja contrário à interpretação ou uso, ou melhor, abuso do CV por certas pessoas e revistas, eu antigamente calculava assim, só para que tenha uma idéia. Sempre em cursos eu argumento contra o CV e meus códigos raramente o fazem aparecer. Isso porque muitas pessoas acham que um CV alto invalidam o experimento ou a análise. Argumento que não. Enfim, segue minha contribuição.
##-----------------------------------------------------------------------------
da <- expand.grid(a=gl(4,1), b=gl(3,1), c=gl(5,1)) da$y <- rnorm(da$a, mean=3)
m0 <- aov(y~a+b*c+Error(a:b), da) summary(m0)
a0 <- summary(m0) str(a0)
qms <- sapply(a0, function(i){ i[[1]]["Residuals",3] }) m <- mean(da$y)
100*qms/m
##-----------------------------------------------------------------------------
À disposição. Walmes.
-------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: < http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140317/0c3f06c6/attac...
------------------------------
Message: 3 Date: Mon, 17 Mar 2014 11:39:37 -0300 From: Jose Claudio Faria <joseclaudio.faria@gmail.com> To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Coeficiente de determinação (cv): experimentos com mais de uma estimativa do erro Message-ID: <CAN+Emd-xv=_jDu0GQKY3x8_gPx8ngw4zfdstE_2T_cEi0= imYw@mail.gmail.com> Content-Type: text/plain; charset=ISO-8859-1
Walmes,
Concordo com seus argumentos, mas considero o cv uma medida útil (comparativamente) da qualidade de um experimento frente a grupos similares. Uso no ensino de graduação em experimentação.
Obrigado pelo código! (Havia encontrado alguns mas não havia gostado do código de nenhum do que vi).
Ab, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9966.9100 - VIVO 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\
2014-03-17 10:29 GMT-03:00 walmes . <walmeszeviani@gmail.com>:
Prezado Faria,
Embora eu seja contrário à interpretação ou uso, ou melhor, abuso do CV por certas pessoas e revistas, eu antigamente calculava assim, só para que tenha uma idéia. Sempre em cursos eu argumento contra o CV e meus códigos raramente o fazem aparecer. Isso porque muitas pessoas acham que um CV alto invalidam o experimento ou a análise. Argumento que não. Enfim, segue minha contribuição.
##-----------------------------------------------------------------------------
da <- expand.grid(a=gl(4,1), b=gl(3,1), c=gl(5,1)) da$y <- rnorm(da$a, mean=3)
m0 <- aov(y~a+b*c+Error(a:b), da) summary(m0)
a0 <- summary(m0) str(a0)
qms <- sapply(a0, function(i){ i[[1]]["Residuals",3] }) m <- mean(da$y)
100*qms/m
##-----------------------------------------------------------------------------
À disposição. Walmes.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
código
mínimo reproduzível.
------------------------------
Subject: Legenda do Digest
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
------------------------------
Fim da Digest R-br, volume 39, assunto 17 *****************************************
-- Gilenio Borges Fernandes Professor Associado IV Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n - Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860