- R-br - DINF UFPR

duvida integers e as.double
by Vitor Aguiar 30 Mar '12

30 Mar '12

Olá Benilton. Por favor, gostaria de tirar uma dúvida a respeito de um código que você me passou há uns dias. É esse código abaixo. grupos <- unique(gsub("(.*)\\.\\d{1}$", "\\1", names(Pop))) tabelaPorGrupo <- function(grp, dat){ cols <- grep(grp, names(dat)) Nomes <- levels(factor(c(dat[,cols[1]], dat[,cols[2]]))) table(factor(Pop[,cols[1]], levels=Nomes), factor(Pop[,cols[2]], levels=Nomes)) } tabelas <- lapply(grupos, tabelaPorGrupo, Pop) names(tabelas) <- grupos tabelas post: http://r-br.2285057.n4.nabble.com/R-br-problema-com-for-loop-td4465197.html Tenho que fazer muitos cálculos com os valores gerados por esse código e, algumas vezes, quando estou trabalhando com valores altos, há warning messages dizendo "NAs produced by integer overflow". Um outro programador me disse que eu poderia resolver isso estocando meus valores como doubles e não como integers, e que eu teria que fazer alguma modificação nesse código para esse fim. Você concorda e poderia me ajudar com isso? obrigado, Vitor Aguiar

2 5

Ajuda - Arima Sazonal
by Vinicius Macedo Magalhães 30 Mar '12

30 Mar '12

Fala, galera. Primeira vez usando este fórum. Tenho tido problemas para rodar meu código para previsão sazonal via ARIMA. Abaixo segue parte do código. >s.ts Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 2004 118 143 169 158 143 135 135 140 135 125 120 120 2005 143 158 180 180 150 150 153 148 150 145 145 143 2006 165 180 190 168 163 163 168 175 180 160 155 150 2007 152 163 185 163 158 175 175 178 168 168 160 160 2008 155 170 195 195 195 185 185 165 165 153 154 159 ##### Início de ARIMA ##### linhaX <- 3 dados <- matrix() nome <- names(db) i=4 get.best.arima <- function (x.ts, maxord=c(1,1,1,1,1,1)) { best.aic <- 1e8 n <- length(x.ts) for (p in 0:maxord[1]) for (d in 0:maxord[2]) for (q in 0:maxord[3]) for (P in 0:maxord[4]) for (D in 0:maxord[5]) for (Q in 0:maxord[6]) { fit <- arima(x.ts, order=c(p,d,q),seasonal=list(order=c(P,D,Q),period=12), method = "CSS") fit.aic <- -2 * fit$loglik + (log(n) + 1) * length(fit$coef) if (fit.aic < best.aic) { best.aic <- fit.aic best.fit <- fit best.model <- c(p,d,q,P,D,Q) } } list(best.aic,best.fit,best.model) } s.ts <- ts(db[,i],start=2004, fre=12) best.arima.ss <- get.best.arima(log(s.ts),maxord = c(2,2,2,2,2,2)) best.fit.ss <- best.arima.ss[[2]] best.arima.ss[[3]] pdf(paste("AR.", i-4, ".pdf", sep=""),width=7,height=5) ts.plot( cbind( window(s.ts, start=2004),exp(predict(best.fit.ss,24)$pred)),col=c("black","green")) #par(new=TRUE) #ts.plot(s.ts,best.fit.ss,col=c("black","red")) legend("bottomright",c("Real","Previsão_Passado","Previsão_Futuro"),col=1:3,lty=1,cex=0.7,merge = TRUE,bg = 'gray90') dev.off() report<-function(Actual,Forecast,Net,R2) { result <- new.env() if (R2==0) R2<-1-(sum((Forecast-Actual)^2)/sum((Actual-mean(Actual))^2)) result$Rsquare<-R2 result$MAPE <- 100*(sum(abs((Actual-Forecast))/Actual))/nrow(Forecast) ## Mean Absolute Percent Error result$CFE <- sum((Actual-Forecast)) ## Cumulative Forecast Error result$ME <- result$CFE/nrow(Forecast) ## Mean Error result$MSE <- sum((Actual-Forecast)^2)/nrow(Forecast) ## Mean Squared Error result$RMSE <- sqrt(result$MSE) ## Root Mean Squared Error result$MAD <- (sum(abs(Actual-Forecast))/nrow(Forecast)) ## Mean Absolute Deviation result$SUMMARY <- summary(Net) return(result) } ## Fim da função report pred.past<- db[,colY] previsao <- predict(best.fit.ss) forecast <- as.matrix(pred.past[1:length(pred.past)]) reportAR <- report(db[,colY],forecast,best.fit.ss,0) dados <- rbind(dados,"Produto:"=nome[colY],R2=reportAR$Rsquare,MAPE=reportAR$MAPE, CFE=reportAR$CFE,ME=reportAR$ME,MSE=reportAR$MSE,RMSE=reportAR$RMSE,MAD=reportAR$MAD) write.table(dados,sep=",",file="ReportAR.csv",) } ## Fim do for return(dados) ##### Fim de ARIMA ##### -- Att., Vinicius Macedo Magalhães (21) 9584-1533

