- R-br - DINF UFPR

substituindo valores (data.frame)
by Rogerio R. Silva 03 Nov '11

03 Nov '11

Caros, Tenho dois arquivos. Um deles, com milhares de linhas devido a replicação temporal, inclui três colunas (campanha, transecto, parcela). O segundo, contém informação sobre variáveis medidas em cada parcela (um único valor). Eu preciso substituir os dados das parcelas (primeiro arquivo) por valores das variáveis (segundo arquivo) condicionando a que transecto a parcela pertença. Um exemplo abaixo para uma variável, mas que não trabalha corretamente. coletas<-data.frame (campanha=c(gl(4,8,64,labels=c("C1","C2","C3","C4"))),transecto=c(gl(8,4,64,labels=c("T1","T2","T3","T4","T5","T6","T7","T8"))),parcela=c(gl(4,1,64,labels=c("P1","P2","P3","P4")))) dados<- data.frame (transecto=c(gl(8,4,32,labels=c("T1","T2","T3","T4","T5","T6","T7","T8"))), parcela=c(gl(4,1,32,labels=c("P1","P2","P3","P4"))), alt=rnorm(32)) x.temp<-coletas for (i in 1:length(levels(coletas$transecto))){ x.temp[which(coletas$parcela[i,]==dados$parcela[i,])] <- dados$Alt[i,] } x.temp Error in coletas$parcela[i, ] : incorrect number of dimensions Obrigado.- rogério

2 2

Carregar muitos arquivos txt
by Eder David Borges da Silva 03 Nov '11

03 Nov '11

Bom dia pessoal, venho mais uma vez pedir ajuda de vocês, O problema é o seguinte, tenho 7 pastas(Cenários), dentro dos mesmo tenho 949 arquivos TXT (1.5 Mb) cada, preciso carregar estes 7*949 arquivos e guardar em um array, tenho que carregar tudo pois tudo isso será plotado. meu código é o seguinte: cenarios <- c(1:7) con <- 1 for (i in cenarios){ setwd(paste('C:\\Users\\l01199\\Desktop\\Modelo\\Resultados\\Modelo_Cenario',i,'_Q1300',sep='')) files <- sort(list.files(pattern = "m000.TXT"))[1:50] tempo <- as.numeric(gsub("m000.TXT","",gsub("CON", "", files))) dim.a <- c(dim(read.table(files[1],head=T)),length(tempo)) if(con==1){ dados <- array(NA,c(length(cenarios),dim.a)) nomes.col <- names(read.table(files[1],head=T)) } con1 <- 1 for (j in files){ dados[con,,,con1] <- as.matrix(read.table(j,head=TRUE)) con1 <- con1+1 } print(con) con <- con+1 } tenho 4 gb de memoria na maquina, porem quando crio o objeto dados que é um array de tamanho c(length(cenarios),dim.a) [1] 7 11069 11 949 Não tenho memoria para tudo isso, alguém tem uma sugestão para guardar isso, dei uma olhada no pacote SOAR, mas como ele guarda e depois tenho que acessar novamente não sei se é a melhor saída(Posso estar equivocado em minha interpretação). Alguém possui um dica para contornar este problema. Att

2 5

Script para atualização de versão
by Marcio Nicolau 02 Nov '11

02 Nov '11

Olá Pessoas! Eu sei que já passaram várias formas para se atualizar o R para uma nova versão, mas gostaria de deixar aqui os meus $0,002 O processo é o seguinte: 1.) Crie uma pasta chamada UpdateR no home do usuário 2.) Copie os arquivos oldVersion.R e newVersion.R 3.) Na versão antiga (antes da atualização) execute o arquivo oldVersion.R no Passo (3) o script copia os arquivos (Rprofile.site e Rconsole) com as "configurações do usuário" - caso tenha outros arquivos, altere o script - para a pasta atual "UpdateR" e salva uma lista com as bibliotecas instaladas. 4.) Desinstale o R 5.) Instale a nova versão 6.) Execute o arquivo newVersion.R no Passo (6) o script copia os arquivos (Rprofile.site e Rconsole) e instala as bibliotecas da lista. Espero que seja util! Abraços e bom feriado! MN (@marcionicolau)

