Juan Gago
11/11/2016
Analyzing Clustered Data: Why and How to Account for Multiple Observations Nested within a Study Participant?
# Assignment 4 - Analyzing Clustered Data
# packages needed for analysis
library(readr)
library(dplyr)
library(gee)
library(lme4)
library(Matrix)
# reading data
mouse <- readr::read_csv("https://raw.githubusercontent.com/elmoen/clustered_data_code/master/PtenAnalysisData.csv")
# Table 3 - Regression Results
# tabulating treatment groups
count(mouse,fa)
# A tibble: 3 × 2
fa n
<int> <int>
1 0 396
2 1 467
3 2 279
# subsetting data for experimental subset
sset1 <- filter(mouse,fa>0)
count(sset1,fa)
# A tibble: 2 × 2
fa n
<int> <int>
1 1 467
2 2 279
sset2<-filter(sset1, pten==0)
table(sset2$fa)
1 2
190 105
# RESULT 1 Neuron-level linear regression
lnNeuron <- lm(somasize ~ fa, data = sset2)
summary(lnNeuron)
Call:
lm(formula = somasize ~ fa, data = sset2)
Residuals:
Min 1Q Median 3Q Max
-32.008 -11.545 -0.926 9.547 66.311
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 93.712 2.739 34.212 <2e-16 ***
fa 3.151 1.905 1.654 0.0992 .
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 15.66 on 293 degrees of freedom
Multiple R-squared: 0.009251, Adjusted R-squared: 0.00587
F-statistic: 2.736 on 1 and 293 DF, p-value: 0.09919
confint(lnNeuron, level = 0.95)
2.5 % 97.5 %
(Intercept) 88.3207189 99.102633
fa -0.5982463 6.899694
# RESULT 2 Mouse-level regression (mean); no weighting
#aggregate the data by mouse
lnmouse<-summarize(group_by(sset2,mouseid), soma_mouse =mean(somasize), fa= mean(fa))
lnmouse
# A tibble: 9 × 3
mouseid soma_mouse fa
<int> <dbl> <dbl>
1 5 99.15844 1
2 6 95.26709 1
3 7 98.00452 1
4 8 95.14335 1
5 9 103.52473 2
6 10 114.39785 2
7 11 80.81386 2
8 12 95.05462 2
9 13 93.73686 2
#model
lnmouseR<- lm(soma_mouse ~ fa , data = lnmouse)
summary(lnmouseR)
Call:
lm(formula = soma_mouse ~ fa, data = lnmouse)
Residuals:
Min 1Q Median 3Q Max
-16.692 -2.451 -1.626 2.265 16.892
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 96.2811 10.4127 9.247 3.58e-05 ***
fa 0.6122 6.3765 0.096 0.926
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 9.505 on 7 degrees of freedom
Multiple R-squared: 0.001315, Adjusted R-squared: -0.1414
F-statistic: 0.009219 on 1 and 7 DF, p-value: 0.9262
confint(lnmouseR, level = 0.95)
2.5 % 97.5 %
(Intercept) 71.65898 120.90324
fa -14.46568 15.69015
# RESULT 3 Mouse-level regression (mean); analytic weights
lnmouseW <-summarize(group_by(sset2, mouseid), soma_mouse2 = mean(somasize), fa2 = mean(fa), weight = n())
lnmouseWR <- lm(soma_mouse2 ~ fa2 , data = lnmouseW, weights = weight)
summary(lnmouseWR)
Call:
lm(formula = soma_mouse2 ~ fa2, data = lnmouseW, weights = weight)
Weighted Residuals:
Min 1Q Median 3Q Max
-50.796 -28.913 -9.164 13.600 74.745
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 93.712 6.804 13.774 2.51e-06 ***
fa2 3.151 4.731 0.666 0.527
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Residual standard error: 38.91 on 7 degrees of freedom
Multiple R-squared: 0.05958, Adjusted R-squared: -0.07477
F-statistic: 0.4435 on 1 and 7 DF, p-value: 0.5268
confint(lnmouseWR, level = 0.95)
2.5 % 97.5 %
(Intercept) 77.623660 109.79969
fa2 -8.037177 14.33862
# RESULT 4 Marginal regression - GEE
neuronMlm <- gee(somasize ~ fa, id = mouseid, family = gaussian, corstr = 'exchangeable', data = sset2)
(Intercept) fa
93.711676 3.150724
summary(neuronMlm)
GEE: GENERALIZED LINEAR MODELS FOR DEPENDENT DATA
gee S-function, version 4.13 modified 98/01/27 (1998)
Model:
Link: Identity
Variance to Mean Relation: Gaussian
Correlation Structure: Exchangeable
Call:
gee(formula = somasize ~ fa, id = mouseid, data = sset2, family = gaussian,
corstr = "exchangeable")
Summary of Residuals:
Min 1Q Median 3Q Max
-31.198808 -11.431308 -0.957507 9.760842 67.120192
Coefficients:
Estimate Naive S.E. Naive z Robust S.E. Robust z
(Intercept) 94.583206 5.046156 18.7436157 4.815801 19.6401828
fa 2.310301 3.237590 0.7135866 4.580108 0.5044207
Estimated Scale Parameter: 245.6234
Number of Iterations: 3
Working Correlation
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
[7,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[13,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[19,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[25,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[31,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[37,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[43,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[49,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[55,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[,7] [,8] [,9] [,10] [,11] [,12]
[1,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[7,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
[13,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[19,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[25,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[31,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[37,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[43,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[49,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[55,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[,13] [,14] [,15] [,16] [,17] [,18]
[1,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[7,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[13,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
[19,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[25,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[31,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[37,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[43,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[49,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[55,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[,19] [,20] [,21] [,22] [,23] [,24]
[1,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[7,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[13,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[19,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
[25,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[31,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[37,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[43,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[49,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[55,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[,25] [,26] [,27] [,28] [,29] [,30]
[1,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[7,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[13,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[19,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[25,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
[31,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[37,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[43,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[49,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[55,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[,31] [,32] [,33] [,34] [,35] [,36]
[1,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[7,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[13,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[19,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[25,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[31,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
[37,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[43,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[49,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[55,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[,37] [,38] [,39] [,40] [,41] [,42]
[1,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[7,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[13,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[19,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[25,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[31,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[37,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
[43,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[49,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[55,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[,43] [,44] [,45] [,46] [,47] [,48]
[1,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[7,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[13,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[19,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[25,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[31,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[37,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[43,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
[49,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[55,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[,49] [,50] [,51] [,52] [,53] [,54]
[1,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[7,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[13,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[19,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[25,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[31,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[37,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[43,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[49,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
[55,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[,55] [,56] [,57] [,58] [,59] [,60]
[1,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[2,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[3,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[4,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[5,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[6,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[7,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[8,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[9,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[10,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[11,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[12,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[13,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[14,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[15,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[16,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[17,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[18,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[19,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[20,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[21,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[22,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[23,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[24,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[25,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[26,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[27,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[28,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[29,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[30,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[31,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[32,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[33,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[34,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[35,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[36,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[37,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[38,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[39,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[40,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[41,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[42,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[43,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[44,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[45,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[46,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[47,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[48,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[49,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[50,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[51,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[52,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[53,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[54,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[55,] 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651
[56,] 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651 0.06057651
[57,] 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651 0.06057651
[58,] 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651 0.06057651
[59,] 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000 0.06057651
[60,] 0.06057651 0.06057651 0.06057651 0.06057651 0.06057651 1.00000000
# RESULT 6 Mixed-effect regressiona
neuron_lmer <- lmer(somasize ~ fa + (1 | mouseid), data = sset2)
summary(neuron_lmer)
Linear mixed model fit by REML ['lmerMod']
Formula: somasize ~ fa + (1 | mouseid)
Data: sset2
REML criterion at convergence: 2430.4
Scaled residuals:
Min 1Q Median 3Q Max
-2.1070 -0.7522 -0.0974 0.6665 3.6621
Random effects:
Groups Name Variance Std.Dev.
mouseid (Intercept) 65.51 8.094
Residual 214.82 14.657
Number of obs: 295, groups: mouseid, 9
Fixed effects:
Estimate Std. Error t value
(Intercept) 95.535 9.274 10.302
fa 1.359 5.768 0.236
Correlation of Fixed Effects:
(Intr)
fa -0.951