Multiple Interval Mapping and Composite Interval Mapping
P. A. Velasquez Vasconez
12 de julio de 2017
QTLs Composite Interval Mapping and Interval Mapping in Mimulus
We need to export the list of maps using the function write_map let us export the 14 linkage groups previously constructed (https://rpubs.com/PAVelasquezVasconez/280676), and to a file named “maps_list.map”
rm(list = ls());ls()## character(0)setwd("C:\\Users\\alex\\Desktop\\bio")
getwd()## [1] "C:/Users/alex/Desktop/bio"library("qtl")
library("onemap")Multiple Interval Mapping for the phenothype pollen viability (pv)
raw_file <- paste(system.file("extdata", package = "onemap"),
"m_feb06.raw", sep = "\\")
#Now we can read the data using the R/qtl function read.cross
fake_f2_qtl <- read.cross("mm", file = "C:\\Users\\alex\\Desktop\\bio\\m_feb06.raw",
mapfile = "maps_list.map")## --Read the following data:
## Type of cross: f2
## Number of individuals: 287
## Number of markers: 418
## Number of phenotypes: 16
## --Cross type: f2#We can to comparate our map with output of R/qtl Genotype simulation given observed marker data
sim_mimulus <- sim.geno(fake_f2_qtl, n.draws = 20,
step = 2, error.prob = 0.001,
map.function= 'kosambi')Conditional genotype probability in each genome position
qtl_prob <- calc.genoprob(sim_mimulus, step=5,
error.prob=0.001, map.function=c("kosambi"),
stepwidth=c("fixed"))Using the function scantwo to perform a two-dimensional scan for corolla width phenotype.
scan_QTL <- scantwo(qtl_prob, pheno.col = 8, method = 'hk')## Warning in checkcovar(cross, pheno.col, addcovar, intcovar, perm.strata, : Dropping 12 individuals with missing phenotypes.## --Running scanone
## --Running scantwo
## (1,1)
## (1,2)
## (1,3)
## (1,4)
## (1,5)
## (1,6)
## (1,7)
## (1,8)
## (1,9)
## (1,10)
## (1,11)
## (1,12)
## (1,13)
## (1,14)
## (2,2)
## (2,3)
## (2,4)
## (2,5)
## (2,6)
## (2,7)
## (2,8)
## (2,9)
## (2,10)
## (2,11)
## (2,12)
## (2,13)
## (2,14)
## (3,3)
## (3,4)
## (3,5)
## (3,6)
## (3,7)
## (3,8)
## (3,9)
## (3,10)
## (3,11)
## (3,12)
## (3,13)
## (3,14)
## (4,4)
## (4,5)
## (4,6)
## (4,7)
## (4,8)
## (4,9)
## (4,10)
## (4,11)
## (4,12)
## (4,13)
## (4,14)
## (5,5)
## (5,6)
## (5,7)
## (5,8)
## (5,9)
## (5,10)
## (5,11)
## (5,12)
## (5,13)
## (5,14)
## (6,6)
## (6,7)
## (6,8)
## (6,9)
## (6,10)
## (6,11)
## (6,12)
## (6,13)
## (6,14)
## (7,7)
## (7,8)
## (7,9)
## (7,10)
## (7,11)
## (7,12)
## (7,13)
## (7,14)
## (8,8)
## (8,9)
## (8,10)
## (8,11)
## (8,12)
## (8,13)
## (8,14)
## (9,9)
## (9,10)
## (9,11)
## (9,12)
## (9,13)
## (9,14)
## (10,10)
## (10,11)
## (10,12)
## (10,13)
## (10,14)
## (11,11)
## (11,12)
## (11,13)
## (11,14)
## (12,12)
## (12,13)
## (12,14)
## (13,13)
## (13,14)
## (14,14)(summary(scan_QTL))## pos1f pos2f lod.full lod.fv1 lod.int pos1a pos2a lod.add
## c1 :c1 135 175 4.45 2.49 1.660 0 175 2.79
## c1 :c2 170 55 4.71 2.75 1.279 175 25 3.43
## c1 :c3 180 180 4.80 2.84 2.423 175 50 2.37
## c1 :c4 105 25 4.54 2.59 1.100 175 140 3.44
## c1 :c5 160 105 5.78 3.83 2.540 175 40 3.24
## c1 :c6 180 160 8.20 4.15 2.519 175 160 5.68
## c1 :c7 175 30 4.75 2.80 1.514 175 25 3.24
## c1 :c8 135 220 5.70 3.68 2.079 175 90 3.62
## c1 :c9 0 55 4.55 2.60 1.300 175 55 3.25
## c1 :c10 175 25 5.18 3.22 1.450 175 30 3.73
## c1 :c11 175 55 7.07 5.02 2.609 175 110 4.46
## c1 :c12 165 65 4.87 2.92 1.686 175 135 3.19
## c1 :c13 175 60 10.72 2.51 0.543 170 60 10.18
## c1 :c14 155 220 6.35 2.39 0.628 175 210 5.72
## c2 :c2 5 165 5.08 3.27 2.648 0 50 2.43
## c2 :c3 90 180 5.45 3.64 3.205 50 110 2.24
## c2 :c4 65 105 4.87 3.06 1.679 25 140 3.19
## c2 :c5 35 90 6.06 4.25 2.892 25 40 3.16
## c2 :c6 55 160 6.50 2.45 0.746 50 160 5.75
## c2 :c7 45 50 4.96 3.15 1.891 50 25 3.06
## c2 :c8 55 175 4.53 2.51 1.033 25 90 3.50
## c2 :c9 25 65 4.19 2.38 0.805 25 65 3.39
## c2 :c10 50 45 5.83 4.02 2.076 50 25 3.76
## c2 :c11 55 45 6.02 3.97 2.151 50 55 3.87
## c2 :c12 55 75 6.71 4.90 3.874 50 135 2.84
## c2 :c13 30 60 10.73 2.51 1.384 25 60 9.34
## c2 :c14 5 210 7.31 3.35 1.771 25 210 5.54
## c3 :c3 105 130 3.69 3.33 2.352 65 75 1.34
## c3 :c4 105 150 4.04 2.62 2.242 50 140 1.80
## c3 :c5 110 40 4.89 3.55 3.007 50 40 1.88
## c3 :c6 145 135 7.19 3.14 2.897 50 160 4.29
## c3 :c7 50 80 3.69 2.33 1.996 75 25 1.70
## c3 :c8 135 15 6.49 4.47 4.032 110 90 2.46
## c3 :c9 160 65 3.84 2.16 1.615 55 55 2.22
## c3 :c10 10 25 5.77 4.06 3.621 110 25 2.14
## c3 :c11 160 60 6.13 4.08 3.607 165 55 2.52
## c3 :c12 15 170 3.77 2.98 2.430 110 135 1.34
## c3 :c13 110 55 12.39 4.18 3.595 50 60 8.80
## c3 :c14 150 210 8.38 4.43 3.774 110 210 4.61
## c4 :c4 145 230 3.97 2.56 1.489 90 140 2.48
## c4 :c5 20 40 5.76 4.35 3.001 140 40 2.76
## c4 :c6 330 130 7.47 3.42 2.205 175 160 5.27
## c4 :c7 325 30 5.54 4.13 2.703 260 25 2.84
## c4 :c8 50 205 5.16 3.14 1.889 265 90 3.27
## c4 :c9 145 50 4.90 3.22 2.035 140 65 2.87
## c4 :c10 260 25 5.77 4.07 2.587 140 30 3.19
## c4 :c11 295 110 5.42 3.36 1.752 140 55 3.66
## c4 :c12 250 175 5.10 3.69 2.767 140 135 2.34
## c4 :c13 140 60 12.08 3.87 2.167 140 60 9.91
## c4 :c14 175 210 8.99 5.03 3.819 140 210 5.17
## c5 :c5 55 100 4.12 2.79 1.743 55 60 2.38
## c5 :c6 80 160 7.74 3.69 2.070 75 160 5.67
## c5 :c7 70 20 5.58 4.22 2.535 85 25 3.05
## c5 :c8 45 65 5.36 3.34 2.179 40 90 3.18
## c5 :c9 75 60 6.78 5.10 4.033 40 55 2.75
## c5 :c10 115 20 5.64 3.93 2.646 40 25 2.99
## c5 :c11 55 5 4.44 2.39 0.716 40 55 3.73
## c5 :c12 115 60 4.16 2.83 2.072 40 135 2.09
## c5 :c13 40 55 11.40 3.19 1.738 40 60 9.67
## c5 :c14 100 215 6.05 2.09 0.685 95 210 5.36
## c6 :c6 160 220 7.23 3.18 1.891 140 160 5.34
## c6 :c7 165 75 6.76 2.71 1.550 160 35 5.21
## c6 :c8 160 240 7.00 2.96 0.778 160 90 6.23
## c6 :c9 45 65 8.15 4.11 2.152 160 55 6.00
## c6 :c10 160 15 7.63 3.59 1.368 160 30 6.27
## c6 :c11 135 105 7.68 3.63 1.594 160 105 6.08
## c6 :c12 110 200 5.71 1.67 0.778 160 135 4.93
## c6 :c13 160 50 20.42 12.21 5.973 160 60 14.45
## c6 :c14 160 240 10.06 6.01 2.051 160 210 8.01
## c7 :c7 15 20 5.11 3.75 2.390 30 35 2.72
## c7 :c8 35 90 4.96 2.94 1.526 25 90 3.44
## c7 :c9 20 65 3.75 2.07 0.592 25 55 3.16
## c7 :c10 25 20 5.21 3.50 2.109 25 30 3.10
## c7 :c11 35 150 4.61 2.56 1.208 25 55 3.41
## c7 :c12 40 10 5.27 3.91 3.229 25 135 2.04
## c7 :c13 50 55 10.57 2.36 0.502 25 60 10.07
## c7 :c14 25 225 7.07 3.11 1.213 25 210 5.86
## c8 :c8 75 90 5.56 3.53 1.593 90 150 3.96
## c8 :c9 90 65 5.76 3.74 1.973 90 65 3.79
## c8 :c10 250 150 5.11 3.09 1.371 90 30 3.74
## c8 :c11 105 55 6.24 4.19 2.069 90 55 4.17
## c8 :c12 215 175 5.66 3.64 2.999 90 135 2.66
## c8 :c13 90 60 12.55 4.34 1.415 90 60 11.14
## c8 :c14 85 210 7.74 3.78 2.025 90 210 5.71
## c9 :c9 15 55 4.14 2.46 1.451 65 115 2.69
## c9 :c10 145 0 5.89 4.18 2.299 55 30 3.59
## c9 :c11 65 105 5.45 3.39 1.686 55 110 3.76
## c9 :c12 45 135 4.42 2.74 2.032 55 135 2.39
## c9 :c13 55 60 10.61 2.40 0.637 65 60 9.98
## c9 :c14 125 225 7.04 3.08 1.332 65 210 5.71
## c10:c10 20 25 3.71 2.00 1.059 25 160 2.65
## c10:c11 25 115 6.21 4.16 2.039 25 55 4.17
## c10:c12 25 10 5.18 3.47 2.933 30 0 2.25
## c10:c13 0 60 11.30 3.09 1.793 25 60 9.51
## c10:c14 55 240 7.37 3.42 2.236 160 210 5.14
## c11:c11 55 110 5.06 3.01 1.535 55 110 3.53
## c11:c12 110 150 4.70 2.65 1.959 55 135 2.74
## c11:c13 70 60 10.37 2.16 0.492 110 60 9.88
## c11:c14 50 240 7.11 3.15 1.022 110 210 6.09
## c12:c12 0 85 4.94 4.15 2.984 0 20 1.96
## c12:c13 70 50 9.58 1.37 0.834 135 60 8.75
## c12:c14 20 60 6.85 2.89 2.228 0 210 4.63
## c13:c13 55 160 11.51 3.29 1.213 60 160 10.29
## c13:c14 60 215 11.54 3.32 0.245 60 215 11.29
## c14:c14 150 210 6.96 3.00 2.065 195 205 4.89
## lod.av1
## c1 :c1 0.833
## c1 :c2 1.472
## c1 :c3 0.420
## c1 :c4 1.489
## c1 :c5 1.288
## c1 :c6 1.632
## c1 :c7 1.284
## c1 :c8 1.597
## c1 :c9 1.300
## c1 :c10 1.772
## c1 :c11 2.411
## c1 :c12 1.230
## c1 :c13 1.963
## c1 :c14 1.765
## c2 :c2 0.621
## c2 :c3 0.435
## c2 :c4 1.381
## c2 :c5 1.355
## c2 :c6 1.708
## c2 :c7 1.256
## c2 :c8 1.476
## c2 :c9 1.576
## c2 :c10 1.948
## c2 :c11 1.822
## c2 :c12 1.026
## c2 :c13 1.127
## c2 :c14 1.576
## c3 :c3 0.979
## c3 :c4 0.381
## c3 :c5 0.543
## c3 :c6 0.246
## c3 :c7 0.336
## c3 :c8 0.442
## c3 :c9 0.543
## c3 :c10 0.436
## c3 :c11 0.474
## c3 :c12 0.545
## c3 :c13 0.584
## c3 :c14 0.651
## c4 :c4 1.067
## c4 :c5 1.349
## c4 :c6 1.219
## c4 :c7 1.424
## c4 :c8 1.249
## c4 :c9 1.186
## c4 :c10 1.479
## c4 :c11 1.613
## c4 :c12 0.922
## c4 :c13 1.700
## c4 :c14 1.207
## c5 :c5 1.045
## c5 :c6 1.621
## c5 :c7 1.686
## c5 :c8 1.162
## c5 :c9 1.071
## c5 :c10 1.286
## c5 :c11 1.677
## c5 :c12 0.753
## c5 :c13 1.452
## c5 :c14 1.404
## c6 :c6 1.292
## c6 :c7 1.164
## c6 :c8 2.180
## c6 :c9 1.953
## c6 :c10 2.218
## c6 :c11 2.035
## c6 :c12 0.887
## c6 :c13 6.236
## c6 :c14 3.959
## c7 :c7 1.359
## c7 :c8 1.415
## c7 :c9 1.479
## c7 :c10 1.393
## c7 :c11 1.354
## c7 :c12 0.678
## c7 :c13 1.857
## c7 :c14 1.896
## c8 :c8 1.942
## c8 :c9 1.770
## c8 :c10 1.714
## c8 :c11 2.121
## c8 :c12 0.643
## c8 :c13 2.924
## c8 :c14 1.753
## c9 :c9 1.012
## c9 :c10 1.878
## c9 :c11 1.709
## c9 :c12 0.711
## c9 :c13 1.764
## c9 :c14 1.752
## c10:c10 0.941
## c10:c11 2.120
## c10:c12 0.537
## c10:c13 1.297
## c10:c14 1.180
## c11:c11 1.476
## c11:c12 0.692
## c11:c13 1.664
## c11:c14 2.129
## c12:c12 1.167
## c12:c13 0.536
## c12:c14 0.666
## c13:c13 2.081
## c13:c14 3.080
## c14:c14 0.932Using permutations to assess significance level
CIMthresold <- scantwo(qtl_prob, method="hk",
n.perm=100, verbose=T,pheno.col=8)## Warning in checkcovar(cross, pheno.col, addcovar, intcovar, perm.strata, : Dropping 12 individuals with missing phenotypes.## Doing permutation in batch mode ...save(CIMthresold, file = "CIMthresold")
load("CIMthresold")
summary(CIMthresold)## pv (100 permutations)
## full fv1 int add av1 one
## 5% 9.12 6.67 5.92 6.09 3.26 3.53
## 10% 8.68 6.43 5.60 5.85 3.06 3.31CIM<-summary(CIMthresold, alpha = 0.05)Permutation distribution plot: significant threshold in percentile 95 for full, fv1, add, av1, one
plot_hk_CIM<-hist(CIMthresold$full, col=c("blue"))
abline(plot_hk_CIM, v=CIM$full,CIMthresold,
col="red", lty=6)
plot_hk_CIM<-hist(CIMthresold$fv1, col=c("blue"))
abline(plot_hk_CIM, v=CIM$fv1,CIMthresold,
col="red", lty=6)
plot_hk_CIM<-hist(CIMthresold$int, col=c("blue"))
abline(plot_hk_CIM, v=CIM$int,CIMthresold,
col="red", lty=6)
plot_hk_CIM<-hist(CIMthresold$add, col=c("blue"))
abline(plot_hk_CIM, v=CIM$add,CIMthresold,
col="red", lty=6)
plot_hk_CIM<-hist(CIMthresold$av1, col=c("blue"))
abline(plot_hk_CIM, v=CIM$av1,CIMthresold,
col="red", lty=6)
plot_hk_CIM<-hist(CIMthresold$one, col=c("blue"))
abline(plot_hk_CIM, v=CIM$one,CIMthresold,
col="red", lty=6)
####QTL interaction
(qtl_inter <- summary(scan_QTL,perms=CIMthresold,
alpha=c(0.1), pvalues = T))## pos1f pos2f lod.full pval lod.fv1 pval lod.int pval pos1a
## c6 :c13 160 50 20.4 0.00 12.21 0.0 5.973 0.04 160
## c6 :c14 160 240 10.1 0.01 6.01 0.3 2.051 1.00 160
## c13:c14 60 215 11.5 0.00 3.32 1.0 0.245 1.00 60
## pos2a lod.add pval lod.av1 pval
## c6 :c13 60 14.45 0 6.24 0.00
## c6 :c14 210 8.01 0 3.96 0.02
## c13:c14 215 11.29 0 3.08 0.09QTL observation on the linkage map
out_qtl_mim <- makeqtl(qtl_prob, chr=c(6,13),
pos=c(160,50), what="prob")
plot(out_qtl_mim)
#### Fitting the multiple-QTL model
out_fit1 <- fitqtl(qtl_prob, pheno.col=8, qtl=out_qtl_mim,
method="hk", get.ests=TRUE)## Warning in fitqtlengine(pheno = pheno, qtl = qtl, covar = covar, formula = formula, : Dropping 12 individuals with missing phenotypes.summary(out_fit1)##
## fitqtl summary
##
## Method: Haley-Knott regression
## Model: normal phenotype
## Number of observations : 275
##
## Full model result
## ----------------------------------
## Model formula: y ~ Q1 + Q2
##
## df SS MS LOD %var Pvalue(Chi2) Pvalue(F)
## Model 4 3.086815 0.77170374 13.5733 20.33179 8.615331e-13 1.341371e-12
## Error 270 12.095394 0.04479776
## Total 274 15.182209
##
##
## Drop one QTL at a time ANOVA table:
## ----------------------------------
## df Type III SS LOD %var F value Pvalue(Chi2) Pvalue(F)
## 6@160.0 2 1.248 5.865 8.221 13.93 0 1.75e-06 ***
## 13@50.0 2 2.092 9.526 13.779 23.35 0 4.44e-10 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
##
## Estimated effects:
## -----------------
## est SE t
## Intercept 0.38368 0.01315 29.187
## 6@160.0a -0.09661 0.01976 -4.889
## 6@160.0d -0.07303 0.02815 -2.594
## 13@50.0a 0.12744 0.01916 6.652
## 13@50.0d 0.08574 0.03011 2.847Calculate penalties
pen <- calc.penalties(CIMthresold, alpha = 0.1)
summary(CIMthresold)## pv (100 permutations)
## full fv1 int add av1 one
## 5% 9.12 6.67 5.92 6.09 3.26 3.53
## 10% 8.68 6.43 5.60 5.85 3.06 3.31print(pen)## main heavy light
## 3.306125 5.603522 3.123193Perform stepwise to confirm the model
out_stp <- stepwiseqtl(qtl_prob,pheno.col=8,
penalties = pen, max.qtl=4, method="hk",
verbose=T, refine.locations = TRUE)## -Initial scan
## initial lod: 8.273739
## ** new best ** (pLOD increased by 4.9676)
## no.qtl = 1 pLOD = 4.967614 formula: y ~ Q1
## -Step 1
## ---Scanning for additive qtl
## plod = 7.955972
## ---Scanning for QTL interacting with Q1
## plod = 10.00113
## ---Refining positions
## --- Moved a bit
## no.qtl = 2 pLOD = 10.68679 formula: y ~ Q1 + Q2 + Q1:Q2
## ** new best ** (pLOD increased by 5.7192)
## -Step 2
## ---Scanning for additive qtl
## plod = 11.15303
## ---Scanning for QTL interacting with Q1
## plod = 6.434318
## ---Scanning for QTL interacting with Q2
## plod = 6.585002
## ---Look for additional interactions
## ---Refining positions
## no.qtl = 3 pLOD = 11.15303 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3
## ** new best ** (pLOD increased by 0.4662)
## -Step 3
## ---Scanning for additive qtl
## plod = 11.33371
## ---Scanning for QTL interacting with Q1
## plod = 7.392262
## ---Scanning for QTL interacting with Q2
## plod = 7.657056
## ---Scanning for QTL interacting with Q3
## plod = 10.71093
## ---Look for additional interactions
## plod = 6.159536
## ---Refining positions
## no.qtl = 4 pLOD = 11.33371 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3 + Q4
## ** new best ** (pLOD increased by 0.1807)
## -Starting backward deletion
## ---Dropping Q4
## no.qtl = 3 pLOD = 11.15303 formula: y ~ Q1 + Q2 + Q3 + Q1:Q2
## ---Refining positions
## ---Dropping Q3
## no.qtl = 2 pLOD = 10.68679 formula: y ~ Q1 + Q2 + Q1:Q2
## ---Refining positions
## ---Dropping Q1:Q2
## no.qtl = 2 pLOD = 6.961054 formula: y ~ Q1 + Q2
## ---Refining positions
## --- Moved a bit
## ---Dropping Q2
## no.qtl = 1 pLOD = 4.967614 formula: y ~ Q1
## ---Refining positions
## ---One last pass through refineqtlout_stp## QTL object containing genotype probabilities.
##
## name chr pos n.gen
## Q1 6@160.0 6 160.0 3
## Q2 8@90.0 8 90.0 3
## Q3 13@50.0 13 50.0 3
## Q4 13@159.6 13 159.6 3
##
## Formula: y ~ Q1 + Q2 + Q3 + Q4 + Q1:Q3
##
## pLOD: 11.334ref<-refineqtl(qtl_prob,
pheno.col=8,chr=c(6,8,13,13),
pos = c(160,90,50,159))## pos: 160 90 50 159
## Iteration 1
## Q2 pos: 90 -> 88
## LOD increase: 0.149
## Q4 pos: 159 -> 156
## LOD increase: 0.281
## Q1 pos: 160 -> 172
## LOD increase: 0.155
## Q3 pos: 50 -> 54
## LOD increase: 0.041
## all pos: 160 90 50 159 -> 172 88 54 156
## LOD increase at this iteration: 0.627
## Iteration 2
## Q3 pos: 54 -> 54
## LOD increase: 0
## Q2 pos: 88 -> 88
## LOD increase: 0
## Q1 pos: 172 -> 162.9964
## LOD increase: 0.717
## Q4 pos: 156 -> 156
## LOD increase: 0
## all pos: 172 88 54 156 -> 162.9964 88 54 156
## LOD increase at this iteration: 0.717
## Iteration 3
## Q4 pos: 156 -> 156
## LOD increase: 0
## Q2 pos: 88 -> 90
## LOD increase: 0.086
## Q3 pos: 54 -> 54
## LOD increase: 0
## Q1 pos: 162.9964 -> 162
## LOD increase: 0.081
## all pos: 162.9964 88 54 156 -> 162 90 54 156
## LOD increase at this iteration: 0.166
## Iteration 4
## Q1 pos: 162 -> 162
## LOD increase: 0
## Q2 pos: 90 -> 90
## LOD increase: 0
## Q4 pos: 156 -> 156
## LOD increase: 0
## Q3 pos: 54 -> 54
## LOD increase: 0
## all pos: 162 90 54 156 -> 162 90 54 156
## LOD increase at this iteration: 0
## overall pos: 160 90 50 159 -> 162 90 54 156
## LOD increase overall: 1.51Comparison between the Composite Interval Mapping (CIM) and Multiple Interval Mapping (MIM)
fake_f2_qtl <- calc.genoprob(fake_f2_qtl, step = 2)
out.hk <- scanone(fake_f2_qtl, pheno.col =8, method = "hk")## Warning in checkcovar(cross, pheno.col, addcovar, intcovar, perm.strata, : Dropping 12 individuals with missing phenotypes.operm.hk <- scanone(fake_f2_qtl, method="hk", n.perm=100) ## Warning in checkcovar(cross, pheno.col, addcovar, intcovar, perm.strata, : Dropping 11 individuals with missing phenotypes.## Doing permutation in batch mode ...MI1<-summary(operm.hk, alpha = 0.05) [1,1]
(QTL_IM<-summary(out.hk, perms=operm.hk, alpha=0.05,
pvalues=T))## chr pos lod pval
## c6.loc162 6 162 4.13 0.04
## c13.loc58 13 58 8.31 0.00
## c14.loc210 14 210 4.05 0.05qc <- c("6", "13")
qp <- c(162, 58)
fake.f2 <- subset(fake_f2_qtl, chr=qc)
fake.f2 <- sim.geno(fake.f2, n.draws=8, step=2, err=0.001)
(qtl_IntervalMapping <- makeqtl(fake.f2, qc, qp, what="draws"))## QTL object containing imputed genotypes, with 8 imputations.
##
## name chr pos n.gen
## Q1 6@162.0 6 162 3
## Q2 13@58.0 13 58 3plot(qtl_IntervalMapping)
#achar markers with QTLs
######## Composite Interval Mapping (CIM)
qtl_probC <- calc.genoprob(fake_f2_qtl, step=1,
error.prob=0.001, map.function=c("kosambi"),
stepwidth=c("fixed"))
#Covariates selection
cov<-(stepwiseqtl(qtl_probC, pheno.col=8, method=c("hk"),
covar = NULL, model=c("normal"),
incl.markers=TRUE, refine.locations=TRUE,
additive.only=FALSE, scan.pairs=FALSE,
keeplodprofile=TRUE, keeptrace=FALSE,
verbose=TRUE, require.fullrank=FALSE))## -Initial scan
## initial lod: 8.279218
## ** new best ** (pLOD increased by 4.7592)
## no.qtl = 1 pLOD = 4.759218 formula: y ~ Q1
## -Step 1
## ---Scanning for additive qtl
## plod = 7.53468
## ---Scanning for QTL interacting with Q1
## plod = 10.15381
## ---Refining positions
## --- Moved a bit
## no.qtl = 2 pLOD = 10.86002 formula: y ~ Q1 + Q2 + Q1:Q2
## ** new best ** (pLOD increased by 6.1008)
## -Step 2
## ---Scanning for additive qtl
## plod = 10.91642
## ---Scanning for QTL interacting with Q1
## plod = 7.826683
## ---Scanning for QTL interacting with Q2
## plod = 7.992514
## ---Look for additional interactions
## ---Refining positions
## --- Moved a bit
## no.qtl = 3 pLOD = 11.11771 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3
## ** new best ** (pLOD increased by 0.2577)
## -Step 3
## ---Scanning for additive qtl
## plod = 11.06593
## ---Scanning for QTL interacting with Q1
## plod = 8.34772
## ---Scanning for QTL interacting with Q2
## plod = 8.803038
## ---Scanning for QTL interacting with Q3
## plod = 9.714898
## ---Look for additional interactions
## plod = 7.816965
## ---Refining positions
## --- Moved a bit
## no.qtl = 4 pLOD = 11.08155 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3 + Q4
## -Step 4
## ---Scanning for additive qtl
## plod = 9.944729
## ---Scanning for QTL interacting with Q1
## plod = 8.59307
## ---Scanning for QTL interacting with Q2
## plod = 7.511044
## ---Scanning for QTL interacting with Q3
## plod = 9.163309
## ---Scanning for QTL interacting with Q4
## plod = 9.307072
## ---Look for additional interactions
## plod = 8.805918
## ---Refining positions
## --- Moved a bit
## no.qtl = 5 pLOD = 10.34722 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3 + Q4 + Q5
## -Step 5
## ---Scanning for additive qtl
## plod = 9.642601
## ---Scanning for QTL interacting with Q1
## plod = 7.570209
## ---Scanning for QTL interacting with Q2
## plod = 6.837968
## ---Scanning for QTL interacting with Q3
## plod = 9.110951
## ---Scanning for QTL interacting with Q4
## plod = 8.663964
## ---Scanning for QTL interacting with Q5
## plod = 10.07941
## ---Look for additional interactions
## plod = 8.539705
## ---Refining positions
## --- Moved a bit
## no.qtl = 6 pLOD = 10.14991 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3 + Q4 + Q5 + Q6 + Q5:Q6
## -Step 6
## ---Scanning for additive qtl
## plod = 9.489171
## ---Scanning for QTL interacting with Q1
## plod = 8.131762
## ---Scanning for QTL interacting with Q2
## plod = 6.871478
## ---Scanning for QTL interacting with Q3
## plod = 8.185784
## ---Scanning for QTL interacting with Q4
## plod = 8.29245
## ---Scanning for QTL interacting with Q5
## plod = 7.015311
## ---Scanning for QTL interacting with Q6
## plod = 7.367391
## ---Look for additional interactions
## plod = 8.683088
## ---Refining positions
## --- Moved a bit
## no.qtl = 7 pLOD = 9.58532 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3 + Q4 + Q5 + Q6 + Q5:Q6 + Q7
## -Step 7
## ---Scanning for additive qtl
## plod = 8.730081
## ---Scanning for QTL interacting with Q1
## plod = 5.105551
## ---Scanning for QTL interacting with Q2
## plod = 7.231426
## ---Scanning for QTL interacting with Q3
## plod = 8.188776
## ---Scanning for QTL interacting with Q4
## plod = 7.933156
## ---Scanning for QTL interacting with Q5
## plod = 6.469789
## ---Scanning for QTL interacting with Q6
## plod = 6.297228
## ---Scanning for QTL interacting with Q7
## plod = 9.137022
## ---Look for additional interactions
## plod = 8.303687
## ---Refining positions
## --- Moved a bit
## no.qtl = 8 pLOD = 9.633914 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3 + Q4 + Q5 + Q6 + Q5:Q6 + Q7 + Q8 + Q7:Q8
## -Step 8
## ---Scanning for additive qtl
## plod = 8.24634
## ---Scanning for QTL interacting with Q1
## plod = 5.464871
## ---Scanning for QTL interacting with Q2
## plod = 7.05731
## ---Scanning for QTL interacting with Q3
## plod = 8.761807
## ---Scanning for QTL interacting with Q4
## plod = 8.77653
## ---Scanning for QTL interacting with Q5
## plod = 5.689705
## ---Scanning for QTL interacting with Q6
## plod = 6.274078
## ---Scanning for QTL interacting with Q7
## plod = 5.498029
## ---Scanning for QTL interacting with Q8
## plod = 6.007742
## ---Look for additional interactions
## plod = 7.370286
## ---Refining positions
## --- Moved a bit
## no.qtl = 9 pLOD = 9.459231 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3 + Q4 + Q5 + Q6 + Q5:Q6 + Q7 + Q8 + Q7:Q8 + Q9 + Q4:Q9
## -Step 9
## ---Scanning for additive qtl
## plod = 9.526092
## ---Scanning for QTL interacting with Q1
## plod = 6.272767
## ---Scanning for QTL interacting with Q2
## plod = 6.982674
## ---Scanning for QTL interacting with Q3
## plod = 9.160758
## ---Scanning for QTL interacting with Q4
## plod = 6.644181
## ---Scanning for QTL interacting with Q5
## plod = 6.541418
## ---Scanning for QTL interacting with Q6
## plod = 6.193493
## ---Scanning for QTL interacting with Q7
## plod = 7.128738
## ---Scanning for QTL interacting with Q8
## plod = 6.162382
## ---Scanning for QTL interacting with Q9
## plod = 6.281188
## ---Look for additional interactions
## plod = 7.64303
## ---Refining positions
## --- Moved a bit
## no.qtl = 10 pLOD = 10.3338 formula: y ~ Q1 + Q2 + Q1:Q2 + Q3 + Q4 + Q5 + Q6 + Q5:Q6 + Q7 + Q8 + Q7:Q8 + Q9 + Q4:Q9 + Q10
## -Starting backward deletion
## ---Dropping Q9
## no.qtl = 9 pLOD = 9.348867 formula: y ~ Q1 + Q2 + Q3 + Q4 + Q5 + Q6 + Q7 + Q8 + Q9 + Q1:Q2 + Q5:Q6 + Q7:Q8
## ---Refining positions
## --- Moved a bit
## ---Dropping Q9
## no.qtl = 8 pLOD = 9.186686 formula: y ~ Q1 + Q2 + Q3 + Q4 + Q5 + Q6 + Q7 + Q8 + Q1:Q2 + Q5:Q6 + Q7:Q8
## ---Refining positions
## --- Moved a bit
## ---Dropping Q6
## no.qtl = 7 pLOD = 9.62124 formula: y ~ Q1 + Q2 + Q3 + Q4 + Q5 + Q6 + Q7 + Q1:Q2 + Q6:Q7
## ---Refining positions
## --- Moved a bit
## ---Dropping Q5
## no.qtl = 6 pLOD = 10.22447 formula: y ~ Q1 + Q2 + Q3 + Q4 + Q5 + Q6 + Q1:Q2 + Q5:Q6
## ---Refining positions
## --- Moved a bit
## ---Dropping Q5
## no.qtl = 5 pLOD = 10.0146 formula: y ~ Q1 + Q2 + Q3 + Q4 + Q5 + Q1:Q2
## ---Refining positions
## --- Moved a bit
## ** new best ** (pLOD increased by 0.1253)
## ---Dropping Q3
## no.qtl = 4 pLOD = 11.45419 formula: y ~ Q1 + Q2 + Q3 + Q4 + Q1:Q2
## ---Refining positions
## --- Moved a bit
## ** new best ** (pLOD increased by 0.2541)
## ---Dropping Q3
## no.qtl = 3 pLOD = 11.37541 formula: y ~ Q1 + Q2 + Q3 + Q1:Q2
## ---Refining positions
## --- Moved a bit
## ---Dropping Q3
## no.qtl = 2 pLOD = 8.372117 formula: y ~ Q1 + Q2 + Q1:Q2
## ---Refining positions
## --- Moved a bit
## ---Dropping Q1:Q2
## no.qtl = 2 pLOD = 6.365183 formula: y ~ Q1 + Q2
## ---Refining positions
## --- Moved a bit
## ---Dropping Q2
## no.qtl = 1 pLOD = 4.759218 formula: y ~ Q1
## ---Refining positions
## ---One last pass through refineqtl(out.cim.com <- cim(qtl_probC, pheno.col=8, n.marcovar= 4,
window=5, method="hk", error.prob=0.001,
map.function=c("kosambi")))## chr pos lod
## AA341 1 0.000000 5.578655e-01
## c1.loc1 1 1.000000 4.976519e-01
## c1.loc2 1 2.000000 4.351178e-01
## c1.loc3 1 3.000000 3.720416e-01
## c1.loc4 1 4.000000 3.105186e-01
## c1.loc5 1 5.000000 2.527121e-01
## c1.loc6 1 6.000000 2.005532e-01
## c1.loc7 1 7.000000 1.554712e-01
## c1.loc8 1 8.000000 1.182319e-01
## c1.loc9 1 9.000000 8.892071e-02
## c1.loc10 1 10.000000 6.705274e-02
## BA374 1 10.809136 5.423725e-02
## c1.loc11 1 11.000000 5.372529e-02
## c1.loc12 1 12.000000 5.412359e-02
## c1.loc13 1 13.000000 5.944333e-02
## c1.loc14 1 14.000000 6.916310e-02
## c1.loc15 1 15.000000 8.262195e-02
## c1.loc16 1 16.000000 9.903787e-02
## c1.loc17 1 17.000000 1.175395e-01
## c1.loc18 1 18.000000 1.372101e-01
## c1.loc19 1 19.000000 1.571393e-01
## c1.loc20 1 20.000000 1.764764e-01
## c1.loc21 1 21.000000 1.944784e-01
## c1.loc22 1 22.000000 2.105460e-01
## c1.loc23 1 23.000000 2.242456e-01
## c1.loc24 1 24.000000 2.353138e-01
## c1.loc25 1 25.000000 2.436489e-01
## AAT225 1 25.490405 2.467310e-01
## c1.loc26 1 26.000000 2.388703e-01
## c1.loc27 1 27.000000 2.237061e-01
## c1.loc28 1 28.000000 2.091097e-01
## c1.loc29 1 29.000000 1.952252e-01
## c1.loc30 1 30.000000 1.821907e-01
## c1.loc31 1 31.000000 1.701314e-01
## c1.loc32 1 32.000000 1.591533e-01
## c1.loc33 1 33.000000 1.493381e-01
## c1.loc34 1 34.000000 1.407400e-01
## c1.loc35 1 35.000000 1.333862e-01
## c1.loc36 1 36.000000 1.272813e-01
## c1.loc37 1 37.000000 1.224159e-01
## c1.loc38 1 38.000000 1.187787e-01
## c1.loc39 1 39.000000 1.163713e-01
## c1.loc40 1 40.000000 1.152216e-01
## c1.loc41 1 41.000000 1.153956e-01
## c1.loc42 1 42.000000 1.170037e-01
## c1.loc43 1 43.000000 1.202004e-01
## c1.loc44 1 44.000000 1.251771e-01
## c1.loc45 1 45.000000 1.321487e-01
## c1.loc46 1 46.000000 1.413344e-01
## c1.loc47 1 47.000000 1.529366e-01
## c1.loc48 1 48.000000 1.671188e-01
## AA461 1 48.930439 1.827121e-01
## c1.loc49 1 49.000000 1.814259e-01
## c1.loc50 1 50.000000 1.666621e-01
## c1.loc51 1 51.000000 1.605114e-01
## c1.loc52 1 52.000000 1.647562e-01
## c1.loc53 1 53.000000 1.808243e-01
## c1.loc54 1 54.000000 2.095601e-01
## c1.loc55 1 55.000000 2.510402e-01
## c1.loc56 1 56.000000 3.044848e-01
## c1.loc57 1 57.000000 3.682937e-01
## c1.loc58 1 58.000000 4.402060e-01
## c1.loc59 1 59.000000 5.175454e-01
## c1.loc60 1 60.000000 5.974980e-01
## c1.loc61 1 61.000000 6.773659e-01
## c1.loc62 1 62.000000 7.547613e-01
## BA416 1 62.518350 7.930600e-01
## c1.loc63 1 63.000000 8.172631e-01
## c1.loc64 1 64.000000 8.656911e-01
## c1.loc65 1 65.000000 9.112483e-01
## c1.loc66 1 66.000000 9.536059e-01
## c1.loc67 1 67.000000 9.921251e-01
## c1.loc68 1 68.000000 1.025902e+00
## c1.loc69 1 69.000000 1.053864e+00
## c1.loc70 1 70.000000 1.074912e+00
## AAT333 1 70.716835 1.085178e+00
## c1.loc71 1 71.000000 1.081211e+00
## c1.loc72 1 72.000000 1.060777e+00
## c1.loc73 1 73.000000 1.029194e+00
## c1.loc74 1 74.000000 9.852087e-01
## c1.loc75 1 75.000000 9.283991e-01
## c1.loc76 1 76.000000 8.594692e-01
## c1.loc77 1 77.000000 7.803849e-01
## c1.loc78 1 78.000000 6.942604e-01
## c1.loc79 1 79.000000 6.049863e-01
## CA174 1 79.005194 6.045262e-01
## c1.loc80 1 80.000000 4.874200e-01
## c1.loc81 1 81.000000 3.757788e-01
## c1.loc82 1 82.000000 2.746710e-01
## c1.loc83 1 83.000000 1.877111e-01
## c1.loc84 1 84.000000 1.173225e-01
## c1.loc85 1 85.000000 6.458114e-02
## c1.loc86 1 86.000000 2.927204e-02
## c1.loc87 1 87.000000 1.011542e-02
## c1.loc88 1 88.000000 5.088286e-03
## c1.loc89 1 89.000000 1.176563e-02
## c1.loc90 1 90.000000 2.762141e-02
## CA389C 1 90.331790 3.444340e-02
## c1.loc91 1 91.000000 3.357311e-02
## c1.loc92 1 92.000000 3.184563e-02
## c1.loc93 1 93.000000 2.959254e-02
## c1.loc94 1 94.000000 2.685896e-02
## c1.loc95 1 95.000000 2.376899e-02
## c1.loc96 1 96.000000 2.050639e-02
## c1.loc97 1 97.000000 1.727467e-02
## BA129 1 97.866105 1.464400e-02
## c1.loc98 1 98.000000 1.545022e-02
## c1.loc99 1 99.000000 2.246817e-02
## c1.loc100 1 100.000000 3.118949e-02
## c1.loc101 1 101.000000 4.156364e-02
## c1.loc102 1 102.000000 5.338693e-02
## c1.loc103 1 103.000000 6.629452e-02
## c1.loc104 1 104.000000 7.978691e-02
## c1.loc105 1 105.000000 9.329100e-02
## c1.loc106 1 106.000000 1.062415e-01
## c1.loc107 1 107.000000 1.181592e-01
## c1.loc108 1 108.000000 1.287046e-01
## c1.loc109 1 109.000000 1.376964e-01
## AAT267 1 109.527878 1.417808e-01
## c1.loc110 1 110.000000 1.287016e-01
## c1.loc111 1 111.000000 1.101224e-01
## c1.loc112 1 112.000000 1.050263e-01
## c1.loc113 1 113.000000 1.134778e-01
## c1.loc114 1 114.000000 1.342050e-01
## c1.loc115 1 115.000000 1.645478e-01
## c1.loc116 1 116.000000 2.006746e-01
## c1.loc117 1 117.000000 2.381026e-01
## c1.loc118 1 118.000000 2.724380e-01
## c1.loc119 1 119.000000 3.001311e-01
## BA214 1 119.520752 3.110909e-01
## c1.loc120 1 120.000000 3.495544e-01
## c1.loc121 1 121.000000 4.399620e-01
## c1.loc122 1 122.000000 5.417086e-01
## c1.loc123 1 123.000000 6.512208e-01
## c1.loc124 1 124.000000 7.635362e-01
## c1.loc125 1 125.000000 8.728119e-01
## c1.loc126 1 126.000000 9.731302e-01
## c1.loc127 1 127.000000 1.059396e+00
## c1.loc128 1 128.000000 1.128058e+00
## c1.loc129 1 129.000000 1.177449e+00
## c1.loc130 1 130.000000 1.207707e+00
## c1.loc131 1 131.000000 1.220388e+00
## BD270 1 131.469461 1.220958e+00
## c1.loc132 1 132.000000 1.240634e+00
## c1.loc133 1 133.000000 1.267364e+00
## c1.loc134 1 134.000000 1.277789e+00
## c1.loc135 1 135.000000 1.269344e+00
## c1.loc136 1 136.000000 1.240295e+00
## c1.loc137 1 137.000000 1.190127e+00
## c1.loc138 1 138.000000 1.119836e+00
## c1.loc139 1 139.000000 1.032016e+00
## c1.loc140 1 140.000000 9.306938e-01
## c1.loc141 1 141.000000 8.208897e-01
## c1.loc142 1 142.000000 7.080117e-01
## c1.loc143 1 143.000000 5.972196e-01
## c1.loc144 1 144.000000 4.929016e-01
## c1.loc145 1 145.000000 3.983554e-01
## c1.loc146 1 146.000000 3.156901e-01
## c1.loc147 1 147.000000 2.459099e-01
## AA173C 1 147.198235 2.337242e-01
## c1.loc148 1 148.000000 2.584451e-01
## c1.loc149 1 149.000000 2.955888e-01
## c1.loc150 1 150.000000 3.403988e-01
## c1.loc151 1 151.000000 3.934083e-01
## c1.loc152 1 152.000000 4.547048e-01
## c1.loc153 1 153.000000 5.237566e-01
## c1.loc154 1 154.000000 5.993274e-01
## c1.loc155 1 155.000000 6.795240e-01
## MgSTS336 1 155.959657 7.586132e-01
## c1.loc156 1 156.000000 7.609673e-01
## c1.loc157 1 157.000000 8.196402e-01
## c1.loc158 1 158.000000 8.773668e-01
## c1.loc159 1 159.000000 9.334014e-01
## c1.loc160 1 160.000000 9.869835e-01
## c1.loc161 1 161.000000 1.037376e+00
## c1.loc162 1 162.000000 1.083919e+00
## c1.loc163 1 163.000000 1.126075e+00
## c1.loc164 1 164.000000 1.163477e+00
## c1.loc165 1 165.000000 1.195943e+00
## c1.loc166 1 166.000000 1.223478e+00
## CC138 1 166.173331 1.227725e+00
## c1.loc167 1 167.000000 1.233613e+00
## c1.loc168 1 168.000000 1.236811e+00
## c1.loc169 1 169.000000 1.234978e+00
## c1.loc170 1 170.000000 1.227515e+00
## c1.loc171 1 171.000000 1.214087e+00
## c1.loc172 1 172.000000 1.194693e+00
## c1.loc173 1 173.000000 1.169675e+00
## c1.loc174 1 174.000000 1.139674e+00
## c1.loc175 1 175.000000 1.105528e+00
## c1.loc176 1 176.000000 1.068156e+00
## c1.loc177 1 177.000000 1.028449e+00
## c1.loc178 1 178.000000 9.871935e-01
## c1.loc179 1 179.000000 9.450396e-01
## BC194C 1 179.053121 9.428067e-01
## c1.loc180 1 180.000000 9.111478e-01
## c1.loc181 1 181.000000 8.749897e-01
## c1.loc182 1 182.000000 8.360976e-01
## c1.loc183 1 183.000000 7.945839e-01
## c1.loc184 1 184.000000 7.506346e-01
## c1.loc185 1 185.000000 7.045143e-01
## c1.loc186 1 186.000000 6.565686e-01
## c1.loc187 1 187.000000 6.072219e-01
## c1.loc188 1 188.000000 5.569722e-01
## c1.loc189 1 189.000000 5.063796e-01
## c1.loc190 1 190.000000 4.560498e-01
## c1.loc191 1 191.000000 4.066123e-01
## c1.loc192 1 192.000000 3.586939e-01
## c1.loc193 1 193.000000 3.128902e-01
## c1.loc194 1 194.000000 2.697371e-01
## c1.loc195 1 195.000000 2.296872e-01
## c1.loc196 1 196.000000 1.930910e-01
## c1.loc197 1 197.000000 1.601881e-01
## c1.loc198 1 198.000000 1.311061e-01
## c1.loc199 1 199.000000 1.058665e-01
## MgSTS98 1 199.550665 9.363795e-02
## c1.loc200 1 200.000000 8.203669e-02
## c1.loc201 1 201.000000 5.999291e-02
## c1.loc202 1 202.000000 4.404538e-02
## c1.loc203 1 203.000000 3.456979e-02
## c1.loc204 1 204.000000 3.139915e-02
## c1.loc205 1 205.000000 3.395872e-02
## MgSTS212 1 205.242039 3.534195e-02
## CC385C 2 0.000000 1.862388e-01
## c2.loc1 2 1.000000 1.821384e-01
## c2.loc2 2 2.000000 1.763397e-01
## c2.loc3 2 3.000000 1.687932e-01
## c2.loc4 2 4.000000 1.596210e-01
## c2.loc5 2 5.000000 1.491492e-01
## c2.loc6 2 6.000000 1.379095e-01
## c2.loc7 2 7.000000 1.266009e-01
## c2.loc8 2 8.000000 1.160098e-01
## c2.loc9 2 9.000000 1.069070e-01
## c2.loc10 2 10.000000 9.994243e-02
## c2.loc11 2 11.000000 9.556783e-02
## c2.loc12 2 12.000000 9.400003e-02
## c2.loc13 2 13.000000 9.522905e-02
## BA125 2 13.959318 9.885441e-02
## c2.loc14 2 14.000000 9.883393e-02
## c2.loc15 2 15.000000 9.863191e-02
## c2.loc16 2 16.000000 9.912529e-02
## c2.loc17 2 17.000000 1.004714e-01
## c2.loc18 2 18.000000 1.028160e-01
## c2.loc19 2 19.000000 1.062733e-01
## c2.loc20 2 20.000000 1.109011e-01
## c2.loc21 2 21.000000 1.166737e-01
## AA420 2 21.281016 1.184734e-01
## c2.loc22 2 22.000000 1.460644e-01
## c2.loc23 2 23.000000 1.877091e-01
## c2.loc24 2 24.000000 2.303533e-01
## c2.loc25 2 25.000000 2.707170e-01
## c2.loc26 2 26.000000 3.055683e-01
## c2.loc27 2 27.000000 3.321801e-01
## c2.loc28 2 28.000000 3.487183e-01
## c2.loc29 2 29.000000 3.544577e-01
## c2.loc30 2 30.000000 3.497725e-01
## c2.loc31 2 31.000000 3.359313e-01
## CC130 2 31.802670 3.194695e-01
## c2.loc32 2 32.000000 3.192062e-01
## c2.loc33 2 33.000000 3.151007e-01
## CA228 2 33.432057 3.118808e-01
## c2.loc34 2 34.000000 3.218578e-01
## c2.loc35 2 35.000000 3.396737e-01
## c2.loc36 2 36.000000 3.569574e-01
## c2.loc37 2 37.000000 3.727372e-01
## c2.loc38 2 38.000000 3.860465e-01
## c2.loc39 2 39.000000 3.960570e-01
## c2.loc40 2 40.000000 4.022041e-01
## MgSTS316 2 40.506117 4.037538e-01
## c2.loc41 2 41.000000 4.028038e-01
## c2.loc42 2 42.000000 4.007928e-01
## c2.loc43 2 43.000000 3.986760e-01
## c2.loc44 2 44.000000 3.964504e-01
## c2.loc45 2 45.000000 3.941108e-01
## c2.loc46 2 46.000000 3.916493e-01
## c2.loc47 2 47.000000 3.890544e-01
## c2.loc48 2 48.000000 3.863107e-01
## c2.loc49 2 49.000000 3.833990e-01
## c2.loc50 2 50.000000 3.802966e-01
## c2.loc51 2 51.000000 3.769787e-01
## c2.loc52 2 52.000000 3.734200e-01
## c2.loc53 2 53.000000 3.695972e-01
## c2.loc54 2 54.000000 3.654924e-01
## c2.loc55 2 55.000000 3.610957e-01
## c2.loc56 2 56.000000 3.564083e-01
## c2.loc57 2 57.000000 3.514450e-01
## c2.loc58 2 58.000000 3.462348e-01
## c2.loc59 2 59.000000 3.408208e-01
## c2.loc60 2 60.000000 3.352588e-01
## c2.loc61 2 61.000000 3.296138e-01
## c2.loc62 2 62.000000 3.239569e-01
## c2.loc63 2 63.000000 3.183606e-01
## c2.loc64 2 64.000000 3.128952e-01
## c2.loc65 2 65.000000 3.076253e-01
## c2.loc66 2 66.000000 3.026071e-01
## c2.loc67 2 67.000000 2.978868e-01
## c2.loc68 2 68.000000 2.934997e-01
## c2.loc69 2 69.000000 2.894703e-01
## c2.loc70 2 70.000000 2.858129e-01
## c2.loc71 2 71.000000 2.825326e-01
## BA301 2 71.284018 2.816755e-01
## c2.loc72 2 72.000000 2.694183e-01
## c2.loc73 2 73.000000 2.528442e-01
## c2.loc74 2 74.000000 2.371583e-01
## c2.loc75 2 75.000000 2.226555e-01
## c2.loc76 2 76.000000 2.096953e-01
## c2.loc77 2 77.000000 1.986645e-01
## c2.loc78 2 78.000000 1.899317e-01
## c2.loc79 2 79.000000 1.838040e-01
## c2.loc80 2 80.000000 1.804972e-01
## c2.loc81 2 81.000000 1.801242e-01
## BC167 2 81.299111 1.805832e-01
## c2.loc82 2 82.000000 1.849296e-01
## c2.loc83 2 83.000000 1.869079e-01
## c2.loc84 2 84.000000 1.844098e-01
## c2.loc85 2 85.000000 1.784104e-01
## c2.loc86 2 86.000000 1.701461e-01
## c2.loc87 2 87.000000 1.608974e-01
## c2.loc88 2 88.000000 1.517960e-01
## BA172 2 88.988905 1.437852e-01
## c2.loc89 2 89.000000 1.436892e-01
## c2.loc90 2 90.000000 1.372819e-01
## CB280 2 90.815691 1.360256e-01
## c2.loc91 2 91.000000 1.255992e-01
## c2.loc92 2 92.000000 8.483225e-02
## c2.loc93 2 93.000000 7.528523e-02
## c2.loc94 2 94.000000 9.871516e-02
## c2.loc95 2 95.000000 1.522824e-01
## c2.loc96 2 96.000000 2.294577e-01
## MgSTS49 2 96.167923 2.440287e-01
## c2.loc97 2 97.000000 2.194227e-01
## MgSTS106 2 97.065335 2.165176e-01
## c2.loc98 2 98.000000 1.889339e-01
## c2.loc99 2 99.000000 1.558922e-01
## c2.loc100 2 100.000000 1.219952e-01
## c2.loc101 2 101.000000 9.044021e-02
## MgSTS457 2 101.685321 7.162999e-02
## c2.loc102 2 102.000000 7.647191e-02
## c2.loc103 2 103.000000 9.255325e-02
## c2.loc104 2 104.000000 1.084871e-01
## c2.loc105 2 105.000000 1.228161e-01
## c2.loc106 2 106.000000 1.342168e-01
## BD175 2 106.013733 1.343468e-01
## c2.loc107 2 107.000000 1.329289e-01
## c2.loc108 2 108.000000 1.324187e-01
## c2.loc109 2 109.000000 1.329100e-01
## c2.loc110 2 110.000000 1.344418e-01
## MgSTS56 2 110.656502 1.360087e-01
## c2.loc111 2 111.000000 1.375622e-01
## MgSTS589 2 111.428017 1.383996e-01
## c2.loc112 2 112.000000 1.404826e-01
## c2.loc113 2 113.000000 1.451354e-01
## c2.loc114 2 114.000000 1.510252e-01
## c2.loc115 2 115.000000 1.580411e-01
## c2.loc116 2 116.000000 1.659905e-01
## c2.loc117 2 117.000000 1.746186e-01
## c2.loc118 2 118.000000 1.836413e-01
## c2.loc119 2 119.000000 1.927863e-01
## c2.loc120 2 120.000000 2.018282e-01
## CC132 2 120.762952 2.085417e-01
## c2.loc121 2 121.000000 2.117427e-01
## c2.loc122 2 122.000000 2.235959e-01
## c2.loc123 2 123.000000 2.318952e-01
## c2.loc124 2 124.000000 2.359461e-01
## c2.loc125 2 125.000000 2.355113e-01
## c2.loc126 2 126.000000 2.308492e-01
## MgSTS565 2 126.910244 2.234887e-01
## c2.loc127 2 127.000000 2.257784e-01
## c2.loc128 2 128.000000 2.438387e-01
## MgSTS513 2 128.340332 2.460059e-01
## c2.loc129 2 129.000000 2.475400e-01
## c2.loc130 2 130.000000 2.492803e-01
## c2.loc131 2 131.000000 2.499573e-01
## c2.loc132 2 132.000000 2.492465e-01
## c2.loc133 2 133.000000 2.469823e-01
## c2.loc134 2 134.000000 2.432250e-01
## c2.loc135 2 135.000000 2.382751e-01
## c2.loc136 2 136.000000 2.326192e-01
## c2.loc137 2 137.000000 2.268269e-01
## BD99 2 137.889020 2.220058e-01
## c2.loc138 2 138.000000 2.201555e-01
## BC126 2 138.504875 2.119460e-01
## c2.loc139 2 139.000000 2.129164e-01
## c2.loc140 2 140.000000 2.153898e-01
## c2.loc141 2 141.000000 2.186050e-01
## c2.loc142 2 142.000000 2.226413e-01
## c2.loc143 2 143.000000 2.275731e-01
## c2.loc144 2 144.000000 2.334654e-01
## c2.loc145 2 145.000000 2.403682e-01
## c2.loc146 2 146.000000 2.483106e-01
## c2.loc147 2 147.000000 2.572939e-01
## c2.loc148 2 148.000000 2.672865e-01
## c2.loc149 2 149.000000 2.782183e-01
## c2.loc150 2 150.000000 2.899785e-01
## c2.loc151 2 151.000000 3.024159e-01
## c2.loc152 2 152.000000 3.153426e-01
## c2.loc153 2 153.000000 3.285417e-01
## c2.loc154 2 154.000000 3.417781e-01
## c2.loc155 2 155.000000 3.548117e-01
## c2.loc156 2 156.000000 3.674117e-01
## c2.loc157 2 157.000000 3.793691e-01
## c2.loc158 2 158.000000 3.905080e-01
## c2.loc159 2 159.000000 4.006923e-01
## c2.loc160 2 160.000000 4.098294e-01
## c2.loc161 2 161.000000 4.178692e-01
## c2.loc162 2 162.000000 4.248001e-01
## c2.loc163 2 163.000000 4.306435e-01
## c2.loc164 2 164.000000 4.354454e-01
## c2.loc165 2 165.000000 4.392697e-01
## c2.loc166 2 166.000000 4.421903e-01
## c2.loc167 2 167.000000 4.442846e-01
## CA279 2 167.720113 4.453213e-01
## CB187 3 0.000000 6.004143e-01
## c3.loc1 3 1.000000 5.241653e-01
## c3.loc2 3 2.000000 4.576377e-01
## c3.loc3 3 3.000000 4.012144e-01
## c3.loc4 3 4.000000 3.552775e-01
## c3.loc5 3 5.000000 3.202067e-01
## c3.loc6 3 6.000000 2.963637e-01
## c3.loc7 3 7.000000 2.840655e-01
## c3.loc8 3 8.000000 2.835432e-01
## c3.loc9 3 9.000000 2.948860e-01
## BC70 3 9.504368 3.050234e-01
## c3.loc10 3 10.000000 2.674240e-01
## c3.loc11 3 11.000000 2.001354e-01
## c3.loc12 3 12.000000 1.487473e-01
## c3.loc13 3 13.000000 1.169548e-01
## c3.loc14 3 14.000000 1.066076e-01
## c3.loc15 3 15.000000 1.172811e-01
## c3.loc16 3 16.000000 1.463872e-01
## c3.loc17 3 17.000000 1.897869e-01
## c3.loc18 3 18.000000 2.426624e-01
## BC586C 3 18.170200 2.521813e-01
## c3.loc19 3 19.000000 2.650233e-01
## c3.loc20 3 20.000000 2.854910e-01
## c3.loc21 3 21.000000 3.106284e-01
## c3.loc22 3 22.000000 3.390943e-01
## c3.loc23 3 23.000000 3.691700e-01
## c3.loc24 3 24.000000 3.990426e-01
## c3.loc25 3 25.000000 4.271010e-01
## c3.loc26 3 26.000000 4.521555e-01
## BD292 3 26.126935 4.550560e-01
## c3.loc27 3 27.000000 4.227789e-01
## c3.loc28 3 28.000000 3.844362e-01
## c3.loc29 3 29.000000 3.467835e-01
## c3.loc30 3 30.000000 3.123532e-01
## c3.loc31 3 31.000000 2.835466e-01
## c3.loc32 3 32.000000 2.621082e-01
## c3.loc33 3 33.000000 2.487792e-01
## AA404 3 33.407539 2.456504e-01
## c3.loc34 3 34.000000 2.454769e-01
## c3.loc35 3 35.000000 2.448755e-01
## c3.loc36 3 36.000000 2.438733e-01
## BD429 3 36.053128 2.438094e-01
## c3.loc37 3 37.000000 2.496071e-01
## c3.loc38 3 38.000000 2.605435e-01
## c3.loc39 3 39.000000 2.754135e-01
## c3.loc40 3 40.000000 2.925253e-01
## CA131 3 40.232537 2.965941e-01
## c3.loc41 3 41.000000 2.824210e-01
## c3.loc42 3 42.000000 2.652522e-01
## c3.loc43 3 43.000000 2.502880e-01
## c3.loc44 3 44.000000 2.381102e-01
## c3.loc45 3 45.000000 2.290343e-01
## c3.loc46 3 46.000000 2.230631e-01
## c3.loc47 3 47.000000 2.199067e-01
## AAT312 3 47.746051 2.190905e-01
## c3.loc48 3 48.000000 2.202195e-01
## c3.loc49 3 49.000000 2.251031e-01
## c3.loc50 3 50.000000 2.305502e-01
## c3.loc51 3 51.000000 2.363835e-01
## c3.loc52 3 52.000000 2.423539e-01
## c3.loc53 3 53.000000 2.481627e-01
## c3.loc54 3 54.000000 2.534947e-01
## c3.loc55 3 55.000000 2.580517e-01
## c3.loc56 3 56.000000 2.615799e-01
## BA75 3 56.907760 2.637271e-01
## c3.loc57 3 57.000000 2.629506e-01
## c3.loc58 3 58.000000 2.536473e-01
## c3.loc59 3 59.000000 2.430145e-01
## c3.loc60 3 60.000000 2.313456e-01
## c3.loc61 3 61.000000 2.190163e-01
## c3.loc62 3 62.000000 2.064535e-01
## c3.loc63 3 63.000000 1.940902e-01
## c3.loc64 3 64.000000 1.823155e-01
## c3.loc65 3 65.000000 1.714319e-01
## c3.loc66 3 66.000000 1.616313e-01
## c3.loc67 3 67.000000 1.529924e-01
## c3.loc68 3 68.000000 1.454969e-01
## c3.loc69 3 69.000000 1.390564e-01
## CB156 3 69.304823 1.372959e-01
## c3.loc70 3 70.000000 1.342143e-01
## c3.loc71 3 71.000000 1.317325e-01
## c3.loc72 3 72.000000 1.311549e-01
## c3.loc73 3 73.000000 1.319010e-01
## c3.loc74 3 74.000000 1.334328e-01
## MgSTS574B 3 74.884218 1.351110e-01
## c3.loc75 3 75.000000 1.348976e-01
## c3.loc76 3 76.000000 1.329359e-01
## c3.loc77 3 77.000000 1.308056e-01
## c3.loc78 3 78.000000 1.285423e-01
## c3.loc79 3 79.000000 1.262045e-01
## c3.loc80 3 80.000000 1.238787e-01
## c3.loc81 3 81.000000 1.216812e-01
## c3.loc82 3 82.000000 1.197584e-01
## c3.loc83 3 83.000000 1.182823e-01
## c3.loc84 3 84.000000 1.174415e-01
## c3.loc85 3 85.000000 1.174274e-01
## c3.loc86 3 86.000000 1.184167e-01
## c3.loc87 3 87.000000 1.205529e-01
## c3.loc88 3 88.000000 1.239291e-01
## c3.loc89 3 89.000000 1.285764e-01
## c3.loc90 3 90.000000 1.344600e-01
## c3.loc91 3 91.000000 1.414821e-01
## c3.loc92 3 92.000000 1.494931e-01
## c3.loc93 3 93.000000 1.583057e-01
## c3.loc94 3 94.000000 1.677123e-01
## c3.loc95 3 95.000000 1.775005e-01
## c3.loc96 3 96.000000 1.874664e-01
## CC286 3 96.768543 1.951116e-01
## c3.loc97 3 97.000000 1.980056e-01
## c3.loc98 3 98.000000 2.111510e-01
## c3.loc99 3 99.000000 2.251027e-01
## c3.loc100 3 100.000000 2.396925e-01
## c3.loc101 3 101.000000 2.546781e-01
## c3.loc102 3 102.000000 2.697500e-01
## c3.loc103 3 103.000000 2.845505e-01
## c3.loc104 3 104.000000 2.987039e-01
## c3.loc105 3 105.000000 3.118534e-01
## c3.loc106 3 106.000000 3.236971e-01
## c3.loc107 3 107.000000 3.340162e-01
## c3.loc108 3 108.000000 3.426893e-01
## c3.loc109 3 109.000000 3.496927e-01
## c3.loc110 3 110.000000 3.550870e-01
## c3.loc111 3 111.000000 3.589973e-01
## c3.loc112 3 112.000000 3.615896e-01
## c3.loc113 3 113.000000 3.630499e-01
## AG19 3 113.246604 3.632565e-01
## c3.loc114 3 114.000000 3.469917e-01
## c3.loc115 3 115.000000 3.255195e-01
## c3.loc116 3 116.000000 3.044351e-01
## c3.loc117 3 117.000000 2.839810e-01
## c3.loc118 3 118.000000 2.644242e-01
## c3.loc119 3 119.000000 2.460418e-01
## c3.loc120 3 120.000000 2.291023e-01
## c3.loc121 3 121.000000 2.138418e-01
## c3.loc122 3 122.000000 2.004436e-01
## c3.loc123 3 123.000000 1.890225e-01
## c3.loc124 3 124.000000 1.796192e-01
## c3.loc125 3 125.000000 1.722051e-01
## c3.loc126 3 126.000000 1.666946e-01
## c3.loc127 3 127.000000 1.629626e-01
## c3.loc128 3 128.000000 1.608616e-01
## BD286 3 128.645982 1.602998e-01
## c3.loc129 3 129.000000 1.694427e-01
## c3.loc130 3 130.000000 1.955886e-01
## c3.loc131 3 131.000000 2.192085e-01
## c3.loc132 3 132.000000 2.371443e-01
## c3.loc133 3 133.000000 2.470798e-01
## MgSTS474 3 133.645130 2.486506e-01
## c3.loc134 3 134.000000 2.605778e-01
## c3.loc135 3 135.000000 2.924459e-01
## c3.loc136 3 136.000000 3.186906e-01
## c3.loc137 3 137.000000 3.369015e-01
## c3.loc138 3 138.000000 3.465423e-01
## MgSTS214 3 138.390567 3.481963e-01
## c3.loc139 3 139.000000 2.844597e-01
## c3.loc140 3 140.000000 1.864398e-01
## c3.loc141 3 141.000000 1.055442e-01
## c3.loc142 3 142.000000 4.961381e-02
## c3.loc143 3 143.000000 2.246284e-02
## CA220 3 143.614625 1.985463e-02
## c3.loc144 3 144.000000 1.981546e-02
## c3.loc145 3 145.000000 1.974698e-02
## c3.loc146 3 146.000000 1.967340e-02
## c3.loc147 3 147.000000 1.952314e-02
## c3.loc148 3 148.000000 1.922168e-02
## c3.loc149 3 149.000000 1.869558e-02
## c3.loc150 3 150.000000 1.788008e-02
## c3.loc151 3 151.000000 1.673041e-02
## c3.loc152 3 152.000000 1.523518e-02
## c3.loc153 3 153.000000 1.342816e-02
## c3.loc154 3 154.000000 1.139391e-02
## c3.loc155 3 155.000000 9.263859e-03
## c3.loc156 3 156.000000 7.202819e-03
## c3.loc157 3 157.000000 5.389338e-03
## MgSTS425 3 157.887462 4.127328e-03
## c3.loc158 3 158.000000 3.789792e-03
## c3.loc159 3 159.000000 1.571633e-03
## c3.loc160 3 160.000000 1.080407e-03
## c3.loc161 3 161.000000 2.642628e-03
## c3.loc162 3 162.000000 6.444024e-03
## c3.loc163 3 163.000000 1.247341e-02
## c3.loc164 3 164.000000 2.050738e-02
## c3.loc165 3 165.000000 3.014395e-02
## c3.loc166 3 166.000000 4.087300e-02
## c3.loc167 3 167.000000 5.215814e-02
## c3.loc168 3 168.000000 6.350596e-02
## c3.loc169 3 169.000000 7.450966e-02
## MgSTS95 3 169.392376 7.863910e-02
## c3.loc170 3 170.000000 7.653852e-02
## c3.loc171 3 171.000000 7.512299e-02
## c3.loc172 3 172.000000 7.666531e-02
## c3.loc173 3 173.000000 8.156874e-02
## MgSTS50 3 173.527693 8.555857e-02
## c3.loc174 3 174.000000 9.299078e-02
## c3.loc175 3 175.000000 1.135151e-01
## c3.loc176 3 176.000000 1.411748e-01
## c3.loc177 3 177.000000 1.768538e-01
## MgSTS535 3 177.362540 1.918144e-01
## c3.loc178 3 178.000000 1.998092e-01
## c3.loc179 3 179.000000 2.105340e-01
## c3.loc180 3 180.000000 2.182790e-01
## c3.loc181 3 181.000000 2.224666e-01
## c3.loc182 3 182.000000 2.228769e-01
## c3.loc183 3 183.000000 2.196870e-01
## CC378 3 183.025773 2.195625e-01
## CB329 4 0.000000 4.961021e-02
## c4.loc1 4 1.000000 5.474728e-02
## c4.loc2 4 2.000000 6.038103e-02
## c4.loc3 4 3.000000 6.629177e-02
## c4.loc4 4 4.000000 7.219513e-02
## MgSTS542B 4 4.183834 7.324563e-02
## c4.loc5 4 5.000000 5.633836e-02
## c4.loc6 4 6.000000 3.831224e-02
## c4.loc7 4 7.000000 2.401701e-02
## c4.loc8 4 8.000000 1.395586e-02
## c4.loc9 4 9.000000 8.227386e-03
## c4.loc10 4 10.000000 6.493122e-03
## BC512 4 10.403851 6.771804e-03
## c4.loc11 4 11.000000 6.752430e-03
## c4.loc12 4 12.000000 1.417047e-02
## c4.loc13 4 13.000000 2.854633e-02
## c4.loc14 4 14.000000 4.642546e-02
## MgSTS347 4 14.489772 5.534531e-02
## c4.loc15 4 15.000000 4.771986e-02
## c4.loc16 4 16.000000 3.424536e-02
## c4.loc17 4 17.000000 2.327100e-02
## c4.loc18 4 18.000000 1.503593e-02
## c4.loc19 4 19.000000 9.413560e-03
## MgSTS132 4 19.992076 6.005071e-03
## c4.loc20 4 20.000000 6.048811e-03
## c4.loc21 4 21.000000 1.273081e-02
## c4.loc22 4 22.000000 2.127405e-02
## c4.loc23 4 23.000000 3.163554e-02
## c4.loc24 4 24.000000 4.443138e-02
## c4.loc25 4 25.000000 6.091371e-02
## c4.loc26 4 26.000000 8.283304e-02
## BA210 4 26.991980 1.118714e-01
## c4.loc27 4 27.000000 1.121016e-01
## CA399 4 27.654986 1.330897e-01
## c4.loc28 4 28.000000 1.293924e-01
## c4.loc29 4 29.000000 1.193685e-01
## c4.loc30 4 30.000000 1.107395e-01
## c4.loc31 4 31.000000 1.038228e-01
## c4.loc32 4 32.000000 9.888750e-02
## c4.loc33 4 33.000000 9.611580e-02
## c4.loc34 4 34.000000 9.557005e-02
## c4.loc35 4 35.000000 9.717328e-02
## c4.loc36 4 36.000000 1.007104e-01
## c4.loc37 4 37.000000 1.058516e-01
## c4.loc38 4 38.000000 1.121931e-01
## c4.loc39 4 39.000000 1.193060e-01
## c4.loc40 4 40.000000 1.267834e-01
## CA415 4 40.302047 1.290433e-01
## c4.loc41 4 41.000000 1.151259e-01
## c4.loc42 4 42.000000 9.566295e-02
## c4.loc43 4 43.000000 7.712093e-02
## c4.loc44 4 44.000000 5.985139e-02
## c4.loc45 4 45.000000 4.422543e-02
## c4.loc46 4 46.000000 3.061004e-02
## c4.loc47 4 47.000000 1.933666e-02
## c4.loc48 4 48.000000 1.066538e-02
## c4.loc49 4 49.000000 4.750521e-03
## c4.loc50 4 50.000000 1.614812e-03
## c4.loc51 4 51.000000 1.138588e-03
## c4.loc52 4 52.000000 3.067466e-03
## c4.loc53 4 53.000000 7.037551e-03
## c4.loc54 4 54.000000 1.261296e-02
## c4.loc55 4 55.000000 1.932807e-02
## c4.loc56 4 56.000000 2.672682e-02
## c4.loc57 4 57.000000 3.439366e-02
## c4.loc58 4 58.000000 4.197367e-02
## c4.loc59 4 59.000000 4.918196e-02
## CA497 4 59.029103 4.938190e-02
## c4.loc60 4 60.000000 4.013444e-02
## c4.loc61 4 61.000000 3.074469e-02
## c4.loc62 4 62.000000 2.222064e-02
## c4.loc63 4 63.000000 1.528855e-02
## c4.loc64 4 64.000000 1.052370e-02
## c4.loc65 4 65.000000 8.238272e-03
## c4.loc66 4 66.000000 8.444025e-03
## BD243 4 66.195272 8.754579e-03
## c4.loc67 4 67.000000 9.087547e-03
## c4.loc68 4 68.000000 1.090598e-02
## c4.loc69 4 69.000000 1.439502e-02
## c4.loc70 4 70.000000 1.955840e-02
## c4.loc71 4 71.000000 2.625345e-02
## c4.loc72 4 72.000000 3.419435e-02
## c4.loc73 4 73.000000 4.298605e-02
## c4.loc74 4 74.000000 5.218070e-02
## c4.loc75 4 75.000000 6.134036e-02
## AA66 4 75.956325 6.971639e-02
## c4.loc76 4 76.000000 6.915972e-02
## c4.loc77 4 77.000000 5.688708e-02
## c4.loc78 4 78.000000 4.615431e-02
## c4.loc79 4 79.000000 3.738853e-02
## c4.loc80 4 80.000000 3.079508e-02
## c4.loc81 4 81.000000 2.637447e-02
## c4.loc82 4 82.000000 2.396671e-02
## CA183C 4 82.575110 2.339494e-02
## c4.loc83 4 83.000000 2.254742e-02
## c4.loc84 4 84.000000 2.230552e-02
## c4.loc85 4 85.000000 2.477410e-02
## c4.loc86 4 86.000000 3.007847e-02
## c4.loc87 4 87.000000 3.812169e-02
## c4.loc88 4 88.000000 4.858052e-02
## c4.loc89 4 89.000000 6.094616e-02
## c4.loc90 4 90.000000 7.459965e-02
## c4.loc91 4 91.000000 8.889979e-02
## c4.loc92 4 92.000000 1.032600e-01
## AA166C 4 92.962317 1.166806e-01
## c4.loc93 4 93.000000 1.142743e-01
## c4.loc94 4 94.000000 7.365955e-02
## c4.loc95 4 95.000000 7.872209e-02
## c4.loc96 4 96.000000 1.247917e-01
## c4.loc97 4 97.000000 2.025244e-01
## c4.loc98 4 98.000000 2.994748e-01
## c4.loc99 4 99.000000 4.019486e-01
## c4.loc100 4 100.000000 4.968928e-01
## MgSTS477 4 100.811402 5.608301e-01
## c4.loc101 4 101.000000 5.605575e-01
## c4.loc102 4 102.000000 5.526463e-01
## c4.loc103 4 103.000000 5.329181e-01
## c4.loc104 4 104.000000 5.014910e-01
## c4.loc105 4 105.000000 4.604268e-01
## c4.loc106 4 106.000000 4.132959e-01
## MgSTS362 4 106.389774 3.942893e-01
## c4.loc107 4 107.000000 3.680492e-01
## c4.loc108 4 108.000000 3.215360e-01
## c4.loc109 4 109.000000 2.721376e-01
## c4.loc110 4 110.000000 2.219203e-01
## c4.loc111 4 111.000000 1.733860e-01
## c4.loc112 4 112.000000 1.290670e-01
## c4.loc113 4 113.000000 9.106777e-02
## c4.loc114 4 114.000000 6.071848e-02
## c4.loc115 4 115.000000 3.844912e-02
## MgSTS262 4 115.872250 2.536113e-02
## c4.loc116 4 116.000000 2.900519e-02
## c4.loc117 4 117.000000 7.151927e-02
## c4.loc118 4 118.000000 1.398832e-01
## c4.loc119 4 119.000000 2.345162e-01
## c4.loc120 4 120.000000 3.531522e-01
## c4.loc121 4 121.000000 4.904856e-01
## c4.loc122 4 122.000000 6.384967e-01
## c4.loc123 4 123.000000 7.875464e-01
## c4.loc124 4 124.000000 9.279885e-01
## c4.loc125 4 125.000000 1.051774e+00
## MgSTS492 4 125.794648 1.134433e+00
## c4.loc126 4 126.000000 1.136483e+00
## c4.loc127 4 127.000000 1.140366e+00
## c4.loc128 4 128.000000 1.133121e+00
## c4.loc129 4 129.000000 1.114253e+00
## c4.loc130 4 130.000000 1.084263e+00
## c4.loc131 4 131.000000 1.044609e+00
## c4.loc132 4 132.000000 9.974903e-01
## c4.loc133 4 133.000000 9.454818e-01
## c4.loc134 4 134.000000 8.911598e-01
## MgSTS228 4 134.700064 8.531024e-01
## c4.loc135 4 135.000000 8.812774e-01
## c4.loc136 4 136.000000 9.481796e-01
## c4.loc137 4 137.000000 9.643921e-01
## BB208 4 137.123795 9.624957e-01
## c4.loc138 4 138.000000 1.041991e+00
## c4.loc139 4 139.000000 1.122525e+00
## CA297 4 139.713603 1.164946e+00
## c4.loc140 4 140.000000 1.160411e+00
## c4.loc141 4 141.000000 1.131772e+00
## c4.loc142 4 142.000000 1.083235e+00
## c4.loc143 4 143.000000 1.016527e+00
## c4.loc144 4 144.000000 9.355177e-01
## c4.loc145 4 145.000000 8.454962e-01
## c4.loc146 4 146.000000 7.521166e-01
## BA279C 4 146.521949 7.039955e-01
## c4.loc147 4 147.000000 7.163077e-01
## c4.loc148 4 148.000000 7.379706e-01
## c4.loc149 4 149.000000 7.522203e-01
## c4.loc150 4 150.000000 7.571883e-01
## c4.loc151 4 151.000000 7.515140e-01
## c4.loc152 4 152.000000 7.346467e-01
## c4.loc153 4 153.000000 7.070063e-01
## c4.loc154 4 154.000000 6.699401e-01
## BD130 4 154.968905 6.269746e-01
## c4.loc155 4 155.000000 6.241479e-01
## c4.loc156 4 156.000000 5.393167e-01
## BA113 4 156.864152 4.800528e-01
## c4.loc157 4 157.000000 4.774264e-01
## c4.loc158 4 158.000000 4.613985e-01
## c4.loc159 4 159.000000 4.529311e-01
## c4.loc160 4 160.000000 4.532975e-01
## AA346C 4 160.486475 4.566847e-01
## c4.loc161 4 161.000000 4.451242e-01
## c4.loc162 4 162.000000 4.379202e-01
## c4.loc163 4 163.000000 4.538980e-01
## c4.loc164 4 164.000000 4.958659e-01
## c4.loc165 4 165.000000 5.656789e-01
## c4.loc166 4 166.000000 6.638225e-01
## c4.loc167 4 167.000000 7.890710e-01
## c4.loc168 4 168.000000 9.383048e-01
## c4.loc169 4 169.000000 1.106557e+00
## c4.loc170 4 170.000000 1.287318e+00
## c4.loc171 4 171.000000 1.473060e+00
## c4.loc172 4 172.000000 1.655906e+00
## BC192 4 172.962577 1.822060e+00
## c4.loc173 4 173.000000 1.823495e+00
## c4.loc174 4 174.000000 1.858815e+00
## c4.loc175 4 175.000000 1.885993e+00
## c4.loc176 4 176.000000 1.902593e+00
## c4.loc177 4 177.000000 1.906290e+00
## c4.loc178 4 178.000000 1.895196e+00
## c4.loc179 4 179.000000 1.868219e+00
## c4.loc180 4 180.000000 1.825368e+00
## c4.loc181 4 181.000000 1.767888e+00
## c4.loc182 4 182.000000 1.698184e+00
## c4.loc183 4 183.000000 1.619513e+00
## c4.loc184 4 184.000000 1.535529e+00
## c4.loc185 4 185.000000 1.449813e+00
## c4.loc186 4 186.000000 1.365485e+00
## c4.loc187 4 187.000000 1.284977e+00
## BC216 4 187.690596 1.232671e+00
## c4.loc188 4 188.000000 1.248058e+00
## c4.loc189 4 189.000000 1.293709e+00
## c4.loc190 4 190.000000 1.332402e+00
## c4.loc191 4 191.000000 1.364004e+00
## c4.loc192 4 192.000000 1.388483e+00
## c4.loc193 4 193.000000 1.405861e+00
## c4.loc194 4 194.000000 1.416162e+00
## c4.loc195 4 195.000000 1.419372e+00
## c4.loc196 4 196.000000 1.415406e+00
## c4.loc197 4 197.000000 1.404085e+00
## c4.loc198 4 198.000000 1.385135e+00
## c4.loc199 4 199.000000 1.358191e+00
## c4.loc200 4 200.000000 1.322829e+00
## c4.loc201 4 201.000000 1.278613e+00
## c4.loc202 4 202.000000 1.225171e+00
## c4.loc203 4 203.000000 1.162286e+00
## c4.loc204 4 204.000000 1.090020e+00
## CA384 4 204.377632 1.060555e+00
## c4.loc205 4 205.000000 1.069990e+00
## c4.loc206 4 206.000000 1.084784e+00
## c4.loc207 4 207.000000 1.098799e+00
## c4.loc208 4 208.000000 1.111687e+00
## c4.loc209 4 209.000000 1.123044e+00
## c4.loc210 4 210.000000 1.132414e+00
## c4.loc211 4 211.000000 1.139292e+00
## c4.loc212 4 212.000000 1.143152e+00
## c4.loc213 4 213.000000 1.143460e+00
## c4.loc214 4 214.000000 1.139719e+00
## c4.loc215 4 215.000000 1.131495e+00
## c4.loc216 4 216.000000 1.118469e+00
## c4.loc217 4 217.000000 1.100469e+00
## c4.loc218 4 218.000000 1.077502e+00
## c4.loc219 4 219.000000 1.049770e+00
## c4.loc220 4 220.000000 1.017667e+00
## c4.loc221 4 221.000000 9.817576e-01
## c4.loc222 4 222.000000 9.427414e-01
## c4.loc223 4 223.000000 9.014014e-01
## c4.loc224 4 224.000000 8.585523e-01
## c4.loc225 4 225.000000 8.149894e-01
## c4.loc226 4 226.000000 7.714463e-01
## c4.loc227 4 227.000000 7.285647e-01
## c4.loc228 4 228.000000 6.868765e-01
## c4.loc229 4 229.000000 6.467980e-01
## c4.loc230 4 230.000000 6.086331e-01
## c4.loc231 4 231.000000 5.725837e-01
## c4.loc232 4 232.000000 5.387636e-01
## BD411 4 232.552539 5.211273e-01
## c4.loc233 4 233.000000 5.521964e-01
## c4.loc234 4 234.000000 6.265197e-01
## c4.loc235 4 235.000000 7.069388e-01
## c4.loc236 4 236.000000 7.931240e-01
## c4.loc237 4 237.000000 8.845530e-01
## c4.loc238 4 238.000000 9.805025e-01
## c4.loc239 4 239.000000 1.080054e+00
## c4.loc240 4 240.000000 1.182108e+00
## c4.loc241 4 241.000000 1.285420e+00
## c4.loc242 4 242.000000 1.388640e+00
## c4.loc243 4 243.000000 1.490369e+00
## c4.loc244 4 244.000000 1.589215e+00
## c4.loc245 4 245.000000 1.683854e+00
## c4.loc246 4 246.000000 1.773086e+00
## c4.loc247 4 247.000000 1.855876e+00
## c4.loc248 4 248.000000 1.931386e+00
## c4.loc249 4 249.000000 1.998998e+00
## c4.loc250 4 250.000000 2.058314e+00
## c4.loc251 4 251.000000 2.109152e+00
## c4.loc252 4 252.000000 2.151527e+00
## c4.loc253 4 253.000000 2.185630e+00
## c4.loc254 4 254.000000 2.211798e+00
## c4.loc255 4 255.000000 2.230486e+00
## c4.loc256 4 256.000000 2.242236e+00
## c4.loc257 4 257.000000 2.247652e+00
## c4.loc258 4 258.000000 2.247372e+00
## MgSTS455 4 258.472195 2.245460e+00
## c4.loc259 4 259.000000 2.301516e+00
## c4.loc260 4 260.000000 2.393507e+00
## c4.loc261 4 261.000000 2.461197e+00
## c4.loc262 4 262.000000 2.499430e+00
## c4.loc263 4 263.000000 2.504676e+00
## c4.loc264 4 264.000000 2.475707e+00
## c4.loc265 4 265.000000 2.413934e+00
## c4.loc266 4 266.000000 2.323295e+00
## c4.loc267 4 267.000000 2.209710e+00
## c4.loc268 4 268.000000 2.080256e+00
## BC321 4 268.310983 2.038213e+00
## c4.loc269 4 269.000000 1.994881e+00
## c4.loc270 4 270.000000 1.926338e+00
## c4.loc271 4 271.000000 1.852147e+00
## c4.loc272 4 272.000000 1.773289e+00
## c4.loc273 4 273.000000 1.690773e+00
## c4.loc274 4 274.000000 1.605604e+00
## c4.loc275 4 275.000000 1.518754e+00
## c4.loc276 4 276.000000 1.431153e+00
## c4.loc277 4 277.000000 1.343675e+00
## c4.loc278 4 278.000000 1.257134e+00
## c4.loc279 4 279.000000 1.172283e+00
## c4.loc280 4 280.000000 1.089801e+00
## c4.loc281 4 281.000000 1.010295e+00
## c4.loc282 4 282.000000 9.342873e-01
## c4.loc283 4 283.000000 8.622110e-01
## CC270 4 283.893499 8.014742e-01
## c4.loc284 4 284.000000 7.956216e-01
## c4.loc285 4 285.000000 7.485042e-01
## c4.loc286 4 286.000000 7.160577e-01
## c4.loc287 4 287.000000 6.969277e-01
## c4.loc288 4 288.000000 6.888754e-01
## c4.loc289 4 289.000000 6.890474e-01
## c4.loc290 4 290.000000 6.943280e-01
## c4.loc291 4 291.000000 7.017077e-01
## c4.loc292 4 292.000000 7.085980e-01
## c4.loc293 4 293.000000 7.130367e-01
## MgSTS234 4 293.638872 7.139801e-01
## c4.loc294 4 294.000000 7.038594e-01
## c4.loc295 4 295.000000 6.727326e-01
## c4.loc296 4 296.000000 6.378627e-01
## c4.loc297 4 297.000000 6.002266e-01
## c4.loc298 4 298.000000 5.610269e-01
## c4.loc299 4 299.000000 5.215375e-01
## c4.loc300 4 300.000000 4.829270e-01
## c4.loc301 4 301.000000 4.461094e-01
## CB257 4 301.727951 4.208436e-01
## c4.loc302 4 302.000000 4.146374e-01
## c4.loc303 4 303.000000 3.910241e-01
## c4.loc304 4 304.000000 3.665206e-01
## c4.loc305 4 305.000000 3.412951e-01
## c4.loc306 4 306.000000 3.155418e-01
## c4.loc307 4 307.000000 2.894774e-01
## c4.loc308 4 308.000000 2.633355e-01
## c4.loc309 4 309.000000 2.373618e-01
## c4.loc310 4 310.000000 2.118076e-01
## c4.loc311 4 311.000000 1.869253e-01
## c4.loc312 4 312.000000 1.629625e-01
## c4.loc313 4 313.000000 1.401579e-01
## c4.loc314 4 314.000000 1.187361e-01
## c4.loc315 4 315.000000 9.890351e-02
## c4.loc316 4 316.000000 8.084300e-02
## c4.loc317 4 317.000000 6.470941e-02
## c4.loc318 4 318.000000 5.062464e-02
## c4.loc319 4 319.000000 3.867355e-02
## c4.loc320 4 320.000000 2.890096e-02
## c4.loc321 4 321.000000 2.131021e-02
## c4.loc322 4 322.000000 1.586366e-02
## c4.loc323 4 323.000000 1.248516e-02
## AAT367 4 323.342539 1.178833e-02
## c4.loc324 4 324.000000 1.192097e-02
## c4.loc325 4 325.000000 1.210754e-02
## c4.loc326 4 326.000000 1.226977e-02
## c4.loc327 4 327.000000 1.240151e-02
## c4.loc328 4 328.000000 1.249642e-02
## c4.loc329 4 329.000000 1.254816e-02
## c4.loc330 4 330.000000 1.255062e-02
## c4.loc331 4 331.000000 1.249822e-02
## c4.loc332 4 332.000000 1.238617e-02
## c4.loc333 4 333.000000 1.221071e-02
## c4.loc334 4 334.000000 1.196943e-02
## c4.loc335 4 335.000000 1.166143e-02
## c4.loc336 4 336.000000 1.128748e-02
## c4.loc337 4 337.000000 1.085013e-02
## c4.loc338 4 338.000000 1.035368e-02
## c4.loc339 4 339.000000 9.804095e-03
## c4.loc340 4 340.000000 9.208878e-03
## c4.loc341 4 341.000000 8.576795e-03
## c4.loc342 4 342.000000 7.917620e-03
## c4.loc343 4 343.000000 7.241838e-03
## c4.loc344 4 344.000000 6.560347e-03
## c4.loc345 4 345.000000 5.884196e-03
## c4.loc346 4 346.000000 5.224363e-03
## c4.loc347 4 347.000000 4.591596e-03
## c4.loc348 4 348.000000 3.996307e-03
## c4.loc349 4 349.000000 3.448521e-03
## AA384 4 349.121610 3.386156e-03
## BB186 5 0.000000 2.965114e-01
## c5.loc1 5 1.000000 2.720590e-01
## c5.loc2 5 2.000000 2.471939e-01
## c5.loc3 5 3.000000 2.225545e-01
## c5.loc4 5 4.000000 1.986826e-01
## c5.loc5 5 5.000000 1.760057e-01
## c5.loc6 5 6.000000 1.548330e-01
## CC124 5 6.843141 1.383227e-01
## c5.loc7 5 7.000000 1.343850e-01
## c5.loc8 5 8.000000 1.113205e-01
## c5.loc9 5 9.000000 9.286705e-02
## c5.loc10 5 10.000000 7.970754e-02
## c5.loc11 5 11.000000 7.209237e-02
## c5.loc12 5 12.000000 6.984029e-02
## c5.loc13 5 13.000000 7.238812e-02
## c5.loc14 5 14.000000 7.887930e-02
## c5.loc15 5 15.000000 8.827634e-02
## c5.loc16 5 16.000000 9.947990e-02
## MgSTS40 5 16.815324 1.091990e-01
## c5.loc17 5 17.000000 1.128067e-01
## c5.loc18 5 18.000000 1.352952e-01
## c5.loc19 5 19.000000 1.628471e-01
## c5.loc20 5 20.000000 1.959212e-01
## c5.loc21 5 21.000000 2.349483e-01
## c5.loc22 5 22.000000 2.803061e-01
## c5.loc23 5 23.000000 3.322861e-01
## c5.loc24 5 24.000000 3.910534e-01
## c5.loc25 5 25.000000 4.566017e-01
## c5.loc26 5 26.000000 5.287091e-01
## c5.loc27 5 27.000000 6.069003e-01
## c5.loc28 5 28.000000 6.904252e-01
## c5.loc29 5 29.000000 7.782618e-01
## c5.loc30 5 30.000000 8.691468e-01
## c5.loc31 5 31.000000 9.616357e-01
## c5.loc32 5 32.000000 1.054186e+00
## c5.loc33 5 33.000000 1.145250e+00
## c5.loc34 5 34.000000 1.233374e+00
## c5.loc35 5 35.000000 1.317273e+00
## c5.loc36 5 36.000000 1.395899e+00
## c5.loc37 5 37.000000 1.468467e+00
## c5.loc38 5 38.000000 1.534467e+00
## MgSTS586 5 38.146257 1.543475e+00
## c5.loc39 5 39.000000 1.549705e+00
## c5.loc40 5 40.000000 1.548528e+00
## c5.loc41 5 41.000000 1.537317e+00
## c5.loc42 5 42.000000 1.515557e+00
## c5.loc43 5 43.000000 1.483214e+00
## c5.loc44 5 44.000000 1.440735e+00
## CB230 5 44.251410 1.428649e+00
## c5.loc45 5 45.000000 1.430593e+00
## BB210 5 45.881845 1.444285e+00
## c5.loc46 5 46.000000 1.436987e+00
## c5.loc47 5 47.000000 1.365401e+00
## BD316C 5 47.696944 1.307019e+00
## c5.loc48 5 48.000000 1.342072e+00
## c5.loc49 5 49.000000 1.453655e+00
## c5.loc50 5 50.000000 1.553957e+00
## c5.loc51 5 51.000000 1.638671e+00
## c5.loc52 5 52.000000 1.704984e+00
## c5.loc53 5 53.000000 1.752248e+00
## c5.loc54 5 54.000000 1.782264e+00
## c5.loc55 5 55.000000 1.799053e+00
## AA378 5 55.824982 1.806861e+00
## c5.loc56 5 56.000000 1.811744e+00
## c5.loc57 5 57.000000 1.831334e+00
## c5.loc58 5 58.000000 1.837081e+00
## c5.loc59 5 59.000000 1.831864e+00
## c5.loc60 5 60.000000 1.820698e+00
## c5.loc61 5 61.000000 1.809724e+00
## c5.loc62 5 62.000000 1.804757e+00
## BD173 5 62.525368 1.806042e+00
## c5.loc63 5 63.000000 1.847136e+00
## c5.loc64 5 64.000000 1.927876e+00
## c5.loc65 5 65.000000 1.999849e+00
## c5.loc66 5 66.000000 2.063670e+00
## c5.loc67 5 67.000000 2.120709e+00
## CA264 5 67.456667 2.144825e+00
## c5.loc68 5 68.000000 2.174289e+00
## c5.loc69 5 69.000000 2.223322e+00
## c5.loc70 5 70.000000 2.264544e+00
## c5.loc71 5 71.000000 2.297136e+00
## c5.loc72 5 72.000000 2.320356e+00
## c5.loc73 5 73.000000 2.333531e+00
## c5.loc74 5 74.000000 2.336049e+00
## c5.loc75 5 75.000000 2.327359e+00
## c5.loc76 5 76.000000 2.307015e+00
## c5.loc77 5 77.000000 2.274748e+00
## c5.loc78 5 78.000000 2.230577e+00
## MgSTS245 5 78.782071 2.188080e+00
## c5.loc79 5 79.000000 2.213919e+00
## c5.loc80 5 80.000000 2.325191e+00
## c5.loc81 5 81.000000 2.418639e+00
## c5.loc82 5 82.000000 2.488680e+00
## c5.loc83 5 83.000000 2.530565e+00
## c5.loc84 5 84.000000 2.541025e+00
## c5.loc85 5 85.000000 2.518854e+00
## c5.loc86 5 86.000000 2.465294e+00
## c5.loc87 5 87.000000 2.384095e+00
## AA268 5 87.891587 2.293265e+00
## c5.loc88 5 88.000000 2.297239e+00
## c5.loc89 5 89.000000 2.331069e+00
## c5.loc90 5 90.000000 2.358469e+00
## c5.loc91 5 91.000000 2.378380e+00
## c5.loc92 5 92.000000 2.389798e+00
## c5.loc93 5 93.000000 2.391854e+00
## c5.loc94 5 94.000000 2.383891e+00
## c5.loc95 5 95.000000 2.365552e+00
## CC342 5 95.546866 2.351188e+00
## c5.loc96 5 96.000000 2.337313e+00
## c5.loc97 5 97.000000 2.298045e+00
## c5.loc98 5 98.000000 2.248084e+00
## c5.loc99 5 99.000000 2.188882e+00
## c5.loc100 5 100.000000 2.122027e+00
## c5.loc101 5 101.000000 2.049114e+00
## c5.loc102 5 102.000000 1.971663e+00
## AA454C 5 102.896157 1.899632e+00
## c5.loc103 5 103.000000 1.900594e+00
## c5.loc104 5 104.000000 1.908428e+00
## c5.loc105 5 105.000000 1.913840e+00
## c5.loc106 5 106.000000 1.917396e+00
## c5.loc107 5 107.000000 1.919734e+00
## c5.loc108 5 108.000000 1.921466e+00
## c5.loc109 5 109.000000 1.923058e+00
## c5.loc110 5 110.000000 1.924719e+00
## c5.loc111 5 111.000000 1.926323e+00
## c5.loc112 5 112.000000 1.927385e+00
## AA163 5 112.220121 1.927468e+00
## c5.loc113 5 113.000000 1.928064e+00
## c5.loc114 5 114.000000 1.927102e+00
## c5.loc115 5 115.000000 1.923897e+00
## c5.loc116 5 116.000000 1.918022e+00
## c5.loc117 5 117.000000 1.908933e+00
## c5.loc118 5 118.000000 1.895970e+00
## c5.loc119 5 119.000000 1.878379e+00
## c5.loc120 5 120.000000 1.855364e+00
## c5.loc121 5 121.000000 1.826155e+00
## c5.loc122 5 122.000000 1.790102e+00
## c5.loc123 5 123.000000 1.746751e+00
## c5.loc124 5 124.000000 1.695920e+00
## c5.loc125 5 125.000000 1.637727e+00
## c5.loc126 5 126.000000 1.572596e+00
## c5.loc127 5 127.000000 1.501219e+00
## c5.loc128 5 128.000000 1.424507e+00
## BC506 5 128.407939 1.392113e+00
## c5.loc129 5 129.000000 1.340544e+00
## c5.loc130 5 130.000000 1.254196e+00
## c5.loc131 5 131.000000 1.169781e+00
## c5.loc132 5 132.000000 1.087883e+00
## c5.loc133 5 133.000000 1.009022e+00
## c5.loc134 5 134.000000 9.336482e-01
## c5.loc135 5 135.000000 8.621366e-01
## c5.loc136 5 136.000000 7.947771e-01
## c5.loc137 5 137.000000 7.317694e-01
## c5.loc138 5 138.000000 6.732206e-01
## c5.loc139 5 139.000000 6.191462e-01
## c5.loc140 5 140.000000 5.694768e-01
## c5.loc141 5 141.000000 5.240694e-01
## c5.loc142 5 142.000000 4.827225e-01
## c5.loc143 5 143.000000 4.451920e-01
## c5.loc144 5 144.000000 4.112081e-01
## c5.loc145 5 145.000000 3.804899e-01
## c5.loc146 5 146.000000 3.527565e-01
## c5.loc147 5 147.000000 3.277364e-01
## MgSTS282A 5 147.593699 3.141133e-01
## c5.loc148 5 148.000000 3.153887e-01
## c5.loc149 5 149.000000 3.057234e-01
## c5.loc150 5 150.000000 2.802125e-01
## c5.loc151 5 151.000000 2.428159e-01
## MgSTS255 5 151.641858 2.148873e-01
## c5.loc152 5 152.000000 2.118613e-01
## c5.loc153 5 153.000000 2.033830e-01
## c5.loc154 5 154.000000 1.950029e-01
## c5.loc155 5 155.000000 1.868062e-01
## c5.loc156 5 156.000000 1.788855e-01
## c5.loc157 5 157.000000 1.713389e-01
## c5.loc158 5 158.000000 1.642673e-01
## c5.loc159 5 159.000000 1.577706e-01
## c5.loc160 5 160.000000 1.519440e-01
## c5.loc161 5 161.000000 1.468740e-01
## c5.loc162 5 162.000000 1.426339e-01
## c5.loc163 5 163.000000 1.392808e-01
## c5.loc164 5 164.000000 1.368529e-01
## c5.loc165 5 165.000000 1.353683e-01
## c5.loc166 5 166.000000 1.348258e-01
## c5.loc167 5 167.000000 1.352057e-01
## c5.loc168 5 168.000000 1.364733e-01
## c5.loc169 5 169.000000 1.385809e-01
## c5.loc170 5 170.000000 1.414716e-01
## c5.loc171 5 171.000000 1.450818e-01
## c5.loc172 5 172.000000 1.493432e-01
## c5.loc173 5 173.000000 1.541848e-01
## c5.loc174 5 174.000000 1.595339e-01
## c5.loc175 5 175.000000 1.653177e-01
## c5.loc176 5 176.000000 1.714639e-01
## c5.loc177 5 177.000000 1.779020e-01
## c5.loc178 5 178.000000 1.845641e-01
## c5.loc179 5 179.000000 1.913860e-01
## AA118C 5 179.683514 1.960950e-01
## CC457 6 0.000000 4.443749e-01
## c6.loc1 6 1.000000 4.766722e-01
## c6.loc2 6 2.000000 4.986671e-01
## c6.loc3 6 3.000000 5.089607e-01
## c6.loc4 6 4.000000 5.070769e-01
## c6.loc5 6 5.000000 4.935170e-01
## c6.loc6 6 6.000000 4.696803e-01
## c6.loc7 6 7.000000 4.376617e-01
## MgSTS508 6 7.932642 4.026583e-01
## c6.loc8 6 8.000000 4.171329e-01
## c6.loc9 6 9.000000 6.459503e-01
## MgSTS529 6 9.125086 6.738032e-01
## c6.loc10 6 10.000000 6.466034e-01
## c6.loc11 6 11.000000 6.146013e-01
## c6.loc12 6 12.000000 5.828628e-01
## c6.loc13 6 13.000000 5.525747e-01
## c6.loc14 6 14.000000 5.246947e-01
## c6.loc15 6 15.000000 4.998225e-01
## c6.loc16 6 16.000000 4.781503e-01
## c6.loc17 6 17.000000 4.595091e-01
## AA311 6 17.033907 4.589312e-01
## c6.loc18 6 18.000000 4.436936e-01
## c6.loc19 6 19.000000 4.275302e-01
## c6.loc20 6 20.000000 4.112428e-01
## c6.loc21 6 21.000000 3.950885e-01
## c6.loc22 6 22.000000 3.792974e-01
## c6.loc23 6 23.000000 3.640641e-01
## c6.loc24 6 24.000000 3.495459e-01
## MgSTS440 6 24.559826 3.418102e-01
## c6.loc25 6 25.000000 3.377333e-01
## c6.loc26 6 26.000000 3.371740e-01
## c6.loc27 6 27.000000 3.507058e-01
## c6.loc28 6 28.000000 3.797388e-01
## c6.loc29 6 29.000000 4.243152e-01
## c6.loc30 6 30.000000 4.830202e-01
## c6.loc31 6 31.000000 5.531567e-01
## c6.loc32 6 32.000000 6.311366e-01
## c6.loc33 6 33.000000 7.129820e-01
## c6.loc34 6 34.000000 7.948157e-01
## MgSTS25 6 34.836948 8.607082e-01
## c6.loc35 6 35.000000 8.502671e-01
## c6.loc36 6 36.000000 7.893470e-01
## c6.loc37 6 37.000000 7.384258e-01
## c6.loc38 6 38.000000 7.025172e-01
## c6.loc39 6 39.000000 6.854616e-01
## c6.loc40 6 40.000000 6.888775e-01
## c6.loc41 6 41.000000 7.117127e-01
## MgSTS542A 6 41.573283 7.322358e-01
## c6.loc42 6 42.000000 7.301894e-01
## c6.loc43 6 43.000000 7.225148e-01
## c6.loc44 6 44.000000 7.100773e-01
## c6.loc45 6 45.000000 6.919160e-01
## c6.loc46 6 46.000000 6.671813e-01
## c6.loc47 6 47.000000 6.353317e-01
## c6.loc48 6 48.000000 5.963323e-01
## c6.loc49 6 49.000000 5.507912e-01
## c6.loc50 6 50.000000 4.999789e-01
## c6.loc51 6 51.000000 4.457035e-01
## c6.loc52 6 52.000000 3.900715e-01
## MgSTS545 6 52.898334 3.407356e-01
## c6.loc53 6 53.000000 3.431059e-01
## c6.loc54 6 54.000000 3.711179e-01
## MgSTS314 6 54.931725 4.053483e-01
## c6.loc55 6 55.000000 4.069409e-01
## c6.loc56 6 56.000000 4.311305e-01
## c6.loc57 6 57.000000 4.556712e-01
## c6.loc58 6 58.000000 4.787743e-01
## MgSTS120 6 58.009326 4.789744e-01
## c6.loc59 6 59.000000 4.090872e-01
## c6.loc60 6 60.000000 3.373681e-01
## BD179 6 60.184324 3.244392e-01
## c6.loc61 6 61.000000 2.578668e-01
## c6.loc62 6 62.000000 1.824741e-01
## c6.loc63 6 63.000000 1.178001e-01
## MgSTS459 6 63.838498 7.398768e-02
## c6.loc64 6 64.000000 7.580324e-02
## c6.loc65 6 65.000000 9.099895e-02
## c6.loc66 6 66.000000 1.130055e-01
## c6.loc67 6 67.000000 1.415511e-01
## c6.loc68 6 68.000000 1.757187e-01
## c6.loc69 6 69.000000 2.139977e-01
## c6.loc70 6 70.000000 2.544255e-01
## c6.loc71 6 71.000000 2.948033e-01
## AA277 6 71.530555 3.153595e-01
## c6.loc72 6 72.000000 3.140224e-01
## c6.loc73 6 73.000000 3.099316e-01
## BA117 6 73.919497 3.046489e-01
## c6.loc74 6 74.000000 3.038373e-01
## c6.loc75 6 75.000000 2.933615e-01
## c6.loc76 6 76.000000 2.824704e-01
## c6.loc77 6 77.000000 2.714018e-01
## c6.loc78 6 78.000000 2.604177e-01
## c6.loc79 6 79.000000 2.497682e-01
## c6.loc80 6 80.000000 2.396496e-01
## c6.loc81 6 81.000000 2.301699e-01
## c6.loc82 6 82.000000 2.213378e-01
## c6.loc83 6 83.000000 2.130760e-01
## AA158 6 83.914642 2.059134e-01
## c6.loc84 6 84.000000 2.047752e-01
## c6.loc85 6 85.000000 1.926409e-01
## c6.loc86 6 86.000000 1.834445e-01
## c6.loc87 6 87.000000 1.777417e-01
## c6.loc88 6 88.000000 1.757130e-01
## c6.loc89 6 89.000000 1.770391e-01
## c6.loc90 6 90.000000 1.809266e-01
## c6.loc91 6 91.000000 1.862996e-01
## c6.loc92 6 92.000000 1.920669e-01
## AA270 6 92.010223 1.921240e-01
## c6.loc93 6 93.000000 1.689646e-01
## c6.loc94 6 94.000000 1.468002e-01
## c6.loc95 6 95.000000 1.274001e-01
## c6.loc96 6 96.000000 1.119894e-01
## c6.loc97 6 97.000000 1.011285e-01
## c6.loc98 6 98.000000 9.461245e-02
## MgSTS105 6 98.240082 9.362264e-02
## c6.loc99 6 99.000000 9.734185e-02
## c6.loc100 6 100.000000 1.014481e-01
## MgSTS21 6 100.322546 1.025492e-01
## c6.loc101 6 101.000000 1.089807e-01
## c6.loc102 6 102.000000 1.217176e-01
## c6.loc103 6 103.000000 1.384886e-01
## c6.loc104 6 104.000000 1.593278e-01
## c6.loc105 6 105.000000 1.839822e-01
## c6.loc106 6 106.000000 2.118231e-01
## c6.loc107 6 107.000000 2.417855e-01
## c6.loc108 6 108.000000 2.723745e-01
## c6.loc109 6 109.000000 3.017746e-01
## c6.loc110 6 110.000000 3.280676e-01
## c6.loc111 6 111.000000 3.495200e-01
## c6.loc112 6 112.000000 3.648571e-01
## c6.loc113 6 113.000000 3.734381e-01
## MgSTS58 6 113.873257 3.754101e-01
## c6.loc114 6 114.000000 3.900228e-01
## c6.loc115 6 115.000000 5.134974e-01
## MgSTS431 6 115.971222 6.232671e-01
## c6.loc116 6 116.000000 6.227346e-01
## c6.loc117 6 117.000000 5.890026e-01
## MgSTS229 6 117.885640 5.357472e-01
## c6.loc118 6 118.000000 5.311855e-01
## c6.loc119 6 119.000000 4.847435e-01
## c6.loc120 6 120.000000 4.267350e-01
## c6.loc121 6 121.000000 3.579394e-01
## c6.loc122 6 122.000000 2.816032e-01
## c6.loc123 6 123.000000 2.036676e-01
## c6.loc124 6 124.000000 1.317417e-01
## c6.loc125 6 125.000000 7.290394e-02
## MgSTS220 6 125.892139 3.509332e-02
## c6.loc126 6 126.000000 3.389141e-02
## c6.loc127 6 127.000000 2.728964e-02
## c6.loc128 6 128.000000 3.181435e-02
## c6.loc129 6 129.000000 5.171145e-02
## c6.loc130 6 130.000000 9.074692e-02
## c6.loc131 6 131.000000 1.507520e-01
## c6.loc132 6 132.000000 2.307804e-01
## BD169 6 132.799623 3.070156e-01
## c6.loc133 6 133.000000 3.225415e-01
## c6.loc134 6 134.000000 4.094977e-01
## c6.loc135 6 135.000000 5.095810e-01
## c6.loc136 6 136.000000 6.168453e-01
## c6.loc137 6 137.000000 7.197031e-01
## CC171C 6 137.316182 7.492132e-01
## c6.loc138 6 138.000000 7.671826e-01
## c6.loc139 6 139.000000 7.781085e-01
## c6.loc140 6 140.000000 7.551918e-01
## c6.loc141 6 141.000000 6.808967e-01
## c6.loc142 6 142.000000 5.529435e-01
## c6.loc143 6 143.000000 3.976201e-01
## c6.loc144 6 144.000000 2.558290e-01
## CC381 6 144.953978 1.575645e-01
## c6.loc145 6 145.000000 1.570825e-01
## c6.loc146 6 146.000000 1.399602e-01
## c6.loc147 6 147.000000 1.085710e-01
## c6.loc148 6 148.000000 7.475320e-02
## CA152C 6 148.757348 6.368290e-02
## c6.loc149 6 149.000000 5.882841e-02
## c6.loc150 6 150.000000 3.588042e-02
## c6.loc151 6 151.000000 5.786752e+00
## c6.loc152 6 152.000000 5.753999e+00
## AAT300 6 152.880103 5.599811e+00
## c6.loc153 6 153.000000 5.594846e+00
## c6.loc154 6 154.000000 5.476417e+00
## c6.loc155 6 155.000000 5.213705e+00
## MgSTS22 6 155.485954 4.421793e-01
## c6.loc156 6 156.000000 3.772191e-01
## c6.loc157 6 157.000000 2.514873e-01
## c6.loc158 6 158.000000 1.796363e-01
## c6.loc159 6 159.000000 1.837482e-01
## c6.loc160 6 160.000000 2.348401e-01
## c6.loc161 6 161.000000 2.996256e-01
## c6.loc162 6 162.000000 3.606141e-01
## CA283 6 162.744227 3.997933e-01
## CC262 6 162.996417 4.013515e-01
## c6.loc163 6 163.000000 4.013188e-01
## c6.loc164 6 164.000000 3.994862e-01
## c6.loc165 6 165.000000 4.144795e-01
## c6.loc166 6 166.000000 4.433016e-01
## c6.loc167 6 167.000000 4.759624e-01
## c6.loc168 6 168.000000 5.019432e-01
## MgSTS504A 6 168.784943 5.138459e-01
## c6.loc169 6 169.000000 5.162786e-01
## c6.loc170 6 170.000000 5.218474e-01
## c6.loc171 6 171.000000 5.154674e-01
## c6.loc172 6 172.000000 4.957103e-01
## c6.loc173 6 173.000000 4.640424e-01
## c6.loc174 6 174.000000 4.243338e-01
## c6.loc175 6 175.000000 3.812646e-01
## c6.loc176 6 176.000000 3.388214e-01
## CC126 6 176.646949 3.130722e-01
## c6.loc177 6 177.000000 3.051683e-01
## c6.loc178 6 178.000000 2.818390e-01
## c6.loc179 6 179.000000 2.572135e-01
## c6.loc180 6 180.000000 2.312745e-01
## c6.loc181 6 181.000000 2.043105e-01
## c6.loc182 6 182.000000 1.770246e-01
## c6.loc183 6 183.000000 1.504805e-01
## c6.loc184 6 184.000000 1.258666e-01
## c6.loc185 6 185.000000 1.041854e-01
## c6.loc186 6 186.000000 8.603155e-02
## c6.loc187 6 187.000000 7.154583e-02
## MgSTS426 6 187.500207 6.564635e-02
## c6.loc188 6 188.000000 6.341210e-02
## c6.loc189 6 189.000000 5.739798e-02
## MgSTS28 6 189.795497 5.154729e-02
## c6.loc190 6 190.000000 4.652741e-02
## c6.loc191 6 191.000000 5.332585e-02
## MgSTS453 6 191.171812 5.885118e-02
## MgSTS504B 6 191.546586 7.829204e-02
## c6.loc192 6 192.000000 6.656821e-02
## c6.loc193 6 193.000000 4.286695e-02
## c6.loc194 6 194.000000 2.498423e-02
## c6.loc195 6 195.000000 1.514573e-02
## c6.loc196 6 196.000000 1.400539e-02
## c6.loc197 6 197.000000 2.071314e-02
## MgSTS606 6 197.433185 2.557193e-02
## c6.loc198 6 198.000000 7.580677e-02
## MgSTS323 6 198.829643 2.030061e-01
## c6.loc199 6 199.000000 1.982288e-01
## c6.loc200 6 200.000000 1.624986e-01
## MgSTS467 6 200.380094 1.467648e-01
## MgSTS430 6 200.676038 1.035490e-01
## c6.loc201 6 201.000000 1.258626e-01
## c6.loc202 6 202.000000 2.098124e-01
## c6.loc203 6 203.000000 3.073845e-01
## c6.loc204 6 204.000000 4.053337e-01
## c6.loc205 6 205.000000 4.906393e-01
## MgSTS456 6 205.979038 5.526602e-01
## c6.loc206 6 206.000000 5.517637e-01
## c6.loc207 6 207.000000 5.125256e-01
## c6.loc208 6 208.000000 4.878077e-01
## c6.loc209 6 209.000000 4.847282e-01
## c6.loc210 6 210.000000 5.052528e-01
## BC392 6 210.589085 5.271893e-01
## c6.loc211 6 211.000000 5.311368e-01
## c6.loc212 6 212.000000 5.389710e-01
## c6.loc213 6 213.000000 5.436395e-01
## c6.loc214 6 214.000000 5.445660e-01
## c6.loc215 6 215.000000 5.413875e-01
## c6.loc216 6 216.000000 5.339971e-01
## c6.loc217 6 217.000000 5.225552e-01
## c6.loc218 6 218.000000 5.074649e-01
## c6.loc219 6 219.000000 4.893149e-01
## c6.loc220 6 220.000000 4.687997e-01
## BC243 6 220.502879 4.578702e-01
## c6.loc221 6 221.000000 4.626499e-01
## c6.loc222 6 222.000000 4.742810e-01
## c6.loc223 6 223.000000 4.885717e-01
## c6.loc224 6 224.000000 5.054680e-01
## c6.loc225 6 225.000000 5.247985e-01
## c6.loc226 6 226.000000 5.462567e-01
## c6.loc227 6 227.000000 5.693888e-01
## c6.loc228 6 228.000000 5.935903e-01
## c6.loc229 6 229.000000 6.181166e-01
## c6.loc230 6 230.000000 6.421082e-01
## c6.loc231 6 231.000000 6.646334e-01
## c6.loc232 6 232.000000 6.847472e-01
## c6.loc233 6 233.000000 7.015643e-01
## c6.loc234 6 234.000000 7.143372e-01
## c6.loc235 6 235.000000 7.225322e-01
## c6.loc236 6 236.000000 7.258890e-01
## c6.loc237 6 237.000000 7.244543e-01
## c6.loc238 6 238.000000 7.185807e-01
## c6.loc239 6 239.000000 7.088883e-01
## c6.loc240 6 240.000000 6.961955e-01
## c6.loc241 6 241.000000 6.814296e-01
## c6.loc242 6 242.000000 6.655344e-01
## c6.loc243 6 243.000000 6.493869e-01
## c6.loc244 6 244.000000 6.337371e-01
## c6.loc245 6 245.000000 6.191722e-01
## c6.loc246 6 246.000000 6.061063e-01
## c6.loc247 6 247.000000 5.947908e-01
## c6.loc248 6 248.000000 5.853361e-01
## c6.loc249 6 249.000000 5.777413e-01
## c6.loc250 6 250.000000 5.719240e-01
## c6.loc251 6 251.000000 5.677473e-01
## c6.loc252 6 252.000000 5.650433e-01
## c6.loc253 6 253.000000 5.636306e-01
## c6.loc254 6 254.000000 5.633273e-01
## MgSTS480 6 254.797106 5.637620e-01
## c6.loc255 6 255.000000 5.623353e-01
## c6.loc256 6 256.000000 5.551505e-01
## c6.loc257 6 257.000000 5.478328e-01
## c6.loc258 6 258.000000 5.404294e-01
## c6.loc259 6 259.000000 5.329991e-01
## c6.loc260 6 260.000000 5.256136e-01
## c6.loc261 6 261.000000 5.183591e-01
## c6.loc262 6 262.000000 5.113374e-01
## c6.loc263 6 263.000000 5.046668e-01
## c6.loc264 6 264.000000 4.984824e-01
## c6.loc265 6 265.000000 4.929362e-01
## c6.loc266 6 266.000000 4.881952e-01
## c6.loc267 6 267.000000 4.844391e-01
## c6.loc268 6 268.000000 4.818563e-01
## c6.loc269 6 269.000000 4.806379e-01
## c6.loc270 6 270.000000 4.809699e-01
## c6.loc271 6 271.000000 4.830247e-01
## c6.loc272 6 272.000000 4.869496e-01
## c6.loc273 6 273.000000 4.928565e-01
## c6.loc274 6 274.000000 5.008110e-01
## c6.loc275 6 275.000000 5.108227e-01
## c6.loc276 6 276.000000 5.228396e-01
## c6.loc277 6 277.000000 5.367453e-01
## c6.loc278 6 278.000000 5.523616e-01
## c6.loc279 6 279.000000 5.694557e-01
## c6.loc280 6 280.000000 5.877519e-01
## c6.loc281 6 281.000000 6.069470e-01
## c6.loc282 6 282.000000 6.267268e-01
## c6.loc283 6 283.000000 6.467826e-01
## c6.loc284 6 284.000000 6.668266e-01
## c6.loc285 6 285.000000 6.866029e-01
## c6.loc286 6 286.000000 7.058960e-01
## c6.loc287 6 287.000000 7.245349e-01
## c6.loc288 6 288.000000 7.423931e-01
## c6.loc289 6 289.000000 7.593862e-01
## c6.loc290 6 290.000000 7.754675e-01
## c6.loc291 6 291.000000 7.906215e-01
## c6.loc292 6 292.000000 8.048578e-01
## c6.loc293 6 293.000000 8.182039e-01
## c6.loc294 6 294.000000 8.306993e-01
## c6.loc295 6 295.000000 8.423896e-01
## c6.loc296 6 296.000000 8.533210e-01
## c6.loc297 6 297.000000 8.635363e-01
## c6.loc298 6 298.000000 8.730708e-01
## c6.loc299 6 299.000000 8.819481e-01
## c6.loc300 6 300.000000 8.901776e-01
## BC125 6 300.005891 8.902237e-01
## CB162 7 0.000000 5.008895e-01
## c7.loc1 7 1.000000 4.942327e-01
## c7.loc2 7 2.000000 4.859058e-01
## c7.loc3 7 3.000000 4.804818e-01
## c7.loc4 7 4.000000 4.820234e-01
## c7.loc5 7 5.000000 4.933668e-01
## BA220 7 5.503593 5.030819e-01
## c7.loc6 7 6.000000 5.132187e-01
## c7.loc7 7 7.000000 5.333187e-01
## c7.loc8 7 8.000000 5.523856e-01
## c7.loc9 7 9.000000 5.694913e-01
## c7.loc10 7 10.000000 5.835050e-01
## c7.loc11 7 11.000000 5.934698e-01
## c7.loc12 7 12.000000 5.988883e-01
## BD209 7 12.614372 5.999360e-01
## c7.loc13 7 13.000000 6.001525e-01
## c7.loc14 7 14.000000 5.961454e-01
## c7.loc15 7 15.000000 5.852591e-01
## c7.loc16 7 16.000000 5.680166e-01
## c7.loc17 7 17.000000 5.459305e-01
## c7.loc18 7 18.000000 5.210467e-01
## CC283 7 18.081364 5.189801e-01
## c7.loc19 7 19.000000 4.630643e-01
## c7.loc20 7 20.000000 4.351693e-01
## c7.loc21 7 21.000000 4.486377e-01
## CA258C 7 21.373625 4.633300e-01
## c7.loc22 7 22.000000 5.314103e-01
## c7.loc23 7 23.000000 6.584141e-01
## c7.loc24 7 24.000000 7.987719e-01
## c7.loc25 7 25.000000 9.361431e-01
## c7.loc26 7 26.000000 1.051076e+00
## CB246 7 26.085247 1.059226e+00
## c7.loc27 7 27.000000 1.098062e+00
## c7.loc28 7 28.000000 1.074898e+00
## c7.loc29 7 29.000000 9.518760e-01
## MgSTS574a 7 29.701254 8.139448e-01
## c7.loc30 7 30.000000 7.636095e-01
## c7.loc31 7 31.000000 5.901362e-01
## c7.loc32 7 32.000000 4.221464e-01
## AA280 7 32.611508 3.285641e-01
## CB173 7 32.612021 3.284702e-01
## c7.loc33 7 33.000000 3.196622e-01
## c7.loc34 7 34.000000 2.975682e-01
## c7.loc35 7 35.000000 2.737378e-01
## c7.loc36 7 36.000000 2.076102e+00
## c7.loc37 7 37.000000 2.020811e+00
## c7.loc38 7 38.000000 1.970583e+00
## AA361 7 38.474226 1.948090e+00
## c7.loc39 7 39.000000 1.899379e+00
## c7.loc40 7 40.000000 1.807933e+00
## c7.loc41 7 41.000000 5.509707e-02
## c7.loc42 7 42.000000 8.347673e-02
## c7.loc43 7 43.000000 1.148240e-01
## c7.loc44 7 44.000000 1.466204e-01
## c7.loc45 7 45.000000 1.767789e-01
## c7.loc46 7 46.000000 2.036184e-01
## c7.loc47 7 47.000000 2.258977e-01
## BD170 7 47.612890 2.369105e-01
## c7.loc48 7 48.000000 2.451915e-01
## c7.loc49 7 49.000000 2.677552e-01
## c7.loc50 7 50.000000 2.920239e-01
## c7.loc51 7 51.000000 3.181298e-01
## c7.loc52 7 52.000000 3.459911e-01
## c7.loc53 7 53.000000 3.752661e-01
## c7.loc54 7 54.000000 4.053605e-01
## c7.loc55 7 55.000000 4.354914e-01
## c7.loc56 7 56.000000 4.647936e-01
## c7.loc57 7 57.000000 4.924391e-01
## c7.loc58 7 58.000000 5.177398e-01
## c7.loc59 7 59.000000 5.402128e-01
## c7.loc60 7 60.000000 5.595994e-01
## c7.loc61 7 61.000000 5.758474e-01
## c7.loc62 7 62.000000 5.890694e-01
## c7.loc63 7 63.000000 5.994926e-01
## c7.loc64 7 64.000000 6.074115e-01
## c7.loc65 7 65.000000 6.131482e-01
## MgSTS251 7 65.961893 6.169064e-01
## c7.loc66 7 66.000000 6.202758e-01
## c7.loc67 7 67.000000 7.115248e-01
## c7.loc68 7 68.000000 8.049999e-01
## c7.loc69 7 69.000000 8.982160e-01
## c7.loc70 7 70.000000 9.886903e-01
## c7.loc71 7 71.000000 1.074189e+00
## c7.loc72 7 72.000000 1.152932e+00
## c7.loc73 7 73.000000 1.223705e+00
## c7.loc74 7 74.000000 1.285883e+00
## BC334 7 74.789453 1.328735e+00
## c7.loc75 7 75.000000 1.329091e+00
## c7.loc76 7 76.000000 1.326218e+00
## c7.loc77 7 77.000000 1.315037e+00
## c7.loc78 7 78.000000 1.294809e+00
## c7.loc79 7 79.000000 1.265244e+00
## c7.loc80 7 80.000000 1.226611e+00
## c7.loc81 7 81.000000 1.179773e+00
## c7.loc82 7 82.000000 1.126124e+00
## c7.loc83 7 83.000000 1.067439e+00
## c7.loc84 7 84.000000 1.005668e+00
## c7.loc85 7 85.000000 9.427253e-01
## AAT261 7 85.857502 8.891903e-01
## BA158 8 0.000000 1.260371e-01
## c8.loc1 8 1.000000 1.289003e-01
## c8.loc2 8 2.000000 1.319397e-01
## c8.loc3 8 3.000000 1.350183e-01
## c8.loc4 8 4.000000 1.379959e-01
## c8.loc5 8 5.000000 1.407357e-01
## c8.loc6 8 6.000000 1.431112e-01
## c8.loc7 8 7.000000 1.450148e-01
## c8.loc8 8 8.000000 1.463652e-01
## c8.loc9 8 9.000000 1.471130e-01
## c8.loc10 8 10.000000 1.472442e-01
## c8.loc11 8 11.000000 1.467791e-01
## c8.loc12 8 12.000000 1.457688e-01
## c8.loc13 8 13.000000 1.442874e-01
## c8.loc14 8 14.000000 1.424244e-01
## c8.loc15 8 15.000000 1.402752e-01
## c8.loc16 8 16.000000 1.379335e-01
## c8.loc17 8 17.000000 1.354853e-01
## c8.loc18 8 18.000000 1.330047e-01
## BB167 8 18.124756 1.326988e-01
## c8.loc19 8 19.000000 1.146555e-01
## c8.loc20 8 20.000000 9.983985e-02
## AA246 8 20.481085 9.528924e-02
## c8.loc21 8 21.000000 7.682595e-02
## c8.loc22 8 22.000000 4.557408e-02
## c8.loc23 8 23.000000 2.167247e-02
## c8.loc24 8 24.000000 7.235226e-03
## MgSTS76 8 24.555188 4.055775e-03
## c8.loc25 8 25.000000 5.196022e-03
## c8.loc26 8 26.000000 1.264844e-02
## c8.loc27 8 27.000000 2.741156e-02
## c8.loc28 8 28.000000 5.005481e-02
## c8.loc29 8 29.000000 8.098355e-02
## c8.loc30 8 30.000000 1.203999e-01
## c8.loc31 8 31.000000 1.682723e-01
## c8.loc32 8 32.000000 2.243166e-01
## c8.loc33 8 33.000000 2.879928e-01
## c8.loc34 8 34.000000 3.585181e-01
## c8.loc35 8 35.000000 4.349001e-01
## c8.loc36 8 36.000000 5.159868e-01
## c8.loc37 8 37.000000 6.005303e-01
## c8.loc38 8 38.000000 6.872557e-01
## c8.loc39 8 39.000000 7.749279e-01
## MgSTS621 8 39.531011 8.212560e-01
## c8.loc40 8 40.000000 8.223567e-01
## c8.loc41 8 41.000000 8.206100e-01
## c8.loc42 8 42.000000 8.146866e-01
## c8.loc43 8 43.000000 8.063616e-01
## CC359 8 43.687407 7.999878e-01
## c8.loc44 8 44.000000 8.171053e-01
## c8.loc45 8 45.000000 8.704462e-01
## c8.loc46 8 46.000000 9.199264e-01
## c8.loc47 8 47.000000 9.641981e-01
## c8.loc48 8 48.000000 1.002160e+00
## c8.loc49 8 49.000000 1.033062e+00
## c8.loc50 8 50.000000 1.056562e+00
## c8.loc51 8 51.000000 1.072714e+00
## AAT217 8 51.684621 1.079714e+00
## c8.loc52 8 52.000000 1.079509e+00
## c8.loc53 8 53.000000 1.075348e+00
## c8.loc54 8 54.000000 1.065776e+00
## c8.loc55 8 55.000000 1.050891e+00
## c8.loc56 8 56.000000 1.031004e+00
## c8.loc57 8 57.000000 1.006620e+00
## c8.loc58 8 58.000000 9.783937e-01
## c8.loc59 8 59.000000 9.470863e-01
## c8.loc60 8 60.000000 9.135147e-01
## c8.loc61 8 61.000000 8.785091e-01
## c8.loc62 8 62.000000 8.428740e-01
## c8.loc63 8 63.000000 8.073535e-01
## c8.loc64 8 64.000000 7.726038e-01
## CC338C 8 64.145125 7.677006e-01
## c8.loc65 8 65.000000 8.442076e-01
## c8.loc66 8 66.000000 9.376536e-01
## c8.loc67 8 67.000000 1.033907e+00
## c8.loc68 8 68.000000 1.131272e+00
## c8.loc69 8 69.000000 1.227812e+00
## c8.loc70 8 70.000000 1.321465e+00
## c8.loc71 8 71.000000 1.410203e+00
## c8.loc72 8 72.000000 1.492177e+00
## c8.loc73 8 73.000000 1.565843e+00
## c8.loc74 8 74.000000 1.630044e+00
## c8.loc75 8 75.000000 1.684056e+00
## MgSTS538 8 75.084128 1.688086e+00
## c8.loc76 8 76.000000 1.853767e+00
## c8.loc77 8 77.000000 2.032249e+00
## c8.loc78 8 78.000000 2.198440e+00
## c8.loc79 8 79.000000 2.342768e+00
## MgSTS330 8 79.928051 2.450930e+00
## c8.loc80 8 80.000000 2.457002e+00
## c8.loc81 8 81.000000 2.538292e+00
## c8.loc82 8 82.000000 2.610743e+00
## c8.loc83 8 83.000000 2.672295e+00
## c8.loc84 8 84.000000 2.721115e+00
## MgSTS381 8 84.007733 2.721436e+00
## c8.loc85 8 85.000000 2.837417e+00
## c8.loc86 8 86.000000 2.949914e+00
## c8.loc87 8 87.000000 3.055177e+00
## c8.loc88 8 88.000000 3.150912e+00
## c8.loc89 8 89.000000 3.235596e+00
## c8.loc90 8 90.000000 3.308527e+00
## c8.loc91 8 91.000000 3.369671e+00
## MgSTS31 8 91.177444 3.379237e+00
## c8.loc92 8 92.000000 3.101980e+00
## c8.loc93 8 93.000000 2.735651e+00
## c8.loc94 8 94.000000 2.368672e+00
## CA122 8 94.985080 2.037681e+00
## c8.loc95 8 95.000000 2.038242e+00
## c8.loc96 8 96.000000 2.071621e+00
## c8.loc97 8 97.000000 2.094947e+00
## c8.loc98 8 98.000000 2.106763e+00
## c8.loc99 8 99.000000 2.105970e+00
## c8.loc100 8 100.000000 2.092001e+00
## c8.loc101 8 101.000000 2.064931e+00
## c8.loc102 8 102.000000 2.025518e+00
## c8.loc103 8 103.000000 1.975133e+00
## c8.loc104 8 104.000000 1.915602e+00
## c8.loc105 8 105.000000 1.849003e+00
## c8.loc106 8 106.000000 1.777453e+00
## BD189C 8 106.468881 1.742955e+00
## c8.loc107 8 107.000000 1.679001e+00
## c8.loc108 8 108.000000 1.560581e+00
## c8.loc109 8 109.000000 1.447536e+00
## c8.loc110 8 110.000000 1.342100e+00
## c8.loc111 8 111.000000 1.246122e+00
## c8.loc112 8 112.000000 1.160871e+00
## c8.loc113 8 113.000000 1.086910e+00
## c8.loc114 8 114.000000 1.024058e+00
## c8.loc115 8 115.000000 9.714491e-01
## c8.loc116 8 116.000000 9.276824e-01
## c8.loc117 8 117.000000 8.910271e-01
## c8.loc118 8 118.000000 8.596509e-01
## BA314 8 118.108715 8.565106e-01
## c8.loc119 8 119.000000 8.009836e-01
## c8.loc120 8 120.000000 7.352002e-01
## c8.loc121 8 121.000000 6.674188e-01
## c8.loc122 8 122.000000 5.993450e-01
## c8.loc123 8 123.000000 5.326317e-01
## c8.loc124 8 124.000000 4.688054e-01
## c8.loc125 8 125.000000 4.092060e-01
## c8.loc126 8 126.000000 3.549412e-01
## c8.loc127 8 127.000000 3.068559e-01
## c8.loc128 8 128.000000 2.655172e-01
## c8.loc129 8 129.000000 2.312133e-01
## c8.loc130 8 130.000000 2.039677e-01
## CB272 8 130.339325 1.963424e-01
## c8.loc131 8 131.000000 1.984906e-01
## c8.loc132 8 132.000000 2.001823e-01
## c8.loc133 8 133.000000 2.027095e-01
## c8.loc134 8 134.000000 2.105094e-01
## c8.loc135 8 135.000000 2.285509e-01
## MgSTS571 8 135.932428 2.583401e-01
## c8.loc136 8 136.000000 2.661373e-01
## c8.loc137 8 137.000000 4.095623e-01
## c8.loc138 8 138.000000 5.918073e-01
## MgSTS504C 8 138.345967 6.588193e-01
## c8.loc139 8 139.000000 6.225107e-01
## c8.loc140 8 140.000000 5.440807e-01
## MgSTS590 8 140.577524 4.905008e-01
## c8.loc141 8 141.000000 5.026022e-01
## c8.loc142 8 142.000000 5.324329e-01
## c8.loc143 8 143.000000 5.634550e-01
## c8.loc144 8 144.000000 5.952378e-01
## c8.loc145 8 145.000000 6.271918e-01
## c8.loc146 8 146.000000 6.585637e-01
## c8.loc147 8 147.000000 6.884506e-01
## c8.loc148 8 148.000000 7.158410e-01
## c8.loc149 8 149.000000 7.396835e-01
## c8.loc150 8 150.000000 7.589798e-01
## c8.loc151 8 151.000000 7.728883e-01
## c8.loc152 8 152.000000 7.808213e-01
## c8.loc153 8 153.000000 7.825134e-01
## c8.loc154 8 154.000000 7.780475e-01
## c8.loc155 8 155.000000 7.678321e-01
## c8.loc156 8 156.000000 7.525382e-01
## c8.loc157 8 157.000000 7.330119e-01
## c8.loc158 8 158.000000 7.101798e-01
## c8.loc159 8 159.000000 6.849655e-01
## c8.loc160 8 160.000000 6.582258e-01
## c8.loc161 8 161.000000 6.307098e-01
## AAT296 8 161.347816 6.211376e-01
## c8.loc162 8 162.000000 5.703306e-01
## c8.loc163 8 163.000000 4.836414e-01
## c8.loc164 8 164.000000 3.909925e-01
## c8.loc165 8 165.000000 2.983690e-01
## c8.loc166 8 166.000000 2.120138e-01
## c8.loc167 8 167.000000 1.373659e-01
## c8.loc168 8 168.000000 7.819144e-02
## c8.loc169 8 169.000000 3.615955e-02
## AAT211 8 169.063157 3.410117e-02
## c8.loc170 8 170.000000 2.856059e-02
## c8.loc171 8 171.000000 2.276607e-02
## c8.loc172 8 172.000000 1.752486e-02
## AAT39 8 172.440820 1.549640e-02
## c8.loc173 8 173.000000 1.723276e-02
## c8.loc174 8 174.000000 2.154873e-02
## c8.loc175 8 175.000000 2.704053e-02
## c8.loc176 8 176.000000 3.319612e-02
## c8.loc177 8 177.000000 3.945710e-02
## c8.loc178 8 178.000000 4.528871e-02
## c8.loc179 8 179.000000 5.025061e-02
## c8.loc180 8 180.000000 5.405144e-02
## c8.loc181 8 181.000000 5.657281e-02
## c8.loc182 8 182.000000 5.785836e-02
## CC450 8 182.893525 5.809633e-02
## c8.loc183 8 183.000000 5.796283e-02
## c8.loc184 8 184.000000 5.598311e-02
## c8.loc185 8 185.000000 5.333472e-02
## c8.loc186 8 186.000000 5.103558e-02
## c8.loc187 8 187.000000 5.005286e-02
## c8.loc188 8 188.000000 5.107828e-02
## CA305 8 188.575316 5.270691e-02
## c8.loc189 8 189.000000 4.899273e-02
## c8.loc190 8 190.000000 6.068528e-02
## c8.loc191 8 191.000000 1.033986e-01
## c8.loc192 8 192.000000 1.751433e-01
## MgSTS59 8 192.411809 2.113506e-01
## c8.loc193 8 193.000000 1.820146e-01
## c8.loc194 8 194.000000 1.325005e-01
## c8.loc195 8 195.000000 8.679288e-02
## c8.loc196 8 196.000000 4.850173e-02
## c8.loc197 8 197.000000 2.053614e-02
## c8.loc198 8 198.000000 4.466716e-03
## c8.loc199 8 199.000000 3.059383e-04
## CYCB 8 199.597527 2.988304e-03
## c8.loc200 8 200.000000 1.860381e-03
## c8.loc201 8 201.000000 1.865671e-04
## c8.loc202 8 202.000000 1.483736e-03
## c8.loc203 8 203.000000 7.545412e-03
## c8.loc204 8 204.000000 2.003380e-02
## c8.loc205 8 205.000000 4.015797e-02
## c8.loc206 8 206.000000 6.844433e-02
## c8.loc207 8 207.000000 1.046698e-01
## AAT242 8 207.535270 1.269318e-01
## c8.loc208 8 208.000000 1.245096e-01
## c8.loc209 8 209.000000 1.184870e-01
## c8.loc210 8 210.000000 1.114154e-01
## c8.loc211 8 211.000000 1.033778e-01
## c8.loc212 8 212.000000 9.452904e-02
## c8.loc213 8 213.000000 8.509606e-02
## c8.loc214 8 214.000000 7.536928e-02
## c8.loc215 8 215.000000 6.568307e-02
## c8.loc216 8 216.000000 5.638797e-02
## c8.loc217 8 217.000000 4.781830e-02
## c8.loc218 8 218.000000 4.026053e-02
## c8.loc219 8 219.000000 3.392771e-02
## c8.loc220 8 220.000000 2.894411e-02
## c8.loc221 8 221.000000 2.534173e-02
## c8.loc222 8 222.000000 2.306814e-02
## c8.loc223 8 223.000000 2.200272e-02
## c8.loc224 8 224.000000 2.197766e-02
## MgSTS563 8 224.498765 2.229229e-02
## c8.loc225 8 225.000000 3.988022e-02
## MgSTS69 8 225.146098 4.730766e-02
## c8.loc226 8 226.000000 5.464583e-02
## c8.loc227 8 227.000000 6.644269e-02
## c8.loc228 8 228.000000 8.248959e-02
## c8.loc229 8 229.000000 1.036215e-01
## c8.loc230 8 230.000000 1.305728e-01
## c8.loc231 8 231.000000 1.638752e-01
## c8.loc232 8 232.000000 2.037591e-01
## c8.loc233 8 233.000000 2.500806e-01
## c8.loc234 8 234.000000 3.022916e-01
## c8.loc235 8 235.000000 3.594644e-01
## c8.loc236 8 236.000000 4.203666e-01
## c8.loc237 8 237.000000 4.835757e-01
## c8.loc238 8 238.000000 5.476069e-01
## c8.loc239 8 239.000000 6.110347e-01
## c8.loc240 8 240.000000 6.725901e-01
## MgSTS537 8 240.866778 7.235465e-01
## c8.loc241 8 241.000000 7.154929e-01
## c8.loc242 8 242.000000 6.475271e-01
## c8.loc243 8 243.000000 5.691200e-01
## c8.loc244 8 244.000000 4.843615e-01
## c8.loc245 8 245.000000 3.983843e-01
## CA210 8 245.396574 3.652841e-01
## c8.loc246 8 246.000000 3.373150e-01
## c8.loc247 8 247.000000 2.913534e-01
## c8.loc248 8 248.000000 2.476340e-01
## c8.loc249 8 249.000000 2.074696e-01
## c8.loc250 8 250.000000 1.717128e-01
## c8.loc251 8 251.000000 1.408429e-01
## c8.loc252 8 252.000000 1.150613e-01
## BA372 8 252.006970 1.149012e-01
## BA445 9 0.000000 3.831922e-01
## c9.loc1 9 1.000000 3.738628e-01
## c9.loc2 9 2.000000 3.641847e-01
## c9.loc3 9 3.000000 3.541747e-01
## c9.loc4 9 4.000000 3.438651e-01
## c9.loc5 9 5.000000 3.333055e-01
## c9.loc6 9 6.000000 3.225630e-01
## c9.loc7 9 7.000000 3.117228e-01
## c9.loc8 9 8.000000 3.008862e-01
## c9.loc9 9 9.000000 2.901683e-01
## c9.loc10 9 10.000000 2.796942e-01
## c9.loc11 9 11.000000 2.695949e-01
## c9.loc12 9 12.000000 2.600012e-01
## c9.loc13 9 13.000000 2.510388e-01
## c9.loc14 9 14.000000 2.428229e-01
## c9.loc15 9 15.000000 2.354530e-01
## c9.loc16 9 16.000000 2.290099e-01
## c9.loc17 9 17.000000 2.235523e-01
## c9.loc18 9 18.000000 2.191159e-01
## c9.loc19 9 19.000000 2.157133e-01
## c9.loc20 9 20.000000 2.133349e-01
## c9.loc21 9 21.000000 2.119511e-01
## c9.loc22 9 22.000000 2.115151e-01
## c9.loc23 9 23.000000 2.119659e-01
## AAT222 9 23.657744 2.127084e-01
## c9.loc24 9 24.000000 2.103098e-01
## c9.loc25 9 25.000000 2.025218e-01
## c9.loc26 9 26.000000 1.937105e-01
## c9.loc27 9 27.000000 1.840217e-01
## c9.loc28 9 28.000000 1.736623e-01
## c9.loc29 9 29.000000 1.628884e-01
## c9.loc30 9 30.000000 1.519854e-01
## c9.loc31 9 31.000000 1.412412e-01
## c9.loc32 9 32.000000 1.309196e-01
## c9.loc33 9 33.000000 1.212370e-01
## c9.loc34 9 34.000000 1.123480e-01
## c9.loc35 9 35.000000 1.043417e-01
## c9.loc36 9 36.000000 9.724633e-02
## MgSTS536 9 36.543697 9.378025e-02
## c9.loc37 9 37.000000 1.077773e-01
## c9.loc38 9 38.000000 1.446326e-01
## c9.loc39 9 39.000000 1.905015e-01
## c9.loc40 9 40.000000 2.460846e-01
## c9.loc41 9 41.000000 3.117354e-01
## c9.loc42 9 42.000000 3.873377e-01
## c9.loc43 9 43.000000 4.722078e-01
## c9.loc44 9 44.000000 5.650446e-01
## c9.loc45 9 45.000000 6.639491e-01
## c9.loc46 9 46.000000 7.665241e-01
## c9.loc47 9 47.000000 8.700461e-01
## c9.loc48 9 48.000000 9.716853e-01
## c9.loc49 9 49.000000 1.068738e+00
## c9.loc50 9 50.000000 1.158833e+00
## c9.loc51 9 51.000000 1.240081e+00
## c9.loc52 9 52.000000 1.311155e+00
## c9.loc53 9 53.000000 1.371299e+00
## CA196 9 53.442274 1.394215e+00
## c9.loc54 9 54.000000 1.449454e+00
## c9.loc55 9 55.000000 1.532178e+00
## c9.loc56 9 56.000000 1.592671e+00
## c9.loc57 9 57.000000 1.630506e+00
## c9.loc58 9 58.000000 1.645344e+00
## c9.loc59 9 59.000000 1.636766e+00
## CC540 9 59.594000 1.620489e+00
## c9.loc60 9 60.000000 1.663146e+00
## c9.loc61 9 61.000000 1.760107e+00
## c9.loc62 9 62.000000 1.839641e+00
## c9.loc63 9 63.000000 1.896038e+00
## c9.loc64 9 64.000000 1.925026e+00
## c9.loc65 9 65.000000 1.924511e+00
## c9.loc66 9 66.000000 1.894893e+00
## c9.loc67 9 67.000000 1.838856e+00
## c9.loc68 9 68.000000 1.760771e+00
## CA261 9 68.488694 1.716400e+00
## c9.loc69 9 69.000000 1.619674e+00
## c9.loc70 9 70.000000 1.452546e+00
## c9.loc71 9 71.000000 1.319352e+00
## AA95 9 71.977726 1.217599e+00
## c9.loc72 9 72.000000 1.215928e+00
## c9.loc73 9 73.000000 1.135034e+00
## c9.loc74 9 74.000000 1.044426e+00
## c9.loc75 9 75.000000 9.459089e-01
## c9.loc76 9 76.000000 8.428663e-01
## c9.loc77 9 77.000000 7.397176e-01
## c9.loc78 9 78.000000 6.409887e-01
## BA311 9 78.278370 6.150099e-01
## c9.loc79 9 79.000000 5.797375e-01
## c9.loc80 9 80.000000 5.128885e-01
## c9.loc81 9 81.000000 4.314879e-01
## c9.loc82 9 82.000000 3.435755e-01
## c9.loc83 9 83.000000 2.574422e-01
## MgSTS611 9 83.010077 2.566181e-01
## c9.loc84 9 84.000000 2.957111e-01
## c9.loc85 9 85.000000 3.418114e-01
## c9.loc86 9 86.000000 3.910414e-01
## c9.loc87 9 87.000000 4.395892e-01
## c9.loc88 9 88.000000 4.844029e-01
## LFY 9 88.280035 4.959401e-01
## c9.loc89 9 89.000000 4.960094e-01
## c9.loc90 9 90.000000 4.891931e-01
## c9.loc91 9 91.000000 4.733504e-01
## c9.loc92 9 92.000000 4.483122e-01
## c9.loc93 9 93.000000 4.151023e-01
## c9.loc94 9 94.000000 3.758333e-01
## c9.loc95 9 95.000000 3.332629e-01
## c9.loc96 9 96.000000 2.901958e-01
## MgSTS638 9 96.076261 2.869916e-01
## c9.loc97 9 97.000000 3.112085e-01
## c9.loc98 9 98.000000 3.408923e-01
## c9.loc99 9 99.000000 3.735705e-01
## c9.loc100 9 100.000000 4.083704e-01
## c9.loc101 9 101.000000 4.442008e-01
## c9.loc102 9 102.000000 4.798572e-01
## c9.loc103 9 103.000000 5.141711e-01
## c9.loc104 9 104.000000 5.461733e-01
## c9.loc105 9 105.000000 5.752319e-01
## c9.loc106 9 106.000000 6.011272e-01
## c9.loc107 9 107.000000 6.240433e-01
## c9.loc108 9 108.000000 6.444843e-01
## c9.loc109 9 109.000000 6.631439e-01
## c9.loc110 9 110.000000 6.807678e-01
## c9.loc111 9 111.000000 6.980384e-01
## c9.loc112 9 112.000000 7.155005e-01
## BA396C 9 112.026802 7.159705e-01
## c9.loc113 9 113.000000 7.126868e-01
## c9.loc114 9 114.000000 7.065424e-01
## c9.loc115 9 115.000000 6.969662e-01
## c9.loc116 9 116.000000 6.834047e-01
## c9.loc117 9 117.000000 6.654736e-01
## c9.loc118 9 118.000000 6.430478e-01
## c9.loc119 9 119.000000 6.163212e-01
## c9.loc120 9 120.000000 5.858084e-01
## c9.loc121 9 121.000000 5.522788e-01
## c9.loc122 9 122.000000 5.166433e-01
## BC108 9 122.011777 5.162190e-01
## c9.loc123 9 123.000000 5.335720e-01
## c9.loc124 9 124.000000 5.438399e-01
## c9.loc125 9 125.000000 5.440723e-01
## BB218 9 125.661648 5.383424e-01
## c9.loc126 9 126.000000 5.394777e-01
## c9.loc127 9 127.000000 5.372368e-01
## c9.loc128 9 128.000000 5.256065e-01
## c9.loc129 9 129.000000 5.046830e-01
## c9.loc130 9 130.000000 4.759022e-01
## c9.loc131 9 131.000000 4.414643e-01
## c9.loc132 9 132.000000 4.037436e-01
## MgSTS600 9 132.375176 3.892830e-01
## c9.loc133 9 133.000000 3.669694e-01
## MgSTS470 9 133.631210 3.367029e-01
## c9.loc134 9 134.000000 3.454309e-01
## c9.loc135 9 135.000000 3.692073e-01
## c9.loc136 9 136.000000 3.952445e-01
## c9.loc137 9 137.000000 4.241590e-01
## MgSTS558 9 137.845176 4.493030e-01
## c9.loc138 9 138.000000 4.436312e-01
## c9.loc139 9 139.000000 4.103781e-01
## MgSTS468 9 139.381880 4.009531e-01
## c9.loc140 9 140.000000 4.288222e-01
## c9.loc141 9 141.000000 4.730798e-01
## c9.loc142 9 142.000000 5.145414e-01
## c9.loc143 9 143.000000 5.519142e-01
## c9.loc144 9 144.000000 5.844742e-01
## CA167 9 144.618001 6.019743e-01
## c9.loc145 9 145.000000 5.982632e-01
## CA217 9 145.062788 5.976635e-01
## AAT372 10 0.000000 5.847386e-01
## c10.loc1 10 1.000000 5.700479e-01
## c10.loc2 10 2.000000 5.630702e-01
## c10.loc3 10 3.000000 5.625068e-01
## c10.loc4 10 4.000000 5.661365e-01
## c10.loc5 10 5.000000 5.715498e-01
## c10.loc6 10 6.000000 5.766944e-01
## c10.loc7 10 7.000000 5.801303e-01
## c10.loc8 10 8.000000 5.810401e-01
## MgSTS609 10 8.521260 5.803937e-01
## c10.loc9 10 9.000000 5.960023e-01
## c10.loc10 10 10.000000 6.270652e-01
## c10.loc11 10 11.000000 6.548091e-01
## c10.loc12 10 12.000000 6.780276e-01
## c10.loc13 10 13.000000 6.956508e-01
## c10.loc14 10 14.000000 7.068929e-01
## c10.loc15 10 15.000000 7.113669e-01
## c10.loc16 10 16.000000 7.091306e-01
## MgSTS308 10 16.467936 7.059197e-01
## c10.loc17 10 17.000000 7.133250e-01
## c10.loc18 10 18.000000 7.184871e-01
## c10.loc19 10 19.000000 7.119632e-01
## MgSTS579 10 19.158859 7.099100e-01
## c10.loc20 10 20.000000 7.709961e-01
## c10.loc21 10 21.000000 8.436300e-01
## c10.loc22 10 22.000000 9.135943e-01
## c10.loc23 10 23.000000 9.780035e-01
## c10.loc24 10 24.000000 1.034156e+00
## MgSTS293 10 24.358964 1.051764e+00
## c10.loc25 10 25.000000 1.092847e+00
## c10.loc26 10 26.000000 1.153166e+00
## c10.loc27 10 27.000000 1.206935e+00
## c10.loc28 10 28.000000 1.252487e+00
## c10.loc29 10 29.000000 1.288513e+00
## c10.loc30 10 30.000000 1.314160e+00
## c10.loc31 10 31.000000 1.329070e+00
## BA497 10 31.136678 1.330265e+00
## c10.loc32 10 32.000000 1.322373e+00
## c10.loc33 10 33.000000 1.307899e+00
## c10.loc34 10 34.000000 1.287306e+00
## c10.loc35 10 35.000000 1.260600e+00
## c10.loc36 10 36.000000 1.228316e+00
## c10.loc37 10 37.000000 1.191497e+00
## c10.loc38 10 38.000000 1.151540e+00
## c10.loc39 10 39.000000 1.109962e+00
## CC392 10 39.981286 1.068952e+00
## c10.loc40 10 40.000000 1.068679e+00
## c10.loc41 10 41.000000 1.052506e+00
## c10.loc42 10 42.000000 1.033485e+00
## c10.loc43 10 43.000000 1.011688e+00
## c10.loc44 10 44.000000 9.872224e-01
## c10.loc45 10 45.000000 9.602299e-01
## c10.loc46 10 46.000000 9.308917e-01
## c10.loc47 10 47.000000 8.994277e-01
## c10.loc48 10 48.000000 8.660965e-01
## c10.loc49 10 49.000000 8.311935e-01
## c10.loc50 10 50.000000 7.950449e-01
## c10.loc51 10 51.000000 7.580006e-01
## c10.loc52 10 52.000000 7.204235e-01
## c10.loc53 10 53.000000 6.826776e-01
## c10.loc54 10 54.000000 6.451157e-01
## c10.loc55 10 55.000000 6.080673e-01
## c10.loc56 10 56.000000 5.718297e-01
## CA96 10 56.399386 5.577263e-01
## c10.loc57 10 57.000000 5.414252e-01
## c10.loc58 10 58.000000 5.106386e-01
## c10.loc59 10 59.000000 4.759617e-01
## c10.loc60 10 60.000000 4.387304e-01
## c10.loc61 10 61.000000 4.006243e-01
## MgSTS48 10 61.417826 3.849267e-01
## c10.loc62 10 62.000000 4.010510e-01
## c10.loc63 10 63.000000 4.275256e-01
## c10.loc64 10 64.000000 4.542009e-01
## MgSTS37 10 64.336832 4.636392e-01
## c10.loc65 10 65.000000 4.683601e-01
## c10.loc66 10 66.000000 4.755054e-01
## c10.loc67 10 67.000000 4.826061e-01
## c10.loc68 10 68.000000 4.895980e-01
## c10.loc69 10 69.000000 4.964038e-01
## c10.loc70 10 70.000000 5.029299e-01
## c10.loc71 10 71.000000 5.090632e-01
## c10.loc72 10 72.000000 5.146687e-01
## c10.loc73 10 73.000000 5.195888e-01
## c10.loc74 10 74.000000 5.236466e-01
## c10.loc75 10 75.000000 5.266527e-01
## c10.loc76 10 76.000000 5.284180e-01
## c10.loc77 10 77.000000 5.287687e-01
## c10.loc78 10 78.000000 5.275638e-01
## c10.loc79 10 79.000000 5.247095e-01
## c10.loc80 10 80.000000 5.201708e-01
## c10.loc81 10 81.000000 5.139747e-01
## c10.loc82 10 82.000000 5.062079e-01
## c10.loc83 10 83.000000 4.970079e-01
## c10.loc84 10 84.000000 4.865510e-01
## c10.loc85 10 85.000000 4.750374e-01
## c10.loc86 10 86.000000 4.626786e-01
## BB281 10 86.717156 4.534388e-01
## BC379 10 86.717640 4.534245e-01
## c10.loc87 10 87.000000 4.487399e-01
## c10.loc88 10 88.000000 4.322209e-01
## c10.loc89 10 89.000000 4.160630e-01
## c10.loc90 10 90.000000 4.002182e-01
## c10.loc91 10 91.000000 3.844752e-01
## c10.loc92 10 92.000000 3.685117e-01
## c10.loc93 10 93.000000 3.519957e-01
## c10.loc94 10 94.000000 3.346989e-01
## c10.loc95 10 95.000000 3.165716e-01
## CA289 10 95.755888 3.024341e-01
## c10.loc96 10 96.000000 2.951801e-01
## c10.loc97 10 97.000000 2.637701e-01
## c10.loc98 10 98.000000 2.303674e-01
## c10.loc99 10 99.000000 1.961133e-01
## BC498 10 99.138056 1.914467e-01
## c10.loc100 10 100.000000 2.050027e-01
## c10.loc101 10 101.000000 2.247500e-01
## c10.loc102 10 102.000000 2.453935e-01
## c10.loc103 10 103.000000 2.635747e-01
## AAT374 10 103.023049 2.639406e-01
## c10.loc104 10 104.000000 3.469164e-01
## c10.loc105 10 105.000000 4.375296e-01
## c10.loc106 10 106.000000 5.277742e-01
## BD263 10 106.384981 5.609515e-01
## c10.loc107 10 107.000000 5.925296e-01
## BB122 10 107.983706 6.300258e-01
## c10.loc108 10 108.000000 6.350970e-01
## c10.loc109 10 109.000000 9.286061e-01
## MgSTS388 10 109.725454 1.094199e+00
## c10.loc110 10 110.000000 1.103409e+00
## MgSTS509 10 110.032951 1.104495e+00
## c10.loc111 10 111.000000 1.036420e+00
## MgSTS27 10 111.894934 9.511441e-01
## c10.loc112 10 112.000000 9.558286e-01
## c10.loc113 10 113.000000 9.985892e-01
## c10.loc114 10 114.000000 1.036426e+00
## c10.loc115 10 115.000000 1.068132e+00
## c10.loc116 10 116.000000 1.092674e+00
## c10.loc117 10 117.000000 1.109283e+00
## c10.loc118 10 118.000000 1.117538e+00
## BC131 10 118.741293 1.118231e+00
## c10.loc119 10 119.000000 1.121696e+00
## c10.loc120 10 120.000000 1.134262e+00
## c10.loc121 10 121.000000 1.144884e+00
## c10.loc122 10 122.000000 1.152227e+00
## MgSTS539 10 122.016472 1.152310e+00
## c10.loc123 10 123.000000 9.878202e-01
## c10.loc124 10 124.000000 8.154151e-01
## c10.loc125 10 125.000000 6.458937e-01
## c10.loc126 10 126.000000 4.886784e-01
## c10.loc127 10 127.000000 3.531422e-01
## CA198 10 127.558137 2.900567e-01
## c10.loc128 10 128.000000 3.308779e-01
## c10.loc129 10 129.000000 4.568563e-01
## c10.loc130 10 130.000000 6.176391e-01
## BB259 10 130.603714 7.226339e-01
## c10.loc131 10 131.000000 8.092882e-01
## c10.loc132 10 132.000000 1.036088e+00
## MgSTS351 10 132.863395 1.215710e+00
## c10.loc133 10 133.000000 1.227255e+00
## MgSTS435 10 133.296773 1.257056e+00
## c10.loc134 10 134.000000 1.275822e+00
## c10.loc135 10 135.000000 1.276538e+00
## c10.loc136 10 136.000000 1.250585e+00
## AA100 10 136.873940 1.212513e+00
## c10.loc137 10 137.000000 1.202386e+00
## c10.loc138 10 138.000000 1.113855e+00
## c10.loc139 10 139.000000 1.019000e+00
## MgSTS43 10 139.827277 9.431211e-01
## c10.loc140 10 140.000000 9.525759e-01
## c10.loc141 10 141.000000 1.011793e+00
## c10.loc142 10 142.000000 1.077526e+00
## c10.loc143 10 143.000000 1.148578e+00
## c10.loc144 10 144.000000 1.222792e+00
## c10.loc145 10 145.000000 1.297168e+00
## c10.loc146 10 146.000000 1.368229e+00
## c10.loc147 10 147.000000 1.432560e+00
## c10.loc148 10 148.000000 1.487350e+00
## c10.loc149 10 149.000000 1.530767e+00
## AA153C 10 149.757332 1.555546e+00
## c10.loc150 10 150.000000 1.562975e+00
## c10.loc151 10 151.000000 1.594300e+00
## c10.loc152 10 152.000000 1.626301e+00
## c10.loc153 10 153.000000 1.658808e+00
## c10.loc154 10 154.000000 1.691606e+00
## c10.loc155 10 155.000000 1.724427e+00
## c10.loc156 10 156.000000 1.756943e+00
## c10.loc157 10 157.000000 1.788760e+00
## c10.loc158 10 158.000000 1.819413e+00
## c10.loc159 10 159.000000 1.848365e+00
## c10.loc160 10 160.000000 1.875007e+00
## c10.loc161 10 161.000000 1.898666e+00
## c10.loc162 10 162.000000 1.918614e+00
## c10.loc163 10 163.000000 1.934091e+00
## c10.loc164 10 164.000000 1.944332e+00
## c10.loc165 10 165.000000 1.948602e+00
## c10.loc166 10 166.000000 1.946239e+00
## c10.loc167 10 167.000000 1.936696e+00
## c10.loc168 10 168.000000 1.919592e+00
## c10.loc169 10 169.000000 1.894747e+00
## c10.loc170 10 170.000000 1.862209e+00
## c10.loc171 10 171.000000 1.822268e+00
## c10.loc172 10 172.000000 1.775452e+00
## c10.loc173 10 173.000000 1.722497e+00
## c10.loc174 10 174.000000 1.664309e+00
## c10.loc175 10 175.000000 1.601908e+00
## c10.loc176 10 176.000000 1.536375e+00
## c10.loc177 10 177.000000 1.468792e+00
## c10.loc178 10 178.000000 1.400190e+00
## c10.loc179 10 179.000000 1.331514e+00
## c10.loc180 10 180.000000 1.263595e+00
## c10.loc181 10 181.000000 1.197130e+00
## c10.loc182 10 182.000000 1.132685e+00
## c10.loc183 10 183.000000 1.070693e+00
## c10.loc184 10 184.000000 1.011468e+00
## c10.loc185 10 185.000000 9.552153e-01
## c10.loc186 10 186.000000 9.020515e-01
## c10.loc187 10 187.000000 8.520172e-01
## c10.loc188 10 188.000000 8.050936e-01
## c10.loc189 10 189.000000 7.612154e-01
## c10.loc190 10 190.000000 7.202834e-01
## MgSTS75 10 190.951541 6.839743e-01
## c10.loc191 10 191.000000 6.887150e-01
## c10.loc192 10 192.000000 7.875261e-01
## c10.loc193 10 193.000000 8.825454e-01
## c10.loc194 10 194.000000 9.677812e-01
## c10.loc195 10 195.000000 1.037367e+00
## c10.loc196 10 196.000000 1.086263e+00
## c10.loc197 10 197.000000 1.110970e+00
## c10.loc198 10 198.000000 1.110099e+00
## MgSTS70 10 198.018805 1.109843e+00
## c10.loc199 10 199.000000 1.020544e+00
## c10.loc200 10 200.000000 9.137320e-01
## c10.loc201 10 201.000000 7.965739e-01
## c10.loc202 10 202.000000 6.759873e-01
## MgSTS332 10 202.825113 5.790436e-01
## c10.loc203 10 203.000000 5.677606e-01
## c10.loc204 10 204.000000 5.008147e-01
## c10.loc205 10 205.000000 4.322002e-01
## c10.loc206 10 206.000000 3.641923e-01
## c10.loc207 10 207.000000 2.990749e-01
## c10.loc208 10 208.000000 2.389445e-01
## c10.loc209 10 209.000000 1.854858e-01
## MgSTS511 10 209.247780 1.734992e-01
## c10.loc210 10 210.000000 1.670875e-01
## c10.loc211 10 211.000000 1.610779e-01
## c10.loc212 10 212.000000 1.586556e-01
## c10.loc213 10 213.000000 1.603065e-01
## c10.loc214 10 214.000000 1.661255e-01
## c10.loc215 10 215.000000 1.757885e-01
## c10.loc216 10 216.000000 1.886217e-01
## CB309 10 216.420391 1.947172e-01
## MgSTS598 11 0.000000 5.326910e-01
## c11.loc1 11 1.000000 5.392998e-01
## c11.loc2 11 2.000000 5.378528e-01
## c11.loc3 11 3.000000 5.280874e-01
## c11.loc4 11 4.000000 5.106882e-01
## MgSTS26 11 4.260334 5.051180e-01
## c11.loc5 11 5.000000 5.107736e-01
## MgSTS93 11 5.470128 5.094059e-01
## c11.loc6 11 6.000000 5.454933e-01
## c11.loc7 11 7.000000 6.124209e-01
## c11.loc8 11 8.000000 6.740422e-01
## c11.loc9 11 9.000000 7.267562e-01
## c11.loc10 11 10.000000 7.681098e-01
## c11.loc11 11 11.000000 7.972693e-01
## c11.loc12 11 12.000000 8.149467e-01
## MgSTS91 11 12.162428 8.168293e-01
## c11.loc13 11 13.000000 7.150469e-01
## c11.loc14 11 14.000000 5.945331e-01
## c11.loc15 11 15.000000 4.796870e-01
## MgSTS20 11 15.429581 4.333816e-01
## c11.loc16 11 16.000000 4.520720e-01
## c11.loc17 11 17.000000 4.821516e-01
## c11.loc18 11 18.000000 5.066138e-01
## c11.loc19 11 19.000000 5.233893e-01
## c11.loc20 11 20.000000 5.311891e-01
## c11.loc21 11 21.000000 5.298186e-01
## BC526C 11 21.871990 5.218151e-01
## c11.loc22 11 22.000000 5.211807e-01
## c11.loc23 11 23.000000 5.139135e-01
## c11.loc24 11 24.000000 5.025874e-01
## c11.loc25 11 25.000000 4.873640e-01
## c11.loc26 11 26.000000 4.687090e-01
## c11.loc27 11 27.000000 4.474085e-01
## c11.loc28 11 28.000000 4.245390e-01
## c11.loc29 11 29.000000 4.013850e-01
## c11.loc30 11 30.000000 3.793108e-01
## c11.loc31 11 31.000000 3.596023e-01
## c11.loc32 11 32.000000 3.433084e-01
## c11.loc33 11 33.000000 3.311126e-01
## BA387 11 33.415683 3.273424e-01
## c11.loc34 11 34.000000 3.214713e-01
## c11.loc35 11 35.000000 3.135021e-01
## c11.loc36 11 36.000000 3.095401e-01
## c11.loc37 11 37.000000 3.111073e-01
## c11.loc38 11 38.000000 3.195189e-01
## c11.loc39 11 39.000000 3.356869e-01
## c11.loc40 11 40.000000 3.599753e-01
## c11.loc41 11 41.000000 3.921500e-01
## c11.loc42 11 42.000000 4.314322e-01
## c11.loc43 11 43.000000 4.766368e-01
## MgSTS19 11 43.201900 4.862807e-01
## c11.loc44 11 44.000000 5.671945e-01
## c11.loc45 11 45.000000 6.797448e-01
## c11.loc46 11 46.000000 8.006725e-01
## c11.loc47 11 47.000000 9.237604e-01
## c11.loc48 11 48.000000 1.041301e+00
## c11.loc49 11 49.000000 1.145501e+00
## c11.loc50 11 50.000000 1.230087e+00
## c11.loc51 11 51.000000 1.291464e+00
## MgSTS620 11 51.110722 1.296734e+00
## c11.loc52 11 52.000000 1.289168e+00
## c11.loc53 11 53.000000 1.271828e+00
## c11.loc54 11 54.000000 1.244512e+00
## c11.loc55 11 55.000000 1.207172e+00
## c11.loc56 11 56.000000 1.160483e+00
## c11.loc57 11 57.000000 1.105833e+00
## c11.loc58 11 58.000000 1.045178e+00
## c11.loc59 11 59.000000 9.808205e-01
## c11.loc60 11 60.000000 9.151254e-01
## c11.loc61 11 61.000000 8.502773e-01
## c11.loc62 11 62.000000 7.880978e-01
## MgSTS360 11 62.390090 7.650146e-01
## c11.loc63 11 63.000000 7.123852e-01
## c11.loc64 11 64.000000 6.203673e-01
## c11.loc65 11 65.000000 5.245842e-01
## c11.loc66 11 66.000000 4.293382e-01
## c11.loc67 11 67.000000 3.393022e-01
## c11.loc68 11 68.000000 2.586394e-01
## c11.loc69 11 69.000000 1.902619e-01
## MgSTS380 11 69.191511 1.787968e-01
## c11.loc70 11 70.000000 1.827859e-01
## c11.loc71 11 71.000000 1.889144e-01
## c11.loc72 11 72.000000 1.962241e-01
## c11.loc73 11 73.000000 2.044577e-01
## c11.loc74 11 74.000000 2.132395e-01
## c11.loc75 11 75.000000 2.221010e-01
## c11.loc76 11 76.000000 2.305275e-01
## c11.loc77 11 77.000000 2.380161e-01
## c11.loc78 11 78.000000 2.441308e-01
## c11.loc79 11 79.000000 2.485433e-01
## c11.loc80 11 80.000000 2.510512e-01
## CC53 11 80.729628 2.516300e-01
## c11.loc81 11 81.000000 2.649211e-01
## c11.loc82 11 82.000000 3.183553e-01
## c11.loc83 11 83.000000 3.747361e-01
## c11.loc84 11 84.000000 4.290081e-01
## c11.loc85 11 85.000000 4.762169e-01
## BD340 11 85.563120 4.981567e-01
## c11.loc86 11 86.000000 5.478975e-01
## c11.loc87 11 87.000000 6.742138e-01
## c11.loc88 11 88.000000 8.115471e-01
## c11.loc89 11 89.000000 9.509947e-01
## c11.loc90 11 90.000000 1.082747e+00
## c11.loc91 11 91.000000 1.198339e+00
## BD100 11 91.900452 1.284021e+00
## c11.loc92 11 92.000000 1.291460e+00
## c11.loc93 11 93.000000 1.356820e+00
## c11.loc94 11 94.000000 1.399900e+00
## c11.loc95 11 95.000000 1.416892e+00
## c11.loc96 11 96.000000 1.407142e+00
## c11.loc97 11 97.000000 1.373215e+00
## c11.loc98 11 98.000000 1.320055e+00
## MgSTS87 11 98.112283 1.313209e+00
## c11.loc99 11 99.000000 1.167168e+00
## CYCA 11 99.500975 1.064201e+00
## c11.loc100 11 100.000000 1.428740e+00
## c11.loc101 11 101.000000 2.259785e+00
## MgSTS471 11 101.045365 2.289253e+00
## c11.loc102 11 102.000000 2.467916e+00
## c11.loc103 11 103.000000 2.629574e+00
## c11.loc104 11 104.000000 2.749517e+00
## c11.loc105 11 105.000000 2.814562e+00
## c11.loc106 11 106.000000 2.817497e+00
## c11.loc107 11 107.000000 2.759025e+00
## BA196 11 107.915465 2.658645e+00
## BB124 11 107.915475 2.658644e+00
## c11.loc108 11 108.000000 2.649386e+00
## c11.loc109 11 109.000000 2.528301e+00
## c11.loc110 11 110.000000 2.391166e+00
## c11.loc111 11 111.000000 2.242511e+00
## c11.loc112 11 112.000000 2.086292e+00
## c11.loc113 11 113.000000 1.925919e+00
## c11.loc114 11 114.000000 1.764374e+00
## c11.loc115 11 115.000000 1.604321e+00
## c11.loc116 11 116.000000 1.448175e+00
## c11.loc117 11 117.000000 1.298103e+00
## c11.loc118 11 118.000000 1.155985e+00
## c11.loc119 11 119.000000 1.023347e+00
## c11.loc120 11 120.000000 9.013100e-01
## c11.loc121 11 121.000000 7.905594e-01
## c11.loc122 11 122.000000 6.913527e-01
## c11.loc123 11 123.000000 6.035629e-01
## MgSTS358 11 123.082997 5.968390e-01
## c11.loc124 11 124.000000 6.197706e-01
## c11.loc125 11 125.000000 6.471415e-01
## c11.loc126 11 126.000000 6.770543e-01
## c11.loc127 11 127.000000 7.095757e-01
## c11.loc128 11 128.000000 7.447290e-01
## c11.loc129 11 129.000000 7.824820e-01
## c11.loc130 11 130.000000 8.227348e-01
## c11.loc131 11 131.000000 8.653076e-01
## c11.loc132 11 132.000000 9.099304e-01
## c11.loc133 11 133.000000 9.562354e-01
## c11.loc134 11 134.000000 1.003753e+00
## c11.loc135 11 135.000000 1.051917e+00
## c11.loc136 11 136.000000 1.100070e+00
## c11.loc137 11 137.000000 1.147487e+00
## c11.loc138 11 138.000000 1.193397e+00
## c11.loc139 11 139.000000 1.237020e+00
## c11.loc140 11 140.000000 1.277598e+00
## c11.loc141 11 141.000000 1.314439e+00
## c11.loc142 11 142.000000 1.346951e+00
## c11.loc143 11 143.000000 1.374666e+00
## c11.loc144 11 144.000000 1.397266e+00
## c11.loc145 11 145.000000 1.414589e+00
## c11.loc146 11 146.000000 1.426627e+00
## c11.loc147 11 147.000000 1.433510e+00
## c11.loc148 11 148.000000 1.435490e+00
## c11.loc149 11 149.000000 1.432915e+00
## c11.loc150 11 150.000000 1.426198e+00
## c11.loc151 11 151.000000 1.415792e+00
## c11.loc152 11 152.000000 1.402170e+00
## c11.loc153 11 153.000000 1.385796e+00
## c11.loc154 11 154.000000 1.367117e+00
## BC374 11 154.036519 1.366404e+00
## BB182 12 0.000000 7.833465e-01
## c12.loc1 12 1.000000 7.375001e-01
## c12.loc2 12 2.000000 6.885143e-01
## c12.loc3 12 3.000000 6.371526e-01
## CA140 12 3.406237 6.159659e-01
## c12.loc4 12 4.000000 5.870725e-01
## c12.loc5 12 5.000000 5.399587e-01
## c12.loc6 12 6.000000 4.953072e-01
## c12.loc7 12 7.000000 4.532044e-01
## c12.loc8 12 8.000000 4.135071e-01
## c12.loc9 12 9.000000 3.759061e-01
## c12.loc10 12 10.000000 3.400180e-01
## c12.loc11 12 11.000000 3.054841e-01
## c12.loc12 12 12.000000 2.720541e-01
## c12.loc13 12 13.000000 2.396364e-01
## CB216 12 13.606611 2.205650e-01
## c12.loc14 12 14.000000 2.101024e-01
## c12.loc15 12 15.000000 1.836191e-01
## c12.loc16 12 16.000000 1.579525e-01
## CA378 12 16.315098 1.501724e-01
## c12.loc17 12 17.000000 1.436316e-01
## c12.loc18 12 18.000000 1.308174e-01
## c12.loc19 12 19.000000 1.148064e-01
## c12.loc20 12 20.000000 9.745309e-02
## c12.loc21 12 21.000000 8.146979e-02
## MgSTS36 12 21.545915 7.438736e-02
## c12.loc22 12 22.000000 5.285849e-02
## MgSTS438 12 22.161045 4.600664e-02
## c12.loc23 12 23.000000 4.740730e-02
## c12.loc24 12 24.000000 4.921351e-02
## c12.loc25 12 25.000000 5.132886e-02
## c12.loc26 12 26.000000 5.397277e-02
## c12.loc27 12 27.000000 5.738976e-02
## c12.loc28 12 28.000000 6.181244e-02
## c12.loc29 12 29.000000 6.741936e-02
## c12.loc30 12 30.000000 7.429884e-02
## c12.loc31 12 31.000000 8.242964e-02
## c12.loc32 12 32.000000 9.168348e-02
## c12.loc33 12 33.000000 1.018476e-01
## c12.loc34 12 34.000000 1.126587e-01
## AA296C 12 34.441203 1.175349e-01
## c12.loc35 12 35.000000 1.050345e-01
## c12.loc36 12 36.000000 8.528491e-02
## c12.loc37 12 37.000000 7.013017e-02
## c12.loc38 12 38.000000 6.058331e-02
## c12.loc39 12 39.000000 5.714557e-02
## c12.loc40 12 40.000000 5.970823e-02
## BC330 12 40.470117 6.279870e-02
## c12.loc41 12 41.000000 5.802183e-02
## c12.loc42 12 42.000000 5.029626e-02
## c12.loc43 12 43.000000 4.449580e-02
## c12.loc44 12 44.000000 4.080237e-02
## c12.loc45 12 45.000000 3.934566e-02
## c12.loc46 12 46.000000 4.018905e-02
## c12.loc47 12 47.000000 4.331920e-02
## c12.loc48 12 48.000000 4.864189e-02
## c12.loc49 12 49.000000 5.598572e-02
## c12.loc50 12 50.000000 6.511464e-02
## c12.loc51 12 51.000000 7.574795e-02
## c12.loc52 12 52.000000 8.758478e-02
## c12.loc53 12 53.000000 1.003288e-01
## c12.loc54 12 54.000000 1.137087e-01
## c12.loc55 12 55.000000 1.274919e-01
## c12.loc56 12 56.000000 1.414891e-01
## c12.loc57 12 57.000000 1.555533e-01
## c12.loc58 12 58.000000 1.695725e-01
## c12.loc59 12 59.000000 1.834621e-01
## c12.loc60 12 60.000000 1.971565e-01
## c12.loc61 12 61.000000 2.106026e-01
## c12.loc62 12 62.000000 2.237557e-01
## c12.loc63 12 63.000000 2.365769e-01
## MgSTS282B 12 63.422599 2.418568e-01
## c12.loc64 12 64.000000 2.398773e-01
## c12.loc65 12 65.000000 2.372757e-01
## c12.loc66 12 66.000000 2.356347e-01
## c12.loc67 12 67.000000 2.348136e-01
## c12.loc68 12 68.000000 2.346697e-01
## c12.loc69 12 69.000000 2.350794e-01
## c12.loc70 12 70.000000 2.359600e-01
## MgSTS548 12 70.777701 2.369484e-01
## c12.loc71 12 71.000000 2.333331e-01
## c12.loc72 12 72.000000 2.156020e-01
## c12.loc73 12 73.000000 1.951977e-01
## c12.loc74 12 74.000000 1.717855e-01
## c12.loc75 12 75.000000 1.455765e-01
## c12.loc76 12 76.000000 1.174304e-01
## c12.loc77 12 77.000000 8.880608e-02
## c12.loc78 12 78.000000 6.155703e-02
## c12.loc79 12 79.000000 3.762658e-02
## c12.loc80 12 80.000000 1.872821e-02
## AAT364 12 80.349896 1.358112e-02
## c12.loc81 12 81.000000 1.902189e-02
## c12.loc82 12 82.000000 3.300298e-02
## c12.loc83 12 83.000000 5.404563e-02
## c12.loc84 12 84.000000 8.226972e-02
## c12.loc85 12 85.000000 1.175361e-01
## c12.loc86 12 86.000000 1.594269e-01
## c12.loc87 12 87.000000 2.072522e-01
## c12.loc88 12 88.000000 2.600836e-01
## c12.loc89 12 89.000000 3.168144e-01
## c12.loc90 12 90.000000 3.762344e-01
## c12.loc91 12 91.000000 4.371121e-01
## c12.loc92 12 92.000000 4.982706e-01
## MgSTS113 12 92.173785 5.087606e-01
## c12.loc93 12 93.000000 4.782236e-01
## c12.loc94 12 94.000000 4.376922e-01
## c12.loc95 12 95.000000 3.943174e-01
## c12.loc96 12 96.000000 3.494162e-01
## AAT308 12 96.988834 3.049076e-01
## c12.loc97 12 97.000000 3.049423e-01
## c12.loc98 12 98.000000 3.078231e-01
## c12.loc99 12 99.000000 3.101154e-01
## c12.loc100 12 100.000000 3.116950e-01
## c12.loc101 12 101.000000 3.124380e-01
## c12.loc102 12 102.000000 3.122274e-01
## c12.loc103 12 103.000000 3.109613e-01
## c12.loc104 12 104.000000 3.085618e-01
## c12.loc105 12 105.000000 3.049818e-01
## c12.loc106 12 106.000000 3.002112e-01
## c12.loc107 12 107.000000 2.942786e-01
## c12.loc108 12 108.000000 2.872497e-01
## c12.loc109 12 109.000000 2.792226e-01
## c12.loc110 12 110.000000 2.703202e-01
## c12.loc111 12 111.000000 2.606820e-01
## c12.loc112 12 112.000000 2.504571e-01
## c12.loc113 12 113.000000 2.397975e-01
## c12.loc114 12 114.000000 2.288544e-01
## c12.loc115 12 115.000000 2.177746e-01
## BD143 12 115.890223 2.079210e-01
## c12.loc116 12 116.000000 2.065811e-01
## c12.loc117 12 117.000000 1.936814e-01
## c12.loc118 12 118.000000 1.797200e-01
## c12.loc119 12 119.000000 1.647726e-01
## c12.loc120 12 120.000000 1.490393e-01
## c12.loc121 12 121.000000 1.328413e-01
## c12.loc122 12 122.000000 1.165905e-01
## c12.loc123 12 123.000000 1.007387e-01
## c12.loc124 12 124.000000 8.571691e-02
## c12.loc125 12 125.000000 7.188343e-02
## MgSTS437 12 125.976156 5.977096e-02
## c12.loc126 12 126.000000 6.040342e-02
## c12.loc127 12 127.000000 9.040828e-02
## c12.loc128 12 128.000000 1.263165e-01
## c12.loc129 12 129.000000 1.671693e-01
## c12.loc130 12 130.000000 2.115796e-01
## c12.loc131 12 131.000000 2.578956e-01
## c12.loc132 12 132.000000 3.044123e-01
## BA245C 12 132.482277 3.263361e-01
## c12.loc133 12 133.000000 3.329365e-01
## c12.loc134 12 134.000000 3.487014e-01
## c12.loc135 12 135.000000 3.669471e-01
## c12.loc136 12 136.000000 3.851585e-01
## BA175 12 136.220001 3.888613e-01
## c12.loc137 12 137.000000 3.685309e-01
## c12.loc138 12 138.000000 3.425153e-01
## c12.loc139 12 139.000000 3.172380e-01
## BC376 12 139.074480 3.154169e-01
## c12.loc140 12 140.000000 3.004717e-01
## c12.loc141 12 141.000000 2.835771e-01
## c12.loc142 12 142.000000 2.660913e-01
## c12.loc143 12 143.000000 2.482249e-01
## c12.loc144 12 144.000000 2.302011e-01
## c12.loc145 12 145.000000 2.122307e-01
## c12.loc146 12 146.000000 1.944913e-01
## BC546 12 146.782045 1.808972e-01
## c12.loc147 12 147.000000 1.790256e-01
## c12.loc148 12 148.000000 1.714698e-01
## c12.loc149 12 149.000000 1.657431e-01
## c12.loc150 12 150.000000 1.617724e-01
## c12.loc151 12 151.000000 1.593451e-01
## c12.loc152 12 152.000000 1.581228e-01
## c12.loc153 12 153.000000 1.576852e-01
## c12.loc154 12 154.000000 1.575943e-01
## BD242 12 154.835825 1.574996e-01
## c12.loc155 12 155.000000 1.572692e-01
## c12.loc156 12 156.000000 1.560945e-01
## c12.loc157 12 157.000000 1.550127e-01
## c12.loc158 12 158.000000 1.535411e-01
## c12.loc159 12 159.000000 1.512077e-01
## c12.loc160 12 160.000000 1.476344e-01
## c12.loc161 12 161.000000 1.426107e-01
## c12.loc162 12 162.000000 1.361343e-01
## c12.loc163 12 163.000000 1.284041e-01
## c12.loc164 12 164.000000 1.197686e-01
## c12.loc165 12 165.000000 1.106517e-01
## BC388 12 165.139049 1.093813e-01
## c12.loc166 12 166.000000 1.123694e-01
## c12.loc167 12 167.000000 1.160895e-01
## c12.loc168 12 168.000000 1.198189e-01
## c12.loc169 12 169.000000 1.232718e-01
## c12.loc170 12 170.000000 1.261629e-01
## c12.loc171 12 171.000000 1.282353e-01
## c12.loc172 12 172.000000 1.292917e-01
## c12.loc173 12 173.000000 1.292227e-01
## c12.loc174 12 174.000000 1.280262e-01
## c12.loc175 12 175.000000 1.258117e-01
## c12.loc176 12 176.000000 1.227859e-01
## c12.loc177 12 177.000000 1.192248e-01
## c12.loc178 12 178.000000 1.154348e-01
## c12.loc179 12 179.000000 1.117147e-01
## c12.loc180 12 180.000000 1.083242e-01
## c12.loc181 12 181.000000 1.054631e-01
## c12.loc182 12 182.000000 1.032632e-01
## c12.loc183 12 183.000000 1.017904e-01
## c12.loc184 12 184.000000 1.010533e-01
## BC266 12 184.850027 1.009789e-01
## c12.loc185 12 185.000000 1.029633e-01
## c12.loc186 12 186.000000 1.166237e-01
## c12.loc187 12 187.000000 1.305811e-01
## c12.loc188 12 188.000000 1.444458e-01
## c12.loc189 12 189.000000 1.578242e-01
## c12.loc190 12 190.000000 1.703697e-01
## c12.loc191 12 191.000000 1.818219e-01
## c12.loc192 12 192.000000 1.920233e-01
## CC114 12 192.509318 1.966934e-01
## c12.loc193 12 193.000000 1.964568e-01
## c12.loc194 12 194.000000 1.970967e-01
## c12.loc195 12 195.000000 1.992515e-01
## c12.loc196 12 196.000000 2.029048e-01
## c12.loc197 12 197.000000 2.080207e-01
## c12.loc198 12 198.000000 2.145467e-01
## c12.loc199 12 199.000000 2.224169e-01
## c12.loc200 12 200.000000 2.315551e-01
## c12.loc201 12 201.000000 2.418770e-01
## c12.loc202 12 202.000000 2.532927e-01
## c12.loc203 12 203.000000 2.657093e-01
## c12.loc204 12 204.000000 2.790327e-01
## c12.loc205 12 205.000000 2.931695e-01
## BA394 12 205.113998 2.948140e-01
## CB172 13 0.000000 5.678224e-02
## c13.loc1 13 1.000000 5.872618e-02
## c13.loc2 13 2.000000 6.110803e-02
## c13.loc3 13 3.000000 6.388865e-02
## c13.loc4 13 4.000000 6.702552e-02
## c13.loc5 13 5.000000 7.047145e-02
## c13.loc6 13 6.000000 7.417357e-02
## c13.loc7 13 7.000000 7.807279e-02
## c13.loc8 13 8.000000 8.210376e-02
## c13.loc9 13 9.000000 8.619548e-02
## c13.loc10 13 10.000000 9.027279e-02
## c13.loc11 13 11.000000 9.425856e-02
## c13.loc12 13 12.000000 9.807661e-02
## c13.loc13 13 13.000000 1.016551e-01
## c13.loc14 13 14.000000 1.049303e-01
## c13.loc15 13 15.000000 1.078495e-01
## c13.loc16 13 16.000000 1.103744e-01
## c13.loc17 13 17.000000 1.124824e-01
## c13.loc18 13 18.000000 1.141674e-01
## c13.loc19 13 19.000000 1.154389e-01
## c13.loc20 13 20.000000 1.163203e-01
## c13.loc21 13 21.000000 1.168465e-01
## c13.loc22 13 22.000000 1.170600e-01
## c13.loc23 13 23.000000 1.170086e-01
## c13.loc24 13 24.000000 1.167416e-01
## BA334 13 24.422639 1.165767e-01
## c13.loc25 13 25.000000 1.254890e-01
## c13.loc26 13 26.000000 1.450391e-01
## c13.loc27 13 27.000000 1.671285e-01
## c13.loc28 13 28.000000 1.891598e-01
## c13.loc29 13 29.000000 2.093009e-01
## MgSTS577 13 29.920751 2.253024e-01
## c13.loc30 13 30.000000 2.215267e-01
## MgSTS45 13 30.842164 1.745494e-01
## c13.loc31 13 31.000000 1.593373e-01
## c13.loc32 13 32.000000 8.972102e-02
## MgSTS11 13 32.075691 8.671321e-02
## c13.loc33 13 33.000000 7.746506e-02
## c13.loc34 13 34.000000 6.951539e-02
## c13.loc35 13 35.000000 6.422161e-02
## c13.loc36 13 36.000000 6.211568e-02
## c13.loc37 13 37.000000 6.366133e-02
## c13.loc38 13 38.000000 6.913770e-02
## c13.loc39 13 39.000000 7.852141e-02
## c13.loc40 13 40.000000 9.144463e-02
## c13.loc41 13 41.000000 1.073083e-01
## c13.loc42 13 42.000000 1.255493e-01
## MgSTS104 13 42.289204 1.311813e-01
## c13.loc43 13 43.000000 1.472576e-01
## c13.loc44 13 44.000000 1.734520e-01
## c13.loc45 13 45.000000 2.045167e-01
## c13.loc46 13 46.000000 2.414521e-01
## c13.loc47 13 47.000000 2.853405e-01
## c13.loc48 13 48.000000 3.369992e-01
## c13.loc49 13 49.000000 3.960464e-01
## c13.loc50 13 50.000000 4.589730e-01
## c13.loc51 13 51.000000 1.015456e+01
## c13.loc52 13 52.000000 1.023067e+01
## c13.loc53 13 53.000000 1.021965e+01
## MgSTS599 13 53.397071 1.019322e+01
## c13.loc54 13 54.000000 1.047255e+01
## c13.loc55 13 55.000000 1.079156e+01
## c13.loc56 13 56.000000 1.260157e+00
## c13.loc57 13 57.000000 1.273688e+00
## c13.loc58 13 58.000000 1.220210e+00
## CB55 13 58.352924 1.189774e+00
## c13.loc59 13 59.000000 1.115863e+00
## c13.loc60 13 60.000000 9.982933e-01
## c13.loc61 13 61.000000 8.821832e-01
## c13.loc62 13 62.000000 7.718536e-01
## c13.loc63 13 63.000000 6.700332e-01
## c13.loc64 13 64.000000 5.780975e-01
## c13.loc65 13 65.000000 4.964364e-01
## c13.loc66 13 66.000000 4.248040e-01
## c13.loc67 13 67.000000 3.625862e-01
## c13.loc68 13 68.000000 3.089805e-01
## c13.loc69 13 69.000000 2.631060e-01
## c13.loc70 13 70.000000 2.240708e-01
## c13.loc71 13 71.000000 1.910125e-01
## c13.loc72 13 72.000000 1.631243e-01
## c13.loc73 13 73.000000 1.396701e-01
## c13.loc74 13 74.000000 1.199924e-01
## c13.loc75 13 75.000000 1.035137e-01
## c13.loc76 13 76.000000 8.973419e-02
## c13.loc77 13 77.000000 7.822599e-02
## c13.loc78 13 78.000000 6.862576e-02
## c13.loc79 13 79.000000 6.062704e-02
## c13.loc80 13 80.000000 5.397243e-02
## c13.loc81 13 81.000000 4.844637e-02
## CA315 13 81.830484 4.459066e-02
## c13.loc82 13 82.000000 4.162778e-02
## c13.loc83 13 83.000000 2.553058e-02
## c13.loc84 13 84.000000 1.351646e-02
## c13.loc85 13 85.000000 7.335982e-03
## c13.loc86 13 86.000000 8.051149e-03
## c13.loc87 13 87.000000 1.583173e-02
## CB263C 13 87.980797 2.969052e-02
## c13.loc88 13 88.000000 2.940420e-02
## c13.loc89 13 89.000000 1.661982e-02
## c13.loc90 13 90.000000 9.513582e-03
## c13.loc91 13 91.000000 9.961597e-03
## c13.loc92 13 92.000000 1.975102e-02
## c13.loc93 13 93.000000 4.030849e-02
## c13.loc94 13 94.000000 7.240773e-02
## c13.loc95 13 95.000000 1.159326e-01
## c13.loc96 13 96.000000 1.697782e-01
## c13.loc97 13 97.000000 2.319410e-01
## c13.loc98 13 98.000000 2.997871e-01
## c13.loc99 13 99.000000 3.704266e-01
## c13.loc100 13 100.000000 4.410883e-01
## c13.loc101 13 101.000000 5.094096e-01
## c13.loc102 13 102.000000 5.735950e-01
## c13.loc103 13 103.000000 6.324508e-01
## MgSTS55 13 103.008316 6.329112e-01
## c13.loc104 13 104.000000 6.576425e-01
## c13.loc105 13 105.000000 6.789842e-01
## c13.loc106 13 106.000000 6.954009e-01
## c13.loc107 13 107.000000 7.055826e-01
## c13.loc108 13 108.000000 7.083680e-01
## c13.loc109 13 109.000000 7.029084e-01
## c13.loc110 13 110.000000 6.888215e-01
## c13.loc111 13 111.000000 6.662935e-01
## c13.loc112 13 112.000000 6.360980e-01
## c13.loc113 13 113.000000 5.995183e-01
## c13.loc114 13 114.000000 5.581889e-01
## c13.loc115 13 115.000000 5.138937e-01
## c13.loc116 13 116.000000 4.683683e-01
## c13.loc117 13 117.000000 4.231442e-01
## c13.loc118 13 118.000000 3.794538e-01
## CB126 13 118.277636 3.678082e-01
## c13.loc119 13 119.000000 3.321216e-01
## c13.loc120 13 120.000000 2.816744e-01
## c13.loc121 13 121.000000 2.328613e-01
## c13.loc122 13 122.000000 1.880130e-01
## c13.loc123 13 123.000000 1.486302e-01
## c13.loc124 13 124.000000 1.154115e-01
## AAT240 13 124.055192 1.137764e-01
## c13.loc125 13 125.000000 1.744793e-01
## c13.loc126 13 126.000000 2.642559e-01
## c13.loc127 13 127.000000 3.803549e-01
## c13.loc128 13 128.000000 5.206810e-01
## c13.loc129 13 129.000000 6.808489e-01
## c13.loc130 13 130.000000 8.546701e-01
## c13.loc131 13 131.000000 1.035094e+00
## c13.loc132 13 132.000000 1.215283e+00
## MgSTS622 13 132.601495 1.320671e+00
## c13.loc133 13 133.000000 1.347849e+00
## c13.loc134 13 134.000000 1.412996e+00
## c13.loc135 13 135.000000 1.471516e+00
## c13.loc136 13 136.000000 1.521386e+00
## c13.loc137 13 137.000000 1.560807e+00
## c13.loc138 13 138.000000 1.588380e+00
## c13.loc139 13 139.000000 1.603251e+00
## c13.loc140 13 140.000000 1.605194e+00
## c13.loc141 13 141.000000 1.594616e+00
## c13.loc142 13 142.000000 1.572480e+00
## c13.loc143 13 143.000000 1.540165e+00
## MgSTS326 13 143.607075 1.516369e+00
## c13.loc144 13 144.000000 1.522413e+00
## c13.loc145 13 145.000000 1.540160e+00
## c13.loc146 13 146.000000 1.560568e+00
## c13.loc147 13 147.000000 1.582418e+00
## c13.loc148 13 148.000000 1.604138e+00
## c13.loc149 13 149.000000 1.624128e+00
## BA145 13 149.464606 1.632395e+00
## c13.loc150 13 150.000000 1.654976e+00
## c13.loc151 13 151.000000 1.704105e+00
## c13.loc152 13 152.000000 1.761624e+00
## c13.loc153 13 153.000000 1.826202e+00
## c13.loc154 13 154.000000 1.895560e+00
## c13.loc155 13 155.000000 1.966669e+00
## c13.loc156 13 156.000000 2.036159e+00
## c13.loc157 13 157.000000 2.100849e+00
## c13.loc158 13 158.000000 2.158219e+00
## c13.loc159 13 159.000000 2.206682e+00
## CC150 13 159.595164 2.230917e+00
## c13.loc160 13 160.000000 2.226644e+00
## c13.loc161 13 161.000000 2.211769e+00
## c13.loc162 13 162.000000 2.190440e+00
## c13.loc163 13 163.000000 2.162545e+00
## c13.loc164 13 164.000000 2.128522e+00
## c13.loc165 13 165.000000 2.089406e+00
## c13.loc166 13 166.000000 2.046762e+00
## c13.loc167 13 167.000000 2.002493e+00
## c13.loc168 13 168.000000 1.958576e+00
## c13.loc169 13 169.000000 1.916798e+00
## c13.loc170 13 170.000000 1.878557e+00
## BC199 13 170.660726 1.855772e+00
## c13.loc171 13 171.000000 1.845282e+00
## c13.loc172 13 172.000000 1.801590e+00
## c13.loc173 13 173.000000 1.742985e+00
## c13.loc174 13 174.000000 1.674804e+00
## c13.loc175 13 175.000000 1.602233e+00
## c13.loc176 13 176.000000 1.529737e+00
## c13.loc177 13 177.000000 1.460719e+00
## c13.loc178 13 178.000000 1.397455e+00
## MgSTS68 13 178.824945 1.350623e+00
## c13.loc179 13 179.000000 1.306224e+00
## MgSTS419 13 179.820822 1.090519e+00
## c13.loc180 13 180.000000 1.096735e+00
## c13.loc181 13 181.000000 1.134121e+00
## c13.loc182 13 182.000000 1.175340e+00
## c13.loc183 13 183.000000 1.219627e+00
## c13.loc184 13 184.000000 1.265654e+00
## c13.loc185 13 185.000000 1.311538e+00
## c13.loc186 13 186.000000 1.354958e+00
## c13.loc187 13 187.000000 1.393331e+00
## c13.loc188 13 188.000000 1.424025e+00
## BB198 13 188.678039 1.439145e+00
## AP3 14 0.000000 1.468683e-01
## c14.loc1 14 1.000000 1.874058e-01
## c14.loc2 14 2.000000 2.413837e-01
## c14.loc3 14 3.000000 3.062429e-01
## CA150 14 3.205217 3.204492e-01
## c14.loc4 14 4.000000 3.720908e-01
## c14.loc5 14 5.000000 4.442116e-01
## c14.loc6 14 6.000000 5.234801e-01
## c14.loc7 14 7.000000 6.086052e-01
## c14.loc8 14 8.000000 6.977360e-01
## c14.loc9 14 9.000000 7.886309e-01
## c14.loc10 14 10.000000 8.789272e-01
## c14.loc11 14 11.000000 9.664358e-01
## c14.loc12 14 12.000000 1.049372e+00
## CC320 14 12.173649 1.063093e+00
## c14.loc13 14 13.000000 1.067197e+00
## c14.loc14 14 14.000000 1.060471e+00
## c14.loc15 14 15.000000 1.043243e+00
## c14.loc16 14 16.000000 1.019096e+00
## MgSTS583 14 16.151945 1.015095e+00
## c14.loc17 14 17.000000 1.080481e+00
## c14.loc18 14 18.000000 1.148233e+00
## c14.loc19 14 19.000000 1.204460e+00
## c14.loc20 14 20.000000 1.248747e+00
## c14.loc21 14 21.000000 1.281215e+00
## MgSTS24 14 21.214324 1.286626e+00
## c14.loc22 14 22.000000 1.334919e+00
## c14.loc23 14 23.000000 1.383507e+00
## c14.loc24 14 24.000000 1.414383e+00
## c14.loc25 14 25.000000 1.425854e+00
## c14.loc26 14 26.000000 1.418089e+00
## MgSTS483 14 26.141772 1.415540e+00
## c14.loc27 14 27.000000 1.446195e+00
## c14.loc28 14 28.000000 1.482087e+00
## c14.loc29 14 29.000000 1.517792e+00
## c14.loc30 14 30.000000 1.552921e+00
## c14.loc31 14 31.000000 1.587025e+00
## c14.loc32 14 32.000000 1.619600e+00
## c14.loc33 14 33.000000 1.650094e+00
## c14.loc34 14 34.000000 1.677907e+00
## c14.loc35 14 35.000000 1.702415e+00
## c14.loc36 14 36.000000 1.722976e+00
## c14.loc37 14 37.000000 1.738958e+00
## c14.loc38 14 38.000000 1.749762e+00
## c14.loc39 14 39.000000 1.754847e+00
## c14.loc40 14 40.000000 1.753765e+00
## c14.loc41 14 41.000000 1.746187e+00
## c14.loc42 14 42.000000 1.731933e+00
## c14.loc43 14 43.000000 1.710994e+00
## c14.loc44 14 44.000000 1.683539e+00
## c14.loc45 14 45.000000 1.649922e+00
## c14.loc46 14 46.000000 1.610662e+00
## c14.loc47 14 47.000000 1.566421e+00
## c14.loc48 14 48.000000 1.517964e+00
## c14.loc49 14 49.000000 1.466123e+00
## c14.loc50 14 50.000000 1.411747e+00
## c14.loc51 14 51.000000 1.355672e+00
## c14.loc52 14 52.000000 1.298683e+00
## c14.loc53 14 53.000000 1.241490e+00
## c14.loc54 14 54.000000 1.184719e+00
## c14.loc55 14 55.000000 1.128896e+00
## c14.loc56 14 56.000000 1.074455e+00
## MgSTS520 14 56.459365 1.050136e+00
## c14.loc57 14 57.000000 1.078198e+00
## c14.loc58 14 58.000000 1.141002e+00
## c14.loc59 14 59.000000 1.214347e+00
## c14.loc60 14 60.000000 1.291106e+00
## c14.loc61 14 61.000000 1.362340e+00
## c14.loc62 14 62.000000 1.420148e+00
## MgSTS631 14 62.223118 1.430588e+00
## c14.loc63 14 63.000000 1.441267e+00
## c14.loc64 14 64.000000 1.433659e+00
## c14.loc65 14 65.000000 1.396683e+00
## c14.loc66 14 66.000000 1.329101e+00
## c14.loc67 14 67.000000 1.235279e+00
## AA137 14 67.946228 1.130191e+00
## c14.loc68 14 68.000000 1.128543e+00
## c14.loc69 14 69.000000 1.081735e+00
## c14.loc70 14 70.000000 1.005325e+00
## c14.loc71 14 71.000000 9.054282e-01
## MgSTS17 14 71.713586 8.259943e-01
## c14.loc72 14 72.000000 8.105313e-01
## c14.loc73 14 73.000000 7.499523e-01
## c14.loc74 14 74.000000 6.803840e-01
## c14.loc75 14 75.000000 6.035399e-01
## c14.loc76 14 76.000000 5.221018e-01
## c14.loc77 14 77.000000 4.394527e-01
## c14.loc78 14 78.000000 3.591944e-01
## c14.loc79 14 79.000000 2.845916e-01
## c14.loc80 14 80.000000 2.181187e-01
## MgSTS29 14 80.312824 1.993783e-01
## c14.loc81 14 81.000000 1.883267e-01
## c14.loc82 14 82.000000 1.718257e-01
## MgSTS263 14 82.805803 1.584706e-01
## c14.loc83 14 83.000000 1.626506e-01
## c14.loc84 14 84.000000 1.878577e-01
## c14.loc85 14 85.000000 2.190169e-01
## c14.loc86 14 86.000000 2.557229e-01
## c14.loc87 14 87.000000 2.969123e-01
## c14.loc88 14 88.000000 3.409555e-01
## c14.loc89 14 89.000000 3.858890e-01
## c14.loc90 14 90.000000 4.297193e-01
## BD371 14 90.499407 4.505350e-01
## CC93 14 90.500065 4.505452e-01
## c14.loc91 14 91.000000 4.911010e-01
## c14.loc92 14 92.000000 5.720137e-01
## c14.loc93 14 93.000000 6.480671e-01
## c14.loc94 14 94.000000 7.144129e-01
## c14.loc95 14 95.000000 7.668491e-01
## c14.loc96 14 96.000000 8.024415e-01
## c14.loc97 14 97.000000 8.198989e-01
## c14.loc98 14 98.000000 8.196214e-01
## MgSTS133 14 98.897556 8.057779e-01
## c14.loc99 14 99.000000 8.025771e-01
## c14.loc100 14 100.000000 7.628151e-01
## c14.loc101 14 101.000000 7.080242e-01
## c14.loc102 14 102.000000 6.399757e-01
## c14.loc103 14 103.000000 5.629040e-01
## c14.loc104 14 104.000000 4.827678e-01
## c14.loc105 14 105.000000 4.058212e-01
## c14.loc106 14 106.000000 3.371581e-01
## MgSTS350 14 106.568567 3.032174e-01
## c14.loc107 14 107.000000 3.018983e-01
## c14.loc108 14 108.000000 2.976981e-01
## c14.loc109 14 109.000000 2.918552e-01
## c14.loc110 14 110.000000 2.843362e-01
## c14.loc111 14 111.000000 2.751788e-01
## c14.loc112 14 112.000000 2.645031e-01
## c14.loc113 14 113.000000 2.525162e-01
## c14.loc114 14 114.000000 2.395062e-01
## c14.loc115 14 115.000000 2.258229e-01
## c14.loc116 14 116.000000 2.118480e-01
## c14.loc117 14 117.000000 1.979579e-01
## c14.loc118 14 118.000000 1.844899e-01
## c14.loc119 14 119.000000 1.717163e-01
## MgSTS18 14 119.129023 1.701433e-01
## c14.loc120 14 120.000000 9.840879e-02
## c14.loc121 14 121.000000 3.660858e-02
## c14.loc122 14 122.000000 4.162960e-03
## c14.loc123 14 123.000000 4.368052e-03
## c14.loc124 14 124.000000 3.574425e-02
## MgSTS632 14 124.524319 6.270895e-02
## c14.loc125 14 125.000000 6.069322e-02
## c14.loc126 14 126.000000 5.685972e-02
## c14.loc127 14 127.000000 5.378387e-02
## c14.loc128 14 128.000000 5.170724e-02
## c14.loc129 14 129.000000 5.086853e-02
## c14.loc130 14 130.000000 5.145322e-02
## c14.loc131 14 131.000000 5.354281e-02
## c14.loc132 14 132.000000 5.708554e-02
## MgSTS578 14 132.471776 5.920223e-02
## c14.loc133 14 133.000000 5.614210e-02
## c14.loc134 14 134.000000 5.052099e-02
## c14.loc135 14 135.000000 4.546611e-02
## c14.loc136 14 136.000000 4.118284e-02
## MgSTS500 14 136.459544 3.949080e-02
## c14.loc137 14 137.000000 4.196832e-02
## c14.loc138 14 138.000000 4.930455e-02
## CA75 14 138.775710 5.788751e-02
## c14.loc139 14 139.000000 6.131946e-02
## c14.loc140 14 140.000000 7.838863e-02
## c14.loc141 14 141.000000 9.821584e-02
## c14.loc142 14 142.000000 1.208228e-01
## c14.loc143 14 143.000000 1.461440e-01
## c14.loc144 14 144.000000 1.740138e-01
## c14.loc145 14 145.000000 2.041574e-01
## c14.loc146 14 146.000000 2.361915e-01
## c14.loc147 14 147.000000 2.696336e-01
## c14.loc148 14 148.000000 3.039241e-01
## c14.loc149 14 149.000000 3.384571e-01
## c14.loc150 14 150.000000 3.726187e-01
## c14.loc151 14 151.000000 4.058254e-01
## c14.loc152 14 152.000000 4.375583e-01
## BD251 14 152.540775 4.538762e-01
## c14.loc153 14 153.000000 4.598608e-01
## c14.loc154 14 154.000000 4.718779e-01
## BB216 14 154.144590 4.734682e-01
## c14.loc155 14 155.000000 4.860606e-01
## c14.loc156 14 156.000000 5.009566e-01
## c14.loc157 14 157.000000 5.158831e-01
## c14.loc158 14 158.000000 5.306681e-01
## c14.loc159 14 159.000000 5.451049e-01
## c14.loc160 14 160.000000 5.589506e-01
## c14.loc161 14 161.000000 5.719261e-01
## c14.loc162 14 162.000000 5.837237e-01
## c14.loc163 14 163.000000 5.940209e-01
## c14.loc164 14 164.000000 6.025027e-01
## c14.loc165 14 165.000000 6.088894e-01
## c14.loc166 14 166.000000 6.129647e-01
## c14.loc167 14 167.000000 6.145988e-01
## c14.loc168 14 168.000000 6.137620e-01
## c14.loc169 14 169.000000 6.105250e-01
## c14.loc170 14 170.000000 6.050480e-01
## c14.loc171 14 171.000000 5.975621e-01
## MgSTS491 14 171.055828 5.970949e-01
## c14.loc172 14 172.000000 5.742894e-01
## c14.loc173 14 173.000000 5.474759e-01
## c14.loc174 14 174.000000 5.180308e-01
## c14.loc175 14 175.000000 4.862303e-01
## c14.loc176 14 176.000000 4.525921e-01
## c14.loc177 14 177.000000 4.179326e-01
## c14.loc178 14 178.000000 3.834153e-01
## c14.loc179 14 179.000000 3.505763e-01
## c14.loc180 14 180.000000 3.213098e-01
## c14.loc181 14 181.000000 2.977929e-01
## c14.loc182 14 182.000000 2.823309e-01
## c14.loc183 14 183.000000 2.771158e-01
## c14.loc184 14 184.000000 2.839075e-01
## c14.loc185 14 185.000000 3.036840e-01
## c14.loc186 14 186.000000 3.363396e-01
## c14.loc187 14 187.000000 3.805431e-01
## c14.loc188 14 188.000000 4.338538e-01
## c14.loc189 14 189.000000 4.931231e-01
## MgSTS494 14 189.954021 5.521937e-01
## c14.loc190 14 190.000000 5.507068e-01
## c14.loc191 14 191.000000 5.156294e-01
## c14.loc192 14 192.000000 4.767429e-01
## c14.loc193 14 193.000000 4.357760e-01
## c14.loc194 14 194.000000 3.949148e-01
## c14.loc195 14 195.000000 3.564325e-01
## c14.loc196 14 196.000000 3.222963e-01
## c14.loc197 14 197.000000 2.938864e-01
## CA267 14 197.829812 2.751972e-01
## c14.loc198 14 198.000000 2.687810e-01
## c14.loc199 14 199.000000 2.340584e-01
## c14.loc200 14 200.000000 2.047370e-01
## c14.loc201 14 201.000000 1.801548e-01
## c14.loc202 14 202.000000 1.593883e-01
## c14.loc203 14 203.000000 1.414140e-01
## c14.loc204 14 204.000000 1.252182e-01
## c14.loc205 14 205.000000 1.098724e-01
## c14.loc206 14 206.000000 9.461086e-02
## c14.loc207 14 207.000000 7.895655e-02
## c14.loc208 14 208.000000 6.296320e-02
## c14.loc209 14 209.000000 4.763344e-02
## c14.loc210 14 210.000000 3.547105e-02
## c14.loc211 14 211.000000 3.076329e-02
## c14.loc212 14 212.000000 3.861781e-02
## c14.loc213 14 213.000000 6.202741e-02
## MgSTS16 14 213.663268 8.517964e-02
## c14.loc214 14 214.000000 1.040179e-01
## c14.loc215 14 215.000000 3.526007e+00
## c14.loc216 14 216.000000 3.522551e+00
## c14.loc217 14 217.000000 3.435393e+00
## MgSTS398 14 217.024733 3.432122e+00
## c14.loc218 14 218.000000 3.339982e+00
## c14.loc219 14 219.000000 3.193899e+00
## c14.loc220 14 220.000000 2.284218e-01
## c14.loc221 14 221.000000 1.833180e-01
## MgSTS562 14 221.907349 1.557340e-01
## c14.loc222 14 222.000000 1.580617e-01
## c14.loc223 14 223.000000 2.182442e-01
## c14.loc224 14 224.000000 3.461685e-01
## c14.loc225 14 225.000000 5.256306e-01
## c14.loc226 14 226.000000 7.183486e-01
## MgSTS527 14 226.151581 7.459890e-01
## c14.loc227 14 227.000000 6.405194e-01
## c14.loc228 14 228.000000 4.884103e-01
## c14.loc229 14 229.000000 3.413100e-01
## c14.loc230 14 230.000000 2.227741e-01
## c14.loc231 14 231.000000 1.377187e-01
## MgSTS341 14 231.379393 1.134264e-01
## c14.loc232 14 232.000000 1.005203e-01
## c14.loc233 14 233.000000 9.257135e-02
## c14.loc234 14 234.000000 1.019284e-01
## c14.loc235 14 235.000000 1.285013e-01
## c14.loc236 14 236.000000 1.703324e-01
## c14.loc237 14 237.000000 2.239247e-01
## c14.loc238 14 238.000000 2.848896e-01
## c14.loc239 14 239.000000 3.486994e-01
## c14.loc240 14 240.000000 4.113451e-01
## AAT265 14 240.214192 4.241979e-01
## c14.loc241 14 241.000000 3.726528e-01
## c14.loc242 14 242.000000 3.102389e-01
## c14.loc243 14 243.000000 2.540540e-01
## c14.loc244 14 244.000000 2.071427e-01
## c14.loc245 14 245.000000 1.727068e-01
## c14.loc246 14 246.000000 1.537834e-01
## c14.loc247 14 247.000000 1.528464e-01
## c14.loc248 14 248.000000 1.714013e-01
## c14.loc249 14 249.000000 2.096778e-01
## c14.loc250 14 250.000000 2.665118e-01
## c14.loc251 14 251.000000 3.394661e-01
## c14.loc252 14 252.000000 4.251648e-01
## c14.loc253 14 253.000000 5.197509e-01
## c14.loc254 14 254.000000 6.193541e-01
## c14.loc255 14 255.000000 7.204700e-01
## BC542 14 255.999585 8.201588e-01(CIMthresold_CIM <- cim(qtl_probC,pheno.col=8,n.marcovar=4,
method="hk",
imp.method="imp", error.prob=0.001,
map.function="kosambi",n.perm=100))## [,1]
## [1,] 3.174940
## [2,] 2.932094
## [3,] 3.349853
## [4,] 3.631475
## [5,] 2.381306
## [6,] 2.775837
## [7,] 4.482875
## [8,] 3.409779
## [9,] 2.243020
## [10,] 4.220763
## [11,] 2.840082
## [12,] 3.080561
## [13,] 2.400415
## [14,] 3.683820
## [15,] 2.349920
## [16,] 2.617907
## [17,] 3.296576
## [18,] 2.833184
## [19,] 2.787924
## [20,] 2.823027
## [21,] 3.375271
## [22,] 1.918973
## [23,] 2.517402
## [24,] 2.419121
## [25,] 4.078331
## [26,] 3.010602
## [27,] 2.199870
## [28,] 3.040662
## [29,] 2.510001
## [30,] 2.179986
## [31,] 2.724275
## [32,] 2.903108
## [33,] 3.755630
## [34,] 3.616531
## [35,] 3.903709
## [36,] 2.584450
## [37,] 4.156390
## [38,] 3.487951
## [39,] 3.001919
## [40,] 2.460094
## [41,] 3.130576
## [42,] 2.606448
## [43,] 3.382082
## [44,] 2.191260
## [45,] 2.625200
## [46,] 2.801326
## [47,] 3.611087
## [48,] 1.954349
## [49,] 4.400486
## [50,] 3.281310
## [51,] 5.653326
## [52,] 2.540173
## [53,] 3.235774
## [54,] 3.455215
## [55,] 3.466388
## [56,] 1.887684
## [57,] 2.756667
## [58,] 6.212978
## [59,] 3.883534
## [60,] 2.484475
## [61,] 3.918943
## [62,] 4.327715
## [63,] 2.822051
## [64,] 2.707548
## [65,] 2.777770
## [66,] 3.103747
## [67,] 2.889306
## [68,] 2.077486
## [69,] 3.563376
## [70,] 2.388249
## [71,] 3.424105
## [72,] 2.232158
## [73,] 2.999710
## [74,] 1.910747
## [75,] 3.201312
## [76,] 2.467614
## [77,] 3.232175
## [78,] 3.114262
## [79,] 2.871584
## [80,] 2.305370
## [81,] 3.176317
## [82,] 3.296092
## [83,] 3.155953
## [84,] 3.161232
## [85,] 2.075910
## [86,] 3.319450
## [87,] 3.307471
## [88,] 3.314211
## [89,] 3.188370
## [90,] 5.325436
## [91,] 3.678701
## [92,] 3.782238
## [93,] 2.764643
## [94,] 4.076146
## [95,] 3.311139
## [96,] 2.145897
## [97,] 2.826894
## [98,] 4.282678
## [99,] 2.431416
## [100,] 3.488064
## attr(,"class")
## [1] "scanoneperm" "matrix"CIM_Composite<-summary(CIMthresold_CIM, alpha = 0.05)[1,1]
plot_hk_CIM<-plot(CIMthresold_CIM, col=c("blue"))
abline(plot_hk_CIM, v=CIM_Composite,CIMthresold_CIM,
col="red", lty=6)
#QTL composite Interval Mapping
plot(out.cim.com, out.hk,
col = c("blue", "gray"), lwd=0.1)
add.cim.covar(out.cim.com)
abline(h=CIM_Composite, col="blue", lty=6, )
abline(h=operm.hk[1,1], col="gray", lty=6)
(QTL_CIM<-summary(out.cim.com, perms=CIMthresold_CIM, alpha=0.05,
pvalues=T))## chr pos lod pval
## c6.loc151 6 151 5.79 0.01
## c13.loc55 13 55 10.79 0.00#QTL composite Interval Mapping chromosome 6
plot(out.cim.com,chr=6, out.hk,
main=substitute("QTL detected in chromossome 6
by CIM and IM"),
col = c("blue", "gray"), lwd=0.5)
abline(h=CIM_Composite, col="blue", lty=6)
abline(h=operm.hk[1,1], col="gray", lty=6)
#QTL Composite Interval Mapping chromossome 13
plot(out.cim.com,chr=c(13), out.hk,
main=substitute("QTL detected in chromossome 13 by CIM and IM"),
col = c("blue", "gray"), lwd=1)
abline(h=CIM_Composite, col="blue", lty=6)
abline(h=operm.hk[1,1], col="gray", lty=6)
#QTLs in genetic map
#Localization do QTL chromossomes
chr1<-(lodint(out.cim.com, chr=6, expandtomarkers=T))
chr3<-(lodint(out.cim.com, chr=13, expandtomarkers=T))
qc <- c("6", "13")
qp <- c(135, 57)
fake.f2_CIM <- subset(fake_f2_qtl, chr=qc)
fake.f2 <- sim.geno(fake.f2_CIM, n.draws=8, step=2, err=0.001)
(qtl_CIM <- makeqtl(fake.f2, qc, qp, what="draws"))## QTL object containing imputed genotypes, with 8 imputations.
##
## name chr pos n.gen
## Q1 6@134.0 6 134 3
## Q2 13@56.0 13 56 3plot(qtl_CIM)
plot(out.cim.com,chr=6, out.hk,
main=substitute("Chromossome 6
by CIM and IM"),
col = c("blue", "black"), lwd=0.5)
abline(h=CIM_Composite, col="blue", lty=6)
abline(h=operm.hk[1,1], col="black", lty=6)
#### Comparison between CIM and MIM
#Composite Interval Mapping chromossome and Interval Mapping
plot(out.cim.com,chr=c(13), out.hk,
main=substitute("Chromossome 13 by CIM and IM"),
col = c("blue", "black"), lwd=1)
abline(h=CIM_Composite, col="blue", lty=6)
abline(h=operm.hk[1,1], col="black", lty=6)
plotLodProfile(ref, main=substitute("Multiple Interval Mapping"))
par(mfrow=c(1,3))
plot(qtl_IntervalMapping, main=substitute("Interval Mapping"))
plot(qtl_CIM, main=substitute("Composite Interval Mapping"))
plot(ref, chr= c(6, 13), main=substitute("Multiple Interval Mapping"))## Warning in thepos - sapply(map[whchr], min): longitud de objeto mayor no es
## múltiplo de la longitud de uno menor