[R-br] interpretação dos resultados

alanarocha em sapo.pt alanarocha em sapo.pt
Quarta Novembro 7 13:32:35 BRST 2012


Boa tarde, eu tenho um script cujo o bjectivo é relacionar as duas  
variaveis com a .variavel de control.. Queria fazer um texto que  
esplicasse os resultados mas não estou a conseguir interpretalos.  
Gostava de uma ajuda. Cumprimentos Ana
-------------- Próxima Parte ----------
setwd("G:/SIN/Users/Ana Rua/Ana_trabalho")

#objectivo:
#avaliar a correlação de pthi com DOSE_Cinacalcet, controlando os efeitos da variável GRUPO, que é do tipo fator com dois estados, H e C. 

dados<-read.table("correlacaocinacalcet.csv",header=TRUE,sep=";")
dados
names(dados)

grupo <- dados$Grupo
pthi <- dados$Pthi
dose <- dados$Dose_Cinacalcet

grupo
pthi 
dose 

#A função utilizada para construir modelos lineares de regressão é a função 'lm' que tem os seguintes argumentos principais: 
#            lm( formula, data, weights, subset, na.action )
#'formula' - é uma fórmula estatística que indica o modelo a ser ajustado. Possui a mesma forma básica que foi vista na funções gráficas.
#'data' - o conjunto de dados (data.frame).
#'weights' - são os pesos para regressão ponderada.
#'subset' - um vetor com as condições que definem um sub-conjunto dos dados.
#'na.acation' - função que especifica o que fazer no caso de observações perdidas (NA). O valor default é 'na.omit' que elimina as linhas (observações) que possuem observações perdidas nas variáveis definidas na fórmula.
modelo1 <- lm(pthi~dose)
modelo1

#'summary' apresenta um resumo do modelo linear com:
#1. 
#estatísticas descritivas dos resíduos;
#2. 
#teste t dos coeficientes de regressão;
#3. 
#erro padrão da estimativa;
#4. 
#coeficiente de determinação e coef. de det. ajustado;
#5. 
#teste F geral do modelo.
summary(modelo1)
plot(pthi~dose)
anova(modelo1)


# 'pthi ~ dose') 
#modele Pthi como função estatística de dose
modelo2 <- lm(pthi~grupo)
modelo2
summary(modelo2)
anova(modelo2)
plot(pthi~grupo)


# 'Pthi ~ dose + grupo') 
#a variável grupo entrou alterando o intercepto da regressão entre (pthi) e (dose). 
modelo3 <- lm(pthi~dose+grupo)
modelo3
summary(modelo3)
anova(modelo3)
plot(pthi~dose+grupo)



O modelo ajustado pela fórmula 'pthi ~ dose + grupo' é: 
''pthi_HC = \beta_0 + \beta_1 * dose_HC + \beta_2 *
* HC_{\textrm{H}} + \beta_3 * HC_{\text{C}} + \varepsilon_HC''
onde: 
- pthi_HC é Pthi em casa e na clinica;
dose_HC é a dosagem em casa e na clinica;
HC_{\textrm{grupo}} é a variável indicadora para grupo,
isto é, ela tem valor H se o grupo for H e valor C se o grupo for C;
HC_{\textrm{grupo}} 
é a variável indicadora para grupo C.
A variável grupo tem 2 
níveis (levels) 
levels(dados$grupo)

##resultados
      Doente Grupo    Pthi Dose_Cinacalcet
