- R-br - DINF UFPR

Organizando os dados (limpando sinais desnecessários)
by Tropidurus Torquatus 20 Mar '14

20 Mar '14

Tenho a matriz a baixo: > str(bmatrix1) 'data.frame': 2189 obs. of 16 variables: $ numMA.ACC : num 15 15 15 15 15 14 14 8 8 8 ... $ numMA.DON : num 11 11 11 11 11 10 10 5 5 5 ... $ numMA.Mol_ID : num 2 2 2 2 2 4 4 31 31 31 ... $ numMA.ROT : num 5 5 5 5 5 4 4 6 6 6 ... $ charMA.energy : Factor w/ 2189 levels "10.745595054136636 __ ",..: 825 840 855 928 957 768 807 54 66 165 ... ..- attr(*, "names")= chr "CMP1" "CMP2" "CMP3" "CMP4" ... A coluna "$ charMA.energy” É a minha coluna de interesse, como retirar os “__” e alterá-la para coluna numérica mantendo todos as casas decimais? No mais agradeço. Thales

2 2

Teste proporção
by Andre Oliveira 19 Mar '14

19 Mar '14

Boa noite, qual função faz o teste para proporções e calcula IC além da função binom.test(12,29, p = 0.5, conf.level = 0.95)? André Oliveira Souza

4 5

Re: [R-br] IC Variancia
by Robert Iquiapaza 17 Mar '14

17 Mar '14

#http://pages.stat.wisc.edu/~yandell/st571/R/append7.pdf var.interval = function(data, conf.level = 0.95) { df = length(data) - 1 chilower = qchisq((1 - conf.level)/2, df) chiupper = qchisq((1 - conf.level)/2, df, lower.tail = FALSE) v = var(data) c(l.inf=df * v/chiupper, l.sup=df * v/chilower) } plot(x=rnorm(48,34,7)) hist(x) var(x) var.interval(x) Em 15 de março de 2014 09:08, Andre Oliveira <andreolsouza(a)yahoo.com.br> escreveu: > Pessoal qual função do R calcula IC para variancia? > obrigado > > > _______________________________________________ > R-br mailing list > R-br(a)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.

2 1

combine_factor
by Dian Andrade 17 Mar '14

17 Mar '14

Como posso combinar NA com outra categoria a <- c("4", "9","1","1", "4",NA,"9",NA) table(a, useNA = "ifany") 1 4 9 <NA> 2 2 2 2 Neste exemplo como incluir o valor do NA em ? Obrigado

3 3

Re: [R-br] Digest R-br, volume 39, assunto 17
by Gilenio Borges Fernandes 17 Mar '14

17 Mar '14

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(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)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(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)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(a)gmail.com> > To: r-br(a)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(a)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(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)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(a)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/atta… > > > > ------------------------------ > > Message: 3 > Date: Mon, 17 Mar 2014 11:39:37 -0300 > From: Jose Claudio Faria <joseclaudio.faria(a)gmail.com> > To: r-br(a)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(a)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(a)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(a)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(a)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

1 0

Significado de graus de liberdade
by Luiz Roberto Martins Pinto 17 Mar '14

17 Mar '14

*Pessoal,* *Perdoem-me apresentar uma pergunta que nada tem a ver com o R, mas não participo de outra lista com tantos entusiastas por análise de dados e estatística.* *Alguém tem referência (livro ou artigo) que trate do 'Significado dos graus de liberdade'?* *Estou especialmente interessado em entender os graus de liberdade em uma interação entre dois fatores A e B, isto é, para AB, por exemplo na ANOVA. * *Sei que os graus de liberdade para esta interação é igual a (a-1)*(b-1), em que 'a' é o número de níveis de A e 'b' é o número de níveis de B. Mas gostaria de encontrar os resultados entendendo o numero total de observações e o número de parâmetros que estou realmente calculando.* *Utilizando o argumento de Chong Ho Yu "df is a function of both the number of* *observations and the number of variables in one's model".* *Então estou desejoso de entender "both the number of* *observations and the number of variables" para a interação AB.* *Se não for do interesse do pessoal da lista, peço àqueles que puderem ajudar que me enviem enviem e-mail off-lista para luizroberto.uesc(a)gmail.com <luizroberto.uesc(a)gmail.com>.* *Agradeço,* *Luiz Roberto* Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831

3 4

Coeficiente de determinação (cv): experimentos com mais de uma estimativa do erro
by Jose Claudio Faria 17 Mar '14

17 Mar '14

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 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 2

IC Variancia
by Andre Oliveira 15 Mar '14

15 Mar '14

Pessoal qual função do R calcula IC para variancia? obrigado

1 0

Quebrar eixo/escala em gráfico de barras agrupadas no R
by Rodolfo Souza 14 Mar '14

14 Mar '14

Olá pessoal, Alguém sabe como quebrar o eixo/escala em gráfico de barras agrupadas no R? Eu preciso fazer um gráfico e tem um grupo bem maior que os demais, e nesse caso passa ter um problema de escala. Segue o CMR cp1=c(6.81, 3.78, 3.16, 0.66, 4.27, 0.75, 5.64, 2.02, 0.99, 31.00) cp2=c(4.57, 0, 0, 0, 3.81, 0, 6.6, 1.27, 0, 140) cp3=c(0.46, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 0.00, 35.00) grups=c(1,2,3,4,5,6,7,8,9,10) dbh=data.frame(grups,cp1,cp2,cp3) dbhm=t(as.matrix(dbh[,-1])) colnames(dbhm)=dbh$grups barplot(dbhm, beside=T, las=1) O último grupo tem valores bem maior que os demais, por isso que fica ruim ver as variações de cada componente nos outros grupos. Eu queria fazer semelhante a esse ( http://uspace.shef.ac.uk/servlet/JiveServlet/showImage/2-12478-3980/gap_plo…) mas não sei como quando utiliza a função "beside=T". Grato pela atenção. Rodolfo Marcondes Silva Souza

3 2

Manipulação de data frame
by Rodrigo Muniz 14 Mar '14

14 Mar '14

Caros membros da lista, Desculpe por não enviar nenhum exemplo ou lista de comando. Tenho um dataframe com 1500 colunas e 700 linhas, eu preciso fazer uma razão entre todas as colunas, ou seja, C1/C2, C1/C3, C1/C4....C1499/C1500. Fazer essa divisão de forma manual apenas por indexação dos vetores é simples, porém, isso daria muito trabalho e acredito que deve haver um jeito mais fácil de fazer isso. Porém não tenho a menor ideia de como começar a fazer-lo. Por isso peço a ajuda de vocês, se alguém conhece alguma forma, algum exemplo, ou mesmo algum material que eu possa consultar, para que eu tenha esse "input" inicial. Desde já agradeço. Att... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

4 6