- R-br - DINF UFPR

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

Problema com o pacote amélia no debian 5.0 amd64:
by GESER 13 Jul '11

13 Jul '11

Gostaria de fazer uma imputação no linux através do pacote Amélia II versão 1.2-18 e não estou conseguindo. O estranho é que o script roda perfeitamente no Windows. O script está abaixo: #### Início do script #### #Carrego o arquivo para ser imputado# umi99 <- read.table("data/umid99NA.txt", header=TRUE)[,-1] #Carrego o pacote Amélia# require(Amelia) nrep <- 3 system.time(a.out <- amelia(umi99, m=nrep)) #Há mais linhas de comando adiante, mas o problema acontece na linha acima# #### Fim do script #### O erro acontece sempre ao final de uma imputação e diz: Erro em dimnames(impdata$covMatrices)[[1]] <- prepped$theta.names : comprimento de 'dimnames' [1] deve ser o mesmo de 'dims' [3] Se necessário, posso enviar uma imagem do prompt anunciando o erro e o arquivo que carrego para os dados serem imputados. Obrigado, Fabio

4 5

Re: [R-br] Ajuda (possivelmente simples) para reestruturar um data.frame
by Jose Claudio Faria 13 Jul '11

13 Jul '11

Walmes e Benilton, Não deu tempo nem de eu ir tomar um café!!! ;) Muito elegantes as soluções: muito obrigado! Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Pleno UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 1

[Dúvida] Como plotar valores
by Pedro Rafael 13 Jul '11

13 Jul '11

Pessoal, como plotar o valor dos pontos em um gráfico de linha no R? Exemplo: x = 1:10 plot(x,type = "o",pch = 12) Isso vai me gerar um gráfico com 10 pequenos desenhos. Como colocar o seu valor próximo desses desenhos? Pedro

3 2