<div dir="ltr">Pequeno detalhe na cálculo do coeficiente de variacao! 100*sqrt(qms)/m<div>Gilenio Fernandes</div><div><br><div><br></div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 17 de março de 2014 12:00,  <span dir="ltr"><<a href="mailto:r-br-request@listas.c3sl.ufpr.br" target="_blank">r-br-request@listas.c3sl.ufpr.br</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Enviar submissões para a lista de discussão R-br para<br>
        <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
<br>
Para se cadastrar ou descadastrar via WWW, visite o endereço<br>
        <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
ou, via email, envie uma mensagem com a palavra 'help' no assunto ou<br>
corpo da mensagem para<br>
        <a href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a><br>
<br>
Você poderá entrar em contato com a pessoa que gerencia a lista pelo<br>
endereço<br>
        <a href="mailto:r-br-owner@listas.c3sl.ufpr.br">r-br-owner@listas.c3sl.ufpr.br</a><br>
<br>
Quando responder, por favor edite sua linha Assunto assim ela será<br>
mais específica que "Re: Contents of R-br digest..."<br>
<br>
<br>
Tópicos de Hoje:<br>
<br>
   1. Coeficiente de determinação (cv): experimentos com mais de<br>
      uma estimativa do erro (Jose Claudio Faria)<br>
   2. Re: Coeficiente de determinação (cv): experimentos com mais<br>
      de uma estimativa do erro (walmes .)<br>
   3. Re: Coeficiente de determinação (cv): experimentos com mais<br>
      de uma estimativa do erro (Jose Claudio Faria)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Mon, 17 Mar 2014 07:41:15 -0300<br>
From: Jose Claudio Faria <<a href="mailto:joseclaudio.faria@gmail.com">joseclaudio.faria@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: [R-br] Coeficiente de determinação (cv): experimentos com<br>
        mais de uma estimativa do erro<br>
Message-ID:<br>
        <<a href="mailto:CAN%2BEmd8zJtNHF1CKM_6oQcmVV2UX1Mmn0iEuXdQjzgZSCcOF9w@mail.gmail.com">CAN+Emd8zJtNHF1CKM_6oQcmVV2UX1Mmn0iEuXdQjzgZSCcOF9w@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
Pessoal,<br>
<br>
Alguém teria um bom algoritmo em R para determinar a estimativa do<br>
coeficiente de variação (cv - medida estatística básica e clássica<br>
sobre a precisão) dos experimentos com mais de uma estimativa do erro<br>
(por exemplo: parcelas subdivididas, subsubdividads, etc).<br>
<br>
Tenho um feito apenas para experimentos com apenas uma estimativa do erro:<br>
#===============================================================================<br>
# Name           : cv<br>
# Author         : José Cláudio Faria/UESC/DCET<br>
# Date           : 26/7/2006 15:25:07<br>
# Version        : v3<br>
# Aim            : Calculate coefficient variation from aov objects<br>
#===============================================================================<br>
<br>
# Arguments:<br>
# av  aov object<br>
<br>
cv <- function(av)<br>
{<br>
  if(is.null(av))<br>
    stop('Please, check the parameter!')<br>
  qmee <- with(av, sum(residuals^2) / df.residual)<br>
  cv   <- 100 * sqrt(qmee) / mean(av$fitted.values)<br>
  return(round(cv, 2))<br>
}<br>
<br>
# Exemplo de uso:<br>
> av <- aov(Sepal.Length ~ Species, data=iris)<br>
<br>
> cv(av)<br>
[1] 8.81<br>
<br>
Ab,<br>
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<br>
Jose Claudio Faria<br>
Estatistica<br>
UESC/DCET/Brasil<br>
joseclaudio.faria at <a href="http://gmail.com" target="_blank">gmail.com</a><br>
Telefones:<br>
55(73)3680.5545 - UESC<br>
55(73)9966.9100 - VIVO<br>
55(73)9100.7351 - TIM<br>
55(73)8817.6159 - OI<br>
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Mon, 17 Mar 2014 10:29:17 -0300<br>
From: "walmes ." <<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Coeficiente de determinação (cv): experimentos com<br>
        mais de uma estimativa do erro<br>
Message-ID:<br>
        <CAFU=Eka7XYyAJZbGpXc7dP2aP7NKzXS=uyPeqcgyRj-Hg=<a href="mailto:i9KQ@mail.gmail.com">i9KQ@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Prezado Faria,<br>
