Testes com Shortest

Tamanho dos caminhos geodesicos entre A e B

Rede integrada

fint = read.csv("spaths_int.csv")
plot(fint$score ~ fint$dicA, xlab = "Shortest path length", ylab = "S-score")
lmint = lm(fint$score ~ fint$dicA)
abline(lmint, col = "red")

summary(lmint)

## 
## Call:
## lm(formula = fint$score ~ fint$dicA)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -20.467  -2.321  -0.144   1.410  20.638 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -2.706      0.737   -3.67  0.00026 ***
## fint$dicA     -0.265      0.175   -1.52  0.12993    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.53 on 1003 degrees of freedom
## Multiple R-squared:  0.00229,    Adjusted R-squared:  0.00129 
## F-statistic:  2.3 on 1 and 1003 DF,  p-value: 0.13

Rede ppi

fppi = read.csv("spaths_ppi.csv")
plot(fppi$score ~ fppi$dicA, xlab = "Shortest path length", ylab = "S-score")
lmppi = lm(fppi$score ~ fppi$dicA)
abline(lmppi, col = "red")

summary(lmppi)

## 
## Call:
## lm(formula = fppi$score ~ fppi$dicA)
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -19.57  -2.17  -0.13   1.52  17.11 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)    0.392      0.861    0.45     0.65    
## fppi$dicA     -0.843      0.168   -5.02  6.7e-07 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.48 on 597 degrees of freedom
## Multiple R-squared:  0.0406, Adjusted R-squared:  0.039 
## F-statistic: 25.2 on 1 and 597 DF,  p-value: 6.68e-07

Rede reg - inviavel, pois poucos genes no experimento de butland estao na rede regulatoria

Rede met

fmet = read.csv("spaths_met.csv")
plot(fmet$score ~ fmet$dicA, xlab = "Shortest path length", ylab = "S-score")
lmmet = lm(fmet$score ~ fmet$dicA)
abline(lmmet, col = "red")

summary(lmmet)

## 
## Call:
## lm(formula = fmet$score ~ fmet$dicA)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -10.171  -1.986   0.073   1.098  15.022 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   -5.326      0.770   -6.92  6.8e-11 ***
## fmet$dicA      0.293      0.149    1.96    0.051 .  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.72 on 190 degrees of freedom
## Multiple R-squared:  0.0199, Adjusted R-squared:  0.0147 
## F-statistic: 3.85 on 1 and 190 DF,  p-value: 0.0511

NUMERO DE CAMINHOS GEODESICOS

Rede integrada

plot(fint$score ~ fint$dicB, xlab = "Number of shortest paths", ylab = "S-score")
lmint = lm(fint$score ~ fint$dicB)
abline(lmint, col = "red")

summary(lmint)

## 
## Call:
## lm(formula = fint$score ~ fint$dicB)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -20.428  -2.295  -0.118   1.431  20.151 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -3.815547   0.145797  -26.17   <2e-16 ***
## fint$dicB    0.000399   0.000843    0.47     0.64    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.54 on 1003 degrees of freedom
## Multiple R-squared:  0.000223,   Adjusted R-squared:  -0.000774 
## F-statistic: 0.224 on 1 and 1003 DF,  p-value: 0.636

Rede ppi

plot(fppi$score ~ fppi$dicB, xlab = "Number of shortest paths", ylab = "S-score")
lmppi = lm(fppi$score ~ fppi$dicB)
abline(lmppi, col = "red")

summary(lmppi)

## 
## Call:
## lm(formula = fppi$score ~ fppi$dicB)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -20.425  -2.006  -0.037   1.269  16.287 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -3.90181    0.20731  -18.82   <2e-16 ***
## fppi$dicB    0.00587    0.00815    0.72     0.47    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.57 on 597 degrees of freedom
## Multiple R-squared:  0.000866,   Adjusted R-squared:  -0.000807 
## F-statistic: 0.518 on 1 and 597 DF,  p-value: 0.472

Rede met

plot(fmet$score ~ fmet$dicB, xlab = "Number of shortest paths", ylab = "S-score")
lmmet = lm(fmet$score ~ fmet$dicB)
abline(lmmet, col = "red")

summary(lmmet)

## 
## Call:
## lm(formula = fmet$score ~ fmet$dicB)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -10.626  -2.106   0.111   1.322  14.587 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -4.01430    0.28573  -14.05   <2e-16 ***
## fmet$dicB    0.00298    0.00264    1.13     0.26    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.74 on 190 degrees of freedom
## Multiple R-squared:  0.00669,    Adjusted R-squared:  0.00146 
## F-statistic: 1.28 on 1 and 190 DF,  p-value: 0.26