Novembro 2013 - R-br

Manual
by Maicon J. Galiazzi 13 Nov '13

13 Nov '13

Olá, estou desenvolvendo um manual de fácil utilização ao usuário, e no momento estou descrevendo a interfacedo software R Stuio 0.97.551, assim como a barra de ferramentas do mesmo. Podem me indicar algum material no qual posso conferir? Desde já agradeço! Att Maicon J. Galiazzi

2 2

ordenar um data.frame
by Aline Lipsky 13 Nov '13

13 Nov '13

Boa tarde a todos, Gostaria de uma ajuda com o seguinte: Tenho um data frame como esse: Animal Ordem Data Especie Sexo Peso 1 47 1 54 1 55 2 52 2 54 2 56 1 53 1 52 1 58 3 57 3 56 3 59 A estrutura é a seguinte: data.frame': 20 obs. of 6 variables: $ Animal : Factor w/ 10 levels "Barcelona","Baronesa",..: 1 4 7 1 4 7 3 2 10 8 ... $ Ordem: int 1 1 1 2 2 2 1 1 1 1 ... $ Data: Factor w/ 19 levels "01/01/2007","01/02/2004",..: 5 4 2 19 13 8 10 11 15 12 ... $ Especie : Factor w/ 5 levels "","",..: 1 2 3 1 4 5 5 4 4 5 ... $ Sexo : Factor w/ 2 levels "F","M": 2 1 1 1 1 2 2 1 2 2 ... $ Peso: num 47 54 55 52 54 56 53 52 58 57 ... Preciso ordenar primeiro todo o data.frame em ordem crescente de acordo com a ordem e depois ordenar todo o data.frame deacordo com o peso em ordem crescente. Agradeço pela ajuda, Aline

2 2

Processar arquivo .rnw
by Andre Oliveira 13 Nov '13

13 Nov '13

Senhores estou meio enferrujado aqui no Sweave. Fiz um print aqui e não deu o esperado. <<>>= setwd("/media/Meus arquivos/Pendrive/KINGSTON/sweave/apostilaR2") require(gdata) dad<-read.xls("dados.xls") str(dad) attach(dad) require(xtable) test=xtable(dad, align = "llllll") print(test, table.placement = "H") @ Saida impressa no pdf. \begin{table}[ht] \centering \begin{tabular}{llllll} \hline & sexo & esc & sal & idad & filhos \\ \hline 1 & masculino & 1 grau & 4.00 & 26 & 0 \\ 2 & masculino & 1 grau & 4.56 & 32 & 1 \\ 3 & feminino & 1 grau & 5.25 & 36 & 2 \\ 4 & masculino & 2 grau & 5.73 & 20 & 0 \\ 5 & feminino & 1 grau & 6.26 & 40 & 1 \\ 6 & masculino & 1 grau & 6.66 & 28 & 1 \\ 7 & feminino & 1 grau & 6.86 & 41 & 0 \\ 8 & masculino & 1 grau & 7.39 & 43 & 0 \\ 9 & feminino & 2 grau & 7.59 & 34 & 1 \\ 10 & masculino & 2 grau & 7.84 & 23 & 0 \\ 11 & feminino & 2 grau & 8.12 & 33 & 2 \\ 12 & masculino & 1 grau & 8.46 & 27 & 0 \\ 13 & feminino & 2 grau & 8.74 & 37 & 0 \\ 14 & masculino & 1 grau & 8.95 & 44 & 3 \\ 15 & feminino & 2 grau & 9.13 & 30 & 0 \\ 16 & masculino & 2 grau & 9.35 & 38 & 3 \\ 17 & feminino & 2 grau & 9.77 & 31 & 1 \\ 18 & masculino & 1 grau & 9.80 & 39 & 2 \\ 19 & feminino & superior & 10.53 & 25 & 1 \\ 20 & masculino & 2 grau & 10.76 & 37 & 0 \\ 21 & feminino & 2 grau & 11.06 & 30 & 1 \\ 22 & masculino & 1 grau & 11.59 & 34 & 2 \\ 23 & feminino & 2 grau & 12.00 & 41 & 4 \\ 24 & masculino & superior & 12.79 & 26 & 0 \\ 25 & feminino & 2 grau & 13.23 & 32 & 2 \\ 26 & masculino & 2 grau & 13.60 & 35 & 2 \\ 27 & feminino & 1 grau & 13.85 & 46 & 0 \\ 28 & masculino & 2 grau & 14.69 & 29 & 2 \\ 29 & masculino & 2 grau & 14.71 & 40 & 5 \\ 30 & masculino & 2 grau & 15.99 & 35 & 2 \\ 31 & feminino & superior & 16.22 & 31 & 3 \\ 32 & masculino & 2 grau & 16.61 & 36 & 1 \\ 33 & feminino & superior & 17.26 & 43 & 3 \\ 34 & masculino & superior & 18.75 & 33 & 4 \\ 35 & feminino & 2 grau & 19.40 & 48 & 2 \\ 36 & masculino & superior & 23.30 & 42 & 3 \\ \hline \end{tabular} \end{table} Não se espera que imprima uma tabela com linhas? O que fiz de errado? André Oliveira Souza

5 6

Customização de gráficos com a função plot() [Resolvido]
by Tiago Souza Marçal 13 Nov '13

13 Nov '13

Consegui aqui. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

3 3

Contrastes em modelos de regressão de poisson inflacionados por zeros
by ASANTOS 13 Nov '13

13 Nov '13