1   54020637     H  221.6            1080
2   54020698     H  317            2790
3   54020681     H  412.6             930
4   54131227     H  462             900
5   54030338     H  474.9             930
6   54220370     H  547             930
7   54110496     H  547.5             930
8   54110160     H  589             930
9   54140258     H  709             930
10  54110468     H  804.9             930
11  54130380     H  842             930
12  54300217     H  855.4             930
13  54110542     H  911.5             930
14  54300021     H  979.3             930
15  54020732     H 1000             930
16  54300160     H 1054             840
17  54300040     H 1178             900
18  54131239     H 1409             930
19  54140376     H 1579             420
20  54020303     H 2794            1860
21  54030135     H  204             420
22  54140295     H  236            2790
23  54140011     H  253             930
24  54210252     H  334             930
25  54130307     H  358             390
26  54300160     H  360.8             930
27  54210085     H  426             930
28  54030451     H  494.5            1860
29  54020681     H  535.9             900
30  54030344     H  594.6            1860
31  54110542     H  615.7             900
32  54110255     H  659            1860
33  54110468     H  666.1             900
34  54020637     H  705.7             480
35  54300040     H  765.2             930
36  54300217     H  875.8             930
37  54110496     H  953.2             900
38  54131239     H  982             930
39  54300021     H 1275             930
40  54110160     H 1297             900
41  54110070     H 1312            1440
42  54130380     H 1330             930
43  54020732     H 1388             900
44  54020303     H 2060            2790
45  54110386     H 2432             900
46  54210217     H   63.5            1860
47  54110255     H   86.19            1710
48  54210104     H  236            2040
49  54210333     H  249            1860
50  54300040     H  256             390
51  54030135     H  265.9              60
52  54140295     H  322.7            2790
53  54030451     H  354.6            1860
54  54110542     H  436             930
55  54060515     H  458.8             510
56  54110468     H  460.4             930
57  54210252     H  480.1             930
58  54300160     H  480.9             900
59  54210085     H  627.8             930
60  54110496     H  700.5             930
61  54110160     H  700.9             930
62  54020732     H  739.9             870
63  54030344     H  760.7            1860
64  54020681     H  808.6             870
65  54300021     H  830.1             900
66  54210233     H  883             930
67  54140326     H 1033             390
68  54300217     H 1074             900
69  54210181     H 1138           2820
70  54130380     H 1145            900
71  54110070     H 1174          1350
72  54020637     H 1308            870
73  54220660     H 1401           1860
74  54300248     H 1462           1860
75  54020303     H 1536            2700
76  54131031     C   76.5             780
77  54210217     C  159.3             390
78  54210252     C  161.1             780
79  54210181     C  268            2070
80  54300021     C  280.9             390
81  54030451     C  364.1             720
82  54300160     C  574.3             390
83  54110468     C  603.9             390
84  54030135     C  651.8             390
85  54300255     C  727.1             360
86  54110542     C  738.6             390
87  54110255     C  745            390
88  54210233     C  764.5            1170
89  54020681     C  771.2             420
90  54300217     C  772.3             840
91  54140011     C  788.7             540
92  54030344     C  807.5             780
93  54210085     C  813.6             900
94  54110160     C  856.3             780
95  54300268     C  865.1             840
96  54110070     C  882.7             810
97  54140326     C  883.9             540
98  54140295     C  916.3             600
99  54110496     C  923.7             780
100 54210333     C  984.7             780
101 54300068     C 1044             420
102 54210104     C 1084             810
103 54020637     C 1287           2970
104 54020732     C 1332            390
105 54130380     C 1334            600
106 54210180     C 1401            420
107 54300248     C 1452            390
108 54110386     C 2664            780
109 54020303     C 2980           1620
110 54300068     C  176.7             390
111 54210217     C  217             390
112 54210104     C  221.6            1080
113 54110160     C  369.7             720
114 54300160     C  373             420
115 54030451     C  484.4             840
116 54300156     C  512.1            390
117 54030135     C  574.5            390
118 54210252     C  613            780
119 54210181     C  613.3            1080
120 54130380     C  626.3             840
121 54060515     C  706.1             840
122 54140326     C  777.3             720
123 54030344     C  793.2             780
124 54110468     C  819            360
125 54300255     C  822.6             390
126 54110255     C  855.2            390
127 54140011     C  861.5             960
128 54300021     C  948.7             960
129 54110496     C  984.5             660
130 54300268     C 1069            780
131 54140295     C 1130            720
132 54110070     C 1177            780
133 54210180     C 1250            1170
134 54300217     C 1318            780
135 54210333     C 1382             870
136 54300248     C 1404            570
137 54110542     C 1446             360
138 54210085     C 2230            1170
139 54110386     C 2577             720
140 54020303     C 2924            1440
141 54210180     C  116.8            1080
142 54210181     C  229.7             960
143 54210217     C  267.9             270
144 54210252     C  285.8             690
145 54300068     C  311.8             390
146 54300156     C  470.2             420
147 54110255     C  485.9             390
148 54030451     C  526.4             780
149 54110374     C  628.1             720
150 54110070     C  687.1             780
151 54210333     C  702.7            1170
152 54030344     C  722.2             840
153 54140326     C  765.3             840
154 54130380     C  811.8             780
155 54300217     C  845.2             720
156 54030135     C  853.8             420
157 54300255     C  949.3             570
158 54210104     C  978.1             780
159 54140295     C 1076            780
160 54300268     C 1240            780
161 54210085     C 1391           1260
162 54300248     C 1738            780
> 
> names(dados)
[1] "Doente"          "Grupo"           "Pthi"            "Dose_Cinacalcet"

