[R-br] Digest R-br, volume 39, assunto 17
Gilenio Borges Fernandes
gilenio em ufba.br
Segunda Março 17 16:21:41 BRT 2014
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 em listas.c3sl.ufpr.br> escreveu:
> Enviar submissões para a lista de discussão R-br para
> r-br em 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 em listas.c3sl.ufpr.br
>
> Você poderá entrar em contato com a pessoa que gerencia a lista pelo
> endereço
> r-br-owner em 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 em gmail.com>
> To: r-br em 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 em 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 em gmail.com>
> To: "r-br em listas.c3sl.ufpr.br" <r-br em 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 em 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 em gmail.com>
> To: r-br em 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 em 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 em 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 em 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 em 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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140317/496d7654/attachment.html>
Mais detalhes sobre a lista de discussão R-br