- R-br - DINF UFPR

[OFF-TOPIC] código R para plotar um coração com animação
by walmes . 11 Jun '13

11 Jun '13

Saudações pessoal, Alguém tem cópia daquele código R, uma vez compartilhado na lista, acho que na época do R-STAT ainda. O código que falo é aquele que tem um laço for(), as cores mudam com a execução. Dou aula sobre funções hoje e quinta e queria motivar os alunos à presentear suas respectivas com um CMR amoroso. Ou provocar logo a separação dependendo do senso de humor e diferença de nerdismo do casal. Versões simples eu encontrei em: http://stackoverflow.com/questions/8082429/plot-a-heart-in-r http://blog.revolutionanalytics.com/2013/02/make-a-valentines-heart-with-r.… http://rud.is/b/2013/02/14/happy-valentines-day-mrshrbrmstr/ Grato. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

3 3

Problema para instalar pacote rgdal no R 3.0.0 em Ubuntu 12.04
by ASANTOS 11 Jun '13

11 Jun '13

Boa noite pessoal, Não consigo instalar o pacote rgdal no R 3.0.0 no Ubuntu, quando utilizo: install.packages("rgdal") ocorre o segundo problema: ** preparing package for lazy loading Error : package 'lattice' was built before R 3.0.0: please re-install it ERROR: lazy loading failed for package 'sp' * removing '/home/asantos/R/x86_64-pc-linux-gnu-library/3.0/sp' Warning in install.packages : installation of package 'sp' had non-zero exit status ERROR: dependency 'sp' is not available for package 'rgdal' * removing '/home/asantos/R/x86_64-pc-linux-gnu-library/3.0/rgdal' Warning in install.packages : installation of package 'rgdal' had non-zero exit status The downloaded source packages are in '/tmp/Rtmp5gpcLF/downloaded_packages' Alguém poderia me ajudar e ver por que esta ocorrendo este erro Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

2 1

lm e NA
by Daniela Recchia 11 Jun '13

11 Jun '13