> grupo <- dados$Grupo
> pthi <- dados$Pthi
> dose <- dados$Dose_Cinacalcet
> 
> grupo
  [1] H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H
 [38] H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H H
 [75] H C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C
[112] C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C C
[149] C C C C C C C C C C C C C C
Levels: C H
> pthi 
  [1]  221.60  317.00  412.60  462.00  474.90  547.00  547.50  589.00  709.00
 [10]  804.90  842.00  855.40  911.50  979.30 1000.00 1054.00 1178.00 1409.00
 [19] 1579.00 2794.00  204.00  236.00  253.00  334.00  358.00  360.80  426.00
 [28]  494.50  535.90  594.60  615.70  659.00  666.10  705.70  765.20  875.80
 [37]  953.20  982.00 1275.00 1297.00 1312.00 1330.00 1388.00 2060.00 2432.00
 [46]   63.50   86.19  236.00  249.00  256.00  265.90  322.70  354.60  436.00
 [55]  458.80  460.40  480.10  480.90  627.80  700.50  700.90  739.90  760.70
 [64]  808.60  830.10  883.00 1033.00 1074.00 1138.00 1145.00 1174.00 1308.00
 [73] 1401.00 1462.00 1536.00   76.50  159.30  161.10  268.00  280.90  364.10
 [82]  574.30  603.90  651.80  727.10  738.60  745.00  764.50  771.20  772.30
 [91]  788.70  807.50  813.60  856.30  865.10  882.70  883.90  916.30  923.70
[100]  984.70 1044.00 1084.00 1287.00 1332.00 1334.00 1401.00 1452.00 2664.00
[109] 2980.00  176.70  217.00  221.60  369.70  373.00  484.40  512.10  574.50
[118]  613.00  613.30  626.30  706.10  777.30  793.20  819.00  822.60  855.20
[127]  861.50  948.70  984.50 1069.00 1130.00 1177.00 1250.00 1318.00 1382.00
[136] 1404.00 1446.00 2230.00 2577.00 2924.00  116.80  229.70  267.90  285.80
[145]  311.80  470.20  485.90  526.40  628.10  687.10  702.70  722.20  765.30
[154]  811.80  845.20  853.80  949.30  978.10 1076.00 1240.00 1391.00 1738.00
> dose 
  [1] 1080 2790  930  900  930  930  930  930  930  930  930  930  930  930  930
 [16]  840  900  930  420 1860  420 2790  930  930  390  930  930 1860  900 1860
 [31]  900 1860  900  480  930  930  900  930  930  900 1440  930  900 2790  900
 [46] 1860 1710 2040 1860  390   60 2790 1860  930  510  930  930  900  930  930
 [61]  930  870 1860  870  900  930  390  900 2820  900 1350  870 1860 1860 2700
 [76]  780  390  780 2070  390  720  390  390  390  360  390  390 1170  420  840
 [91]  540  780  900  780  840  810  540  600  780  780  420  810 2970  390  600