2 1

loops no R
by Humberto Hazin 30 Mar '12

30 Mar '12

Prezados, Venho mais uma vez pedir ajuda a voces para resolver um problema! baixei uma série de images de satelite, por ano e mes (2004 a 2010) totalizando 163 arquivos ate ai tudo bem! Por termos de comodidade gostaria de criar um loop para nao ter que digitar essas 163 linhas ou mais para extrair latitude, longitude e SST. conforme o codigo abaixo: hdf4to5("S20040012004031.hdf");m0401 <- satinGet("S20040012004031.hdf5",lats=c(10, -30), lons=c(-20, -50), itype="oceancolor") m0401<-(simplifySat(m0401, extent=1));m0401<-extract(m0401);m0401<-data.frame(m0401);m0401$mes<-1; m0401$year<-2004 hdf4to5("S20040322004060.hdf");m0402 <- satinGet("S20040322004060.hdf5",lats=c(10, -30), lons=c(-20, -50), itype="oceancolor") m0402<-(simplifySat(m0402, extent=1));m0402<-extract(m0402);m0402<-data.frame(m0402);m0402$mes<-2; m0402$year<-2004 Alguem poderia me ajudar! Nao coloquei o arquivo pois é muito grande quase 36gigas porém caso precisem posso colocar o parte do codigo para baixarem um ou mais arquivos e testarem Humberto

3 9

Tobit espacial
by Diego andre 30 Mar '12

30 Mar '12

Olá pessoal, alguém sabe como eu posso estimar um tobit espacial com o R? Obrigado. -- Diego André Doutorando em Economia - CAEN/UFC cel: 85-87576600 lattes: http://lattes.cnpq.br/6480130040427049<http://www.orkut.com.br/Interstitial?u=http://lattes.cnpq.br/64801300404270…>

2 1

BRugs_windows 7
by Daniele Portela de oliveira 29 Mar '12

29 Mar '12

Caros colegas eu troquei de computador e o novo veio com windows 7 64 bits. Depois disso não consegui mas usar o pacote BRugs, sempre dá erro quando peço a library. Alguém já teve problema semelhante e poderia me ajudar a solucioná-lo? Sem mai, agradeço > >Daniele Portela de Oliveira > >Doutoranda em Zootecnia >UNESP - Jaboticabal-SP >Cel: 16 8221 6815 > >

2 1

Problemas com fixed breaks
by Julianna Trindade 29 Mar '12

29 Mar '12

Boa Tarde, Gostaria de uma solução para o FixedBreaks obedecer os limites superiores e inferiores das minhas classes de intervalos. No entanto, ele deve associar algum método estatístico e burla minha configuração. Eis o código que eu utilizo: brks <- classIntervals(dados$ano_2006, n = 7,style="fixed", intervalClosure="right", fixedBreaks=c(-2, -0.5, 0, 20, 40, 60, 80, 10000), main="Fixed") brks<- brks$brks Eis o breaks gerado pelo R: c(-11.1, -0.5, 16.665, 41.665, 58.335, 83.335, 110.1) Alguém pode me ajudar, por favor? -- Julianna Trindade "Exalta-te, Senhor, na Tua força! Nós cantaremos e louvaremos o Teu poder." Salmos 21; 13

2 1

Problemas para modificar vetor
by Augusto Ribas 29 Mar '12

29 Mar '12

