Outubro 2013 - R-br

Aggregate
by Fátima Lima Paula 17 Out '13

17 Out '13

Pessoal, estou querendo agregar umas variáveis, mas não sei. Imaginemos o seguinte data frame: sexo=c("1","1","2","2","1","2","1","1","1","1") tcont=c("100","320","24","256","134","290","18","34","15","121") read=c("s","n","n","n","s","s","s","n","n","s") df=data.frame(sexo,idade,read) df Quero agregar essas variáveis de forma que fique a primeira coluna sexo, a segunda coluna a soma do tcont e a terceira coluna a contagem dos read = sim. Alguém pode me ajudar, por favor? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

5 5

auxílio: conjunto de dados e servidor
by Maurício Lordêlo 17 Out '13

17 Out '13

Caro usuário do R, Estou desenvolvendo meu trabalho de doutorado e pretendo obter resultados através de simulações. Consegui fazer todo o procedimento no R através do auxílio dos membros desta lista. No entanto, gostaria de também fazer análises através de dados reais. Como estou com dificuldade em obter estes dados e sei que existem muitos aqui que são de várias áreas diferentes, resolvi solicitar a ajuda nesta lista. O trabalho consiste em um delineamento experimental inteiramente casualizado com resposta binária medida no tempo onde se tem interesse em avaliar o efeito de três ou mais tratamentos. Em último caso, mesmo que a resposta não seja binária, usarei alguma recurso que possa transformá-la. Caso tenha ou conheça alguém que possa ter este conjunto de dados (e lógico, que tenha interesse em disponibilizá-lo) ou até mesmo uma fonte onde possa encontrar ficarei muito agradecido. Aproveitando a oportunidade, como tenho muitos cenários para simular, alguém conhece alguma instituição que disponibilize o servidor para que estas simulações sejam rodadas? Desde já agradeço a quem possa contribuir. Abraço, Maurício Lordêlo

3 3

Re: [R-br] Distribuição de Bernoulli
by Ivan Bezerra Allaman 14 Out '13

14 Out '13

Se foi lhe proposto trabalhar com a função dbern, então se trata de outro software e não o R. Caso exista esta função em outro software, ela é extramente restrita, uma vez que um variável aleatória binomial são "n" realizações de uma variável bernoulli. Outra coisa, sua regressão não é multivariada pelo amor de Deus!!!!!!!! O que caracteriza uma análise uni ou multivariada não são as variáveis preditoras!! Caso ainda queira que alguém lhe ajude, nos informe qual a mensagem de erro ou ainda melhor, nos envie um CMR. Cordialmente, (S,f,P) Allaman

3 3

Correlação policórica
by Felipe Buchbinder 14 Out '13

14 Out '13

Professores, bom dia, Existe alguma função para calcular uma matriz de correlação policórica a partir de uma matriz de respostas a um questionário de escala Likert? A função polychor, do pacote polycor, me retorna apenas um número... Obrigado, Felipe

4 9

Remover NULL de um output da função lapply
by Tiago Souza Marçal 14 Out '13

14 Out '13

Bom dia pessoal, gostaria de saber como remover os NULL do CMR abaixo. x=c("b","a","c") u=lapply(1:3,function(i){ if(x[i]=="a"){ 1} } ) 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 #################################################################

2 1

[OFF-Topic] Dados climáticos de Paris durante maratona de 20 km
by walmes . 12 Out '13

12 Out '13

Saudações à todos, Nos meus cursos para motivar os alunos estou considerando dados reais e maiores que os dados disponíveis nos livros. Dessa vez pedi para que trabalhassem com os dados da maratona de paris cujas planilhas dos últimos anos estão disponíveis em http://www.20kmparis.com/web/resultats_precedents.asp Muito conhecimento se adquire da análise desses dados. É muito gostoso trabalhar esses dados em sala de aula. São infinitas possibilidades de gráficos. Para ver alguns gráficos baixe descompacte o arquivo e rode a aula10.R http://dl.dropboxusercontent.com/u/48140237/ce083.zip No entanto eu gostaria de verificar o efeito do clima sobre os resultados das corridas (no tempo médio, no tempo recorde, nos homens e mulheres, etc). Para isso eu preciso de séries climáticas (ou valores pontuais) de Paris no dia da prova. As datas das provas estão nos panfletos de divulgação da mesma em http://www.20kmparis.com/web/affiches.asp Alguém saberia indicar um site de onde eu pudesse pegar dados climáticos (temperatura e umidade relativa) nessas datas? A idéia e fazer merge de bases de dados, e inclusive testes de hipótese simples envolvendo às variáveis. Muitos maratonistas participam vários anos e tem como ver o desempenho individual. Fiz uma busca no google mas nada que me desse séries históricas. Talvez seja do conhecimento dos que trabalham com séries temporais a origem de dados assim. Quem souber fontes de dados assim, por gentileza me indique. A maratona acontece novamente agora dia 13 de outubro. Quero que os alunos façam uma previsão do tempo recorde por meio da temperatura prevista para Paris ou coisas do tipo. É uma forma de estimulá-los a conhecer os dados. Grato à todos. 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 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

6 8

Script de Bernouilli
by Gisele Correia 12 Out '13

12 Out '13

