Maio 2011 - R-br - DINF UFPR

Usando AIC para seleção de modelo com o comando drop1(). Alguém poderia dar uma olhadinha :)
by Augusto Ribas 18 Mai '11

18 Mai '11

Ola colegas de lista. Eu estava lendo sobre seleção de modelos e fiquei na duvida se estou indo pelo caminho certo no R. Outra coisa, eu não achei exatamente nos livros que li o que eu deveria reportar em um artigo, ao longo destes passos (se estes estiverem fazendo sentido). Eu tentei fazer um exemplo com uma com uma uma variavel resposta e duas variaveis explicativas. Usei o dataset do anscombe que vem no R, usando o primeiro exemplo dele mais uma variavel ao acaso para exemplificar. Aqui vai um codigo de como estava pensando: #Abrindo os dados anscombe x2<-runif(11) #fazendo um plot para olhar panel.hist <- function(x, ...) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(usr[1:2], 0, 1.5) ) h <- hist(x, plot = FALSE) breaks <- h$breaks; nB <- length(breaks) y <- h$counts; y <- y/max(y) rect(breaks[-nB], 0, breaks[-1], y, col="cyan", ...) } panel.cor <- function(x, y, digits=2, prefix="", cex.cor, ...) { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1)) r <- abs(cor(x, y)) txt <- format(c(r, 0.123456789), digits=digits)[1] txt <- paste(prefix, txt, sep="") if(missing(cex.cor)) cex.cor <- 0.8/strwidth(txt) text(0.5, 0.5, txt, cex = cex.cor * r+1) } panel.lm = function(x, y, ...) { tmp <- lm(y ~ x, na.action = na.omit) abline(tmp) points(x, y, ...)} #legal, y1 é relacionado com x1 mas nao com x2 olhando esse grafico pairs(cbind(y1=anscombe$y1,x1=anscombe$x1,x2),diag.panel=panel.hist, lower.panel=panel.cor,upper.panel=panel.lm) #agora fazendo um modelo modelo1<-lm(anscombe$y1~anscombe$x1+x2) modelo1 #vendo premissas par(mfrow=c(2,2)) plot(modelo1) #vendo o quanto os parametros ajudam drop1(modelo1) #testando quem ajuda significativamente drop1(modelo1,test="F") #novo modelo que eu vou ficar ja que vi no drop1 que posso tirar o x2 modelo2<-lm(anscombe$y1~anscombe$x1) modelo2 #vendo premissas pro novo modelo par(mfrow=c(2,2)) plot(modelo2) #o que reportar num artigo de todo o processo? #todo esse processo pode ser inverso tb, começando de um modelo simples e indo pra um mais complexo add1(modelo2, ~anscombe$x1+x2,test="F") É isso. Agradeço a paciencia de te -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 1

coordenadas valor máximo num grafico xy
by jose coelho 18 Mai '11

18 Mai '11

Vivam , ainda sou aprendiz no R... mas gostaria de saber como calcular as coordenada do valor máximo de y em função de x para colocar num gráfico xy. Existe alguma função ? max(y) eu sei qual valor máximo do y mas também quero saber qual o valor de x no qual o y atinge o seu máximo. Muito obrigado jmmcoelho

3 3

OFF TOPIC - Backup do histórico do firefox e decodificação da data de acesso
by Walmes Zeviani 17 Mai '11

17 Mai '11

