- R-br - DINF UFPR

Re: [R-br] Alteração em função de pacote
by Junior Beleti 31 Mai '12

31 Mai '12

Olá Leonard. Estou criando uma função no pacote base, função esta semelhante à função eigen.f já existente. A nova função teria o nome eigem.f Entendi Benilton. Vou tentar colocar um código reproduzivel aqui. Obrigado, Carlos.

1 0

Criar um loop
by Humberto Hazin 31 Mai '12

31 Mai '12

Caros, Venho mais uma vez solicitar ajuda a vocês tenho um banco de dados com cerca de 80 mil linhas conforme o exemplo abaixo. gostaria de atribuir um id para cada barco, exemplo toda a coluna que tivesse o Barco Vieirasa a coluna ID = 1 e assim por diante. Alguem pode me ajudar? Humberto boat Id VIEIRASA COSTA.DE.SAO.MIGUEL VIEIRASA GUARISTE.PRIMEIRO COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS COSMOS GAVILAN.DEL.MAR VIEIRASA GUARISTE GUARISTE RAYMI MAR.DE.LOS.SARGAZOS RAYMI GUARISTE RAYMI COSMOS COSMOS COSMOS

2 2

função image( ) - plotar valores selecionados
by Vitor Aguiar 31 Mai '12

31 Mai '12

Boa tarde, alguém poderia me ajudar com a função image( )? Me desculpe por abusar do tempo de vocês, o texto abaixo é grande porém simples. Tenho uma matrix de valores com que estou fazendo um heatmap usando a função image( ). Gostaria de plotar não apenas os valores observados na matrix, mas também valores que seriam possíveis mas não são observados. Isso significa o seguinte: Tenho uma matrix 16x16. Valores possíveis, mas não observados, correspondem a "0" (zeros) onde a soma (índice da linha + índice da coluna) são iguais ou menores que 15. Exemplo: Se eu tenho um valor "0" na posição [14, 3], eu NÃO quero plotar, pois 14 + 3 = 17. Se eu tenho um valor "0" na posição [11, 3], eu quero plotar, pois 11 + 3 = 14. Eu queria que esses valores fossem plotados em cor diferente do meu heatmap, por exemplo, em azul. Matrix: x = c(3045, 893, 692, 830, 617, 155, 246, 657, 105, 60, 18, 7, 7, 4, 2, 11234, 2985, 2242, 2471, 1575, 366, 503, 1283, 170, 79, 32, 6, 4, 1, 3, 19475, 4756, 3233, 3251, 1810, 409, 575, 1210, 139, 41, 11, 4, 2, 0, 0, 20830, 4739, 2990, 2531, 1346, 298, 325, 612, 60, 17, 1, 0, 1, 0, 0, 15304, 3196, 1885, 1440, 610, 117, 115, 185, 14, 2, 0, 0, 0, 0, 0, 8026, 1535, 806, 539, 223, 33, 37, 39, 0, 0, 0, 0, 0, 0, 0, 3300, 562, 286, 141, 45, 14, 5, 12, 0, 0, 0, 0, 0, 0, 0, 1067, 160, 65, 40, 14, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 277, 47, 6, 2, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 72, 6, 1, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 11, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) xmat = matrix(x, ncol = 12) row15 = rep(0,16) col12 = rep(0,15) col13 = rep(0,15) col14 = rep(0,15) col15 = rep(0,15) xmat = cbind(xmat, col12, col13, col14, col15) xmat = rbind(xmat, row15) dimnames(xmat) = list(0:15, 0:15) xmat # minha matrix original é a primeira xmat (15x12), mas estou criando essa outra com linhas e colunas de zeros pois, por exemplo, 0 + 15 = 15, e preciso plotar esse valor também. Mas se alguém tiver solução mais simples... Desde já agradeço, Vitor Vitor Rezende da Costa Aguiar --------------------------------------------------- PhD candidate in biotechnology Rede Nordeste de Biotecnologia UFES, Brazil --------------------------------------------------- Current Adress: Department of Integrative Biology University of California, Berkeley 2033 Valley Life Sciences Building #4134 Berkeley, CA 94720 Phone: 1 (510) 473-2070

1 1

Re: [R-br] Alteração em função de pacote
by Junior Beleti 31 Mai '12

31 Mai '12

Entendido Benilton. Mudando o foco, você pode me sugerir alguma biblioteca paralela para trabalhar com a função eigen? Fiz experimentos com a multicore, mas não obtive grandes resultados. Att, Carlos

6 7

rcorr do Hmisc
by Edson Lira 31 Mai '12

31 Mai '12