Ola Pessoal, venho aqui pedir uma ajudinha. Eu tenho um conjunto de dados que são varios nomes de especies. Ai eu resumo eles em siglas para fazer um grafico. Com uma dica anterior do Benilton e do Paulo eu aprendi a resumir os nomes em 3 letras, a primeiro do genero e as duas primeira da especie. Meu problema é o seguinte eu tenho especies como: [1] "Aplectana crossodactyli" "Aplectana crucifer" Que as duas viram "Acr" O que eu gostaria de fazer é qd isso acontece, que são poucos casos dar o nome de Acr1 e Acr2, somente nesses casos. Seguem embaixo como estão os dados e como eu estava fazendo as coisas #dados nomes.lin<-c("Anuracanthorhynchus tritaxisentis", "Aplectana crossodactyli", "Aplectana crucifer", "Aplectana delirae", "Aplectana lopesi", "Aplectana membranosa", "Aplectana meridionalis", "Aplectana micropenis", "Aplectana pintoi", "Aplectana vellardi", "Capillaria recondita", "Catadiscus cohni", "Catadiscus cordeiroi", "Catadiscus inopinatus", "Catadiscus marinholutzi", "Catadiscus mirandai", "Catadiscus propinquus", "Catadiscus pygmaeus", "Catadiscus uruguayensis", "Centrorhynchus giganteus", "Centrorhynchus tumidulus", "cestoda não-identificado", "Choledocystus elegans", "Cosmocerca brasiliense", "Cosmocerca cruzi", "Cosmocerca freitasi", "Cosmocerca ornata", "Cosmocerca parva", "Cosmocerca podicipinus", "Cosmocerca rara", "Cosmocerca travassosi", "Cruzia tentaculata", "Cylindrotaenia americana", "Diplostomum compactum", "Falcaustra mascula", "Glypthelmins linguatula", "Glypthelmins palmipedis", "Glypthelmins parva", "Glypthelmins proximus", "Glypthelmins pseudium", "Glypthelmins repandum", "Glypthelmins simulans", "Glypthelmins vetelliniphilum", "Gorgoderina carioca", "Gorgoderina cedroi", "Gorgoderina cryptorchis", "Gorgoderina diaster", "Gorgoderina festoni", "Gorgoderina pigulevskyi", "Gorgodernina parvicava", "Gyrinicola chabaudi", "Haematoloechus freitasi", "Haematoloechus fuellborni", "Haematoloechus legrandi", "Haematoloechus lutzi", "Haematoloechus ozoroi", "Halipegus dubius", "Infindum infindum", "Maicuru solitarium", "Mesocoelium monas", "Mesocoelium travassosi", "Monostoma sulcatum", "Neohaematoloechus iturbei", "Neohaematoloechus neivai", "Ochoterenella convoluta", "Ochoterenella digiticauda", "Opistogliphe amplicavus", "Oswaldocruzia belenensis", "Oswaldocruzia lopesi", "Oswaldocruzia mazzai", "Oswaldocruzia proencai", "Oswaldocruzia subauricularis", "Oswaldocruzia vaucheri", "Oxyascaris oxyascaris", "Oxysomatium baylisi", "Paraoxyascaris travassosi", "Parapharyngodon alvarengai", "Physaloptera retusa", "Physalopteroides venancioi", "Plagiorchis lenti", "Plagiorchis rangeli", "Polystoma lopezromani", "Pseudoacanthocephalus lutzi", "Pteroxyascaris caudacutus", "Pteroxyascaris similis", "Raillietnema minor", "Raillietnema simples", "Raillietnema spectans", "Rhabdias androgyna", "Rhabdias elegans", "Rhabdias fuelleborni", "Rhabdias hermafrodita", "Rhabdias sphaerocephala", "Rudolphitrema rudolphi", "Schrankiana formosula", "Schrankiana freitasi", "Schrankiana fuscus", "Schrankiana inconspicata", "Schrankiana larvata", "Schrankiana schranki", "Schrankianella brasili", "Schulzia subventricosa", "Schulzia travassosi", "Strongyloides carinii", "Strongyloides pereirai", "Subulascaris falcaustriformis", "Thelandros owaldocruzia", "Travtrema stenocotyle") #fazendo siglas nomes.lin<-rownames(as.data.frame(matriz.final2)) for(i in 1:length(rownames(matriz.final2))) { nomes.lin[i]<-paste(substring(strsplit(nomes.lin[i]," ", fixed="TRUE")[[1]][1], 1,1), substring(strsplit(nomes.lin[i]," ", fixed="TRUE")[[1]][2], 1,2), sep="") } #ai o problema é assim: #no caso do Acr length(nomes.lin[which(nomes.lin==unique(nomes.lin)[2])])>1 #ele ocorre mais de uma vez, então qd isso for verdadeiro #tinha que virar Acr1 e Acr2 #então esses objetos which(nomes.lin==unique(nomes.lin)[2]) #que são iguais tinha que virar o Acr1 e Acr2 Eu tava tentando fazer algo tipo um loop pra testar cada situação unica, igual acima, se ela fosse verdadeira, ai eu tentei fazer um outro loop dentro pra ver qts casos era e ir adicionando 1 2 3 do lado da sigla Mas entrou em loop dentro de loop e me embananei todo :( Alguém poderia dar uma mão? -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 2

