Junho 2012 - R-br - DINF UFPR

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 20 Jun '12

20 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 20 Jun '12

20 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

Extrair intervalo de uma coluna...
by andrebvs 19 Jun '12

19 Jun '12

2 2

Extrair intervalo de uma coluna...
by andrebvs 18 Jun '12

18 Jun '12

1 0