- R-br - DINF UFPR

Ajuda modelo hierárquico de poisson com variância robusta
by Sérgio Henrique almeida da silva ju 26 Jun '12

26 Jun '12

Olá Estou enviando uma dúvida de um amigo: Ele está rodando um modelo hierárquico de poisson com variância robusta structure(list(AFDIC = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), F6rec = structure(c(1L, 2L, 1L, 5L, 4L, 5L, 5L, 5L, 2L, 5L, 1L, 4L, 5L, 5L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 4L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 4L, 2L, 2L, 2L, 3L, 2L, 3L, 2L, 2L, 5L, 1L, 1L, 2L, 3L, 2L, 2L, 1L, 2L, 4L, 3L, 1L, 3L, 1L, 3L, 1L, 2L, 1L, 1L, 2L, 2L, 3L, 1L, 5L, 1L, 4L, 2L, 1L, 2L, 2L, 1L, 2L, 4L, 2L, 1L, 3L, 5L, 2L, 1L, 3L, 2L, 4L, 1L, 2L, 4L, 4L, 2L, 2L, 1L, 3L, 2L, 4L, 2L, 2L), .Label = c("1", "2", "3", "4", "5"), class = "factor")), .Names = c("AFDIC", "F6rec"), row.names = c(NA, -99L), class = "data.frame") # Quando ele roda o modelo o resultado do p-valor é para cada categoria da variável, porém ele gostaria também de obter o valor da variável como um todo: banco$F6rec=as.factor(banco$F6rec) library(epicalc) library(lmtest) library(sandwich) model2=glm(AFDIC~F6rec,family=poisson,data=banco) coeftest(model2, vcov=sandwich) z test of coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.77259 0.68465 -4.0496 5.13e-05 *** F6rec2 1.49165 0.73551 2.0281 0.04255 * F6rec3 1.16315 0.93207 1.2479 0.21206 F6rec4 -15.53000 0.74810 -20.7592 < 2.2e-16 *** F6rec5 0.47000 1.16994 0.4017 0.68788 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #Observe que o resultado é para cada categoria de F6rec (F6rec2, F6rec3, F6rec4, F6rec,5), porém ele gostaria de obter também o p-valor para o F6rec como um todo, quando ele rodar o modelo multivariado ele que decidir se mantem essa variável ou não. Vocês têm alguma ideia? Abraços -- 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) 94429486/78101651 id: 123*20942

4 7

Solução snow+Rmpi
by Junior Beleti 26 Jun '12

26 Jun '12

Boa noite Benilton. Buscando alternativas para a solução do problema apresentado no post: "[R-br] Alteração em função de pacote (novo tópico)", vi que em um dos tópicos de discussão você mencionou que utilizava os pacotes snow + Rmpi. Para a resolução deste meu "problema" esses pacotes poderiam ser utilizados? Att, -- Junior Beleti

2 1

Criar dependência de pacote na criação de pacote
by Emmanuel Arnhold 25 Jun '12

25 Jun '12

Como faço para criar uma dependência de um pacote na criação de um pacote. Vou utilizar funções de outros pacotes e gostaria que o pacote criado e suas dependências (outros pacotes) fossem instaladas e requeridas automaticamente. Como fazer?? Agradeço desde já qualquer dica. Emmanuel

2 1

Relacionamento entre duas variáveis
by Edson Lira 25 Jun '12

25 Jun '12

Gente tenho em um banco uma variável bairro, quero relacioná-la com a variável zona. Cada zona é composta por um número de bairros, exemplo: zona Norte: Cidade Nova, Cidade de Deus, Nova Cidade etc. zona Sul: Petrópolis, São Francisco, Mauazinho, Japiim, etc zona Centro Sul: Aleixo, Vieira Alves, Novo Aleixo zona Leste: Armando Mendes, Jorge Teixeira, etc Como fazer esta relação? []'s Edson Lira Estatístico Manaus-Amazonas

2 1

Integração do R com o SPSS
by Edmar Caldas 25 Jun '12

25 Jun '12