4 6

teste de médias para interação
by Daniele Portela de oliveira 01 Nov '11

01 Nov '11

Caros colegas gostaria de fazer um teste de médias para interação pelo pacote multcomp. Onde eu poderia encontrar um exemplo para distribuição de poisson ou binomial? Daniele Portela de Oliveira Mestranda Produção Animal - UEM Cel: 44 9949 3525 67 9617 0965

4 6

Duvida com métodos de plot pra função stripchart
by Augusto Ribas 01 Nov '11

01 Nov '11

Ola pessoal. Estava usando o comando stripchart para fazer alguns gráficos e gostaria de saber se alguém tem alguma forma simples de fazer um grafico similar a esse da figura do systat. http://www.intesoft.com/produits/tech/systat/images/density/pic7.jpg aqui tem outra imagem mas com um boxplot do lado, mas eu quero so as bolinhas alinhadas centralizadas mesmo. http://www.solutions4u-asia.com/PDT/SYSTAT/sigmaplot/SP1110/DotdensityGraph… Acredito que seria alguma coisa tipo o method="stack" mas que ele não alinha-se os pontos a esquerda mas centralizado. Não encontrei mais métodos fora o stack e o jitter. Espero ter conseguido me explicar. Bem é isso, segue um cmr de exemplo: #dados Densidade<-round(c(rnorm(30,2),rnorm(30,5))) Local<-rep(c("A","B"),each=30) #grafico stripchart(Densidade~Local,vertical=T,method="stack",offset=0.6,at=c(1.2,2.2),frame.plot=F,pch=16,cex=1.4) -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

4 5

Fwd: Windows binary of 2.14.0
by Leonard de Assis 01 Nov '11

01 Nov '11

-------- Mensagem original -------- Assunto: Windows binary of 2.14.0 Data: Tue, 01 Nov 2011 16:03:48 -0400 De: Duncan Murdoch <murdoch.duncan(a)gmail.com> Para: r-announce(a)lists.r-project.org <r-announce(a)r-project.org> To all Windows users: The binary build of 2.14.0 that was uploaded yesterday was missing Cairo support. I have rebuilt it, and uploaded a new copy. You can tell which one you have by running "svg()", which works on the new one, but not the old one. You can tell which one is on your CRAN mirror by looking at the last changed date: it is today (November 1) for the corrected build. Sorry for the inconvenience.... Duncan Murdoch _______________________________________________ R-announce(a)r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-announce

1 0

Problemas com split+lapply
by Rubem Kaipper Ceratti 01 Nov '11

01 Nov '11

Caros listeiros, estou com uma base de dados para a qual preciso calcular médias móveis para dois atributos com relação a certos níveis de fatores. Para isso, pensei em utilizar uma combinação de split+lapply, mas mesmo com apenas pouco mais de 126.000 observações o split do banco de dados leva uma eternidade (um dos fatores possui 1242 níveis). Posteriormente tentei utilizar a função ddply (pacote plyr), mas retorna um erro. Alguma sugestão? Segue um exemplo abaixo (base com menos observações que a original): P<-2000:2010 Rep<-letters[1:5] TF<-1:2 CC<-paste("H-",101:299,sep="") a<-expand.grid(P=P,R=Rep,TF=TF,CC=CC) a$y<-rnorm(nrow(a)) a$z<-runif(nrow(a)) library(zoo) ma.fun<-function(x,k){ mm.y<-rollmean(x$y,k) mm.z<-rollmean(x$z,k) cbind(x[-c(1:2),-c(5:6)],ma.y=c(mm.y),ma.z=c(mm.z)) } # Abordagem 1: Nem chega a completar o split a.1<-split(a,list(a$R,a$TF,a$CC),drop=T) a.2<-lapply(a.1,ma.fun,k=3) # Abodagem 2: Mensagem de erro: # Erro em do.ply(i) : # task 1 failed - "adicionando classe "factor" a um objeto inválido" # Além disso: Mensagens de aviso perdidas: # 1: <anonymous>: ... may be used in an incorrect context: ‘.fun(piece, ...)’ # # 2: <anonymous>: ... may be used in an incorrect context: ‘.fun(piece, ...)’ library(plyr) library(doSNOW) library(foreach) ncpu<-4 cl<-makeSOCKcluster(rep("localhost",ncpu)) registerDoSNOW(cl) clusterEvalQ(cl,c(library(plyr),library(zoo))) a.2<-ddply(a,c('R','CC','TF'),ma.fun,k=3,.parallel=T) Grato, Rubem