Bom dia, estou com dificuldades em rodar um modelo no R que acredito ser devido aos NA. Apesar do software aparentemente remover os NA's na hora de rodar a funcao lm, quando ajusto meu modelo model<-lm(cpstotalnorm~indconfident+indpride+indenergy+indgrowth+indseethingsnew+indinputstrategy +indtrainingavailable+grppointsofview+grpconflictresol+grprightpeople+grpdifferentviews+grprethink +grplessonss) recebo um erro "error contrasts can be applied only to factors with 2 or more levels", pesquisei sobre o erro e vi que poderiam ser dois casos: ou as variáveis poderiam estar definidas como caracteres e nao como fatores - o que nao é o caso, ou há NA's no conjunto de dados - é o caso. Tentei retirar entao os NA's dos dados para rodr o modelo, porém usando a funcao: datanew<-na.omit(data) o R retira todas as observacoes, e nao só o NA, por exemplo: > datanew$indseethingsnew factor(0) Levels: 1 2 3 4 5 6 7 enquanto se retiro somente os NA's dessa variável, no conjunto de dados original, ele mantém todas as outras observacoes: > c<-na.omit(data$indseethingsnew) > c [1] 2 2 2 3 1 3 1 2 1 1 2 3 3 3 3 2 2 1 3 4 1 1 4 3 2 3 1 2 2 1 5 2 1 3 1 1 1 6 3 1 2 1 1 1 1 2 1 2 1 3 1 [52] 3 2 5 1 5 1 1 3 2 1 2 2 1 2 2 1 1 2 2 1 2 1 2 2 1 1 2 2 3 2 2 1 2 2 2 2 3 1 3 2 2 1 3 1 2 2 1 2 3 2 1 [103] 2 1 2 4 2 1 1 2 2 1 3 3 2 2 1 1 6 2 1 2 4 2 3 3 2 2 2 2 1 3 6 1 1 3 2 3 3 1 2 3 2 2 7 2 3 2 1 3 4 1 2 [154] 3 2 1 3 attr(,"na.action") [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 [26] 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 [51] 85 86 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 [76] 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 [101] 211 212 213 214 215 216 217 218 219 attr(,"class") [1] "omit" Levels: 1 2 3 4 5 6 7 Entao meu novo conjunto de dados datanew contém apenas zeros! Alguém teria alguma sugestao? Obrigada, Daniela -- Daniela Rodrigues Recchia M.Sc. in Statistics Tel. (Mobil): +49 01578 7588382 Kißlegg - Germany “The standard of success in life is absolutely the amount of joy you feel” Esther & Jerry Hick

5 7

Re: [R-br] Remover NA: PCA
by Helena Turon 09 Jun '13

09 Jun '13

No meu caso era um vetor, que eu queria plotar o boxplot sem o NA. Assim poderia aparecer os pontos discrepantes. Mas, você tem razão se é uma matriz afetará o resultado. Aí deu certo. Helena

2 1

Remover NA: PCA
by Vinícius Lionel Mateus 09 Jun '13

09 Jun '13

Prezados, Estou executando o seguinte CRM: read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> TSF > TSF > names(TSF) [1] "Date" "Label" "PMF" "BC" "Acet" "Form" "Cl." "NO3." "PO43." "SO42." "Na." "K." [13] "Mg2." "Ca2." "NH4." > TSFa <- TSF[,2:15] > model <- prcomp(TSFa, scale = TRUE) Erro em colMeans(x, na.rm = TRUE) : 'x' deve ser numérico > model <- prcomp(TSFa, na.rm = TRUE, scale = TRUE) Erro em colMeans(x, na.rm = TRUE) : 'x' deve ser numérico Tentei diferentes métodos para remover o "NAs", mas não obtive sucesso. Alguém tem alguma dica? -- Vinícius -- -- Atenciosamente, VINÍCIUS LIONEL MATEUS, M.Sc (http://lattes.cnpq.br/6501001637020665) Bacharel em Química - Doutorando em Química Analítica Laboratório de Química Atmosférica - Departamento de Química Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900 Telefone: (+55) (21) 3527-1327 (+55) (21) 9358-8051 www.puc-rio.br

4 10

Remover NA: PCA
by Helena Turon 08 Jun '13

08 Jun '13

Depois de transformar a variável para numérica eu ordenei, sem NA, deu certo. Tsfa$x = as.numeric(Tsfa$x) Tsfa$x <- sort(Tsfa$x[!is.na(Tsfa$x)]) Boa sorte. Helena

2 1

[Dúvida] Estimação por Máxima Verossimilhança.
by Pedro Rafael 07 Jun '13

07 Jun '13

Senhores tenho uma dúvida. Na verdade não é dúvida, apenas quero sugestões. Tenho uma função densidade de probabilidade "complicada". Trate-se de uma distribuição chamada Kwmarashwamy Weibull Poisson. Tenho alguns bancos de dados e gostaria de verificar o ajustamento dessa distribuição à estes bancos de dados. Estou estimando os parâmetros pelo método de máxima verossimilhança. Essa distribuição tem suporte nos reais positivos (x>0) e todos os seus parâmetros são positivos. Optei em utilizar o método L-BFGS-G para restringir a busca nos reais positivos. Segue abaixo o comando. Nesse exemplo não houve convergência. Percebi que os chutes iniciais influenciam muito as estimativas dos parâmetros no caso em que há convergência. Usando métodos de maximização diferentes em muitos casos há grandes diferenças nas estimativas. Existe alguma forma mais tranquila e direta para encontrar as estimativas pelo método de máxima verossimilhança em R? O código segue abaixo: vero <- function(par,x){ a = par[1] b = par[2] c = par[3] lambda = par[4] beta = par[5] -sum(log((a*b*c*lambda*(beta^c)*(x^(c-1))*((1-exp(-(x*beta)^c))^(a-1)) * ((1-(1-exp(-(beta*x)^c))^a)^(b-1)) * exp(-lambda*(1-(1-(1-exp(-(beta*x)^c))^a)^b) - (beta*x)^c))/(1-exp(-lambda)))) } dados = c(17.23, 28.92, 33.00, 41.52, 42.12, 45.60, 48.80, 51.84, 51.96, 54.12, 55.56, 67.80, 68.64, 68.64,68.88, 84.12, 93.12, 98.64, 105.12, 105.84, 127.92, 128.04, 173.40) optim(par=c(1,1,1,1,1),fn=vero, method="L-BFGS-B",x=dados/1000, lower=c(0.001,0.001,0.001,0.001,0.001), upper=c(Inf,Inf,Inf,Inf,Inf)) O engraçado é que essa distribuição é encaixada com a distribuição Weibull. Quando tento ajustar a Weibull à esses dados há convergência. O que vocês acham que devo fazer? [ ], Pedro Rafael Diniz Marinho.

6 12

Inserir NA em objeto do tipo lista
by Aline Lipsky 07 Jun '13

07 Jun '13

Bom dia, Eu tenho um objeto do tipo lista com mais de 6 mil vetores com tamanhos diferentes.Como faço para acrescentar NA nos vetores de tamanho menor e transformar esse objeto em um data frame? Grata, Aline

2 1

encoding de páginas web readLines XML
by walmes . 07 Jun '13

07 Jun '13

Saudações, Estou fazendo um estudo com vídeos do youtube. A partir dos links eu retiro informações deles, como o título. No entanto, não estou conseguindo acertar o encoding, seja para ler corretamente a página ou para tranformar as strings após a leitura. Com o CMR abaixo eu leio o código de um vídeo do youtube e pego seu título, o problema é que os acentos não são corretamente representados, e eu preciso que sejam. Tentei atribuir um encoding na readLines() mas nenhum dos usados (latin1 e utf-8) funcionaram. Tentei usar a iconv() mas sem sucesso. Uso ubuntu 12.04. Sessão R está no CMR. library(XML) # lê a página do vídeo com esse endereço lines <- readLines("http://www.youtube.com/watch?v=P3m6rJR2yDI") # lines <- readLines("http://www.youtube.com/watch?v=P3m6rJR2yDI", encoding=?) # organiza h <- htmlTreeParse(lines, asText=TRUE, useInternalNodes=TRUE) summary(h) # pega o título do vídeo n <- unlist(getNodeSet(doc=h, path="//title", fun=xmlValue)) n Encoding(n) # retira os caracteres estranhos substituindo por vazio iconv(nap, from="UTF-8", to="ASCII", sub="") # iconv(n, from=?, to=?) *# Como sai # "MatemÃ¡tica - EstatÃstica: Conceito e frequÃªncias - YouTube" # Como eu gostaria # "Matemática - Estatística: Conceito e frequências - YouTube"* ## > sessionInfo() ## R version 3.0.1 (2013-05-16) ## Platform: i686-pc-linux-gnu (32-bit) ## locale: ## [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C ## [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8 ## [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=pt_BR.UTF-8 ## [7] LC_PAPER=C LC_NAME=C ## [9] LC_ADDRESS=C LC_TELEPHONE=C ## [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C ## attached base packages: ## [1] stats graphics grDevices utils datasets methods base ## other attached packages: ## [1] XML_3.96-1.1 ## loaded via a namespace (and not attached): ## [1] compiler_3.0.1 tools_3.0.1 ## > Alguém saberia solucionar esse problema? Desde já grato. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

4 4

Adicionar Na em uma lista
by Aline Lipsky 06 Jun '13

06 Jun '13

Boa tarde, Eu tenho um objeto do tipo lista com mais de 6 mil vetores com tamanhos diferentes. Eu consigo acrescentar o NA em um elemento da lista: especies[[1]][[7]] <- NA Se alguém puder me ajudar, como faço para repetir esse passo em todos os elementos da minha lista? Tentei com a função sapply, mas não dá certo. Grata, Aline

5 5