Agosto 2011 - R-br

[OFF- Dúvida] Linguagem de programação para trabalhar com estatística.
by Pedro Rafael 30 Ago '11

30 Ago '11

Pessoal qual seria uma outra boa linguagem de programação para se trabalhar com estatística? Estudo C já faz algum tempo e vi ser bastante interessante pelo fato de puder unir C com o R. Estava pensando em migrar para o C++, será que é mais interessante? Existe outras linguagens de programação que vocês indicam? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

9 18

Dúvida na função glht multcomp
by Marcelo Cardoso mello 30 Ago '11

30 Ago '11

Na função glht quando é solicitado comparações múltiplas com o argumento Tukey, o summary apresenta a saída abaixo. É feito um teste de tukey ou t? Porque a saída apresenta um valor de t? Multiple Comparisons of Means: Tukey Contrasts Linear Hypotheses: Estimate Std. Error t value Pr(>|t|) B - A == 0 4.667 1.477 3.160 0.10735 C - A == 0 20.583 1.732 11.887 < 0.001 *** D - A == 0 35.083 1.732 20.261 < 0.001 *** C - B == 0 15.917 1.732 9.192 0.00313 ** D - B == 0 30.417 1.732 17.566 < 0.001 *** D - C == 0 14.500 1.809 8.017 0.00471 ** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Adjusted p values reported -- single-step method)

2 1

Função não encontrada
by . . 30 Ago '11

30 Ago '11

Olá pessoal, ao rodar a função abaixo, obtenho o erro Error in f(x, ...) : could not find function "f2" ...mas f2 está definida logo acima de f... Alguma dica? Obrigado. func <- Vectorize(FUN= function(y, frac, rate, sad, samp="Poisson", trunc=0, ...){ f0 <- function(y,frac,n) { f1 <- function(y,frac,n){ dpois(y,frac*n) } dcom <- paste("d",deparse(substitute(sad)),sep="") dots <- c(as.name("n"),list(...)) f2 <- call(dcom,dots) f <- function(n){ f1(y,frac,n)*f2(n) } myintegrate <- function() { r <- 0 r1 <- 1 x1 <- 0 dx <- 20 while(r1 > 10e-500) { r1 <- integrate(f,x1,x1+dx)$value r <- r + r1 x1 <- x1 + dx } integrate(f,x1,Inf)$valu } myintegrate() } f0(y,frac,n)/(1-f0(trunc,frac,n)) },"y") func(10, 0.1, 0.1, exp)

3 4

Paste em data.frame
by Leandro Marino 30 Ago '11

30 Ago '11

Caros, Certa vez perguntei aqui na lista como separar uma variável em diversas colunas. Consegui uma solução bem interessante. Agora preciso de um caminho reverso. Alguma sugestão? ### separando as colunas x <- sprintf("%010d",1:30) taux <- max(nchar(x)) temp <- t(data.frame(lapply(x,substring,1:taux,1:taux))) rownames(temp) <- 1:nrow(temp) ### desejo algo do tipo: temp2 <- paste(temp[,1],temp[,2],temp[,3],temp[,4],temp[,5],temp[,6],temp[,7],temp[,8],temp[,9],temp[,10],sep='') No entanto tenho um numero de colunas variavel, alguma sugestão?! Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 R (21) 2103-9600 R.:236 ( leandro(a)cesgranrio.org.br  Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

4 3

legenda, renomear eixo x no boxplot
by Daniele Portela de oliveira 30 Ago '11

30 Ago '11

Caros colegas, mais uma vez eu recorro a vcs para resolver um problema com gráfico. Gostaria que aparecessem os valores de média e intervalos de credibilidade em cada caixa; segundo, no eixo x, aparecesse em forma de letra grega sigma^2 para representar os componentes de variancia σ2a, σ2p, σ2ep, σ2e. Abaixo tem um exemplo da minha tentativa mas não deu certo. Aceito sugestões!!! #exemplo require(coda) mat.pn=read.table("cadeia.txt", header=F) mat.pn=as.matrix(mat.pn) # falar que é uma matriz y=mcmc(mat.pn) vara= y[,1] varp= y[,3] varep= y[,4] vare= y[,2] library(grDevices) x.1=cbind(vara,varp,varep,vare) y.1=boxplot(x.1,range=0,ylab='Observações',font.main=3, font.lab=1, main = "PN", cex.main = 1.2, font.main = 3) axis(1, at = c(sigma^2(a), sigma^2(p),sigma^2(ep),sigma^2(e)), labels = expression(sigma^2(a), sigma^2(p),sigma^2(ep),sigma^2(e))) Abraços Daniele Portela de Oliveira Mestranda Produção Animal - UEM Cel: 44 9949 3525 67 9617 0965

2 1

OFF TOPIC - artigo antigo
by Pedro Emmanuel Alvarenga Americano do Brasil 29 Ago '11

29 Ago '11

Amigos de R, Achei um artigo antigo que gostaria de ler, no entanto nao tenho acesso ao texto completo. Alguem teria acesso ao texto completo que poderia mandra pra mim? http://onlinelibrary.wiley.com.ez68.periodicos.capes.gov.br/doi/10.1002/sim… ou http://www.ncbi.nlm.nih.gov/pubmed/2692111 , Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365 Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 3

Função eigen
by Junior Beleti 29 Ago '11

29 Ago '11

Boa Noite.. Estou analisando um script que realiza um loop e chama a função eigen. Com o pacote debug constatei que a execução desse script gasta muito tempo nessa função. Vi também que ela faz chamada para código fortran. Gostaria de saber se existe outra função que realiza o mesmo que a eigen. Att, -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

distribuições de probabilidade contínuas
by Guilherme Moraes Ferraudo 28 Ago '11

28 Ago '11

Boa tarde, por favor, alguém sabe me indicar alguma referência on-line que mostre as distribuições contínuas no R??? Muito obrigado. Abs -- Guilherme Moraes Ferraudo Jaboticabal/Campinas - SP http://lattes.cnpq.br/2096118558794430

2 2

Re: [R-br] Correlação Canônica!!
by Thiago De paula protásio 28 Ago '11

28 Ago '11

Ok Gustavo;Obrigado! Dei uma olhada no pacote vegan e não consegui exatamente o que eu preciso.... O meu interesse é apenas correlacionar dois grupos de várias da madeira e explicar as correlações existentes e fazer um agrupamento preliminar dos materiais genéticos. De qualquer forma, ainda não consegui descobrir o que a função plt.indiv () do pacote CCA está plotando. Não sei o que está sendo considerado como eixo X (Dimensão 1) e eixo Y (Dimensão 2)...Essa função plota os tratamentos....e é isso que estou precisando. Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

parcelas subdividas
by Gilson Sanchez 27 Ago '11

27 Ago '11

Pessoal mais uma duvida, estou rodando dados em parcelas subdivididas, mas como eu faço comparação das médias pelo teste Tukey? obrigado! -- Eng. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.

2 1