[106]  420  390  780 1620  390  390 1080  720  420  840  390  390  780 1080  840
[121]  840  720  780  360  390  390  960  960  660  780  720  780 1170  780  870
[136]  570  360 1170  720 1440 1080  960  270  690  390  420  390  780  720  780
[151] 1170  840  840  780  720  420  570  780  780  780 1260  780


modelo1 <- lm(pthi~dose)
modelo1

Call:
lm(formula = pthi ~ dose)

Coefficients:
(Intercept)         dose  
   738.3511       0.1177  

summary(modelo1)
Call:
lm(formula = pthi ~ dose)

#estatísticas descritivas dos resíduos;
Residuals:
    Min      1Q  Median      3Q     Max 
-893.70 -363.16  -44.26  236.59 2051.04 

#teste t dos coeficientes de regressão;
Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept) 738.35105   84.15885   8.773 2.47e-15 ***
dose          0.11766    0.07688   1.530    0.128    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

#erro padrão da estimativa;
Residual standard error: 545.3 on 160 degrees of freedom

#coeficiente de determinação e coef. de det. ajustado;
Multiple R-squared: 0.01443,    Adjusted R-squared: 0.008268 

#teste F geral do modelo.
F-statistic: 2.342 on 1 and 160 DF,  p-value: 0.1279 

anova(modelo1)
Analysis of Variance Table

Response: pthi
           Df   Sum Sq Mean Sq F value Pr(>F)
dose        1   696390  696390  2.3422 0.1279
Residuals 160 47571107  297319               
plot(pthi~dose)

modelo2 <- lm(pthi~grupo)
modelo2

Call:
lm(formula = pthi ~ grupo)

Coefficients:
(Intercept)       grupoH  
     888.49       -84.84  

summary(modelo2)

Call:
lm(formula = pthi ~ grupo)

Residuals:
    Min      1Q  Median      3Q     Max 
-811.99 -344.45  -94.97  221.10 2091.51 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept)   888.49      58.71  15.134   <2e-16 ***
grupoH        -84.84      86.28  -0.983    0.327    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 547.6 on 160 degrees of freedom
Multiple R-squared: 0.006007,   Adjusted R-squared: -0.0002057 
F-statistic: 0.9669 on 1 and 160 DF,  p-value: 0.3269 

anova(modelo2)
Analysis of Variance Table

Response: pthi
           Df   Sum Sq Mean Sq F value Pr(>F)
grupo       1   289930  289930  0.9669 0.3269
Residuals 160 47977567  299860               

plot(pthi~grupo)

 
modelo3 <- lm(pthi~dose+grupo)
modelo3

Call:
lm(formula = pthi ~ dose + grupo)

Coefficients:
(Intercept)         dose       grupoH  
   760.1821       0.1729    -159.5969  

summary(modelo3)

Call:
lm(formula = pthi ~ dose + grupo)

Residuals:
    Min      1Q  Median      3Q     Max 
-858.70 -322.92  -64.76  219.93 1939.70 

Coefficients:
              Estimate Std. Error t value Pr(>|t|)    
(Intercept)  760.18206   84.59910   8.986 7.13e-16 ***
dose           0.17291    0.08286   2.087   0.0385 *  
grupoH      -159.59687   92.60311  -1.723   0.0868 .  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 541.9 on 159 degrees of freedom
Multiple R-squared: 0.0325,     Adjusted R-squared: 0.02033 
F-statistic: 2.671 on 2 and 159 DF,  p-value: 0.07231 

anova(modelo3)
Analysis of Variance Table

Response: pthi
           Df   Sum Sq Mean Sq F value  Pr(>F)  
dose        1   696390  696390  2.3711 0.12559  
grupo       1   872380  872380  2.9703 0.08675 .
Residuals 159 46698728  293703                  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 



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