Dezembro 2014 - R-br

parallel no R
by Fernanda De Bastiani 17 Dez '14

17 Dez '14

Caros, alguém pode me indicar algum material (por ex. tutorial) que tenha exemplos de como fazer paralelização de função no R, além do manual do pacote parallel? Desde já muito obrigada,Fernanda

4 3

Mestrado em Modelagem Matemática e Computacional
by Wagner Tassinari 16 Dez '14

16 Dez '14

Gostaria de anunciar, que o Programa em Modelagem Matemática e Computacional (PPGMMC) está com as inscrições abertas para o Processo Seletivo de ingresso para o ano de 2015 . As inscrições são gratuitas e estão abertas do dia 19 de novembro de 2014 a 13 de janeiro de 2015. Conforme consta no edital, o endereço de inscrição é: Secretaria do PROPPG - UFRRJ – Campus Seropédica, sala 115, Prédio P1 BR 465 Km 7, Seropédica, RJ CEP.: 23897-000 Tel./Fax: (21) 26814707 As inscrições podem ser efetuadas pessoalmente (das 9:00 às 11:30 horas e das 13:30 às 16:00 horas) ou através de postagem. Quaisquer dúvidas ou informações os candidatos podem entrar em contato através dos emails: ppgmmc.ufrrj(a)gmail.com ppgmmc(a)ufrrj.br Site: Home Telefones: (21) 2681-4707 (21) 2681-4600 | | | | | | | | | HomeNotícias do MMC-UFRRJ | | | | Visualizar em r1.ufrrj.br | Visualizado por Yahoo | | | | |

1 0

Cálculo do deviance das variáveis explicativas em modelo de Poisson inflacionado de zeros
by ASANTOS 15 Dez '14

15 Dez '14

Bom dia Pessoal, Estou querendo saber como realizar o cálculo do deviance das variáveis explicativas em um modelo de Poisson inflacionado de zeros com a função zeroinfl() do pacote pscl, pois para um modelo do tipo glm eu faria uma anova(modelo, test="Chi"), mas um objeto criado pela função zeroinfl() não permite anova() e o teste de Vuong serve apenas para comparação entre modelos, existe alguma opção? Segue CRM: #------------------------------------------------------------------ # Definições da sessão. rm(list=ls()) require(pscl) #------------------------------------------------------------------ # Dados artificiais. da <- expand.grid(trat=gl(4,1), tempo=1:10) X <- model.matrix(~trat+tempo, da); ncol(X) betas <- c(0.1,0.9,0.6,0.3,0.7) eta <- X%*%betas y1 <- rpois(da$trat, lambda=exp(eta)) y2 <- rbinom(y1, size=1, prob=0.7) da$y <- y1*y2 str(da) #------------------------------------------------------------------ # Ajuste do modelo. m0 <- zeroinfl(y~trat+tempo|trat, data=da) summary(m0) #------------------------------------------------------------------ Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

2 2

Numeração de figuras Sweave na classe artigo
by Andre Oliveira 14 Dez '14

14 Dez '14

Pessoal,bom dia. Minhas figuras estão saindo assim. Figura 1Figura 2 Gostaria que fosse impressa assim. Figura 1.1Figura 1.2 Estou trabalhando assim. \SweaveOpts{width=15,height=12}\setkeys{Gin}{width=0.8\textwidth}\begin{figure}[H]\centering<<fig=T,keep.source=TRUE>>= par(mfrow=c(1,2))curve(dnorm(x,mean=1,sd=sqrt(2)),lwd=2,from=-6,to=17,ylab="Densidade")curve(dnorm(x,mean=10,sd=sqrt(2)),col=2,lwd=2,add=T)legend('topright',legend=c(expression(mu==1),expression(mu==10)),text.col=c(1,2),cex=1.5,bty="n")title("Comparando a distribuição normal com médias diferentes")curve(dnorm(x,mean=4,sd=sqrt(3)),lwd=2,from=-15,to=25,ylab="Densidade")curve(dnorm(x,mean=4,sd=sqrt(20)),col=2,lwd=2,add=T)legend('topright',legend=c(expression(sigma^2==3),expression(sigma^2==20)),text.col=c(1,2),cex=1.5,bty="n")title("Comparando a distribuição normal com variâncias diferentes")@\caption{Curva normal}\label{curvas}\end{figure} Tem como acertar isto? obrigado André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