Bom Dia! Preciso de uma ajuda. Estou fazendo um teste no R dentro da sintaxe do SPSS através do programa abaixo, o resultado saí no outuput (saida do spss) mas o resultado também deveria aparecer no editor de dados do SPSS. veja o exemplo: http://people.cst.cmich.edu/lee1c/spss/Integrate_R.htm begin program R. library(MASS) y=spssdata.GetDataFromSPSS() x <-c(1,2,3,4,5,6,7,8,9) x mean(x) end program. segue a resposta do output so SPSS. begin program R. library(MASS) y=spssdata.GetDataFromSPSS() x <-c(1,2,3,4,5,6,7,8,9) x mean(x) end program. [1] 1 2 3 4 5 6 7 8 9 [1] 5 Edmar

3 2

Livrom de R
by Edmar Caldas 25 Jun '12

25 Jun '12

Alguém conhece um livro de programação em R? Bom. Edmar

4 4

Ajuda na análise de seleção recorrente
by Mauricio Couto 25 Jun '12

25 Jun '12

Pessoal. gostaria da ajuda de vocês para analisar um modelo de seleção recorrente no R, o problema são os efeitos aninhados com interação (R/AS). Necessito estimar os componentes de variância para estimar por exemplo a herdabilidade entre outros. Modelo estatístico: Yijkl = µ + Ai + Sj + ASij+R/ASijk +F/Sjl + AS/Sijl + ε ijkl μ é a média Ai é o efeito do ambiente Sj é o efeito de set ASij é o efeito da interação ambiente e set R/ASijk é o efeito da repetição dentro de ambiente e set F/Sjl é o efeito da família dentro de set AF/Sijl é o efeito da interação ambiente e família dentro de set ɛijkl é o erro experimental Obs? Vejam como eu tentei fazer: > couto<-read.table("C:/Users/Mauricio/Desktop/SR.txt", h=T) >couto<-aov(CE~AMB*SET+REP/AMB*SET+Famílias/SET+AMB*Famílias/SET,data=couto) > summary(couto) Obs: só que a ANOVA saiu errada. Desde de já agradeço. Mauricio Farias Couto Doutorando em Genética e Melhoramento de Plantas (UENF-RJ) Mestre em Estatística Aplicada e Biometria (UFV-MG) Engº Agronomo (UESC-BA) Cel: (22) 8117-3725

2 3

Variaveis Dummy
by Julia Maria 25 Jun '12

25 Jun '12

Walmes, segue o CMR de um exemplo que consta no livro Analise de Modelos de Regressao Linear do Charnet. Pag 152. Depois, ja na pagina 243, é refeito com variaveis dummy, usando "0" para mulher e "1"para homem. #####Exemplo Variaveis Dummy ####### #Os dados referem-se a a salarios e tempo de experiencia no cargo #de gerente de agencias bancarias de um grande banco. #Alem disso, gostaria de saber se há diferenças quando #o sexo é levado em conta. y<-c(1.93,3.17,2.27,3.13,2.77,3.09,2.65,2.22,2.85,3.23,2.82,1.90,2.53,2.56,4.22,4.09,3.60,4.70,3.14,2.99,4.74,4.11,2.36,4.09,4.50,2.90,4.48) x<-c(0,17,5,15,9,15,8,5,13,20,11,1,6,7,23,20,18,27,11,10,29,23,4,22,25,9,25) reg1<-lm(y~x) ## Depois fizemos incorporando as dummies, do seguinte modo: z<-c(rep(0,12),rep(1,15)) #até 12 feminino - 15 restante masculino xz<-x*z #incluindo a interação xz reg<-lm(y~x+z+xz) summary(reg) ### Um outro modo: regm<-lm(y[z=="0"]~x[z=="0"]) #modelo para mulheres summary(regm) regh<-lm(y[z=="1"]~x[z=="1"]) #modelo para homens plot(x[z=="1"],regh$residuals) summary(regh) E se eu quiser entao, associar um segundo fator de variaveis dummy, alem de sexo, poderia colocar "pós graduação", 0 se nao possui, 1 se possui.

2 1

Coeficiente de correlação de concordância
by Fernando Antonio de souza 24 Jun '12

24 Jun '12

Caros amigos, Tenho uma dúvida sobre o uso do coeficiente de correlação de concordância (CCC) e gostaria de uma opinião Eu montei um modelo estatístico utilizando meta análise. Fiz algumas análises de ajuste do modelo como avaliação da regressão dos valores observados versus preditos, avaliação dos resíduos. Desejo realizar a análise do CCC e é aí que tenho uma dúvida. Faz sentido utilizá-lo com o mesmo banco de dados utilizado para gerar o modelo? Ou ele só deve ser utilizado para avaliar o modelo com valores preditos com outro banco de dados independente? Gratos pela atenção