Gente, quero fazer uma correlacção entre n variáveis, no Hmisc tem a função rcorr,vejam o exemplo abaixo: require(HmisC) x <- c(-2, -1, 0, 1, 2) y <- c(4, 1, 0, 1, 4) z <- c(1, 2, 3, 4, NA) v <- c(1, 2, 3, 4, 5) rcorr(cbind(x,y,z,v)) No meu banco são 90 variáveis como fzer o teste acima sem ter que declarar uma por uma? Algo do tipo rcorr(cbind(x:v)) Bom dia! Edson Lira Estatístico Manaus-Amazonas

2 1

Retirar dados de uma função no R
by Ana Lúcia Souza Silva Mateus 29 Mai '12

29 Mai '12

Saudações, Tenho a seguinte dúvida, se algúem poder me responder ficaria grata. Estou trabalhando com a seguinte função no R (stdiagn), apresentada abaixo: *> stdiagn function (pts, stkh, stse, stmc = 0, Dzero = FALSE) { oldpar <- par(mfrow = c(2, 2)) pointmap(as.points(pts)) title("Data map") oprod <- outer(stkh$ks, stkh$kt) st.D <- stkh$kst - oprod st.R <- st.D/(stse) if (!Dzero) { persp(stkh$s, stkh$t, st.D, xlab = "Distance", ylab = "Time", zlab = "D", expand = 0.5, ticktype = "detailed", theta = -30, shade = 0.4, cex = 0.7) title("D plot") } else { persp(stkh$s, stkh$t, st.D/oprod, xlab = "Distance", ylab = "Time", zlab = "Dzero", expand = 0.5, ticktype = "detailed", theta = -30, shade = 0.4, cex = 0.7) title("Dzero plot") } plot(outer(stkh$ks, stkh$kt), st.R, xlab = "K(s)K(t)", ylab = "R") abline(h = 0) title("Residual Plot") if (length(stmc) > 1) { hist(c(stmc$t0, stmc$t), nclass = 50, include.lowest = TRUE, xlab = "Test statistic", main = "") y.75 <- par()$usr[4] * 0.75 y.8 <- par()$usr[4] * 0.8 lines(rep(stmc$t0, 2), c(0, y.75), lwd = 2) text(stmc$t0, y.8, "Data Statistic") title("MC results") } else { plot(c(0, 1), c(0, 1), type = "n", axes = FALSE, xlab = "", ylab = "") ad <- par()$adj par(adj = 0.5) text(0.5, 0.5, "No MC available") par(adj = ad) } par(oldpar) } <environment: namespace:splancs> > * Gostaria de obter os valores de * st.D e st.R . Desde já agradeço pelo auxilio. *

5 8

R com php
by Vinicius Brito Rocha 29 Mai '12

29 Mai '12

Senhores, Estou precisando usar o R como uma engine por trás de um portal em desenvolvido em php. Estou precisando de orientação/material. Alguém já fez algum tipo de integração R + web? abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 1

frequência de vetores
by Márcio Augusto Ferreira Rodrigues 29 Mai '12

29 Mai '12

Prezados, tenho um arquivo da seguinte forma z10z11z12z13z20z21z22z23z30z31z32z33 000100100010 010000010010 010000010001 001000100001 010000100100 010000010010 010001000100 cada linha eu considero como um vetor e gostaria de saber a frequência de cada um, pois tenho 1000 linhas e sei que so possui 64 tipos diferentes de vetores. Assim gostaria de saber a frenquência de cada um dos 64 vetores diferentes que tenho . Grato, Márcio

5 6

Survey
by Emerson Cotta Bodevan 28 Mai '12

28 Mai '12

Prezados, preciso criar uma máscara de um questionário para aplicá-lo numa comunidade. Aquela máscara que o Epi-info nos permite criar. O R tem recurso para isso? Senão, alguém teria indicação de um outro software que não seja o Epi-info? É que gostaria de algum que fosse multiplataforma ou que rodasse pelo menos no LINUX Ubuntu. Agradeço qualquer ajuda, -- *Emerson*

6 6

Fwd: Ajuda com vetores
by Fernando Colugnati 28 Mai '12

28 Mai '12

Prezados, tenho a seguinte matriz de estrutura Id a1 p1 p2 p3 p4 p5 p6 p7 1 2 3 4 0 0 0 0 0 1 6 3 0 0 0 0 0 0 2 7 5 3 2 0 0 0 0 2 3 6 0 0 0 0 0 0 Preciso de uma matriz a partir desta da seguinte forma: Id a1 p 1 2 3 1 2 4 1 6 3 2 7 5 2 7 3 2 7 2 2 3 6 Ou seja, preciso gerar quantas linhas forem diferentes de 0 os p1 - p7. É possível fazendo via um apply, sem loop? Abs e obrigado -- Fernando A.B. Colugnati -- Fernando A.B. Colugnati

2 2