Versão do R 2.14.2 mac os travando
by Edson Lira 29 Mar '12

29 Mar '12

Estou rodando um banco do excell xlsx com 1,6 mb. Consigo ler o arquivo, porém, quando uso a rotina names(), o computador está travando. Achei estranho, visto que minha máquina é um mac 2.7 ghz de clok com 4 gb de memória sendo um i7. Alguém tem alguma sugestão? Estou pensando em baixar novamente o R e instalar novamente. Ah! não desisntalei a versão anterior, pode ser isso? []'s. Edson Lira Estatístico Manaus-Amazonas

3 3

Comparar modelos com gnls
by Jonas Brito 29 Mar '12

29 Mar '12

Quero comparar modelos para duas espécies de peixes. Consigo fazer usando nls mas não gnls. Alguém tem alguma idéia de como fazer? dados<-structure(list(Peixes = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("1", "2"), class = "factor"), Tempo = c(0L, 0L, 0L, 30L, 30L, 30L, 60L, 60L, 60L, 90L, 90L, 90L, 240L, 240L, 240L, 270L, 270L, 270L, 0L, 0L, 0L, 30L, 30L, 30L, 60L, 60L, 60L, 90L, 90L, 90L, 240L, 240L, 240L, 270L, 270L, 270L), Peso = c(6.09, 6.09, 6.09, 65.73, 52, 53.13, 102.88, 115.8, 129.31, 161.06, 174.94, 213.63, 233.67, 291.93, 301.33, 258.47, 326.84, 337.04, 8.92, 8.92, 8.92, 78.8, 72.07, 85.8, 167.82, 146.13, 179.31, 255.75, 230.73, 207.06, 320.67, 280.2, 243.93, 376.55, 441.96, 349.84)), .Names = c("Peixes", "Tempo", "Peso"), row.names = c(NA, -36L), class = "data.frame") library(nlme) m1 <-gnls(Peso~a*(1-b*(exp(-c*Tempo)))^3, weights=varPower(), data=dados, start=c(a=350,b=1, c=0.01)) m2<- nls(Peso~a*(1-b*(exp(-c*Tempo)))^3,data=dados, start=c(a=300,b=10, c=0.05)) m3<-gnls(Peso~a[Peixes]*(1-b[Peixes]*(exp(-c[Peixes]*Tempo)))^3,, weights=varPower(), data=dados, start=list(a=c(300,300),b=c(10,10), c=c(0.05,0.05))) m4<-nls(Peso~a[Peixes]*(1-b[Peixes]*(exp(-c[Peixes]*Tempo)))^3,data=dados, start=list(a=c(300,300),b=c(10,10), c=c(0.05,0.05))) anova(m1,m3) anova(m2,m4)

3 2

ajuda contar subperiodos no R
by Lucas Fernandes de Souza 29 Mar '12

29 Mar '12

Boa tarde. Tenho a seguinte planilha, onde as colunas são anos. Cada linha, represente um determinado dia do ciclo da cultura do milho. Os valores "0" representam que nesse dia, não houve chuva. Os valores "1" representam os dias que choveram. Pois bem, preciso contar os períodos de chuva e seca em cada ano. Porém, não é o perído total, e sim os subperíodos. Exemplo: No ano de 1979, tenho 1 dia de chuva (1° subperiodo), depois tenho 6 dias de seca (2° subperiodo), depois 1 dia de chuva (3 subperiodo), depois 2 dias de seca (4 subperiodo) e assim vai..... Existe um jeito de fazer isso no R automaticamente? Tenho vários anos para fazer isso..se for ficar contando, vou demorar muito... Desde já agradeço a ajuda. obrigado Exemplo: 1979 1980 1 0 0 0 0 1 0 1 0 1 0 0 0 0 1 1 0 0 0 0 1 1 0 1 1 1 0 1 1 0 1 0 1 0 1 0 1 1 0 1 0 1 1 1 1 0 1 0 1 0 1 0 1 0 -- Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

4 5