Outubro 2011 - R-br

comando weight lme
by Fernando Antonio de souza 11 Out '11

11 Out '11

Caros amigos Estou analisando dados em um modelo misto linear utilizando o pacote nlme contendo uma matrix de correlação do tipo compSymm e necessito utilizar como argumento "weights" o* inverso do quadrado do erro padrão* de cada média. Alguem sabe como faço isto. Sei que o argumento weigue utiliza algus objetos varfunction mas não compreendo muito bem seus objetos e não sei se algum deles fazem o que desejo. Alguém pode me ajudar? consumo<-lme(Fixed=Y~X, data=stpierre,random=X|Estudo) consumo<-update(consumo, random=compSymm(~X), Weight=?) Preciso usar o W2 como peso sendo: w2 =w1/ w sendo w1= o inverso do quadrado médio do erro w= média dos quadrados médios

2 1

Modelo de Efeitos Mistos: implementação
by Rodrigo Fadini 11 Out '11

11 Out '11

Prezados membros do R-br Estou auxiliando um aluno com as análises do seu mestrado, mas não tenho experiência com os Modelos de Efeitos Mistos. O aluno tenta avaliar mudanças na biomassa da floresta em função do tipo de floresta (manejada vs. não manejada) ao longo de 4 anos. Dentro de cada área (manejada e não manejada) ele possui várias parcelas de 1 hectare distribuídas sistematicamente, onde ele faz medidas nas mesmas árvores ao longo dos anos. Em cada parcela, ele divide as árvores em categorias (1-10 cm, 11 a 20 cm, e por aí vai). Eu sugeri a ele utilizar um Modelo Misto, por conta dessa estrutura hierárquica dos dados (parcelas dentro dos tipos de floresta) e também por conta das medidas repetidas às mesmas parcelas. O fato é que estou com problemas na implementação do modelo. Enfim, alguém poderia nos ajudar? Att., Rodrigo

2 2

Como salvar uma função sem precisar carregar
by Marcelo Cardoso mello 11 Out '11

11 Out '11

Como salvar a função abaixo de forma que não seja necessário carregar toda vez que quiser usar? Usar ela como se usa qualquer função do base. Obs.: No Windows. Obrigado por qualquer dica. x<-2 y<-1 f<-function(x,y){x*y+x} f(x,y)

2 1

[OFF-topic] muitos seminários sobre DOE, modelos mistos e relacionados (com vídeo)
by Walmes Zeviani 10 Out '11

10 Out '11

Pessoal, Acidentalmente tive a felicidade de encontrar seminários sobre diversos tópicos associados à delineamento experimental e modelos mistos. Quase todos possuem vídeo e slides. Vale muito a pena conferir e salvar aí na sua barra de favoritos. Aproveitem! http://www.newton.ac.uk/programmes/DAE/seminars/ À disposição. 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

Valores iniciais em hwe.cc - teste de verossimilhança para equilíbrio de HW
by Pedro Emmanuel Alvarenga Americano do Brasil 10 Out '11

10 Out '11

Amigos de R, Estou tentando entender valores iniciais que essa função pede. http://finzi.psych.upenn.edu/R/library/gap/html/hwe.cc.html Infelizmente o autor parece nao ter explicado bem o que significa esses valores iniciais e os exemplos tambem não ajudam muito. Eu tentei modificar progressivamente esses valores (de 0 a 100) e parece que os resultados se modificam mas nao modificaram a interpretação inicial. Mandei um email pro mesmo na semana passada que não retornou. Então estou apelando para amigos da lista que poderiam ja ter utilizado essa função e dar uma luz. Estou com receio de estar incorrendo em erro nao entendendo exatamente como esse valor modifica o resultado final, e quais valores seriam mais apropriados na minha situação. Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 1

Ajuste de regressão não-linear com bloco!
by Ivan Bezerra Allaman 10 Out '11

10 Out '11

Bom noite senhores! Estou analisando uma base de dados que se refere a crescimento de peixes avaliado em duas densidades. Em cada densidade, utilizou-se 2 tanques como bloco, sendo amostrado em torno de 50 peixes por taque. Os animais foram avaliados durante 12 meses, sendo amostrados mensalmente. Andei pesquisando na internet, mais não encontrei nenhum trabalho que tivesse ajustado um modelo não-linear em um delineamento em blocos. Alguém já viu algum trabalho, ou poderia me indicar algum caminho para tal tipo de análise? O ideal é que os parâmetros estimados para cada densidade fossem livre dos efeitos de bloco. Desde já muito obrigado pela atenção! \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

