Fevereiro 2013 - R-br

Step GLM - remover valores ausentes?
by Elisa Seyboth 20 Fev '13

20 Fev '13

Olá! Estou tentando fazer uma seleção de GLMs através da função step(). Porém, ocorre o seguinte erro: Erro em step(M1) : Número de linhas em uso mudou: remover valores ausentes? Acho que ele poderia estar relacionado aos NAs presentes na tabela, mas eles já haviam sido eliminados da tabela a ser utilizada no primeiro modelo (M1) com a função na.omit(). Teriam alguma sugestão? Desde já agradeço, Elisa. Oc. Elisa Seyboth Mestranda em Oceanografia Biológica Laboratório de Tartarugas e Mamíferos Marinhos Universidade Federal do Rio Grande - FURG Tel.: (53) 32336537 / 81093713

2 2

função sumary do R
by alanarocha＠sapo.pt 20 Fev '13

20 Fev '13

Boa tarde, eu estou a trabalhar com as estatisticas descritivas no R. a leitura do ficheiro não deu erros mas a função sumary não faz nada. código: setwd("G:/SIN/Users/Ana Rua/indicadores_ano") dados<-read.table("indicadores_08to12.csv",header=TRUE,sep=";") head(dados) Clinica X2008 X2009 X2010 X2011 X2012 1 PBarca 70% 72% 72% 76% 77% 2 Braga 74% 75% 78% 76% 73% 3 Fafe 74% 75% 74% 72% 76% 4 Maia 70% 75% 71% 76% 80% 5 VNGaia 80% 83% 84% 84% 81% 6 SMFeira 73% 75% 74% 78% 75% summary(dados) Clinica X2008 X2009 X2010 X2011 X2012 Abrantes: 1 74% : 5 75% : 8 74% : 4 78% : 6 81% : 7 Almada : 1 : 4 76% : 5 79% : 4 76% : 5 76% : 4 Alverca : 1 69% : 3 71% : 4 : 3 75% : 4 80% : 4 Amadora : 1 70% : 3 : 3 72% : 3 : 3 75% : 3 APDP : 1 71% : 3 68% : 2 73% : 3 77% : 3 77% : 3 Barreiro: 1 76% : 3 70% : 2 77% : 3 72% : 2 79% : 3 (Other) :31 (Other):16 (Other):13 (Other):17 (Other):14 (Other):13 porque aparece o X nas no nome das variaveis? e não percebo porque o sumary apresenta isto. Cumprimentos Ana Rocha

4 3

Testando frequências de distribuição normal usando o pacote MASS
by ASANTOS 19 Fev '13

19 Fev '13

Bom tarde Pessoal, Estou tentando testar frequências de distribuição normal para um vetor com 100,000 valores com o pacote MASS, o problema é que quando faço o teste de Qui quadrado entre a distribuição empírica o teste não funciona: Erro em chisq.test(x = table(emp), y = table(est)) : 'x' e 'y' devem ter ao menos 2 níveis Quando faço para binomial negativa funciona, mas para distribuição normal não. Alguém teria uma outra solução, segue CRM: require(MASS) res<-rnorm(100000)### Com 100,000 valores não posso usar shapiro.test ( k <- fitdistr(res,"normal")) par <- k$estimate (md <- par[1])#Média (ep <- par[2])#sd (N <- length(res)) est <-N*dnorm(res,mean=md,sd=ep) ## Estimativas ##Teste Chi - Quadrado para valiar a significância da distribuição fecdf <- ecdf(res) ###ecdf- Função de distribuição cumulativa empírica knotsX <- knots(fecdf) emp <- fecdf(c(knotsX,Inf)) chisq.test(x=table(emp),y=table(est))## Teste de Chi quadrado 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 3

Conexão em postgres
by Edson Lira 19 Fev '13

19 Fev '13

Boa Tarde, Gostaria de saber como faço para criar um conexão com o postgres. segue o exemplo que fiz require(RPostgreSQL) drv <- dbDriver("PostgreSQL") con <- dbConnect(drv,user="consulta",password="idoctor@)!)",dbname="iDoctor_DW") mas essa conexão é local, gostaria de acessar o proprio servidor de uma outra máquina, como faço? Edson Lira Estatístico Manaus-Amazonas

1 0

tabela muito grande
by Edson Lira 18 Fev '13

18 Fev '13

Bom dia a todos. Estou analisando um banco que necessito fazer a seguinte rotina: GEN<-table(m$LOCALDENAS,m$MUNICIPIOR,m$GENOTIPAGE,m$REG) ftable(addmargins(GEN, c(1,2,3,4), FUN=list(Total = sum))) Que dá essa saída: D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 CAREIRO A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 CASTANHO A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 COARI A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 CODAJAS A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 EIRUNEPE A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 IPIXUNA A 0 0 0 0 0 0 0 0 0 0 0 0 D 0 0 0 0 0 0 0 0 0 0 0 0 F 0 0 0 0 0 0 0 0 0 0 0 0 Total 0 0 0 0 0 0 0 0 0 0 0 0 ITAITUBA A 0 0 0 0 0 0 0 0 0 0 0 0 A sáida é muito grande e não consigo visualizar na tela de saída do R. Como eliminar os zeros. Edson Lira Estatístico Manaus-Amazonas