5 19

Linguagens para simulações Estatísticas e Matemáticas
by Pedro Rafael 01 Nov '11

01 Nov '11

Pessoal qual seria algumas linguagens interessantes fora o R boa para simulações estatística e matemáticas para se estudar? Qual dessas é melhor para trabalhar com banco de dados. Estou estudando o R mas queria estudar outra linguagem em paralelo. A linguagem Ox, O-Matrix, Scilab ou Octave cumpre bem esse papel? -- Saudações, Pedro Rafael Diniz Marinho. Estatístico - Secretaria de Estado da Saúde - PB.

1 0

Tabela de contingẽncia com freq. relativa: dúvida sobre validade teste qui-quadrado
by Jose Claudio Faria 01 Nov '11

01 Nov '11

Pessoal, Recorro aos mais especializados sobre a validade do teste qui-quadrado aplicados a uma tabela de contingência (dados e script abaixo). Os dados são de freq. relativa (cada linha fecha em 100%) de 5 estágios vegetativos (Seedling, Infant, Juvenile, Imature e Adult) do palmito juçara em 4 ambientes (MAR, ECO, MEF e MAZ). dad <- structure(list(Seedling = c(2.1664766249, 0.4830917874, 2.1466905188, 11.8236472946), Infant = c(53.2497149373, 55.0724637681, 65.831842576, 51.9038076152), Juvenile = c(33.7514253136, 33.8164251208, 26.4758497317, 28.0561122244), Imature = c(8.3238312429, 8.2125603865, 3.7567084079, 5.8116232465), Adult = c(2.5085518814, 2.4154589372, 1.7889087657, 2.4048096192)), .Names = c("Seedling", "Infant", "Juvenile", "Imature", "Adult"), class = "data.frame", row.names = c("MAR", "ECO", "MEF", "MAZ")) chisq.test(dad, correct=F)$observed chisq.test(dad, correct=F)$expected chisq.test(dad, correct=F) O objetivo da pesquisadora é saber se os ambientes podem ser considerados homogêneos (similares) ou não em relação as freq. dos estágios vegetativos. Havia aplicado o teste (foi significativo) e a pesquisadora comentou que ela havia lido que o teste do qui-quadrado não poderia ser aplicado em tabelas de contingência desse tipo. Sugestões e indicação de testes alternativos são bem vindas. Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - Prof. Pleno UESC/DCET/Brasil joseclaudio.faria at gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

5 11

Qual distribuição segue meus dados?
by Daniel Dantas 31 Out '11

31 Out '11

Boa noite pessoal, mais uma vez recorro a ajuda de vocês. Abaixo tenho frequencias de furtos e roubos que aconteceram mês a mês.A pergunta que tenho que resolver é a seguinte: Aplicar o teste de igualdade de distribuições para identificar qual é o tipo de distribuição que melhor descreve a distribuição de frequências. Alguém tem uma luz? ###### conjunto de dados (Esses dados são frequencias de furtos e roubos que aconteceram por mês, ou seja, temos 18 meses observados )x <- c(20,19,28,44,23,52,34,36,20,43,62,69,49,85,65,27,34,21)hist(x) Muito obrigado desde já e desejo a todos um excelente final de semana,Daniel

4 7