2 1

Re: [R-br] [R_STAT] Teste de esfericidade
by Fernando Henrique Ferraz Pereira da Rosa 10 Out '11

10 Out '11

Michele, Acho que você estava se referindo a função sphericiy.test que eu implementei e fiz referẽncia em um material sobre ANOVA com medidas repetidas. Deixei o código disponível aqui: http://www.fernandohrosa.com.br/en/P/sphericity-test-for-covariance-matrice… Fernando, 2008/2/25 michele barbosa <mixxele(a)yahoo.com.br> > ** > > Como encontrar no help conteúdo para teste de esfericidade? > Algo referente ao comando: > *sphericity.test* > ** > obrigada > Michele > > ------------------------------ > Abra sua conta no Yahoo! Mail<http://br.rd.yahoo.com/mail/taglines/mail/*http://br.mail.yahoo.com/>, > o único sem limite de espaço para armazenamento! > > __._,_.___ > Mensagens neste tópico > <http://br.groups.yahoo.com/group/R_STAT/message/2683;_ylc=X3oDMTM1ZG1uYWw5B…>( > 1) Responder (através da web) > <http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJxdDB1bmdlBF9TAzk3N…>| Adicionar > um novo tópico > <http://br.groups.yahoo.com/group/R_STAT/post;_ylc=X3oDMTJmOXFwdmd2BF9TAzk3N…> > Mensagens<http://br.groups.yahoo.com/group/R_STAT/messages;_ylc=X3oDMTJmNW9qdjR2BF9TA…>| > Arquivos<http://br.groups.yahoo.com/group/R_STAT/files;_ylc=X3oDMTJnaG5zYzJzBF9TAzk3…>| > Fotos<http://br.groups.yahoo.com/group/R_STAT/photos;_ylc=X3oDMTJmbWRzZ3BjBF9TAzk…>| > Links<http://br.groups.yahoo.com/group/R_STAT/links;_ylc=X3oDMTJnamVkZWw1BF9TAzk3…>| Banco > de dados<http://br.groups.yahoo.com/group/R_STAT/database;_ylc=X3oDMTJkamg0ZWVqBF9TA…>| > Enquetes<http://br.groups.yahoo.com/group/R_STAT/polls;_ylc=X3oDMTJncjBzczVxBF9TAzk3…> > [image: Yahoo! Grupos]<http://br.groups.yahoo.com/;_ylc=X3oDMTJlZmFiNHFxBF9TAzk3NDkwNDM1BGdycElkAz…> > Alterar configurações via web<http://br.groups.yahoo.com/group/R_STAT/join;_ylc=X3oDMTJnaGRiMnBvBF9TAzk3N…>(Requer Yahoo! ID) > Alterar configurações via e-mail: Alterar recebimento para lista diária de > mensagens<R_STAT-digest@yahoogrupos.com.br?subject=+Recebimento+de+e-mail:+Lista+de+mensagens>| Alterar > formato para o tradicional<R_STAT-traditional@yahoogrupos.com.br?subject=Alterar+formato+de+distribui%C3%A7%C3%A3o:+Tradicional> > Visite seu Grupo > <http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJldjFsajNkBF9TAzk3NDkwND…>| Termos > de uso do Yahoo! Grupos <http://br.yahoo.com/info/utos.html> | Sair do > grupo <R_STAT-unsubscribe(a)yahoogrupos.com.br?subject=> > Atividade nos últimos dias > > - 1 > Novos usuários<http://br.groups.yahoo.com/group/R_STAT/members;_ylc=X3oDMTJnNzczN3RxBF9TAz…> > > Visite seu Grupo > <http://br.groups.yahoo.com/group/R_STAT;_ylc=X3oDMTJmanAxcGFzBF9TAzk3NDkwND…> > Yahoo! Mail > > Conecte-se ao mundo<http://us.lrd.yahoo.com/_ylc=X3oDMTJtcW92NGJmBF9TAzk3NDkwNDM3BF9wAzEEZ3JwSW…> > > Proteção anti-spam > > Muito mais espaço > Yahoo! Barra > > Instale grátis<http://us.lrd.yahoo.com/_ylc=X3oDMTJwczc4ZjY2BF9TAzk3NDkwNDM3BF9wAzIEZ3JwSW…> > > Buscar sites na web > > Checar seus e-mails . > Yahoo! Grupos > > Crie seu próprio grupo<http://br.groups.yahoo.com/;_ylc=X3oDMTJwcWI3bGZvBF9TAzk3NDkwNDM3BF9wAzMEZ3…> > > A melhor forma de comunicação > > . > > __,_._,___ > -- "Though this be randomness, yet there is structure in't." Fernando H Rosa - Statistician http://www.fernandohrosa.com.br / http://www.feferraz.net - Estatística, Matemática e Computação BankReview.com.br <http://www.bankreview.com.br/> - Escolha melhor seus serviços financeiros! AprendaAlemao.net <http://aprendaalemao.net/> - Seu ponto de partida para melhorar seu Alemão! @fhrosa

