Março 2013 - R-br

Re: [R-br] teste de Bartlett
by Ivan Bezerra Allaman 04 Mar '13

04 Mar '13

Lembrando que todos os testes implementados em softwares que existem são one-way. Em 2007 se não me engano propuseram uma extensão do teste de Levene para análises multi-way. (S,f,P) Allaman

1 0

Gerar mapa de distribuição por UF
by Kaue Veras 04 Mar '13

04 Mar '13

Bom dia Senhores (as), Eu tenho uma base com três colunas e mais de 100000 linhas que contém a UF o CEP e a Quantidade de clientes por cada CEP/UF, eu gostaria de saber se alguém conhece algum pacote ou meio que eu possa inserir esses dados numa distribuição por UF no mapa do Brasil. Desde já agradeço. Atenciosamente, Kauê P. Veras da Cunha Estatística - UERJ Cel: (21) 8254-9601

3 12

Usando o termo Error() com model.matrix para simular dados de Medidas Repetidas
by Augusto Ribas 04 Mar '13

04 Mar '13

Ola, normalmente para tentar entender melhor algumas analises, eu tento gerar exemplo e ver como ficam os resultados e gráficos. Normalmente eu uso o seguinte método: ################################################################### #gerando os tratamentos example<-data.frame(time=factor(rep(c(0,1,2),each=60)), trat=rep(rep(c("controle","dose2","dose3"),each=20),each=60), id=factor(rep(1:60,3))) str(example) #Usando model.matrix para fazer a matriz que multiplica os parâmetros, se quiser simular uma "two way" matrix.1<-model.matrix(~time*trat,data=example) head(matrix.1) #Gerando valores dos parâmetros param.1<-c(1:9) #Gerando os dados de resposta, que é a multiplicação da matriz pelos valores dos parâmetros response.1<-rnorm(nrow(matrix.1),matrix.1%*%param.1,1) #Assim a analise vai estimar os parâmetros usados, e eu posso ver vários parâmetros diferentes e o que acontece. summary(lm(response.1~example$time*example$trat)) #Mas eu gostaria de simular algo com medidas repetidas, tipo usando o termo Error() que da para usar com a função AOV. #Mas eu fiquei meio perdido como proceder, para fazer o mesmo esquema acima mas com medidas repetidas, tendo o error term. matrix.2<-model.matrix(~time*trat+Error(id),data=example) #Olhando como se faz a mesma coisa com a função gls do pacote nlme, suspeito que teria que fazer uma matriz de correlação, mas ai eu não sei como faria aquela multiplicação la em cima. #Alguém saberia como proceder? ##################################################################################### Os comando para fazer a medidas repetidas com nlme que me refiro são desse tutorial: http://www.ats.ucla.edu/stat/r/seminars/Repeated_Measures/repeated_measures… Olhando isso e a função aov eu vejo que tenho que separar o Error, mas não estou sabendo como. library(nlme) longg <- groupedData(pulse ~ as.numeric(exertype) * as.numeric(time) | id, data = exer) fit.cs <- gls(pulse ~ exertype * time, data = longg, corr = corCompSymm(, form= ~ 1 | id) ) summary(fit.cs) Agradeço a atenção Abraços Augusto Ribas -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

1 0

[OT] Um bom livro de estatística
by Rodolfo Timoteo da Silva 03 Mar '13

03 Mar '13

Senhores, sou graduando em Engenharia de Produção e estou interessado em comprar um bom livro impresso de estatística//que tenha os detalhes e "por quês" matemáticos. Se existir alguma publicação com a utilização do R (pt_BR ou en), mais interessante ainda. Qualquer opinião é bem vinda. Att

3 2

Dados com distribuição espaço-temporal no R
by Gilbert Queiroz 03 Mar '13

03 Mar '13

Amigos, Como eu entro com os dados no R para obter uma modelagem de dados espaço-temporais???? Tenho um arquivo onde nas colunas estão as variáveis que quero trabalhar e nas linhas os municípios. Só que eu repito estes municípios para os anos de 2002 a 2011. É assim mesmo???? Abs. Gilbert.

1 0

simulação de dados
by Sérgio Henrique almeida da silva ju 02 Mar '13

02 Mar '13

Olá Como posso simular 2000 amostras aleatórias com 2000 sujeitos cada uma? Quero modelar um desfecho binário (óbitos e não óbitos). A distribuição dos meus dados é: ### Meus Dados id <- rep(1:2000) servidor <- c(rep('Apoio',292),rep('Tecnico',776), rep('Superior',932)) participacao <- c(rep(0,600),rep(1,1400)) # 0 = Não participou do estudo e 1 = Participou do estudo obitos <- ????? Como posso distribuir meus óbitos em relação a variável participacao, sendo uma chande de 2,5 de óbitos de quem não participou em comparação a quem participou? Obrigado -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 / 94429486

2 2

summary de GLMM
by Renata Carvalho 01 Mar '13

01 Mar '13

Bom dia a todos, estou fazendo um glmm com efeito randômico usando a função lmer do pacote lme4, mas quando dou o summary do modelo não consigo verificar se há "overdispersion" nos dados. Por acaso alguém sabe como verifico isso? obrigada Renata

1 0

Quadrado médio de duas características QM(y1+y2)
by Tiago Souza Marçal 01 Mar '13

01 Mar '13

Boa noite pessoal, estou querendo estimar a correlação genética entre duas variáveis (y1 e y2). Entretanto, para isso eu preciso do quadrado médio de a mais b QM (a+b). Alguém sabe como estimar este quadrado médio para as duas variáveis no R. Att. Tiago.

5 11

Re: [R-br] Resolvendo uma distribuição Uniforme SIMPLES.
by Ivan Bezerra Allaman 01 Mar '13

01 Mar '13

Eu não entendi muito bem o que vc quer realmente, mais você pode fazer algo deste tipo: http://nbcgib.uesc.br/lec/prof/ivan/sim/fdpu Para saber como fazer o que está no link abaixo acesse: http://www.rstudio.com/shiny/ (S,f,P) Allaman

2 1