Abril 2012 - R-br - DINF UFPR

Python para R
by Daniel Marcelino 04 Jun '12

04 Jun '12

Caros, alguem da lista conhece programação python? Estou há dias tentando desenvolver uma rotina em R para explorar dados eleitorais, mas não tenho conseguido ir muito longe. Eu tenho um modelo do que mais ou menos eu preciso, mas a rotina está em python. Alguem encara me ajudar nisso? Daniel ----------------------------------------------------- My contacts: Phone 1+514 343 6111 Ext 3799 Skype: d.marcelino 3200 Jean Brillant, Office C5075 Montreal, QC; H3T 1N8 Canada -----------------------------------------------------

3 2

windows.options() - Linux
by Luís Gustavo 07 Mai '12

07 Mai '12

Boa tarde pessoal, Estou com a seguinte dúvida: qual o comando que é similar ao windows.options() no Linux, pois este não existe para o Linux. O objetivo final é criar uma função que me gera a figura da legenda, porém esta função tem rodar tanto no Windows quanto no Linux. Pois gostaria que a figura da legenda se ajuste para os mais variados números de faixas. Abaixo segue um exemplo. (Acredito que função não a mais eficiente, qualquer sugestão é bem-vinda) Obrigado. brks <- c(0,5,10,15,20) r <- 3 cols <- brewer.pal(4,"Blues") gera_figura_legenda=function(brks, cols,r) gera_figura_legenda=function(brks, cols,r) { brks = round(brks,r) dest.fig = tempfile(pattern="legenda",fileext=".png",tmpdir=getwd()) windows.options(width=2.5, height=1.7, pointsize=12, reset = TRUE) CairoPNG(filename = dest.fig, width = 290, height = 290, bg="#FFFFFF00", pointsize = 10) par(bg="transparent", mar=c(0, 0, 0, 0)) brks = round(brks,r) plot(1,1,col="white",axes=F,col.axis="white",xlab="",ylab="") aux1 = 2:(length(brks)-1) aux2 = expand.grid(aux1,aux1) d = aux2[,1]-aux2[,2] aux3 = aux2[d==-1,] legend( "center", title = "Legenda", c(paste("<",brks[2]), paste(brks[aux3[,1]],brks[aux3[,2]],sep="-"), paste(">",brks[(length(brks)-1)])), fill=c(cols), cex=1.3, bg="#FFFFFF") dev.off() return (dest.fig) } -- Luís Gustavo Silva e Silva

5 5

Ajuda com Correlação de Matrizes com tamanhos diferentes
by Débora Martins 02 Mai '12

02 Mai '12

Olá pessoal, Estou analisando alguns dados para um artigo, entretanto não consigo rodar a correlação do jeito que desejo. Meu problema é um pouco parecido com o descrito nesse link http://stats.stackexchange.com/questions/24980/correlation-between-matrices…, pois eu desejo obter um único valor que corresponda a correlação, mas como eu divido os valores do arquivo em algumas dimensões que me interesso, não consigo obter esse valor único. Eu li o arquivo,separo cada coluna na dimensão que desejo analisar e tento fazer a correlação entre as dimensões mas só obtenho a correlação dos itens. teste<-read.table("1_analise.txt", header=T) a<-teste[1] b<-teste[2] c<-teste[3] d<-teste[4] e<-teste[5] estrutura<-cbind(a,b,c,d,e) a<-teste[6] b<-teste[7] c<-teste[8] d<-teste[9] e<-teste[10] f<-teste[11] g<-teste[12] cultura<-cbind(a,b,c,d,e,f,g) correlacao1<-cor(estrutura,cultura) Com a última linha do código, gostaria de obter algo como a tabela 3 deste link http://www.tlainc.com/articl292.htm. -- *Débora de Vasconcelos Martins** *Graduanda em Ciência da Computação Universidade Estadual do Ceará Bolsista do Programa de Educação Tutorial - PET http://petcomp.uece.br

2 3

Duvida Regressão
by Sérgio Henrique almeida da silva ju 01 Mai '12

01 Mai '12