1 0

Contastes não lineares
by Leonardo Ferreira Fontenelle 13 Dez '14

13 Dez '14

Boa tarde a todos! Procurei na web, mas ainda não achei a resposta para esta dúvida. Tenho um objeto *d* criado pela função *svydesign()*, do pacote *survey*. Duas variáveis desse objeto são *x* (uma variável contínua de 0 a 1) e *y* (uma variável dicotômica, 0 ou 1). Depois de ajustar o seguinte modelo: m <- svyglm(formula = y ~ x, design = subset(d, complete.cases(x, y)), family = quasibinomial) Eu gostaria de estimar a proporção (não as chances) de *y* quando *x* é igual a 0; o mesmo para quando *x* é igual a 1; e então calcular a diferença ou a razão entre essas proporções. Pelo que li, *svycontrast()* é a função mais flexível para estimar contrastes, inclusive não lineares. No entanto, ainda não consegui descobrir qual é a sintaxe para eu convencer o *svycontrast()* a estimar o que estou procurando. Alguém poderia me ajudar com isso? Obrigado! Leonardo Ferreira Fontenelle[1] Links: 1. http://lattes.cnpq.br/9234772336296638

1 1

Re: [R-br] Barplot múltiplo com dados de proporção de variáveis categóricas
by Ivan Bezerra Allaman 11 Dez '14

11 Dez '14

Há relação entre as variáveis? Se houver você pode fazer: tab <- ftable(d) mosaicplot(tab) O resto da confecção você tem capacidade de fazer! (s,f,p) Allaman

1 0

Barplot múltiplo com dados de proporção de variáveis categóricas
by ASANTOS 10 Dez '14

10 Dez '14

Bom dia Pessoal, Tenho 3 variáveis categóricas a, b e c, sendo a com três fatores, b com dois e c com quatro. Eu consigo representar a proporção delas em um barplot individualmente, porém gostaria de colocar as três variáveis em um mesmo barplot, isso é possível? Segue CRM: d <- data.frame(a=factor(sample(c('doente','sadio','indeterminado'), 100, rep=TRUE)), b=factor(sample(c('sim','nao'), 100, rep=TRUE)), c=factor(sample(c('ruim','regular','bom','otimo'), 100, rep=TRUE))) a.tb <- prop.table(table(d$a)) barplot(a.tb) b.tb <- prop.table(table(d$b)) barplot(b.tb) c.tb <- prop.table(table(d$c)) barplot(c.tb) Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

1 0

Teste
by Fabiana Goncalves Barbosa 07 Dez '14

07 Dez '14

Tenho os seguintes dados indep inter 288 1887 187 2195 377 111 191 154 165 379 350 156 189 153 100 250 173 186 176 122 431 108 249 1841 213 426 314 252 122 184 635 138 168 169 151 166 154 117 109 141 100 313 1030 125 350 116 329 726 253 420 218 245 123 111 108 135 313 124 275 118 192 100 197 246 222 247 103 567 557 558 299 211 102 184 116 180 119 609 259 155 151 634 343 472 299 454 242 381 167 249 162 248 145 124 126 169 125 151 157 156 208 122 161 106 139 104 101 434 NA 190 NA 171 NA 122 NA 117 NA 2016 NA 349 NA 309 NA 271 NA 255 NA 259 NA 239 NA 158 NA 199 NA 200 NA 168 NA 1368 NA 622 NA 530 NA 373 NA 271 NA 121 NA 188 NA 123 NA 194 NA 125 NA 1811 NA 403 NA 322 NA 229 NA 260 NA 215 NA 151 NA 129 NA 135 NA 104 NA 543 NA 460 NA 286 NA 265 NA 248 NA 164 NA 130 NA 125 NA 110 NA 192 NA 122 NA 103 NA 361 NA 314 NA 262 NA 193 NA 155 NA 513 NA 301 NA 243 NA 177 NA 155 NA 188 NA 193 NA 122 NA 108 NA 109 NA 191 NA 324 NA 163 Que teste uso para verificar se existe diferença entre as duas amostras?

