- R-br - DINF UFPR

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

Produção de tabelas no R - LaTeX
by Andre Oliveira 10 Nov '13

10 Nov '13

Pessoas, existe algum pacote melhor e mais eficiente que o xtable para imprimir tabelas no R? André Oliveira Souza

2 1

Versão de um pacote instalado no R
by Andre Oliveira 10 Nov '13

10 Nov '13

Pessoas, como verificar a versão de um dado pacote instalado no R? obrigado André Oliveira Souza

3 2

geoestatística - analise dos resíduos em área irregular
by Hélio Gallo Rocha 09 Nov '13

09 Nov '13

Caros da lista, bom dia. Situação: 1. Área de estudo tem forma irregular, aproximadamente a forma de um feijão com seu eixo principal com azimute de +/- 45º. 2. o grid para fazer a krigagem tem de ter uma forma quadrada ou retangular. 3. posso dizer que tenho duas áreas 3.1 área interna efetiva 3.2 área externa, que, no caso, não interessa. Problema: Ao analisar os residuos da krigagem, percebi que os mesmos (resíduos) da área externa possuem valores superiores. No meu modo de ver eles estariam influenciando na análise, tornando a análise não confiável Pergunto: 1. Teria como isolar estas áreas? Agradeço desde já sugestões -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 0

geoestatística - AIC
by Hélio Gallo Rocha 08 Nov '13

08 Nov '13

Caros listeiros. Como já haviam citado, Prof. Paulo: *"Pensar se a validacao cruzada é de fato nbecessária. Embora seja usada é bem menos informativa e util que outras medidadas de avaliacao do **ajuste baseadas na verossimilhanã e análise dos residuos"* Éder citando Davis, B. M., 1987. Uses and abuses of cross-validation in geostatistics. Math. Geol., v. 17, p. 563–586. *"Cross-validation cannot confirm that a particular model is or is not the optimum. It is a method to better examine and understand the phenomenon under study using the available data"* situação: Caculando o AIC, tive dois resultados diferentes, como segue CRM: require(geoR);data(s100) vModel <- likfit(s100, ini=c(1,0.5), fix.nugget=T) xv.ml<-xvalid(s100,model=vModel) AIC(vModel) #173.1391 reglin_vModel=glm(xv.ml$data~xv.ml$pred) AIC(reglin_vModel) #142.4229 Qual a razão destes valores diferentes? Desde já agradeço qualquer elucidação -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Substituir valores > 0 por 1 em data.frame
by ASANTOS 08 Nov '13

08 Nov '13

Boa tarde pessoal, No seguinte exemplo: da <- expand.grid(trat=gl(2,1), tempo=1:10) X <- model.matrix(~trat+tempo, da); ncol(X) betas <- c(0.1,0.1,0.3) eta <- X%*%betas y1 <- rpois(da$trat, lambda=exp(eta)) y2 <- rbinom(y1, size=1, prob=0.7) da$y <- y1*y2 head(da) Gostaria de criar um novo vetor ybinom, através da substituição de valores em da$y >0 por 1, e os zeros continuam zeros mesmo. Alguém poderia me dar um help? 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 4