Olá Estou trabalhando com a taxa de uma doença e gostaria de rodar um modelo de regressão. Meus dados não apresentaram distribuição normal e por teste de aderência deu uma distribuição lognomal -> fitdistr(tx_esqui, "lognormal") meanlog sdlog 6.25725753 1.90754231 (0.11630491) (0.08223999) ->ks.test(tx_esqui, "plnorm", sdlog=1.90754231, meanlog=6.25725753,alternative=c("two.sided")) One-sample Kolmogorov-Smirnov test data: tx_esqui D = 0.0701, p-value = 0.1419 alternative hypothesis: two-sided Tem como eu rodar um modelo de regressão com essa variável? Como fazê-lo? Qual a melhor forma de trabalhar com essa variável? Abraços Segue a minha variável dependente: structure(list(tx_esqui = c(3L, 3L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 9L, 10L, 10L, 19L, 22L, 22L, 23L, 23L, 25L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L, 34L, 36L, 39L, 41L, 42L, 45L, 46L, 49L, 49L, 56L, 57L, 58L, 58L, 59L, 61L, 62L, 62L, 68L, 70L, 74L, 76L, 80L, 82L, 82L, 83L, 83L, 86L, 101L, 102L, 107L, 107L, 112L, 113L, 118L, 120L, 121L, 123L, 126L, 135L, 138L, 140L, 153L, 156L, 158L, 179L, 191L, 200L, 201L, 204L, 207L, 209L, 212L, 215L, 223L, 227L, 228L, 233L, 247L, 248L, 249L, 255L, 257L, 258L, 258L, 259L, 264L, 274L, 277L, 279L, 287L, 297L, 309L, 310L, 313L, 314L, 326L, 328L, 329L, 331L, 332L, 332L, 340L, 354L, 358L, 365L, 372L, 379L, 383L, 395L, 396L, 418L, 423L, 447L, 449L, 451L, 453L, 493L, 532L, 535L, 553L, 554L, 565L, 573L, 581L, 584L, 599L, 603L, 604L, 611L, 672L, 703L, 708L, 727L, 756L, 799L, 809L, 812L, 815L, 816L, 862L, 864L, 866L, 932L, 946L, 951L, 964L, 969L, 983L, 985L, 989L, 1032L, 1102L, 1102L, 1112L, 1112L, 1199L, 1218L, 1245L, 1252L, 1283L, 1302L, 1327L, 1327L, 1328L, 1340L, 1367L, 1382L, 1413L, 1441L, 1468L, 1487L, 1493L, 1512L, 1523L, 1567L, 1576L, 1585L, 1598L, 1620L, 1623L, 1634L, 1664L, 1704L, 1813L, 1816L, 1824L, 1842L, 1892L, 1939L, 1977L, 2010L, 2035L, 2060L, 2100L, 2152L, 2184L, 2210L, 2248L, 2282L, 2288L, 2358L, 2462L, 2473L, 2527L, 2541L, 2555L, 2576L, 2603L, 2705L, 2711L, 2895L, 2916L, 2958L, 2963L, 3088L, 3110L, 3304L, 3427L, 3489L, 3535L, 3629L, 3638L, 3682L, 3704L, 3796L, 3932L, 4160L, 4293L, 4294L, 4366L, 4516L, 4715L, 4746L, 5152L, 5263L, 5307L, 5823L, 6025L, 6306L, 6313L, 6330L, 6649L, 7114L, 7158L, 7340L, 7488L, 7665L, 8055L, 8394L, 8773L, 8850L, 8918L, 9273L, 9356L, 10110L, 10281L, 10310L, 10519L, 11111L, 13525L, 14465L, 14955L, 16629L)), .Names = "tx_esqui", class = "data.frame", row.names = c(NA, -269L)) -- 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) 94429486/78101651 id: 123*20942

7 16

Contagem das linhas de uma matriz
by Diego Leal 30 Abr '12

30 Abr '12

Boa tarde pessoal, a questão é a seguinte: queria saber se existe algum comando equivalente ao comando table() para tabelar os valores das linhas de uma matriz. Por exemplo: matriz=matrix(rep(1:20),10,10,byrow=TRUE) funcaoDesconhecida(mat) # essa funcao acima deveria me retornar falando que a linha (1,2,...,10) aparece 5 vezes e a linha (11,12,13,...,20) aparece 5 vezes. alguma ideia? Att, Diego Leal Togni

2 3

Tutorial do WinBUGS no R
by Gilbert Queiroz 30 Abr '12

30 Abr '12

Alguém tem, ou sabe onde tem, tutorial de uso do WinBUGS no R? Abs. Gilbert.

2 1

Pacote GPCM
by Felipe Buchbinder 30 Abr '12

30 Abr '12

Prezados, Alguém sabe onde eu posso encontrar o pacote gpcm para instalá-lo? Eu li sobre ele em um artigo, mas não estou conseguindo encontrá-lo nem no site do CRAN nem através do comando "Install packages" do R. Desde já, obrigado, Felipe

2 4

Ajuda em função para criação do Filtro De Kalman
by Victor Eduardo 30 Abr '12

30 Abr '12

Galera, estou criando uma função para realizar o Filtro de Kalman sendo que em certo momento, preciso transpor uma matriz, só que essa matriz está dentro de um array, logo não consigo realizar esta operação. Está dentro de um for, vou colar aqui a parte que está dando erro: ar_F[,,t]<-ar_Z[,,t]%*%ar_P[,,t]%*%t(ar_Z[,,t])+ar_H[,,t] Acredito que o erro esteja justamente em t(ar_Z[,,t]), pois é um array, e não uma matriz. Como poderia fazer isso? Abraços

2 3

sobre distribuções tempered stable
by Marcelo Zeuli 30 Abr '12

30 Abr '12

Bom dia, Gostaria de saber se a pacote do R que implemente distribuições do tipo tempered stable. Grato, Marcelo Zeuli

2 1

Perda de Parcela.
by Carlos Pereira 30 Abr '12

30 Abr '12

Boa tarde queria saber se poderiam me ajudar, estou com uma duvida de como fazer ANOVA, quando existe a perda de parcela, como posso fazer para que o R estime essa parcela perdida para realizar ANOVA. Em um delineamento em Blocos ao Acaso. Desde já obrigado pela atenção! Carlos Pereira da Silva Bacharel em Estatística - UEPB Campina Grande Cel: 083 - 9672- 1966 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

4 3