Caros membros, Gostaria de fazer a junção de termos qualitativos não significativos, através de análises de contraste de modelos, para verificar a semelhança entre tratamentos, no entanto, usando modelos de regressão de poisson inflacionados por zeros, aproveitando um modelo do Walmes fiz: #------------------------------------------------------------------ # Definições da sessão. rm(list=ls()) require(pscl) require(multcomp) require(lattice) require(latticeExtra) #------------------------------------------------------------------ # Dados artificiaiscom 3 tratamentos e 10 tempos da <- expand.grid(trat=gl(3,1), tempo=1:10) X <- model.matrix(~trat+tempo, da); ncol(X) betas <- c(0.1,0.1,0.3,0.5) eta <- X%*%betas y1 <- rpois(da$trat, lambda=exp(eta)) y2 <- rbinom(y1, size=1, prob=0.7) da$y <- y1*y2 str(da) xyplot(y~tempo|trat, data=da, jitter.x=TRUE) #------------------------------------------------------------------ # Ajuste do modelo. compl.mod <- zeroinfl(y~trat+tempo|trat, data=da) summary(compl.mod) #------------------------------------------------------------------ # Contrate de modelos de poisson inflacionado com zeros sort(tapply(da$y,da$trat,mean))#Ordena as médias trat2<-da$trat levels(trat2) levels(trat2)[1]<-"T1eT3" levels(trat2)[3]<-"T1eT3" levels(trat2) reduc.mod<-zeroinfl(y~trat2+tempo|trat, data=da) # Comparando o modelo completo e o ajunte com junção dos tratamentos 1 e 3 pchisq(-2*(logLik(compl.mod)-logLik(reduc.mod)),df=3,lower.tail=FALSE) # Na minha simulação: #[1] 1 #attr(,"df") #[1] 7 #attr(,"class") #[1] "logLik" Então pergunto: T1 e T3 são iguais e podem ser representados por um único modelo, isto esta correto? Ou existe alguma outra abordagem? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção 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 Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

3 9

Análise discriminante
by Gustavo Faria 13 Nov '13

13 Nov '13

Pessoal, estou precisando fazer uma análise discriminante, qual pacote posso usar? Abraços -- --- Gustavo de Faria Silva Biólogo herpetólogo Crbio 87421/04-P Pós-Graduado em Consultoria e Licenciamento Ambiental Mestrando em Zoologia de Vertebrados Telefone: (31) 9579-8473 (Vivo)

3 3

[OFF-TOPIC] Apostila Análise multivaria para os cursos de administração
by ASANTOS 13 Nov '13

13 Nov '13

Prezados membros, Gostari de saber se alguém dispõem em pdf da Apostila Análise multivaria para os cursos de administração de autoria de Marcos F. Silva, Obrigado, -- ====================================================================== Alexandre dos Santos Proteção 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 Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

4 3

Customização de gráficos com a função plot()
by Tiago Souza Marçal 12 Nov '13

12 Nov '13

Boa noite pessoal, gostaria de fazer um gráfico de dispersão com o CMR abaixo, entretanto, gostaria de variar o símbolo no parâmetro pch a cada tratamento e se possível a cor do símbolo. dados<-structure(list(trat = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L), .Label = c("1", "2", "3", "4"), class = "factor"), y = c(11.8807183390169, 12.9694790587891, 13.5719432329935, 14.1495411288071, 13.7188844205091, 13.7085726729894, 14.188188727139, 14.4207062289035, 14.0558767590019, 14.0546150322635, 13.8372107056864, 13.1017273131374, 13.7745678963343, 14.223215897062, 14.341261216336, 13.8407839797003, 11.9323668975876, 13.9694536897252, 15.1330345462558, 14.7502465943722), x = c(5.15285150048511, 5.23088470939349, 6.20787065204823, 6.51270372604554, 6.65529156385608, 6.52655333928202, 6.17338302680638, 5.49968725261633, 5.86263253631495, 6.16181924187809, 6.59939463078051, 4.83421448767032, 4.38760463552042, 5.77272300776806, 5.558880334431, 6.53341974503353, 4.87311768540645, 5.77134033889711, 4.37253448319016, 7.29343464773757)), .Names = c("trat", "y", "x"), row.names = c(NA, -20L), class = "data.frame") Att. Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

1 0

Problema para calcular erro padrão da média usando with(tapply())
by ASANTOS 12 Nov '13

12 Nov '13

Caros membros, Tenho variáveis respostas com diferentes tamanhos de vetores e na hora de calcular o erro padrão da média usando o with() e tapply(), em função dos diferentes comprimentos de vetores , não estou conseguindo calcular, alguém poderia me ajudar, segue CRM: y1<-rnorm(100,5)### Variável resposta y2<-rnorm(80,20) y3<-rnorm(100,15) y<-c(y1,y2,y3) trat1<-rep(1,100)### Variável explicativa tratamento trat2<-rep(2,80) trat3<-rep(3,100) trat<-c(trat1,trat2,trat3) id1<-rep(60,100)### Variável explicativa idade id2<-rep(84,180) id<-c(id1,id2) d<-cbind(trat,id,y) ep.b3<-with(d, tapply(y, list(trat,id), sd, na.rm = T)) /sqrt(with(d, tapply(y, list(trat,id), length, na.rm = T))) ### Calculando o erro padrão ep.b3 -- ====================================================================== Alexandre dos Santos Proteção 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 Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

3 3

Esconder avisos ao carregar pacotes
by Andre Oliveira 11 Nov '13

11 Nov '13

Ao carregar o pacote tem estes avisos, já dei uma pesquisada e não achei o que queria. Alguém sabe como carrgar e esconder estes avisos? require(gdata,quietly=TRUE) gdata: read.xls support for 'XLS' (Excel 97-2004) files ENABLED. gdata: read.xls support for 'XLSX' (Excel 2007+) files ENABLED. Attaching package: ‘gdata’ The following object is masked from ‘package:stats’: nobs The following object is masked from ‘package:utils’: object.size obrigado

4 6