Novembro 2013 - R-br

[OFF -Tópic] Cotas em concursos para negros e índios!
by Andre Oliveira 12 Jun '49

12 Jun '49

Pessoal, tem um colegar meu que vou fazer um concurso ele é descendente de negro ele é pardo. Ele veio com a pergunta como provar que sou negro, pois, tem vagas reservada para negros e índios? Eu sinceramente não sei e não tenho ideia pois nunca tive esta curiosidade. O que de fato comprova sua raça? Algum documento? O que são usados em concursos para provar? obrigado André Oliveira Souza

5 6

Legenda eixo X
by Luciano F. Sgarbi 10 Jun '49

10 Jun '49

Prezados, Estou tentando centralizar a legenda do eixo X em uma figura com vários gráficos. Estou conseguindo centralizar a legenda da linha inferior como no exemplo, mas as linhas do meio não. par(mfrow=c(2,2)) plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 1",outer=TRUE) ### o problema está aqui, pois a legenda desce para a parte inferior da figura. plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 2",outer=TRUE)## Esse dá certo Alguém podia dar uma luz de como eu posso fazer isso? Abraços -- *Luciano F. Sgarbi* Mestrando em Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás, campus II Goiânia-GO 74001-970 Brazil

2 2

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

Transformação Box Cox
by ana paula coelho madeira 19 Dez '13

19 Dez '13

Boa tarde pessoal, Tenho uma variável resposta constituída de comprimento de raízes, porém tenho muitos valores zero (a planta não desenvolveu a raiz). Tentei fazer uma transformação de box-cox para tentar normalizar os dados, porém quando fiz: boxcox(MVR ~ FTR, data = dados, lambda = seq(-1, 1, 0.1)) Error in boxcox.default(MVR ~ FTR, data = dados, lambda = seq(-1, 1,0.1)) : response variable must be positive Isso ocorre devido ao excesso de zeros? Alguma sugestão de como modelar os dados?Desde já agradeço Ana Paula

6 7

Plotar intervalo de confiança
by Fátima Lima Paula 04 Dez '13

04 Dez '13

Prezados, tenho a taxa de mortalidade por idade e queria plotar o intervalo de confiança. Alguém pode me ajudar? pop=c(37374, 54783, 86049, 126048, 287384) obfaF=c(6,15,27,55,239) txobfaF=(obfaF*365/popfaF)*100 txobfaF idade=c("1","2","3","4","5") plot(idade, txobfaF,type="b",main="Sexo Feminino",xlab="Idade",ylab="Taxa de Mortalidade",lty=2,ylim=c(1,40),xaxt='n') Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

6 14

Intervalos de confiança plot
by Gustavo Faria 01 Dez '13

01 Dez '13

Olá, Fiz uma análise discriminante usando o comando lda no pacote MASS e para construção do gráfico utilizei o comando plot. Gostaria de saber como coloco as elipses dos intervalos de confiança no gráfico. Muito obrigado -- --- Gustavo de Faria Silva Biólogo herpetólogo Crbio 87421/04-P Pós-Graduado em Consultoria e Licenciamento Ambiental Mestrando em Zoologia de Vertebrados Telefone: (31) 9579-8473 (Vivo)

2 1

Ajuda em uma função
by Sérgio Henrique almeida da silva ju 30 Nov '13

30 Nov '13

Amigos Eu tenho um conjunto de dados com a média da idade dos participantes e o desvio padrão da mesma e gostaria de gerar a distribuição dessas idades para que eu possa obter a proporção de uma determinada faixa etária. Eu fiz a função abaixo: prop <- function(x){ for (i in 1:nrow(x)) z[i]=rnorm(x$npart[i], mean = x$idade[i],sd =x$sd[i]) z y = ifelse(z < 45,"menor de 45 anos","maior de 45 anos") tabela <- table(y) tabela <- round(prop.table(tabela)*100,2) resultado <- tabela[2] print(resultado)} Onde npart é o número total de participantes. O que tenho que modificar nessa função para que eu possa obter os resultados? Att -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

3 3

Amostragem hierárquica
by André Lucas de Oliveira Moreira 30 Nov '13

30 Nov '13

Olá pessoal, boa noite. Quero avaliar o efeito da fragmentação sobre a diversidade de aves em duas paisagens (uma com mais habitat outra com mennos e outras coisas que são medidas na paisagem) também vou usar algumas métricas de distúrbio crônico que são medidas dentro dos fragmentos. Então eu pesei em usar uma amostragem hierárquica com paisagem/fragmento/unidades amostrais. Mas a minha dúvida é: como poderia usar este modelo com variáveis contínuas em ambos os eixos? A única variável categórica seria paisagem (uma mais fragmentada e outra menos), as demais seriam contínuas, na verdade algumas com dados de contagem e outras com dados contínuos. Alguma sugestão pessoal? Obrigado, André -- *MSc. André Lucas de O. Moreira* http://lattes.cnpq.br/7258065668864153 <http://www.wikiaves.com.br/perfil_andrelukinhas> http://www.wikiaves.com.br/perfil_andrelukinhas 79 8837-3562 79 9153-5552

2 1

Para submeter um pacote ao CRAN
by Freddy Hernandez Barajas 29 Nov '13

29 Nov '13

