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/attachment-0001.html>

------------------------------

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
*****************************************