- R-br - DINF UFPR

Re: [R-br] [R_STAT] Seleção de Quadrados Médios
by Greice Laureano 23 Abr '11

23 Abr '11

Fernando, (agora mandei para a lista nova :P) O problema é que o comando aov() não aceita esse procedimento... Se eu faço uma Anova do fit ele da a seguinte mensagem de erro: > anova(fit) Erro em UseMethod("anova") : método não aplicável para 'anova' aplicado a um objeto de classe "c('aovlist', 'listof')" 2011/4/23 Fernando Henrique Toledo <fernandohtoledo(a)gmail.com> > > > Greice, bom dia. > > Acesse: > > anova(fit)[,3] > > Isso resultará no vetor dos quadrados médios! a 1ª coluna é dos graus de > liberdade, a segunda das somas de quadrados, 3ª dos quadrados médios ... > > PS: Mudamos de lista, mude p/ a "r-br", esta em breve será desativada! > > att, > FH > > 2011/4/23 greice.laureano <greice.laureano(a)gmail.com> > > >> >> Bom dia, >> >> Gostaria de ajuda para extrair os quadrados médios de uma Anova. >> O banco de dados está abaixo assim como o código para fazer a Anova. >> Muito Obrigado; >> >> Greice. >> >> Dados: >> >> sujeito aferidor tempo valor >> 1 1 1 -2 >> 1 1 2 0 >> 1 1 3 1 >> 1 2 1 2 >> 1 2 2 1 >> 1 2 3 1 >> 2 1 1 16 >> 2 1 2 16 >> 2 1 3 15 >> 2 2 1 12 >> 2 2 2 14 >> 2 2 3 13 >> 3 1 1 5 >> 3 1 2 6 >> 3 1 3 6 >> 3 2 1 4 >> 3 2 2 4 >> 3 2 3 4 >> 4 1 1 11 >> 4 1 2 10 >> 4 1 3 10 >> 4 2 1 9 >> 4 2 2 7 >> 4 2 3 8 >> 5 1 1 7 >> 5 1 2 8 >> 5 1 3 6 >> 5 2 1 5 >> 5 2 2 6 >> 5 2 3 6 >> 6 1 1 -7 >> 6 1 2 -8 >> 6 1 3 -8 >> 6 2 1 -9 >> 6 2 2 -10 >> 6 2 3 -9 >> 7 1 1 18 >> 7 1 2 19 >> 7 1 3 19 >> 7 2 1 17 >> 7 2 2 17 >> 7 2 3 17 >> 8 1 1 4 >> 8 1 2 5 >> 8 1 3 5 >> 8 2 1 5 >> 8 2 2 5 >> 8 2 3 5 >> 9 1 1 0 >> 9 1 2 -3 >> 9 1 3 -2 >> 9 2 1 -7 >> 9 2 2 -6 >> 9 2 3 -5 >> 10 1 1 0 >> 10 1 2 0 >> 10 1 3 -2 >> 10 2 1 1 >> 10 2 2 2 >> 10 2 3 1 >> 11 1 1 -3 >> 11 1 2 -2 >> 11 1 3 -2 >> 11 2 1 -4 >> 11 2 2 -3 >> 11 2 3 -3 >> 12 1 1 3 >> 12 1 2 -1 >> 12 1 3 1 >> 12 2 1 -1 >> 12 2 2 -2 >> 12 2 3 1 >> 13 1 1 7 >> 13 1 2 9 >> 13 1 3 9 >> 13 2 1 4 >> 13 2 2 4 >> 13 2 3 2 >> 14 1 1 -6 >> 14 1 2 -7 >> 14 1 3 -6 >> 14 2 1 -8 >> 14 2 2 -10 >> 14 2 3 -9 >> 15 1 1 1 >> 15 1 2 1 >> 15 1 3 0 >> 15 2 1 -2 >> 15 2 2 -2 >> 15 2 3 -3 >> 16 1 1 -13 >> 16 1 2 -14 >> 16 1 3 -14 >> 16 2 1 -12 >> 16 2 2 -12 >> 16 2 3 -12 >> 17 1 1 2 >> 17 1 2 1 >> 17 1 3 0 >> 17 2 1 -1 >> 17 2 2 0 >> 17 2 3 0 >> 18 1 1 4 >> 18 1 2 4 >> 18 1 3 3 >> 18 2 1 7 >> 18 2 2 6 >> 18 2 3 4 >> 19 1 1 -10 >> 19 1 2 -9 >> 19 1 3 -10 >> 19 2 1 -10 >> 19 2 2 -11 >> 19 2 3 -10 >> 20 1 1 8 >> 20 1 2 9 >> 20 1 3 8 >> 20 2 1 2 >> 20 2 2 8 >> 20 2 3 8 >> 21 1 1 7 >> 21 1 2 6 >> 21 1 3 7 >> 21 2 1 8 >> 21 2 2 7 >> 21 2 3 7 >> 22 1 1 -3 >> 22 1 2 -2 >> 22 1 3 -4 >> 22 2 1 -5 >> 22 2 2 -5 >> 22 2 3 -5 >> 23 1 1 -5 >> 23 1 2 -5 >> 23 1 3 -7 >> 23 2 1 -6 >> 23 2 2 -8 >> 23 2 3 -7 >> 24 1 1 5 >> 24 1 2 5 >> 24 1 3 5 >> 24 2 1 3 >> 24 2 2 4 >> 24 2 3 4 >> 25 1 1 0 >> 25 1 2 -1 >> 25 1 3 -1 >> 25 2 1 -4 >> 25 2 2 -3 >> 25 2 3 -4 >> 26 1 1 7 >> 26 1 2 6 >> 26 1 3 6 >> 26 2 1 4 >> 26 2 2 4 >> 26 2 3 4 >> 27 1 1 -8 >> 27 1 2 -8 >> 27 1 3 -8 >> 27 2 1 -10 >> 27 2 2 -11 >> 27 2 3 -10 >> 28 1 1 1 >> 28 1 2 1 >> 28 1 3 2 >> 28 2 1 1 >> 28 2 2 -1 >> 28 2 3 0 >> 29 1 1 -3 >> 29 1 2 -3 >> 29 1 3 -3 >> 29 2 1 -5 >> 29 2 2 -4 >> 29 2 3 -5 >> >> sj=as.factor(dados$sujeito) >> af=as.factor(dados$aferidor) >> tempo=as.factor(dados$tempo) >> >> fit <- aov(dados$valor ~ sj*af + sj:tempo:af + Error(sj/af)) >> summary(fit) >> >> > __._,_.___ > > <fernandohtoledo@gmail.com?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| através > de email<R_STAT@yahoogrupos.com.br?subject=Res%3A%20Re%3A%20%5BR_STAT%5D%20Sele%E7%E3o%20de%20Quadrados%20M%E9dios>| Responder > através da web<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJyNWprYWVzBF9TAzk3N…>| Adicionar > um novo tópico<http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmNDRncmJ1BF9TAzk3N…> > Mensagens neste tópico<http://br.groups.yahoo.com/group/R_STAT/message/16039;_ylc=X3oDMTM3N2Q4Y2Vi…>( > 2) > Atividade nos últimos dias: > > - Novos usuários<http://br.groups.yahoo.com/group/R_STAT/members;_ylc=X3oDMTJnNmk2bWJoBF9TAz…> > 2 > > Visite seu Grupo<http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmbDdmMWM3BF9TAzk3NDkwND…> > Não deixe de ir atrás da pessoa dos seus sonhos!<http://global.ard.yahoo.com/SIG=15me74fii/M=758712.14532720.14422460.129601…> > [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlc2RmczJyBF9TAzk3NDkwNDM1BGdycElkAz…> > Trocar para: SÃ³ Texto<R_STAT-traditional@yahoogrupos.com.br?subject=Mudar+Formato+de+Envio:+Tradicional>, > Resenha DiÃ¡ria<R_STAT-digest@yahoogrupos.com.br?subject=Envio+de+email:+Resenha>• Sair > do grupo <R_STAT-unsubscribe(a)yahoogrupos.com.br?subject=Sair+do+grupo> • Termos > de uso <http://br.yahoo.com/info/utos.html> > . > > __,_._,___ >