Ola Eu fiz um pacote chamado AA que depende do pacote BB, essa dependência foi incluída no campo Depends do arquivo DESCRIPTION. Eu consigo aplicar o "check" e o "build" sem erros, notas e avisos e assim eu obtenho o arquivo AA_1.0.tar. Quando eu tento instalar o pacote AA num outro computador o qual não tem o pacote BB instalado eu obtenho a seguinte mensagem: ERROR: dependency 'BB' is not available for package 'AA' A questão é: posso submeter meu pacote ao CRAN para conseguir o arquivo AA_1.0.zip ? O erro vai desaparecer? Atenciosamente Freddy Hernández Barajas

2 1

Ajustamento de dados à curva exponencial
by Cícero C Nunes [rh] 28 Nov '13

28 Nov '13

Caro Walmes, boa tarde. Li o comentário sobre curva exponencial no e-mail abaixo endereçado ao grupo. Pois bem, estou com um probleminha para resolver, a saber: Tenho estas duas séries: x<-c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23) y<-c(100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950 ,1000,1050,1100,1150,1200) Utilizei um script apresentado no Livro Conhecendo o R de Luiz A Pertenelli e Márcio P de Mello pág. 165 e não consegui resolver no R. Bom! Ajustei a curva y para exponencial Y no EXCEL conforme abaixo: Y<-c(187,206,227,250,275,303,334,368,405,446,491,541,596,657,723,797,877,966 ,1064,1172,1291,1422,1566) Por favor, como ajusto à exponencial a série y no R? Como chegar ao resultado Y no R? #Script utilizado que não funcionou: x<-c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23) y<-c(100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950 ,1000,1050,1100,1150,1200) dados<-data.frame(x,y) plot(dados) funcao<-y~a*exp(b*x) exponencial<-nls(funcao, dados, start=c(a=1,b=1)) exponencial =====Erro no R ================== 10,11,12,13,14,15,16,17,18,19,20,21,22,23) > y<-c(100,150,200,250,300,350,400,450,500,550,600,650,700,750,800,850,900,950 ,1000,1050,1100,1150,1200) > dados<-data.frame(x,y) > plot(dados) > funcao<-y~a*exp(b*x) > exponencial<-nls(funcao, dados, start=c(a=1,b=1)) Erro em nls(funcao, dados, start = c(a = 1, b = 1)) : gradiente singular > exponencial Erro: objeto 'exponencial' não encontrado > Onde estou errando? Abs. _____ <http://www.posicional.com/> assinatura_email_sites Cícero C Nunes Diretor Técnico Operacional Divisão Remuneração <mailto:c2n@posicional.com> c2n(a)posicional.com +55(16)3397.0226 +55(11)3280.0226 +55(11)9.9978.0026 SKYPE: consultormp Rua Padre Duarte, 1360 salas 6/7 Centro Araraquara, SP CEP.: 14.801-310 Antes de imprimir pense em sua responsabilidade com o meio ambiente Quem ama cuida! De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de walmes . Enviada em: sexta-feira, 8 de novembro de 2013 18:33 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Contrastes em modelos de regressão de poisson inflacionados por zeros Você não deve esperar ver uma curva exponencial para a média das contagens em função do tempo e uma sigmoide para a proporção de zeros não Poisson para os zeros se não simulou os dados com efeito de tempo e não declarou efeito de tempo na porção binomial. Para que isso aconteça a simulação tem que estar condizente, tem que ter os efeitos para que as curvas não sejam as "curvas nulas". #------------------------------------------------------------------ # Definições da sessão. rm(list=ls()) require(pscl) require(VGAM) require(multcomp) require(lattice) require(latticeExtra) #------------------------------------------------------------------ # Dados artificiais. da <- expand.grid(tempo=rep(1:10), trat=gl(3,10)) ## simula porção Poisson X <- model.matrix(~trat+tempo, da) colnames(X) beta.pois <- c(-1,0.25,0.5,0.25) eta.pois <- X%*%beta.pois da$y.pois <- rpois(nrow(X), lambda=exp(eta.pois)) xyplot(y.pois~tempo|trat, da) ## simula porção binomial X <- model.matrix(~tempo, da) beta.bern <- c(0,0.3) eta.bern <- X%*%beta.bern da$y.bern <- rbinom(nrow(X), size=1, prob=1/(1+exp(-eta.bern))) xyplot(y.bern~tempo|trat, da) ## monta a contagem com inflação de zeros da$y <- with(da, ifelse(y.bern==0, 0, y.pois)) xyplot(y~tempo|trat, da) #------------------------------------------------------------------ # Modelo completo. compl.mod <- zeroinfl(y~trat+tempo|trat+tempo, data=da) coef(compl.mod) #------------------------------------------------------------------ # Predição do modelo considerando as duas porções. pred <- expand.grid(tempo=rep(1:10), trat=gl(3,1)) X <- model.matrix(~trat+tempo, pred) i <- grep("^count\\_", names(coef(compl.mod))) eta <- X%*%coef(compl.mod)[i] pred$m.pois <- exp(eta) i <- grep("^zero\\_", names(coef(compl.mod))) eta <- X%*%coef(compl.mod)[i] pred$p.zero <- exp(eta)/(1+exp(eta)) xyplot(y~tempo|trat, data=da, jitter.x=TRUE)+ as.layer(xyplot(m.pois~tempo|trat, data=pred, type="l"))+ as.layer(xyplot(p.zero~tempo|trat, data=pred, type="l", lty=2, lwd=2))+ layer(panel.abline(h=1, lty=2)) xyplot(p.zero~tempo|trat, data=pred, type="l", lty=2, lwd=2) #------------------------------------------------------------------ À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218 ==========================================================================

2 2