- R-br - DINF UFPR

Gráficos em um única janela com um único titulo de x
by Edson Lira 13 Jun '49

13 Jun '49

Caros amigos, a rotina abaixo posiciona 5 gráficos em um janela. como faço para mudar o nome do eixo x de x para "TESTE", por exemplo. No help "curve" mostra como em cada gráfico mudar o nome do eixo x, mas não consegui implementar. Alguém tem alguma alternativa. require(graphics) curve(dnorm(x,mean=0.338,sd=0.163), lty=1,lwd=2,from=0,to=2.5) curve(dnorm(x,mean=0.631,sd=0.348),lty=2,col=2,lwd=2,add=T) curve(dnorm(x,mean=1.2,sd=0.292),lty=3,col=3,lwd=2,add=T) curve(dnorm(x,mean=1.5,sd=0.364),lty=4,col=4,lwd=2,add=T) curve(dnorm(x,mean=1.370,sd=0.224),lty=5,col=5,lwd=2,add=T) legend('topright',legend=c(expression(mu==0.338),expression(mu==0.631),expression(mu==1.2),expression(mu==1.5) ),text.col=c(1,2),cex=1) title("Canetas Bic") [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 2

dúvida função offset para testar o coeficiente angular (regressão linear)
by Rodrigo Plei 13 Jun '49

13 Jun '49

Prezados, Estou utilizando a função *offset *para tesar se o coeficiente angular (b) de uma regressão linear simples: mod1 <- lm(y ~ x+offset(x)) Ocorre que fiz esse mesmo teste no Excel utilizando a fórmula disponível em Zar (5th Edition, 2010, pág. 342), e estou obtendo um valor diferente de p. Estou confuso em qual resultado confiar. Imagino que a função *offset *faça o mesmo tipo de teste t descrito em Zar (2010), pois não? Agradeço desde logo pelos eventuais esclarecimentos! Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777 Google Scholar: http://scholar.google.com.br/citations?user=zX_EezEAAAAJ&hl=pt-BR

2 2

readShapePoly
by Hélio Gallo Rocha 13 Jun '49

13 Jun '49

Prezados da lista Com o intuito de selecionar parte do atributo de um shape executei o exemplo do readShapePoly: library(maptools) xx <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1], IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66")) xxx <- xx[xx$SID74 < 2,] plot(xxx, border="red") #end xxx <- xx[xx$SID74 < 2,] irá selecionar atributos menores que 2, ele é numeric. Fiz um teste com um shape do IBGE, não dava certo, verifiquei que a classe dos atributos por mim escolhidos eram do tipo factor. Duas dúvidas: 1. Existem atributos numéricos no shape, mas teria como selecionar um atributo de classe factor ou character? 2. Como selecionar três atributos em intervalos diferentes? Agradeço desde já. -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 0

Dúvida KS regressão logistica e iteração de variáveis em glm
by Vinicius Brito Rocha 13 Jun '49

13 Jun '49

Pessoal, estou com duas dúvidas diferentes, mas vou tentar aproveitar o mesmo email para saná-las. 1º alguém conhece algum pacote com teste ks para regressão logistica? 2º ajustei um modelo com iterações entre SEXO e UF e outras variáveis. Minha dúvida está no seguinte para o meu problema em específico a interação RJ e M não apresentou um p-valor significativo. É possível controlar em específico essa interação meu modelo glm(cbind(MARCADO,I(1-MARCADO))~UF+SEXO+UF:SEXO-1,data=base,family=binomial(link="logit")) Abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário ** M.Sc. Engenharia de Produção/PO*

2 1

helio rocha wants to share "krigagem" with you
by Helio via Dropbox 13 Jun '49

13 Jun '49

Hi there, Helio (heliogallorocha(a)gmail.com) wants to share some files in a folder called "krigagem" with you via Dropbox. View "krigagem" https://www.dropbox.com/l/2NdOOv8EukCllzjZ9UVRD9/invite Enjoy! - The Dropbox Team

2 1

Re: [R-br] R-br] Tempo de processamento e regressão quantílica
by Roney Fraga Souza 12 Jun '49

12 Jun '49

Alberto Gaidys, Funcionou. Muito obrigado pela ajuda! Roney

1 0

[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

Tempo de processamento e regressão quantílica
by Roney Fraga Souza 12 Jun '49

12 Jun '49

Caros, Estou rodando um modelo de regressão quantílica com o pacote 'quantreg', onde os dados estão salvos como objeto 'ffdf', com 7035146 linhas e aproximadamente 13 colunas. Para realizar um teste quanto a possibilidade computacional dessa regressão fiz o seguinte, rodei uma regressão com apenas 3 variáveis explicativas e todas as linhas, o modelo rodou, mas demorou 15 horas de processamento. Questão: Como posso melhorar esse tempo de processameto? Meu computador é um macbook pro i7, 10 gb ram com ssd. > class(c2000) [1] "ffdf" > str(c2000[,]) 'data.frame': 7035146 obs. of 13 variables: $ V0102: int 12 12 12 12 12 12 12 12 12 12 ... $ V0300: int 1339 1782 1782 2934 5534 5534 11305 11305 19520 19520 ... $ V1001: int 1 1 1 1 1 1 1 1 1 1 ... $ V1004: int NA NA NA NA NA NA NA NA NA NA ... $ V0401: int 1 1 2 1 1 2 1 2 1 2 ... $ V0402: int 1 1 3 1 1 2 1 2 1 3 ... $ V4752: int 25 45 18 23 20 16 43 38 60 20 ... $ V0408: int 4 1 1 4 4 1 4 1 1 4 ... $ V0418: int NA NA NA NA NA NA NA NA NA NA ... $ V4210: int NA NA NA NA NA NA NA NA NA NA ... $ V4300: int 0 4 11 4 5 5 11 15 0 13 ... $ V4513: int 182 182 136 182 182 36 733 865 316 380 ... $ V0453: int 20 48 48 45 49 24 40 25 44 25 ... > system.time( r <- rq(V4513 ~ V0453 + V4752 + V0401, tau=c(0.1,0.5,0.9), data=c2000) ) user system elapsed 56258.453 165.855 56123.961 Atenciosamente Roney

2 2

R-br] Tempo de processamento e regressão quantílica
by Alberto Gaidys 11 Jun '49

11 Jun '49

Roney: Tente incluir method="fn" (opção descrita na "vignette" do quantreg para situações como a sua). Outra alternativa é você extrair as regressões com base em várias amostras e, a partir delas, construir a distribuição dos coeficientes. Espero que ajude. Alberto Gaidys

1 0

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