<br>
Embora eu seja contrário à interpretação ou uso, ou melhor, abuso do CV por<br>
certas pessoas e revistas, eu antigamente calculava assim, só para que<br>
tenha uma idéia. Sempre em cursos eu argumento contra o CV e meus códigos<br>
raramente o fazem aparecer. Isso porque muitas pessoas acham que um CV alto<br>
invalidam o experimento ou a análise. Argumento que não. Enfim, segue minha<br>
contribuição.<br>
<br>
##-----------------------------------------------------------------------------<br>
<br>
da <- expand.grid(a=gl(4,1), b=gl(3,1), c=gl(5,1))<br>
da$y <- rnorm(da$a, mean=3)<br>
<br>
m0 <- aov(y~a+b*c+Error(a:b), da)<br>
summary(m0)<br>
<br>
a0 <- summary(m0)<br>
str(a0)<br>
<br>
qms <- sapply(a0,<br>
              function(i){<br>
                  i[[1]]["Residuals",3]<br>
              })<br>
m <- mean(da$y)<br>
<br>
100*qms/m<br>
<br>
##-----------------------------------------------------------------------------<br>
<br>
À disposição.<br>
Walmes.<br>
<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140317/0c3f06c6/attachment-0001.html" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140317/0c3f06c6/attachment-0001.html</a>><br>

<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Mon, 17 Mar 2014 11:39:37 -0300<br>
From: Jose Claudio Faria <<a href="mailto:joseclaudio.faria@gmail.com">joseclaudio.faria@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] Coeficiente de determinação (cv): experimentos com<br>
        mais de uma estimativa do erro<br>
Message-ID:<br>
        <CAN+Emd-xv=_jDu0GQKY3x8_gPx8ngw4zfdstE_2T_cEi0=<a href="mailto:imYw@mail.gmail.com">imYw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset=ISO-8859-1<br>
<br>
Walmes,<br>
<br>
Concordo com seus argumentos, mas considero o cv uma medida útil<br>
(comparativamente) da qualidade de um experimento frente a grupos<br>
similares. Uso no ensino de graduação em experimentação.<br>
<br>
Obrigado pelo código!<br>
(Havia encontrado alguns mas não havia gostado do código de nenhum do que vi).<br>
<br>
Ab,<br>
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<br>
Jose Claudio Faria<br>
Estatistica<br>
UESC/DCET/Brasil<br>
joseclaudio.faria at <a href="http://gmail.com" target="_blank">gmail.com</a><br>
Telefones:<br>
55(73)3680.5545 - UESC<br>
55(73)9966.9100 - VIVO<br>
55(73)9100.7351 - TIM<br>
55(73)8817.6159 - OI<br>
///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\<br>
<br>
<br>
2014-03-17 10:29 GMT-03:00 walmes . <<a href="mailto:walmeszeviani@gmail.com">walmeszeviani@gmail.com</a>>:<br>
> Prezado Faria,<br>
><br>
> Embora eu seja contrário à interpretação ou uso, ou melhor, abuso do CV por<br>
> certas pessoas e revistas, eu antigamente calculava assim, só para que tenha<br>
> uma idéia. Sempre em cursos eu argumento contra o CV e meus códigos<br>
> raramente o fazem aparecer. Isso porque muitas pessoas acham que um CV alto<br>
> invalidam o experimento ou a análise. Argumento que não. Enfim, segue minha<br>
> contribuição.<br>
><br>
> ##-----------------------------------------------------------------------------<br>
><br>
> da <- expand.grid(a=gl(4,1), b=gl(3,1), c=gl(5,1))<br>
> da$y <- rnorm(da$a, mean=3)<br>
><br>
> m0 <- aov(y~a+b*c+Error(a:b), da)<br>
> summary(m0)<br>
><br>
> a0 <- summary(m0)<br>
> str(a0)<br>
><br>
> qms <- sapply(a0,<br>
>               function(i){<br>
>                   i[[1]]["Residuals",3]<br>
>               })<br>
> m <- mean(da$y)<br>
><br>
> 100*qms/m<br>
><br>
> ##-----------------------------------------------------------------------------<br>
><br>
> À disposição.<br>
> Walmes.<br>
><br>
> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código<br>
> mínimo reproduzível.<br>
<br>
<br>
------------------------------<br>
<br>
Subject: Legenda do Digest<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br>
<br>
------------------------------<br>
<br>
Fim da Digest R-br, volume 39, assunto 17<br>
*****************************************<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr">Gilenio Borges Fernandes<div>Professor Associado IV<br>Universidade Federal da Bahia<br>Instituto de Matemática<br>Departamento de Estatística<br>
Av. Adhemar de Barros, s/n – Ondina.<br>40.170-110 - Salvador - BA, Brasil<br>Tel.: (071)3283-6280/6336  Fax:  (071)3283-6276<br>URL: <a href="http://lattes.cnpq.br/6764860618464860" target="_blank">http://lattes.cnpq.br/6764860618464860</a><br>
</div></div>
</div>