1 0

Problema com odfWeaver
by Diogo Braga Mendes 24 Jun '12

24 Jun '12

Obrigado pela atenção Walmes. Estou usando o seguinte código no R: * * *library(odfWeave)* *odfWeave("Fonte.odt","Saida.odt")* onde o arquivo Fonte.odt tem o seguinte código: *<<cemNumeros, echo=FALSE>>=* *cemNumeros<-round(rnorm(10,74,2.34))* *summary(cemNumeros)* *@* Não usei aspas neste código. o log da função aplicada é: Copying Fonte.odt Setting wd to C:\Users\DIOGO\AppData\Local\Temp\Rtmp4qATRq/odfWeave24193047992 Unzipping ODF file using unzip -o "Fonte.odt" Archive: Fonte.odt extracting: mimetype creating: Configurations2/statusbar/ inflating: Configurations2/accelerator/current.xml creating: Configurations2/floater/ creating: Configurations2/popupmenu/ creating: Configurations2/progressbar/ creating: Configurations2/toolpanel/ creating: Configurations2/menubar/ creating: Configurations2/toolbar/ creating: Configurations2/images/Bitmaps/ inflating: content.xml inflating: manifest.rdf inflating: styles.xml extracting: meta.xml extracting: Thumbnails/thumbnail.png inflating: settings.xml inflating: META-INF/manifest.xml Removing Fonte.odt Creating a Pictures directory Pre-processing the contents Sweaving content.Rnw Writing to file content_1.xml Processing code chunks ... 1 : term verbatim(label=cemNumeros) 'content_1.xml' has been Sweaved Removing content.xml Post-processing the contents Removing content.Rnw Removing styles.xml Renaming styles_2.xml to styles.xml Removing manifest.xml Renaming manifest_2.xml to manifest.xml Removing extra files Packaging file using zip -r "Fonte.odt" . adding: Configurations2/ (164 bytes security) (stored 0%) adding: Configurations2/accelerator/ (164 bytes security) (stored 0%) adding: Configurations2/accelerator/current.xml (164 bytes security) (stored 0%) adding: Configurations2/floater/ (164 bytes security) (stored 0%) adding: Configurations2/images/ (164 bytes security) (stored 0%) adding: Configurations2/images/Bitmaps/ (164 bytes security) (stored 0%) adding: Configurations2/menubar/ (164 bytes security) (stored 0%) adding: Configurations2/popupmenu/ (164 bytes security) (stored 0%) adding: Configurations2/progressbar/ (164 bytes security) (stored 0%) adding: Configurations2/statusbar/ (164 bytes security) (stored 0%) adding: Configurations2/toolbar/ (164 bytes security) (stored 0%) adding: Configurations2/toolpanel/ (164 bytes security) (stored 0%) adding: content.xml (164 bytes security) (deflated 97%) adding: manifest.rdf (164 bytes security) (deflated 54%) adding: META-INF/ (164 bytes security) (stored 0%) adding: META-INF/manifest.xml (164 bytes security) (deflated 83%) adding: meta.xml (164 bytes security) (deflated 56%) adding: mimetype (164 bytes security) (stored 0%) adding: Pictures/ (164 bytes security) (stored 0%) adding: settings.xml (164 bytes security) (deflated 84%) adding: styles.xml (164 bytes security) (deflated 92%) adding: Thumbnails/ (164 bytes security) (stored 0%) adding: Thumbnails/thumbnail.png (164 bytes security) (deflated 31%) Copying Fonte.odt Resetting wd Removing C:\Users\DIOGO\AppData\Local\Temp\Rtmp4qATRq/odfWeave24193047992 Fiz as alterações que sugeriste mas infelizmente não funfou. As versões são BrOffice 3.3 e LibreOffice 3.5. -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Diogo Braga Mendes Bacharel em Estatística, CONRE 9558 Fone: (91) 8239-7799 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ “O âmago da ciência, que sempre se apresenta como a mais importante no que refere a resultados práticos, é a pesquisa altamente teórica e abstrata, nascida da infatigável curiosidade, flexibilidade e força da razão humana” Andrei Sakharov In Análise Harmônica de Processos Estocásticos

2 1