Olá!! Alguém, por gentileza, poderia me dizer por que o modelo está dando erro? Segue: # Gerando dados tabela<-read.csv("morfo.csv",sep=";",head=T) tabela$SX<-as.factor(tabela$SX) #criar a lista de dados para o jags lddados <- list(y=tabela$SX, N=length(tabela$SX), CT=tabela$CT, CF=tabela$CF, CP=tabela$CP, CC=tabela$CC, FO=tabela$FO, AC=tabela$AC, PD=tabela$PD, CD=tabela$CD, BD=tabela$BD, PP=tabela$PP, CPT=tabela$CPT, BP=tabela$BP, PV=tabela$PV, BV=tabela$BV, LC=tabela$LC, CV=tabela$CV, PA=tabela$PA, BA=tabela$BA, CA=tabela$CA, AP=tabela$AP, LP=tabela$LP, CM=tabela$CM, DOH=tabela$DOH, DOV=tabela$DOV) #criar modelo sink('morfo.txt') cat(' model{ for(i in 1:N){ # loop sobre as observações y[i] ~ dbern(p[i]) # saída binária logit(p[i]) <- beta0 # regressão multivariada + beta1*CT[i] + beta2*CF[i] + beta3*CP[i] + beta4*CC[i] + beta5*FO[i] + beta6*AC[i] + beta7*PD[i] + beta8*CD[i] + beta9*BD[i] + beta10*PP[i] + beta11*CPT[i] + beta12*BP[i] + beta13*PV[i] + beta14*BV[i] + beta15*LC[i] + beta16*CV[i] + beta17*PA[i] + beta18*BA[i] + beta19*CA[i] + beta20*AP[i] + beta21*LP[i] + beta22*CM[i] + beta23*DOH[i] + beta24*DOV[i] } # Prioris beta0 ~ dnorm(0,1.0E-6) beta1 ~ dnorm(0,1.0E-6) beta2 ~ dnorm(0,1.0E-6) beta3 ~ dnorm(0,1.0E-6) beta4 ~ dnorm(0,1.0E-6) beta5 ~ dnorm(0,1.0E-6) beta6 ~ dnorm(0,1.0E-6) beta7 ~ dnorm(0,1.0E-6) beta8 ~ dnorm(0,1.0E-6) beta9 ~ dnorm(0,1.0E-6) beta10 ~ dnorm(0,1.0E-6) beta11 ~ dnorm(0,1.0E-6) beta12 ~ dnorm(0,1.0E-6) beta13 ~ dnorm(0,1.0E-6) beta14 ~ dnorm(0,1.0E-6) beta15 ~ dnorm(0,1.0E-6) beta16 ~ dnorm(0,1.0E-6) beta17 ~ dnorm(0,1.0E-6) beta18 ~ dnorm(0,1.0E-6) beta19 ~ dnorm(0,1.0E-6) beta20 ~ dnorm(0,1.0E-6) beta21 ~ dnorm(0,1.0E-6) beta22 ~ dnorm(0,1.0E-6) beta23 ~ dnorm(0,1.0E-6) beta24 ~ dnorm(0,1.0E-6) } ') sink() #carregar a biblioteca jags library(rjags) #definir os parâmetros para as posteriores params <- c("beta0", "beta1", "beta2", "beta3", "beta4", "beta5", "beta6", "beta7", "beta8", "beta9", "beta10", "beta11", "beta12", "beta13", "beta14", "beta15", "beta16", "beta17", "beta18", "beta19", "beta20", "beta21", "beta22", "beta23", "beta24") #inicializar cadeia inicio <- list(beta0=0, beta1=0, beta2=0, beta3=0, beta4=0, beta5=0, beta6=0, beta7=0, beta8=0, beta9=0, beta10=0, beta11=0, beta12=0, beta13=0, beta14=0, beta15=0, beta16=0, beta17=0, beta18=0, beta19=0, beta20=0, beta21=0, beta22=0, beta23=0, beta24=0) #chamar o modelo no jags morfo.m <- jags.model("morfo.txt",data=lddados,inicio,n.chains=1,n.adapt=1000) Abraço, Gisele

1 0

Distribuição de Bernoulli
by Gisele Correia 12 Out '13

12 Out '13

Olá!! Alguém teria um exemplo de distribuição de Bernoulli no R? Até agora não encontrei nenhum trabalho/problema utilizando a função dbern(p). Pretendo resolver o problema utilizando cadeias de Markov. Desde já agradeço, Gisele

2 3

Função-escalonar
by Maria de Fatima Ferreira Almeida 11 Out '13

11 Out '13

Gostaria de saber se existe alguma função do R que faz escalonamento de matrizes. Alguém poderia me ajudar?

4 6

Fw: [STAT-MATH] Oportunidade Profissional - Consultor Estatístico
by pjobenil＠gmail.com 11 Out '13

11 Out '13

Repassando aos colegas. Enviado do Email do Windows De: nati.0401 Enviado: ‎quinta-feira‎, ‎10‎ de ‎outubro‎ de ‎2013 ‎15‎:‎29 Para: STAT-MATH(a)yahoogrupos.com.br Consultor Estatístico – Foco no software R Será responsável por atuar com pesquisas baseadas em dados estatísticos, mapeamento de concorrência, mensuração de análises de mercado, entre outras atividades. Profissional com experiência e profundos conhecimentos no software estatístico R. Ensino Superior completo. Local: zona Oeste de São Paulo Remuneração: a combinar Interessados, encaminhar CV atualizado para nmarques(a)fitrh.com.br ou então, retornar no tel.: (11) 3149-9128. __._,_.___ | através de email | Responder através da web | Adicionar um novo tópico Mensagens neste tópico (1) Atividade nos últimos dias: Novos usuários 13 Visite seu Grupo Yahoo! Grupos Trocar para: SÃ³ Texto, Resenha DiÃ¡ria • Sair do grupo • Termos de uso . __,_._,___

1 0