Julho 2011 - R-br - DINF UFPR

problema rJava - xlsx
by Vinicius Brito Rocha 15 Jul '11

15 Jul '11

Srs. estou tentando instalar o pacote xlxs, mas recebo uma mensagem de erro no rJava. O que pode estar acontecendo? erro: Error : não é possível alocar vetor de tamanho 914.6 Mb Erro: carregamento do pacote/espaço de nomes falhou para 'rJava' o pacote tem apenas 650 kb -- 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 2

Re: [R-br] AJUDA - amova (package ade4)
by Ivan Bezerra Allaman 14 Jul '11

14 Jul '11

Olá Ricardo! Forneça-nos um CMR com uma pequena amostra dos seus dados para que possamos ajudá-lo. O erro me parece que está na forma de como vc está declarando no modelo. Pelo erro, é provável que não exista variação de 'samples' dentro de 'regions'. Veja melhor o seu modelo, e verifique se realmente há a necessidade de se colocar esta fonte de variação (se ela existe mesmo)! Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Re: [R-br] Algoritmo do boot.ci do livro do Davison, A.C. and Hinkley, D.V. (1997)!
by Ivan Bezerra Allaman 14 Jul '11

14 Jul '11

Rubem, A função corr não é uma função específica para ser utilizada com o pacote boot. Pode ser utilizada em qualquer outra situação. Vide ?corr. A diferença da corr para cor está apenas num vetor de pesos que pode ser utilizado na função corr. Caso este vetor não seja declarado, o resultado da corr = cor. Quanto as processos de reamostragem, e lendo com mais atenção o artigo do DiCicci & Efron 1996, você tem toda a razão. O modo como estava sendo feito, estava equivocado, pois a distribuição bivariada ocorre para o par XiYi. Valeu Rubem, são essas chamadas de atenção que enriquecem nossos conhecimentos. Muito obrigado!! Abraço! Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

AJUDA - amova (package ade4)
by Ricardo Fernandes Campos Junior 14 Jul '11

14 Jul '11

Olá, pessoal. Estou tentando usar a função "amova" do pacote ade4, mas não consigo. Comparei a estrutura dos meus objetos com a dos objetos do HELP da função, mas me parece tudo ok. A função roda tudo bem, mas recebo o seguinte resultado: amova(samples,sqrt(distance),structure) $call amova(samples = samp3, distances = sqrt(distance), structures = structurss) $results Df Sum Sq Mean Sq Between regions 2 0.1371102 0.06855508 Between samples Within regions 0 0.0000000 NaN Within samples 41 1.7349423 0.04231567 Total 43 1.8720525 0.04353610 $componentsofcovariance Sigma % Variations Between regions NaN NaN Variations Between samples Within regions NaN NaN Variations Within samples 0.04231567 NaN Total variations NaN NaN $statphi Phi Phi-samples-total NaN Phi-samples-regions NA Phi-regions-total NaN Abaixo o resultado de 'str' dos meus objetos. str(samples) 'data.frame': 19 obs. of 3 variables: $ AF : int 0 7 1 4 0 1 2 4 1 1 ... $ Caatinga: int 1 0 0 0 0 0 0 0 0 0 ... $ Cerrado : int 0 0 0 0 1 0 0 0 0 0 ... str(distance) Class 'dist' atomic [1:171] 0.0411 0.1007 0.0919 0.1007 0.0919 ... ..- attr(*, "Labels")= chr [1:19] "4" "12" "17" "19" ... ..- attr(*, "Size")= int 19 ..- attr(*, "call")= language as.dist.default(m = mat) ..- attr(*, "Diag")= logi FALSE ..- attr(*, "Upper")= logi FALSE str(structure) 'data.frame': 3 obs. of 1 variable: $ regions: Factor w/ 3 levels "AF","Caatinga",..: 1 2 3 -- Abraços, Ricardo Fernandes Campos Junior

1 0

conexao teradata
by Vinicius Brito Rocha 14 Jul '11

14 Jul '11

Srs. alguem sabe como funciona a conexao R ao banco de dados teradata? pelo que vi nao é via odbc. 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 3

Re: [R-br] Algoritmo do boot.ci do livro do Davison, A.C. and Hinkley, D.V. (1997)!
by Ivan Bezerra Allaman 14 Jul '11

14 Jul '11