3 4

[Dúvida] Paralelizar código R.
by Pedro Rafael 17 Fev '13

17 Fev '13

Pessoal não entendo muito de computação paralela. Quais os pacotes mais utilizados em R com esse fim? Só devo paralelizar o código quando tenho um cluster de computadores conectados ou faz sentido também distribuir parte do processo com processos escravos em meu notebook por exemplo? Alguém poderia me dar uma ideia sobre esse assunto? Quem tiver algum material falando sobre computação paralela em R poderia enviar para meu email? Obrigado desde já. [ ], Pedro Rafael Diniz Marinho.

2 3

Letras gregas: como utilizar menos linhas de comandos....
by andrebvs＠bol.com.br 16 Fev '13

16 Fev '13

1 0

Concurso Substituto - Mat. Aplicada - FEIS/UNESP
by Alan R Panosso 16 Fev '13

16 Fev '13

Caros amigos, Encaminho para a divulgação dessa oportunidade. Abraços e a todos e muito obrigado. Alan De: Berenice Damasceno [mailto:berenice@mat.feis.unesp.br] Enviada em: sexta-feira, 15 de fevereiro de 2013 16:22 Para: Berenice Damasceno Assunto: Confirmando período de inscrição: Concurso Substituto - Mat. Aplicada Boa tarde! Peço sua colaboração na divulgação de concurso para professores emergenciais, na Área de Matemática Aplicada, Computação e Estatística, no Departamento de Matemática da Faculdade de Engenharia de Ilha Solteira (FEIS) - UNESP. Inscrições de 18 a 22/02/2013 na Seção de Comunicações da FEIS 02 (duas) vagas para o conjunto de disciplinas "Bioestatística"; "Técnicas Experimentais com Animais"; "Técnicas e Análises Experimentais"; "Estatística e Probabilidade"; "Probabilidade e Estatística". Informações sobre o concurso poderão ser obtidas acessando: http://www.feis.unesp.br/#!/concursos/departamento-de-matematica/ (a partir de 18/02/13) (o edital de abertura de inscrições será publicado no DOE de SP amanhã, dia 16/02/13) Grata pela colaboração, cordialmente, Berenice. ************************************************************************* Profª. Berenice C. Damasceno Departamento de Matemática - Faculdade de Engenharia de Ilha Solteira - FEIS - UNESP Alameda Rio de Janeiro, 266 15 385 - 000 - Ilha Solteira - SP.-Brasil E_mail: berenice(a)mat.feis.unesp.br ou damascenoberenice(a)yahoo.com.br http://www.mat.feis.unesp.br/berenice_damasceno.php Tel.: +55 (18) 37 43 10 91 Fax : +55 (18) 37 43 11 55 Cel.: +55 (15) 97 15 67 88 *************************************************************************

1 0

Contrastes multivariados
by Tiago Souza Marçal 16 Fev '13

16 Fev '13

Boa noite pessoal, estou utilizando a MANOVA para analisar um experimento fatorial com duas características: f1<-c(rep("A",20),rep("B",20)) f2<-rep(c(gl(4,5)),2) res1<-c(0,0, 20, 0, 20, 20, 0, 0, 20, 40, 30, 40, 40, 20, 40, 30, 20, 30, 30, 50, 40, 40, 40, 70, 60, 40, 40, 40, 60, 20, 30, 30, 50, 40, 40, 50, 30, 40, 30, 30)res2<-c(0,0,1,0,1,1,0,0,1,1.75,1.333333,1,1.5,1,1.75,1.333333,1.5,1.666667,1.333333,1.666667,1,1,1.25,1,1.166667,1.25,1.25,1.25,1.333333,1,2.333333,1,1.2,1,1,2.2,1.333333,1.5,1.333333,1)f1<-factor(f1)f2<-factor(f2) y<-cbind(res1,res2) require(mvoutlier) aq.plot(y, alpha=0.05) # O que devo fazer com os outlines multivariados? mav<-manova(y~f1*f2) summary(mav) # Como testar os contrastes multivariados de interesse no R? Ex: Nível A de f1 vs Nível B de f1 dentro do nível um de f2 # Com fazer a análise Stepdown no R? O R possui outros testes multivariados implementados? Desde já eu agradeço pela ajuda. Att. Tiago.

1 0

parcela subdividida
by Ze Henrique 15 Fev '13

15 Fev '13

Boa noite a todos. Gostaria de saber se alguem saberia como faço para calcular os teste de bartlett, Lilliefors para parcela subdividida. Desde já obrigado a todos. -- *José Henrique Soler Guilhen* Graduando em Agronomia Universidade Federal do Espírito Santo - CCAUFES

2 2