3 4

extracao da Longitude correta para dados de SST
by Humberto 07 Dez '14

07 Dez '14

Boa noite pessoal, Estou extraindo de um arquivo ncdf os valores de latitude (lat) longitude( lon) e temperatura da superfície da agua do mar( SST), entretanto a longitude (lon) ela está na projeção 0 a 360 graus. Minha pergunta é como faço para converter ela para -180 a 180 graus? Segue o script que estou utilizando juntamente com o link dos dados. Desde já agradeço a ajuda Humberto lon lat data sst Lon como deveria ser 299.5 10.5 11/30/1981 27.85 -60.000 300.5 10.5 11/30/1981 27.96 -59.500 301.5 10.5 11/30/1981 28.02 -58.500 302.5 10.5 11/30/1981 28.02 -57.500 303.5 10.5 11/30/1981 27.97 -56.500 304.5 10.5 11/30/1981 27.84 -55.500 305.5 10.5 11/30/1981 27.81 -54.500 306.5 10.5 11/30/1981 27.78 -53.500 307.5 10.5 11/30/1981 27.75 -52.500 Link para acesso dos dados https://www.dropbox.com/s/xnab2g365mduykx/SST_81_14.nc?dl=0 Script: require(ncdf) #'ncdf' versao 1.6.6 require(reshape) require(chron) ###sst nc <- open.ncdf("SST_81_14.nc") nc lat <- get.var.ncdf(nc, "lat") lon <- get.var.ncdf(nc, "lon") tunits <- att.get.ncdf(nc,"time","units") # print tunits tunits$value # criar tempo "real" Times time<- get.var.ncdf(nc,"time") Time<-chron(time,origin=c(month=1, day=1, year=1900)) ## array da variavel sst wm <- get.var.ncdf(nc, "sst") str(wm) ## atribuindo dimnames conforme lon, lat e tempo dimnames(wm)[[1]] <- as.character(lon) dimnames(wm)[[2]] <- as.character(lat) dimnames(wm)[[3]] <- as.character(Time) str(wm) ## melt para reestruturação dos dados mwm <- melt(wm) head(mwm)

6 6

RCurl: passando parâmetros e executando função javascript
by Éder Comunello 04 Dez '14

04 Dez '14

Colegas, boa tarde! Estou precisando de ajuda com relação ao pacote *RCurl*. Tomando por exemplo o endereço < http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001>, ao acessá-lo via *browser *pode-se notar dois botões que executam trechos de código em *javascript *e geram links para o *download *de arquivos de dados: * Botão [Arquivo Access] remete ao comando "javascript:criarArq(2254001,2)" * Botão [Arquivo Texto] remete ao comando "javascript:criarArq(2254001,1)" Estou tentando usar o *RCurl *para passar diretamente os parâmetros para criação do *link *para o "Arquivo Texto". Isto feito, posteriormente pretendo usar o *link *gerado para *download* dos dados. Contudo, não estou tendo sucesso e ficarei muito grato por qualquer ajuda. Abaixo segue o código que tentei aplicar: ### <BEGIN> ### library(RCurl) myURL <- 'http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001' ### URL resultante da operação via browser: ### myURL2 <- ' http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001&CriaArq=true&TipoArq=1 ' myCurl <- getCurlHandle(cookiefile = "myCookies.txt", useragent="Mozilla/5.0", verbose = T, autoreferer = T, followlocation = T) myParams <- list(Codigo=2254001, CriaArq='true', TipoArq=1) ### Codigo: ID da estação ### CriaArq: true (executa função?), ### TipoArq: 1 (tipo texto) e 2 (Access) postForm(myURL, .params=myParams, curl=myCurl) getURL(myURL, curl=myCurl) ### Operação falha, retornando página original! ### Em caso de sucesso embutiria no código um link no formato: ### <a href="ARQ/A20131008-164720-14/CHUVAS.ZIP"> ### <END> ### -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

4 3