Verdade!! Valeu. Agora uma coisa que estou quebrando a cabeça é a seguinte. Porque que o boot gerado pela função boot está extremamente diferente do boot gerado na 'unha'?. Vejamos um CMR. x <- c(0.75,2.55,1.40,0.80,2.60) y <- c(0.8249,0.797,0.7762,0.7797,0.7483) da <- data.frame(x,y) corr(da) mean(boot(da,corr,10000)$t) #Esta função eu fiz para evitar NAs, já que algumas amostragens são todas iguais, gerando um desvio padrão zero. sa <- function(x){ y <- sample(x,replace=T) if(y[1:1] == y[2:2] && y[1:1] == y[3:3] && y[1:1] == y[4:4] && y[1:1] == y[5:5]){ z <- sample(x,replace=T)}else{ y} } boot_da <- numeric(10000) n <- 10000) for(i in 1:n){ boot_da[[i]] <- corr(matrix(c(sa(da[,1]),sa(da[,2])),ncol=2)) } boot_da mean(boot_da) Se alguém puder ajudar de novo, agradeço. Abraços! Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Ajuda: tratamento de erros no R
by Jose Claudio Faria 14 Jul '11

14 Jul '11

Pessoal, Quem tiver experiência com tratamento de erros no R poderia fazer o favor de me colocar na veia boa: funções, pacotes, material didático, links, etc? Um pesquisador (aqui da UESC) está precisando proteger um trecho de seu script contra erros. Ele já detectou (e mapeou) os possíveis erros, mas não sabe como tratar. Ou seja, registrar o erro e continuar (quando possível) com suas simulações (bem demoradas relecionadas ao estudo de estabilidade de galáxias - Astrofísica). Minha experiência com tratamento de erros no R é bem simples (detectar erro, dar mensagem e parar) e insuficiente para sua necessidade. Nesse caso é necessário uma abordagem mais consistente e efetiva: 1. Proteger o bloco onde ocorrem os erros; 2. Uma vez que ocorrer verificar qual foi o erro; 3. Registrar sua ocorrência (por exemplo num arquivo_erros.txt) e prosseguir (quando/se ppossível) Toda boa informação é bem vinda! Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Pleno UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

3 2

Algoritmo do boot.ci do livro do Davison, A.C. and Hinkley, D.V. (1997)!
by Ivan Bezerra Allaman 14 Jul '11

14 Jul '11

Bom dia senhores!! Alguém que tem o livro do Davison, A.C. and Hinkley, D.V. (1997) poderia me disponibilizar o algorítimo para cálculo do intervalo de confianção BCa. A função boot.ci da library boot calcula estes intervalos com o parâmetro viesado, e eu não quero. Não consegui modificar a função também, pois ela está amarrada ao output da função boot. Desde já obrigado! Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

4 3

Erro ao ler banco no mac os
by Edson Lira 14 Jul '11

14 Jul '11

Bom dia, lista. Alguém sabe como corrigir esse erro? Ou o que é esse erro? > ag<-read.xls("/Users/edsonlira/IBGE/Censo_2010/Agregado/Agre2010_sinopse_AM.xls",sheet=1) Erro em make.names(col.names, unique = TRUE) : string multibyte inválida 2 > ag<-read.csv2("/Users/edsonlira/IBGE/Censo_2010/Agregado/Agre2010_sinopse_AM.csv") Erro em make.names(col.names, unique = TRUE) : string multibyte inválida 2 Edson Lira Estatístico Manaus-Amazonas

1 1

Dúvida para montar um agrupamento
by Clayton Santos Delfino 14 Jul '11

14 Jul '11

Olá pessoal, boa noite. Fico feliz por ter encontrado uma lista de R em português pois estou começando a aprender e tenho várias dúvidas e espero contar com a preciosa ajuda de todos. Eu já consegui realizar algumas análises com R mas agora empaquei com o seguinte: tenho um data.frame com a seguinte estrutura: identificador municipio regiao principal 1 Sao Paulo 1 Sim 1 Sorocaba 1 Não 1 Curitiba 2 Não 1 Sao Paulo 1 Sim 1 Porto Alegre 3 Não 1 Cuiaba 4 Não 2 Sao Paulo 1 Não 2 Sorocaba 1 Não 2 Curitiba 2 Sim 2 Curitiba 2 Sim 3 Porto Alegre 3 Sim 3 Porto Alegre 3 Sim 3 Porto Alegre 3 Sim 3 Sorocaba 1 Não Eu o que eu estou tentando fazer e ainda não encontrei um meio seria gerar uma tabela com uma linha para cada identificador onde somente o municipio com maior frequencia, considerando-se "principal = Sim" fosse retornado, exemplo: identificador Municipio Principal (frequencia) 1 São Paulo 2 2 Curitiba 2 3 Porto Alegre 3 Se alguem puder me ajudar agredeço desde já. Um abraço a todos. Clayton

7 10