Junho 2013 - R-br - DINF UFPR

filtrando dados
by Hélio Gallo Rocha 12 Jun '13

12 Jun '13

Bom dia a todos, Estou com uma dúvida simples, mas não achei na lista, se puderem ajudar, agradeço desde já: ponto=rep(1:10,3) quadra=c(rep(1,10),rep(2,10),rep(3,10)) y=rnorm(1:30) dados=cbind(quadra,ponto,y);dados Dúvida: Teria como selecionar os dados de "y" retirando os dados da quadra 1? -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

3 7

Selecionar valores de um objeto
by Thiago de Andrade 12 Jun '13

12 Jun '13

Boa noite, amigos. Estou trabalhando com o pacote ltm com modelagens de TRI. Tenho a seguinte saída > EscoresAP<-factor.scores(fit_AfetosPositivosR, resp.patterns=DadosAfetosPositivosReduzido) > EscoresAP Call: grm(data = DadosAfetosPositivosReduzido, Hessian = T) Scoring Method: Empirical Bayes Factor-Scores for specified response patterns: Item4 Item7 Item17 Item18 Item19 Obs Exp z1 se.z1 1 1 1 1 1 1 101 89.462 -1.663 0.408 2 6 6 4 2 2 1 0.002 0.181 0.338 3 1 2 2 2 2 4 2.297 -0.809 0.210 4 1 1 1 1 1 101 89.462 -1.663 0.408 5 1 1 1 1 1 101 89.462 -1.663 0.408 6 2 4 2 3 1 1 0.350 -0.552 0.231 Preciso manusear os valores da coluna "z1" (escores dos sujeitos 1 a 6 no teste), realizando uma série de análises com esses valores. Gostaria de saber se há alguma função por meio da qual eu possa extrair ou selecionar os dados dessa coluna, já que eu tenho mais de 1300 casos e não é muito prático copiar um a um. Agradeço desde já. -- Thiago de Andrade Mestrando em Psicologia Social, PPGP-UFPB Núcleo de Pesquisas em Desenvolvimento Sócio-Moral (NPDSM)

4 3

Separar uma matriz por trechos completos de uma sequência
by Cristiano Santos 12 Jun '13

12 Jun '13

Olá pessoal, Alguém conhece uma forma fácil de implementar isso: Eu tenho essa matriz: cbind(c(1,2,3,7,8,9,10,15,16,17), c(1,1,1,1,1,1,1,1,1,1)) E quero separar a matriz pela sequência da primeira coluna, ou seja, obter uma matriz para cada trecho completo da sequência: list(cbind(c(1,2,3), c(1,1,1)), cbind(c(7,8,9,10), c(1,1,1,1)), cbind(c(15,16,17), c(1,1,1))) Ats, Cristiano Santos.

1 0

RES: código R para plotar um coração com animação
by Alan R Panosso 12 Jun '13

12 Jun '13

Caro Walmes, Genial... Amanhã começo funções e o R com o pessoal da graduação, posso utilizar esse seu exemplo para motiva-los? Atenciosamente Alan De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de walmes . Enviada em: terça-feira, 11 de junho de 2013 19:50 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] [OFF-TOPIC] código R para plotar um coração com animação Ok pessoal, valeu aí por disponibilizar. Ainda bem que a resposta veio a tempo de eu terminar a aula (17:30-19:30) então deu para mostrar a função. Os alunos gostaram. Vai fazer sucesso. Segue o script da referida aula. http://www.leg.ufpr.br/~walmes/ensino/ce083-2013-01/ce083-2013-01-aula15i.R 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 <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218 ==========================================================================

1 0

[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