3 5

Sweave Mostrar Codigo R
by eder＠leg.ufpr.br 23 Abr '11

23 Abr '11

Pessoal, Preciso que no documento final seja mantido o codigo R bem como os resultados processados por ele, por exemplo o arquivo .Rnw abaixo: \documentclass[a4paper]{article} \usepackage{Sweave} \title{Estimativa via verossimilhança} \author{Éder David Borges da Silva} \begin{document} \maketitle \section{Distribuição} <<echo=TRUE>>= ###---------------------------------------------------------------------------### dbetabinom1 <- function(x, size, prob, theta, log=FALSE) { alpha <- prob / theta beta <- (1 - prob) / theta dens <- lchoose(size,x) + lbeta(x+alpha,size-x+beta) - lbeta(alpha,beta) if (!log) dens <- exp(dens) return(dens) } rbetabinom1 <- function(n, size, prob, theta) { # n é o número de elementos da amostra # size, prob e theta os parâmetros da densidade r <- function(x) sample(0:x[1], size=1, replace=TRUE, prob=dbetabinom1(0:x[1], size=x[1], prob=x[2], theta=x[3], log=FALSE)) if (length(n)>1) n <- length(n) return(new("numeric", sapply(data.frame(rbind(size, prob, theta, rep(1,n))), r))) } ###---------------------------------------------------------------------------### @ \end{document} Em minha cabeça eu lembrava que o echo=TRUE fazia isso, ou seja mantinha o codigo, porem isso não esta acontecendo, quando compilo não aparece nada do codigo R. Onde será que estou errando? Att

3 3

Removendo pontos dentro de um grid
by Alexandre dos Santos 22 Abr '11

22 Abr '11

