- R-br - DINF UFPR

[OFF-TOPIC] Indicação de livro
by Elton Cavalcanti 20 Jun '12

20 Jun '12

Bom dia, pessoal Vocês poderiam me indicar um bom livro de estatística que um estudante da área de saúde poderia usar, mas que não fosse tão superficial na abordagem de alguns conceitos estatísticos? Abraços, Elton Cavalcanti

2 2

Saída intergen e avaliação com pacote coda
by Dalinne Chrystian Carvalho dos Santos 20 Jun '12

20 Jun '12

Olá pessoal, tudo bem? Meu nome é Dalinne e eu tenho uma saída do programa intergen com análise Bayesiana no formato "txt" da seguinte maneira: 1 27.90086641691324 -42.05430078214862 6609.247759721497 5.547673239405648 -109.7603421548792 6665.595598994243 2 24.47414461350570 -37.16170819212519 8015.513063087009 5.783134528755414 -121.6600727394209 7681.170349762655...10000 5.741506313680731 -78.70398130578000 5932.380770110504 3.514573296713659 -109.4777360152619 6552.098731256713 Ou seja, as minhas 10.000 interações estão em apenas uma coluna.Queria saber como entrar com esse banco de dados no R mas como um objeto mcmc e depois disso como eu faço paraachar o ponto de burn-in, por meio do pacote coda do R? Obrigada pela atenção, Dalinne Chrystian Carvalho dos Santos

1 0

Pacote DiagnosisMed
by Edson Lira 20 Jun '12

20 Jun '12

Ao tentar instalar o pacote DiagnosisMed veriquei que não consta mais no CRAN. Alguém sabe me explicar o que houve? Qual pacote que pode ser utilizado para as mesmas análises que o mesmo fazia? Bom dia a todos os R-lista []'s. Edson Lira Estatístico Manaus-Amazonas

3 3

visualização cluster
by Vinicius Brito Rocha 20 Jun '12

20 Jun '12

Pessoal, estou plotando resultados de cluster kmeans em duas componentes e então colorindo os grupos, mas estava querendo fazer uma coisa diferente. Gostaria de plotar as fronteiras resultantes dos clusters. Mais ou menos como desenhar as regiões de fronteira de cada cluster. Alguém tem alguma idéia? Abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* 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 2

Testes Estacionariedade Séries Temporais
by Josenilson Gomes 20 Jun '12

20 Jun '12

Olá, Gostaria de saber se alguém já usou algum teste de estacionariedade no R. Tenho interesse em usar os seguintes testes, mas não sei se já estão disponíveis no R: -Dickey-Fuller -Phillips-Perron -KPSS Todos testam se a série é estacionária... Agradeço desde já. Att., Josenílson

5 5

Simulação de distribuição bivariada
by Rodrigo Coster 19 Jun '12

19 Jun '12

Caros, o R tem alguma função para simulação de distribuições bivariadas? Seguem as funções densidades e distribuição da distribuição que tenho intesse em simular (u e v são as variaveis, tauU e tauL os parametros, todos no intervalo (0,1)) dsjc <- function(u,v,tauU,tauL) { ## Adaptado do Código para Matlab de Andrew Patton - http://publicecondukeedu/~ap172/codehtml k1 = 1/log2(2-tauU) k2 = -1/log2(tauL) CL1 = ((1 - (1 - u)^k1)^(k2 - 1)* (1 - u)^(k1 - 1)*(-1 + k1*(k2* (-1 + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1))) + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1))))* (1 - (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(-k2^(-1)))^(k1^(-1))* (1 - (1 - v)^k1)^(k2 - 1)* (1 - v)^(k1 - 1)) CL2 = (((-1 + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1)))^2) * ((1 - (1 - u)^k1)^k2 + (1 - (1 - v)^k1)^k2 - (1 - (1 - u)^k1)^k2* (1 - (1 - v)^k1)^k2)^2) CL1 = CL1/CL2 k1 = 1/log2(2-tauL) k2 = -1/log2(tauU) u = 1-u v = 1-v CL3 = ((1 - (1 - u)^k1)^(k2 - 1)* (1 - u)^(k1 - 1)*(-1 + k1*(k2* (-1 + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1))) + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1))))* (1 - (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(-k2^(-1)))^(k1^(-1))* (1 - (1 - v)^k1)^(k2 - 1)* (1 - v)^(k1 - 1)) CL4 = (((-1 + (-1 + (1 - (1 - u)^k1)^(-k2) + (1 - (1 - v)^k1)^(-k2))^(k2^(-1)))^2) * ((1 - (1 - u)^k1)^k2 + (1 - (1 - v)^k1)^k2 - (1 - (1 - u)^k1)^k2* (1 - (1 - v)^k1)^k2)^2) CL3 = CL3/CL4 CL = 0.5*(CL1+CL3) return(CL) } psjc <- function(U,V,tauU,tauL) { ## Adaptado do Código para Matlab de Andrew Patton - http://publicecondukeedu/~ap172/codehtml K = 1/log2(2-tauU); G = -1/log2(tauL); out1 = 1-((1-(((1-((1-U)^K))^(-G))+((1-((1-V)^K))^(-G))-1)^(-1/G))^(1/K)); K = 1/log2(2-tauL); # switching the upper and lower measures G = -1/log2(tauU); U = 1-U; V = 1-V; out2 = (1-U) + (1-V) - 1 + 1-((1-(((1-((1-U)^K))^(-G))+((1-((1-V)^K))^(-G))-1)^(-1/G))^(1/K)); out1 = 0.5*(out1+out2); return(out1) }

3 4

comando no R
by Eliana Silva 19 Jun '12

19 Jun '12

Pessoal, alguem me pode dizer se existe um comando no R uma t_student p_variada? e para a gama inversa tb gostaria. obrigada

2 1

Ajuda em uma função
by Sérgio Henrique almeida da silva ju 19 Jun '12

19 Jun '12

Olá Estou com mais uma dúvida no desenvolvimento de uma função Preciso analisar um conjunto de dados, porém essa análise é feita a cada duas colunas: a = c(1,3,5,4,3,2,1,4,5,3,4,5,2,1,1,3,4,5,5,3,2) a1 = c(4,5,5,4,4,3,2,1,4,3,4,5,3,2,3,4,5,5,4,4,1) d = c(3,3,4,4,5,4,3,2,4,5,5,3,4,4,5,4,3,2,3,4,5) d1 = c(2,2,3,3,4,3,2,3,4,5,5,5,4,3,2,3,4,5,4,3,2) dados = data.frame(a,a1,d,d1) library(psy) x=kappa(dados[,c(1,2)]) x y=kappa(dados[,c(3,4)]) y z=kappa(dados[,c(n,n+1)]) z Porém gostaria de rodar uma função que calcule isso automaticamente, dando todos os resultados, para isso eu daria a coluna de inicio e a coluna final. Como posso fazer? -- 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

3 7

Duvida unzip no R
by Humberto Hazin 19 Jun '12

19 Jun '12

Olá Pessoal Tenho uma pasta que contem vários arquivos com extensão bz2! Gostaria de saber como descomprimir no R. Alguém pode me ajudar com um exemplo? Humberto

3 4

Script Matriz de Vizinhança
by Isabela Zardo 19 Jun '12

19 Jun '12

Gostaria de scripts pra matriz de vizinhança, alguém pode eme ajudar? Isabela Zardo

3 4