Saudações, Peço desculpas pelo tema fora de tópico. Tenho o backup do histórico do navegador, mozilla firefox, usando o complemento FEBE ( https://addons.mozilla.org/pt-BR/firefox/addon/febe/) O arquivo com o backup do histórico de visitação é um arquivo de extensão *.json. Os registros de páginas visitadas apresentam essa estrutura (essa é apenas uma linha das 15000 que tenho) {"aURI":"https://addons.mozilla.org/","aTitle":"addons.mozilla.org ","aLastVisited":1305653472464928} Imagino que o último termo, o numérico com 16 dígitos seja uma forma de armazenar a data de acesso dessa página. Alguém saberia como decodificar esse número em data, no caso a data é 2011-05-17 14:31. Grato. 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

2 2

eixo log contour plot
by eder＠leg.ufpr.br 17 Mai '11

17 Mai '11

Pessoal, to me batendo aqui, Como pode ser feito para deixar apenas o eixo y em escala logaritima (base 10), e que os valores das variáveis continuem aparecendo no grafico. ## Exemplo x <- seq(pi/4, 5 * pi, length.out = 100) y <- seq(pi/4, 5 * pi, length.out = 100) r <- as.vector(sqrt(outer(x^2, y^2, "+"))) grid <- expand.grid(x=x, y=y) grid$z <- cos(r^2) * exp(-r/(pi^3)) contourplot(z~x*y, grid,col.regions=colorRampPalette(c("green","yellow","red")),region =TRUE) att

3 3

Dúvida básica
by Rodrigo Pereira 17 Mai '11

17 Mai '11

Caros colegas, Estou com um problema básico, mas que não estou conseguindo resolver (falta de prática com o R...) Estou com o seguinte banco de dados: resposta fator sujeito 1 0.00 R 1 2 2.96 R 1 3 10.28 R 1 4 12.42 R 1 5 15.60 R 1 6 14.23 R 1 7 10.55 R 1 8 8.00 R 1 9 6.12 R 1 10 3.24 R 1 11 2.82 R 1 12 1.18 R 1 13 0.42 R 1 14 0.00 R 1 15 0.00 R 1 16 0.00 T 2 17 4.01 T 2 18 19.82 T 2 19 17.83 T 2 20 12.74 T 2 21 7.37 T 2 22 7.73 T 2 23 4.70 T 2 24 4.11 T 2 25 3.11 T 2 26 3.19 T 2 27 2.25 T 2 28 0.44 T 2 29 0.16 T 2 30 0.00 T 2 Estou tentando utilizar a função unstack() para transformar este banco de dados longo em largo da seguinte forma: R 0 T 0 2,96 4,01 10,28 19,82 12,42 17,83 15,6 12,74 14,23 7,37 10,55 7,73 8 4,7 6,12 4,11 3,24 3,11 2,82 3,19 1,18 2,25 0,42 0,44 0 0,16 0 0 Alguém sabe como consigo fazer essa mudança com ou sem a função unstack()? Obrigado e abraço! Rodrigo Pereira (31)9326 3139

5 16

?Erros Correlacionados?
by Marcelo Cardoso mello 17 Mai '11

17 Mai '11

Estou com dificuldades em entender a idéia de erros correlacionados. Erros correlacionados com o que? Correlação entre quais variáveis? Entre os erros “puros” e as medias dos tratamentos? Entre os erros “puros” e os valores preditos dos tratamentos? Ou seriam os erros padronizados? Alguém poderia me mostrar um exemplo onde isto ocorre? A idéia é de correlação linear ou de qualquer tipo de dependência, ou seja, uma parabola me dá correlação linear nula mas ocorre relação e existe dependência (não linear). Abaixo segue um exemplo. Onde posso observer a dependência? t<-as.factor(c(rep(1,3),rep(2,3),rep(3,3), rep(4,3))) b<-as.factor(c(seq(1:3), seq(1:3), seq(1:3), seq(1:3))) r<-c(1,2,3,10,18,20,5,25,26,6,30,32) m<-aov(r~t+b) anova(m) e<-resid(m) p<-predict(m) a<-tapply(r,t, mean) b<-tapply(r,t, sd) d<-ave(r, t, FUN = mean) c<-ave(r, t, FUN = sd) ep<-e/(anova(m)$ "Mean Sq"[2]^0.5) ep1<-e/c plot(e~p, xlab= "Valores Preditos", ylab= "Resíduos") cor(e,p) plot(ep~p, xlab= "Valores Preditos", ylab= "Resíduos Padronizados") cor(ep,p) plot(ep1~p, xlab= "Valores Preditos", ylab= "Resíduos Padronizados") cor(ep1,p) plot(c~d, xlab= "Médias", ylab= "Desvios Padrões") cor(c,d) plot(t, e, xlab= "Tratamentos", ylab= "Resíduos") cor(a,b) cor(d,c)

2 1

Capturar dados via scanner
by Carlos Andrade 17 Mai '11

17 Mai '11

Prezados amigos da lista, Estamos realizando uma enquete com os alunos de Economia sobre a qualidade do ensino. Estamos testando um questionário. Observei que a impressora HP C4480 tudo-em-um, que tem um scanner, quando se troca o cartucho pede para imprimir uma página, com manchas quadradas e listras, e que a coloque no vidro do scanner e de lá a impressora lê e alinha os cartuchos. Então pesquisei na internet como programar o scanner para capturar dados de questionários. Assim pouparíamos o tempo e o aborrecimento da digitação. Nada encontrei de útil. Alguém da lista tem alguma dica ou solução para adaptarmos? De forma, também, a serem lidos os dados para o R. Agradeço qualquer ajuda. -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

7 8

ajuste curva potencial
by Rodrigo Plei 16 Mai '11

16 Mai '11

Prezados, Como ajustar uma curva potencial (power curve) em uma dispersão x, y, mostrando a equação e o R2 no gráfico? Desde já, grato. 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

2 1

Dados de contagem com análise de experimentos
by Edson Lira 16 Mai '11

16 Mai '11

Bom dia a todos. Estou tentando analisar o banco de dados anexo com a rotina abaixo: jat1<-read.csv2("...\\jg_09.csv") names(jat1) stopifnot(complete.cases(jat1$ngalhos) !=is.na(jat1$ngalhos))#eliminar os NA's dos dados jat<-jat1[complete.cases(jat1$ngalhos),] Gostaria de verificar o efeito da aração e tplan na contagem de galhos(ngalhos) de uma árvore. Tentei fazer a transformação para "normalidade" a transormação não foi possível. Tentei fazer um teste para verificar se os dados tem distribuição de "poisson" com o pacote vcd, não sei se fiz o correto. Alguém poderia dar algumas sugestões? []'s Edson Lira Estatístico Manaus-Amazonas

3 2

Colocar pontos no mapa
by Victor Eduardo 16 Mai '11

16 Mai '11

Alguém aqui já colocou pontos em um mapa através da latitude e longitude? Estou tentando fazer e preciso até então utilizar alguns pacotes, como por exemplo, o pacote rgdal, sendo que o mesmo vive dando erro ao instalar, pois não consigo usá-lo. Será que existe uma outra forma para fazer isso? Abraços, Victor Eduardo

4 6