Bom dia pessoal, Estou criando uma malha de pontos da seguinte forma: p0<- c(15,15)#Ponto inicial malha10<- expand.grid(xp=seq(p0[1], by=10, length.out=20), yp=seq(p0[2], by=10, length.out=20))#Criando um grid de pontos a cada 10 metros plot(malha10)# Visualizando a malha de postos Agora, gostaria das coordenadas dos pontos a cada 20m, à partir de p0, mas gostaria de manipular o objeto já criado malha10, sem ter que fazer um novo objeto expand.grid(... by=20), na verdade, gostaria de criar um novo objeto malha20, mas fazendo malha20=malha10 e após remover as coordenadas dos pontos do objeto malha10 que não satisfaçam a condição , tentei um outer() para descobrir as distancias entre todos os pontos e após impor uma condição =20, mas não estou conseguindo não, alguem poderia dar um help, Obrigado, Alexandre dos Santos Ingenieur forestier, Msc. INRA- Biostatistique et Processus Spatiaux (BioSP) Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29

3 5

Transformar image{graphics} em plot{graphics}
by Augusto Ribas 22 Abr '11

22 Abr '11

Ola. Estou com uma duvidade aqui. Vamos supor que eu tenha uma matriz assim com 0 e 1(de presença e ausencia) exemplo<-matrix(rep(0,100),nrow=10) while(sum(exemplo)/length(exemplo)<0.1){ x<-sample(x=1:10,1,replace =T) y<-sample(x=1:10,1,replace =T) exemplo[x,y]<-1 } exemplo e represente ela graficamente com o comando image: image(1:10,1:10,exemplo,col=(c(0,1)),xlab="",ylab="") Mas eu gostaria de pegar essa matriz de exemplo e fazer um grafico com o comando plot. Então eu preciso transformar uma matriz de 0 e 1 em coordenadas x e y pra usar com o comando plot para por exemplo preencher com pontinhos o seguinte plot: Mais ou menos assim: Pense nesse grafico: plot(0,0,xlab="",ylab="",xlim=c(0,10),ylim=c(0,10),type="n") abline(h=0.5:9.5,v=0.5:9.5) Em cada local da matriz que tiver um 1 eu queria que tivesse um pontinho no grafico. Ai eu uso o comando points para representar a matriz dentro daquele plot com pontinhos. Mas eu nao estou conseguindo visualizar como transformar uma matriz de 0 e 1 em coordenadas de pontos X e Y. Eu queria arranja um jeito de transformar qualquer matriz de presença e ausencia (0 e 1) em coordenadas para plotar pontinhos com o comando plot. Será que alguem tem uma ideia simples que resolvera todos meus problemas e que eu não consigo visualizar? Bom feriado a todos. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org/ http://dhost.info/augusto/ Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Re: [R-br] Removendo pontos dentro de um grid
by Ivan Bezerra Allaman 22 Abr '11

22 Abr '11

Alexandre!! Não sei se entendi direito, mais vai uma primeira tentativa. p0<- c(15,15)#Ponto inicial c0 <- seq(10,100,by=10) malha <- list() for(i in 1:length(c0)){ malha[[i]] <- expand.grid(xp=seq(p0[1], by=c0[[i]], length.out=20), yp=seq(p0[2], by=c0[[i]], length.out=20))#Criando um grid de pontos a cada 10 metros } malha[[2]] plot(malha[[2]]) plot(malha[[3]]) Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] Removendo pontos dentro de um grid
by Ivan Bezerra Allaman 22 Abr '11

22 Abr '11

Sacola 2!! M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Removendo pontos dentro de um grid
by Ivan Bezerra Allaman 22 Abr '11

22 Abr '11

sacola 1!! M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Removendo pontos dentro de um grid
by Ivan Bezerra Allaman 22 Abr '11

22 Abr '11

Alexandre! Não sei se entendi direito, mais vai uma primeira tentativa aí. p0<- c(15,15)#Ponto inicial c0 <- seq(10,100,by=10) malha <- list() for(i in 1:length(c0)){ malha[[i]] <- expand.grid(xp=seq(p0[1], by=c0[[i]], length.out=20), yp=seq(p0[2], by=c0[[i]], length.out=20))#Criando um grid de pontos a cada 10 metros } malha[[2]] plot(malha[[2]]) plot(malha[[3]]) Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Atualização do R no Linux
by Andre Oliveira 21 Abr '11

21 Abr '11

Prezados, atualizei o R aqui hoje, porém, estou tendo que instalar todos os pacotes novamente. Tem como eu resolver isto de forma mais suave? André Oliveira Souza

5 11

É possível aumentar a precisão numerica do R?
by Carolina Pena 21 Abr '11

21 Abr '11

Boa tarde pessoal, Estou trabalhando com operações que exigem uma precisão numerica muito grande. Por exemplo: (1.5692524257254742661/(10^41)) + 6.125125476214521152) No entanto, o R está gerando resultados inconsistentes. Alguém sabe se é possível aumentar a precisão do R? Se ninguém souber como fazer no R, mas souber de algum pacote que posso utilizar em C++ também seria muito útil. Obrigada, Carolina

3 2