[R-br] Fwd: Função estimate para o R

Fernando Antonio de souza nandodesouza em gmail.com
Quinta Setembro 22 19:58:04 BRT 2011


---------- Mensagem encaminhada ----------
De: Fernando Antonio de souza <nandodesouza em gmail.com>
Data: 22 de setembro de 2011 13:33
Assunto: Re: [R-br] Função estimate para o R
Para: r-br em listas.c3sl.ufpr.br


Caro Walmes, Os comando utilizados até agora foram os seguintes



RegFixo<-glm(Y~-1+as.factor(estudo)/X,family=gaussian(),data=Pierre3)#considera
os efietos fixos porém considere o efeito do estudo
library("car")
Anova(RegFixo,type="III",test.statistic="F") #produz a soma de quadrados
tipo III
anova(RegFixo,test="F")#produz a soma de quadrados tipo I
summary(RegFixo)


Acontece que necessito estimar o modelo baseado nos erros padrões estimados
pela função glm e assim obter o intercepto total e a inclinação total do
modelo. E não sei qual função ou como fazer isto.


Eu anexei o artigo para q você olhe caso seja necessário. trata-se de um
passo a passo apresentado pelo autor para a realização de meta-analises.Ele
fornece o banco de dados utilizado bem como apresenta os comando utilizados
no SAS por ele. Eu estou repetindo os passo a passos para compreender o
pensamento por trás deste tipo de análise, porém quero realizá-los no R!


Abaixo segue o banco de dados (também coloquei em anexo em formato csv)



Base de dados

  OBS Experimento X Y  1 1 1.16391 -2.92296  2 1 1.94273 -1.03365  3 1
2.50229 -1.00383  4 1 3.98627 1.33325  5 1 4.36177 1.13752  6 1 4.98954
1.59278  7 2 2.1091 -1.33751  8 2 2.28976 -1.14376  9 2 2.45645 -1.09304  10
2 4.95084 1.4274  11 2 5.92572 2.99572  12 3 2.09605 -0.70073  13 3 2.30059
-1.11007  14 3 3.90388 0.17334  15 3 5.15574 2.38882  16 3 5.57553 1.78167
17 3 5.60523 3.09353  18 4 3.1955 0.18097  19 4 4.32691 0.7431  20 4 5.04092
1.947  21 4 5.17521 1.07029  22 4 5.57254 2.71423  23 4 5.59532 2.43982  24
5 2.31878 0.26487  25 5 4.06027 1.99979  26 6 2.48122 1.31151  27 6 2.73224
1.32218  28 6 3.66628 2.52286  29 6 4.15568 2.83423  30 6 4.57694 2.99086
31 6 4.86356 3.90049  32 7 2.767 1.18358  33 7 2.79859 1.70206  34 7 2.97074
1.84639  35 7 3.85289 3.13155  36 7 6.24981 5.27271  37 7 6.80389 5.85844
38 8 2.74107 -0.26691  39 8 3.72445 1.32567  40 8 4.51721 1.04086  41 8
5.18536 2.40105  42 9 3.24046 3.12334  43 9 3.2864 3.52558  44 9 4.52824
4.29043  45 9 5.10468 4.04917  46 9 5.40578 5.55135  47 9 7.25461 7.2311  48
10 3.09531 3.31464  49 10 3.13951 2.5785  50 10 3.86928 4.0819  51 10
6.40856 6.6731  52 10 7.35976 6.8001  53 10 7.95716 8.2098  54 11 3.24974
2.8825  55 11 3.25369 4.8511  56 11 4.4005 4.3538  57 11 4.74897 6.076  58
11 4.96576 5.9062  59 12 3.00403 2.5733  60 12 3.1542 3.9324  61 12 4.95652
5.8456  62 12 5.35746 5.2066  63 12 5.62651 6.3312  64 12 6.53297 6.6019  65
12 6.6958 6.9498  66 12 7.83844 8.1031  67 13 5.44713 6.4838  68 13 5.82682
6.2341  69 13 6.76986 7.8937  70 13 6.87949 6.9477  71 13 8.47819 8.1336  72
14 4.31493 5.228  73 14 5.70501 6.0941  74 14 6.70113 7.6858  75 15 4.11223
5.2129  76 15 4.55329 5.2169  77 15 4.75641 6.1755  78 15 5.0742 5.7492  79
15 5.90125 6.5407  80 15 6.84886 6.7004  81 15 8.87813 10.2425  82 16
4.25114 5.3359  83 16 6.50762 9.1897  84 16 6.52857 8.4259  85 16 6.67694
8.2782  86 16 6.79539 8.3206  87 16 8.62833 10.4385  88 17 6.47237 7.6132
89 17 6.73826 7.7381  90 17 7.57947 8.5328  91 17 8.11267 8.9045  92 17
9.18169 9.7732  93 18 5.49319 6.6492  94 18 5.62711 6.4464  95 18 6.67182
8.5909  96 18 7.4494 9.2077  97 18 8.13165 9.365  98 18 8.31896 9.5025  99
18 8.74606 10.3598  100 18 8.92038 9.4361  101 18 9.43227 11.0051  102 19
6.06067 7.9306  103 19 7.98874 10.1887  104 19 8.35784 10.9856  105 19
9.52789 11.507  106 19 9.68213 12.6051  107 20 7.15668 9.9928  108 20
9.61256
12.2149


grato pela atenção


Em 22 de setembro de 2011 10:16, Walmes Zeviani
<walmeszeviani em gmail.com>escreveu:

> Fernando,
>
> Use alguma das opções
>
> contrast::contrast()
> gmodels::fit.contrast()
> gmodels::estimable()
> multcomp::glht()
>
> Consulte a documentação das funções. Em caso de dúvidas com uso, forneça
> código mínimo reproduzível.
>
> À 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 em ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110922/348775a2/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: Pierre3.csv
Tipo: text/csv
Tamanho: 2391 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110922/348775a2/attachment.csv>


Mais detalhes sobre a lista de discussão R-br