2 2

Função não reconhecida
by Mauro Sznelwar 09 Out '11

09 Out '11

Alguém sabe dizer o que está errado por aqui? set.seed(123) b<-5 sigma2<-2 n<-20 x<-runif(n,10,20) y<-b*x+rnorm(n,0,sqrt(sigma2)) plot(x,y) ############################################### #Función para estimar la distribución posterior f1<-function(beta0,nu,a,b,m,x,y) { A<-0 for(i in 1:n) { for(j in 1:n) { A<-A+x[j]*y[i]*(x[j]*y[i]-x[i]*y[j])/(sum(x^2)+nu) } } B0<-(nu*sum((y-beta0*x)^2))/(sum(x)+nu) med<-(sum(x*y)+nu*beta0)/(sum(x^2)+nu) dp<-sqrt((A+B0+b)/((n+a)*(sum(x^2)+nu))) beta<-rt(m,(n+a)/2)*dp+med gamma1<-(n+a)/2; delta1<-(A+B0+b)/2 sigma2<-1/rgamma(m,gamma1,delta1) w<-matrix(c(beta,sigma2),ncol=2) colnames(z)<-c("beta","sigma2") return(w) } ############################################### #Asignamos los valores para las ditribuciones # apriori y el número de simulaciones beta0<-4 nu<-0.00001 a<-1 b<-1 m<-1000 > simula<-f1(beta0,nu,a,b,m,x,y) Erro em colnames(z) <- c("beta", "sigma2") : objeto 'z' não encontrado > simula Erro: objeto 'simula' não encontrado

3 3

Solução
by Daniel Marcelino 08 Out '11

08 Out '11

Pessoal, tô com o seguinte problema: tenho eleições nos DISTRITO x, y, w, z, e vários eleitores para cada distrito desses. Na minha base tenho uma coluna que é o peso de cada voto em relação ao PARTIDO (A, B, C, D), e outras colunas com as opções de preferência de voto de cada eleitor: TURN1, TURN2, TURN3 etc. Eu preciso saber a soma total desse PESO em cada distrito, dividir isso por dois para obter meu cutting point. Depois, eu preciso calcular a soma da variável PESO para cada partido na coluna TURN1, ou seja, a soma do peso de quem votou para o PARTIDO: A, B, C, D, etc. Agora vem a parte difícil; eu preciso verificar se algum dos partidos obteve o valor do cutting point em cada DISTRITO. Se algum deles conseguiu, o distrito já é salvo no objeto final com a indicação de quem venceu. Caso nenhum partido tenha alcançado o valor do cutting point, é preciso eliminar o partido que obteve menos votos e verificar na segunda coluna TURN2 qual foi a segunda opção de voto do eleitor que votou no partido que acabo de eliminar; o o valor na linha da variável PESO do partido que acabo de eliminar é adicionado ao partido de segunda preferência. Agora, verificamos novamente se algum dos partidos restantes no distrito obteve o valor do cutting point, se sim, salvamos ele no objeto final, caso contrário, procedemos da mesma forma anterior até algum partido obter a maioria dos votos. Alguém tem algum ideia de como começar um script para rodar isso.

2 1

Re: [R-br] Análise de MSE, PRESS e SSE
by Thiago De paula protásio 07 Out '11

07 Out '11

Marcelo; Dê uma olhada na função influence.measures().Ela permite explorar mais algumas coisas da análise de regressão. Talvez isso tb te ajude.Att., Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0