Experiment with simple fine mapping approach
Matthew Stephens
2017-10-17
Last updated: 2017-10-24
Code version: a32416d
Simulate some data
We will simulate data with some “LD structure”. b=3 blocks of p=10 variables, each correlated through some latent Z. The phenotype model is \(Y = X_1+...+X_b + e\) where \(X_1,\dots,X_b\) are each a member of a block, and \(e\) is \(N(0,sd)\).
set.seed(1)
simulate = function(n=100,p=10,b=3,sd=1){
Z = list()
X = list()
Y = list()
for(i in 1:b) Z[[i]] = rnorm(n)
for(i in 1:b) X[[i]] = Z[[i]] + matrix(rnorm(n*p),nrow=n)
for(i in 1:b) Y[[i]] = X[[i]][,1]
X = do.call(cbind,X) # bind columns of X and Y
Y = do.call(cbind,Y)
Y = rowSums(Y) # each of the betas is 1
Y = Y + rnorm(n,sd = sd)
return(list(X=X,Y=Y))
}
d = simulate()Note how the first (causal) variable in each block is more correlated than others.
plot(cor(d$Y,d$X)[1,],main= "correlation of Y with each variable")
#to run this chunk had to install varbvs
#library(devtools)
#install_github("pcarbo/varbvs",subdir = "varbvs-R")
library("varbvs")
sd.grid <- c(0,10^seq(-2,1,length.out = 19))
res = varbvsmix(d$X,Z=NULL,d$Y,sd.grid)Fitting variational approximation for linear regression model with
mixture-of-normals priors.
samples: 100 mixture component sd's: 0.1..3.2
variables: 30 mixture component drop thresh.: 1.0e-08
covariates: 0 fit mixture weights: yes
mixture size: 20 fit residual var. (sigma): yes
intercept: yes convergence tolerance 1.0e-04
variational max. --------- hyperparameters ---------
iter lower bound change sigma mixture sd's mix. weights (drop)
0001 -2.494589e+02 1.9e-01 3.1e+00 [0.1,3] [0.007,0.142] (0)
0002 -2.090072e+02 2.2e-01 1.9e+00 [0.1,3] [0.001,0.243] (0)
0003 -1.908186e+02 2.7e-01 1.4e+00 [0.1,3] [0.000,0.319] (0)
0004 -1.822522e+02 2.8e-01 1.3e+00 [0.1,3] [0.000,0.389] (0)
0005 -1.769232e+02 2.9e-01 1.1e+00 [0.1,3] [0.000,0.437] (0)
0006 -1.728738e+02 3.2e-01 1.1e+00 [0.1,3] [0.000,0.501] (0)
0007 -1.695833e+02 2.9e-01 1.0e+00 [0.1,3] [0.000,0.544] (0)
0008 -1.673229e+02 2.1e-01 9.8e-01 [0.1,3] [0.000,0.597] (0)
0009 -1.657857e+02 1.6e-01 9.6e-01 [0.1,3] [0.000,0.634] (0)
0010 -1.648616e+02 1.2e-01 9.5e-01 [0.1,3] [0.000,0.674] (0)
0011 -1.642329e+02 1.2e-01 9.4e-01 [0.1,3] [0.000,0.702] (0)
0012 -1.638182e+02 9.9e-02 9.4e-01 [0.1,3] [0.000,0.730] (0)
0013 -1.635206e+02 6.6e-02 9.4e-01 [0.1,3] [0.000,0.750] (0)
0014 -1.633047e+02 5.9e-02 9.4e-01 [0.1,3] [0.000,0.769] (0)
0015 -1.631438e+02 3.3e-02 9.4e-01 [0.1,3] [0.000,0.783] (0)
0016 -1.630186e+02 3.3e-02 9.4e-01 [0.1,3] [0.000,0.794] (0)
0017 -1.629230e+02 2.3e-02 9.5e-01 [0.1,3] [0.000,0.804] (0)
0018 -1.628462e+02 2.0e-02 9.5e-01 [0.1,3] [0.000,0.811] (0)
0019 -1.627867e+02 1.6e-02 9.5e-01 [0.1,3] [0.000,0.817] (0)
0020 -1.627383e+02 1.2e-02 9.5e-01 [0.1,3] [0.000,0.822] (0)
0021 -1.627002e+02 9.9e-03 9.5e-01 [0.1,3] [0.000,0.826] (1)
0022 -1.626691e+02 8.2e-03 9.5e-01 [0.1,3] [0.000,0.829] (1)
0023 -1.626442e+02 7.6e-03 9.5e-01 [0.1,3] [0.000,0.831] (1)
0024 -1.626237e+02 7.0e-03 9.5e-01 [0.1,3] [0.000,0.833] (1)
0025 -1.626070e+02 6.5e-03 9.5e-01 [0.1,2] [0.000,0.835] (2)
0026 -1.625932e+02 6.0e-03 9.5e-01 [0.1,2] [0.000,0.836] (2)
0027 -1.625818e+02 5.6e-03 9.5e-01 [0.1,2] [0.000,0.837] (2)
0028 -1.625722e+02 5.3e-03 9.5e-01 [0.1,2] [0.000,0.838] (2)
0029 -1.625642e+02 5.0e-03 9.5e-01 [0.1,2] [0.000,0.839] (2)
0030 -1.625574e+02 4.7e-03 9.5e-01 [0.1,2] [0.000,0.839] (2)
0031 -1.625515e+02 4.5e-03 9.5e-01 [0.1,2] [0.000,0.840] (2)
0032 -1.625465e+02 4.3e-03 9.5e-01 [0.1,2] [0.000,0.840] (2)
0033 -1.625422e+02 4.1e-03 9.5e-01 [0.1,2] [0.000,0.840] (3)
0034 -1.625384e+02 4.0e-03 9.5e-01 [0.1,2] [0.000,0.840] (3)
0035 -1.625350e+02 3.8e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0036 -1.625320e+02 3.6e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0037 -1.625294e+02 3.5e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0038 -1.625270e+02 3.3e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0039 -1.625249e+02 3.2e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0040 -1.625229e+02 3.0e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0041 -1.625212e+02 2.9e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0042 -1.625195e+02 2.8e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0043 -1.625180e+02 2.7e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0044 -1.625166e+02 2.6e-03 9.5e-01 [0.1,2] [0.000,0.841] (3)
0045 -1.625153e+02 2.5e-03 9.5e-01 [0.1,1] [0.000,0.841] (4)
0046 -1.625141e+02 2.4e-03 9.5e-01 [0.1,1] [0.000,0.841] (4)
0047 -1.625130e+02 2.3e-03 9.5e-01 [0.1,1] [0.000,0.841] (4)
0048 -1.625120e+02 2.2e-03 9.5e-01 [0.1,1] [0.000,0.841] (4)
0049 -1.625110e+02 2.1e-03 9.5e-01 [0.1,1] [0.000,0.841] (4)
0050 -1.625101e+02 2.0e-03 9.5e-01 [0.1,1] [0.000,0.841] (4)
0051 -1.625092e+02 1.9e-03 9.5e-01 [0.1,1] [0.000,0.841] (4)
0052 -1.625084e+02 1.8e-03 9.5e-01 [0.1,1] [0.000,0.841] (5)
0053 -1.625076e+02 1.8e-03 9.5e-01 [0.1,1] [0.000,0.841] (6)
0054 -1.625069e+02 1.7e-03 9.5e-01 [0.1,1] [0.000,0.841] (6)
0055 -1.625062e+02 1.6e-03 9.5e-01 [0.1,1] [0.000,0.841] (7)
0056 -1.625055e+02 1.6e-03 9.5e-01 [0.1,1] [0.000,0.841] (7)
0057 -1.625049e+02 1.5e-03 9.5e-01 [0.1,1] [0.000,0.841] (7)
0058 -1.625044e+02 1.4e-03 9.5e-01 [0.1,1] [0.000,0.841] (7)
0059 -1.625038e+02 1.4e-03 9.5e-01 [0.1,1] [0.000,0.841] (7)
0060 -1.625033e+02 1.3e-03 9.5e-01 [0.1,1] [0.000,0.841] (8)
0061 -1.625028e+02 1.3e-03 9.5e-01 [0.1,1] [0.000,0.841] (8)
0062 -1.625023e+02 1.2e-03 9.5e-01 [0.1,1] [0.000,0.841] (8)
0063 -1.625019e+02 1.2e-03 9.5e-01 [0.1,1] [0.000,0.841] (8)
0064 -1.625015e+02 1.1e-03 9.5e-01 [0.1,1] [0.000,0.841] (8)
0065 -1.625011e+02 1.1e-03 9.5e-01 [0.1,1] [0.000,0.841] (9)
0066 -1.625007e+02 1.1e-03 9.5e-01 [0.1,1] [0.000,0.841] (10)
0067 -1.625003e+02 1.0e-03 9.5e-01 [0.1,1] [0.000,0.841] (10)
0068 -1.625000e+02 9.8e-04 9.5e-01 [0.1,1] [0.000,0.841] (10)
0069 -1.624997e+02 9.5e-04 9.5e-01 [0.1,1] [0.000,0.841] (11)
0070 -1.624993e+02 9.1e-04 9.5e-01 [0.1,1] [0.000,0.841] (11)
0071 -1.624991e+02 8.8e-04 9.5e-01 [0.1,1] [0.000,0.841] (11)
0072 -1.624988e+02 8.5e-04 9.5e-01 [0.1,1] [0.000,0.841] (11)
0073 -1.624985e+02 8.2e-04 9.5e-01 [0.1,1] [0.000,0.841] (11)
0074 -1.624983e+02 7.9e-04 9.5e-01 [0.1,1] [0.000,0.841] (11)
0075 -1.624980e+02 7.7e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0076 -1.624978e+02 7.4e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0077 -1.624976e+02 7.2e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0078 -1.624974e+02 7.0e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0079 -1.624972e+02 6.8e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0080 -1.624970e+02 6.5e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0081 -1.624968e+02 6.4e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0082 -1.624967e+02 6.2e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0083 -1.624965e+02 6.0e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0084 -1.624964e+02 5.8e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0085 -1.624962e+02 5.7e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0086 -1.624961e+02 5.5e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0087 -1.624960e+02 5.4e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0088 -1.624959e+02 5.2e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0089 -1.624957e+02 5.1e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0090 -1.624956e+02 5.0e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0091 -1.624955e+02 4.8e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0092 -1.624954e+02 4.7e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0093 -1.624953e+02 4.6e-04 9.5e-01 [0.1,1] [0.000,0.841] (12)
0094 -1.624953e+02 4.5e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0095 -1.624952e+02 4.4e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0096 -1.624951e+02 4.3e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0097 -1.624950e+02 4.2e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0098 -1.624950e+02 4.1e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0099 -1.624949e+02 4.1e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0100 -1.624948e+02 4.0e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0101 -1.624948e+02 3.9e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0102 -1.624947e+02 3.8e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0103 -1.624946e+02 3.8e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0104 -1.624946e+02 3.7e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0105 -1.624945e+02 3.6e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0106 -1.624945e+02 3.6e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0107 -1.624945e+02 3.5e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0108 -1.624944e+02 3.5e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0109 -1.624944e+02 3.4e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0110 -1.624943e+02 3.4e-04 9.5e-01 [0.5,1] [0.000,0.841] (13)
0111 -1.624943e+02 3.3e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0112 -1.624943e+02 3.3e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0113 -1.624942e+02 3.2e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0114 -1.624942e+02 3.2e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0115 -1.624942e+02 3.2e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0116 -1.624941e+02 3.1e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0117 -1.624941e+02 3.1e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0118 -1.624941e+02 3.1e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0119 -1.624941e+02 3.0e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0120 -1.624940e+02 3.0e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0121 -1.624940e+02 3.0e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0122 -1.624940e+02 2.9e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0123 -1.624940e+02 2.9e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0124 -1.624940e+02 2.9e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0125 -1.624939e+02 2.9e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0126 -1.624939e+02 2.8e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0127 -1.624939e+02 2.8e-04 9.5e-01 [0.6,1] [0.000,0.841] (14)
0128 -1.624939e+02 2.8e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0129 -1.624939e+02 2.8e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0130 -1.624939e+02 2.7e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0131 -1.624938e+02 2.7e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0132 -1.624938e+02 2.7e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0133 -1.624938e+02 2.7e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0134 -1.624938e+02 2.7e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0135 -1.624938e+02 2.7e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0136 -1.624938e+02 2.6e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0137 -1.624938e+02 2.6e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0138 -1.624938e+02 2.6e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0139 -1.624938e+02 2.6e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0140 -1.624937e+02 2.6e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0141 -1.624937e+02 2.6e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0142 -1.624937e+02 2.6e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0143 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0144 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0145 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0146 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0147 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0148 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0149 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0150 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0151 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0152 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0153 -1.624937e+02 2.5e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0154 -1.624937e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0155 -1.624937e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0156 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0157 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0158 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0159 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0160 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.841] (15)
0161 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0162 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0163 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0164 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0165 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0166 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0167 -1.624936e+02 2.4e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0168 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0169 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0170 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0171 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0172 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0173 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0174 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0175 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0176 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0177 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0178 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0179 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0180 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0181 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0182 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0183 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0184 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0185 -1.624936e+02 2.3e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0186 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0187 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0188 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0189 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0190 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0191 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0192 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0193 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0194 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0195 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0196 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0197 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0198 -1.624936e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0199 -1.624935e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0200 -1.624935e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0201 -1.624935e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0202 -1.624935e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0203 -1.624935e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0204 -1.624935e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0205 -1.624935e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0206 -1.624935e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0207 -1.624935e+02 2.2e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0208 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0209 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0210 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0211 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0212 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0213 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0214 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0215 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0216 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0217 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0218 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0219 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0220 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0221 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0222 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0223 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0224 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0225 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0226 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0227 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0228 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0229 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0230 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0231 -1.624935e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0232 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0233 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0234 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0235 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0236 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0237 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0238 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0239 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0240 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0241 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0242 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0243 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0244 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0245 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0246 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0247 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0248 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0249 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0250 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0251 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0252 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0253 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0254 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0255 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0256 -1.624935e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0257 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0258 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0259 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0260 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0261 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0262 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0263 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0264 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0265 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0266 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0267 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0268 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0269 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0270 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0271 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0272 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0273 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0274 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0275 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0276 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0277 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0278 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0279 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0280 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0281 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0282 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0283 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0284 -1.624935e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0285 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0286 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0287 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0288 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0289 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0290 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0291 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0292 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0293 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0294 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0295 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0296 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0297 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0298 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0299 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0300 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0301 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0302 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0303 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0304 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0305 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0306 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0307 -1.624935e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.842] (15)
0308 -1.624935e+02 1.8e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0309 -1.624935e+02 1.8e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0310 -1.624935e+02 1.8e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0311 -1.624935e+02 1.8e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0312 -1.624935e+02 1.8e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0313 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0314 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0315 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0316 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0317 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0318 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0319 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0320 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0321 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0322 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0323 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0324 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0325 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0326 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0327 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0328 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0329 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0330 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0331 -1.624935e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0332 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0333 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0334 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0335 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0336 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0337 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0338 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0339 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0340 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0341 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0342 -1.624934e+02 1.7e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0343 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0344 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0345 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0346 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0347 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0348 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0349 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0350 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0351 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0352 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0353 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0354 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0355 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0356 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0357 -1.624934e+02 1.6e-04 9.5e-01 [0.7,1] [0.000,0.842] (16)
0358 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0359 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0360 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0361 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0362 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0363 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0364 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0365 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0366 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0367 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0368 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0369 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0370 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0371 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0372 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0373 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0374 -1.624934e+02 1.6e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0375 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0376 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0377 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0378 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0379 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0380 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0381 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0382 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0383 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0384 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0385 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0386 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0387 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0388 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0389 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0390 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0391 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0392 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0393 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0394 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0395 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0396 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.067,0.842] (17)
0397 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0398 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0399 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0400 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0401 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0402 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0403 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0404 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0405 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0406 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0407 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0408 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0409 -1.624934e+02 1.5e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0410 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0411 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0412 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0413 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0414 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0415 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0416 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0417 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0418 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0419 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0420 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0421 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0422 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0423 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0424 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0425 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0426 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0427 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0428 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0429 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0430 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0431 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0432 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0433 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0434 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0435 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0436 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0437 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0438 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0439 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0440 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0441 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.066,0.842] (17)
0442 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0443 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0444 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0445 -1.624934e+02 1.4e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0446 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0447 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0448 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0449 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0450 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0451 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0452 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0453 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0454 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0455 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0456 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0457 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0458 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0459 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0460 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0461 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0462 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0463 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0464 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0465 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0466 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0467 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0468 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0469 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0470 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0471 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0472 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0473 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0474 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0475 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0476 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0477 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0478 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0479 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0480 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0481 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0482 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0483 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0484 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0485 -1.624934e+02 1.3e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0486 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0487 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0488 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0489 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0490 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.065,0.842] (17)
0491 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0492 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0493 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0494 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0495 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0496 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0497 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0498 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0499 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0500 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0501 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0502 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0503 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0504 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0505 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0506 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0507 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0508 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0509 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0510 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0511 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0512 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0513 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0514 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0515 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0516 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0517 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0518 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0519 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0520 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0521 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0522 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0523 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0524 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0525 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0526 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0527 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0528 -1.624934e+02 1.2e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0529 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0530 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0531 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0532 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0533 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0534 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0535 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0536 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0537 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0538 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0539 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0540 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0541 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0542 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0543 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0544 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0545 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.064,0.842] (17)
0546 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0547 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0548 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0549 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0550 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0551 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0552 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0553 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0554 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0555 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0556 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0557 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0558 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0559 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0560 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0561 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0562 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0563 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0564 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0565 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0566 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0567 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0568 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0569 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0570 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0571 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0572 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0573 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0574 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0575 -1.624934e+02 1.1e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0576 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0577 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0578 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0579 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0580 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0581 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0582 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0583 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0584 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0585 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0586 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0587 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0588 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0589 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0590 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0591 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0592 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0593 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0594 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0595 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0596 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0597 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0598 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0599 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0600 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
0601 -1.624934e+02 1.0e-04 9.5e-01 [0.7,0.8] [0.063,0.842] (17)
plot(res$lfsr)
Try adding three variables that are exact duplicates of the causal ones. Notice it picks just one of the causal ones and zeros out the others.
res2 = varbvsmix(cbind(d$X[,1],d$X[,11],d$X[,21],d$X),Z=NULL,d$Y,sd.grid)Fitting variational approximation for linear regression model with
mixture-of-normals priors.
samples: 100 mixture component sd's: 0.1..3.2
variables: 33 mixture component drop thresh.: 1.0e-08
covariates: 0 fit mixture weights: yes
mixture size: 20 fit residual var. (sigma): yes
intercept: yes convergence tolerance 1.0e-04
variational max. --------- hyperparameters ---------
iter lower bound change sigma mixture sd's mix. weights (drop)
0001 -2.537229e+02 2.1e-01 3.0e+00 [0.1,3] [0.008,0.124] (0)
0002 -2.171322e+02 1.9e-01 2.0e+00 [0.1,3] [0.002,0.167] (0)
0003 -2.033323e+02 1.6e-01 1.6e+00 [0.1,3] [0.001,0.196] (0)
0004 -1.974837e+02 1.7e-01 1.4e+00 [0.1,3] [0.000,0.207] (0)
0005 -1.938046e+02 1.5e-01 1.3e+00 [0.1,3] [0.000,0.230] (0)
0006 -1.910621e+02 1.3e-01 1.3e+00 [0.1,3] [0.000,0.242] (0)
0007 -1.887412e+02 1.5e-01 1.2e+00 [0.1,3] [0.000,0.271] (0)
0008 -1.867797e+02 1.6e-01 1.1e+00 [0.1,3] [0.000,0.289] (0)
0009 -1.850343e+02 1.7e-01 1.1e+00 [0.1,3] [0.000,0.326] (0)
0010 -1.834845e+02 1.6e-01 1.1e+00 [0.1,3] [0.000,0.349] (0)
0011 -1.820895e+02 1.6e-01 1.1e+00 [0.1,3] [0.000,0.389] (0)
0012 -1.807757e+02 1.3e-01 1.0e+00 [0.1,3] [0.000,0.414] (0)
0013 -1.795720e+02 1.9e-01 1.0e+00 [0.1,3] [0.000,0.456] (0)
0014 -1.783348e+02 1.1e-01 1.0e+00 [0.1,3] [0.000,0.482] (0)
0015 -1.772114e+02 2.0e-01 9.9e-01 [0.1,3] [0.000,0.525] (0)
0016 -1.760322e+02 1.6e-01 9.8e-01 [0.1,3] [0.000,0.553] (0)
0017 -1.750030e+02 1.6e-01 9.8e-01 [0.1,3] [0.000,0.596] (0)
0018 -1.738424e+02 2.0e-01 9.7e-01 [0.1,3] [0.000,0.626] (0)
0019 -1.726888e+02 3.8e-01 9.8e-01 [0.1,3] [0.000,0.671] (0)
0020 -1.712877e+02 2.0e-01 9.6e-01 [0.1,3] [0.000,0.702] (0)
0021 -1.703693e+02 2.2e-01 9.6e-01 [0.1,3] [0.000,0.734] (0)
0022 -1.697588e+02 1.0e-01 9.6e-01 [0.1,3] [0.000,0.756] (0)
0023 -1.694008e+02 7.3e-02 9.5e-01 [0.1,3] [0.000,0.772] (0)
0024 -1.691751e+02 4.1e-02 9.5e-01 [0.1,3] [0.000,0.785] (0)
0025 -1.690128e+02 3.6e-02 9.5e-01 [0.1,3] [0.000,0.794] (0)
0026 -1.688948e+02 1.9e-02 9.5e-01 [0.1,3] [0.000,0.801] (0)
0027 -1.688018e+02 1.6e-02 9.5e-01 [0.1,3] [0.000,0.806] (0)
0028 -1.687272e+02 1.4e-02 9.5e-01 [0.1,3] [0.000,0.810] (0)
0029 -1.686653e+02 1.2e-02 9.5e-01 [0.1,3] [0.000,0.813] (0)
0030 -1.686130e+02 1.0e-02 9.5e-01 [0.1,3] [0.000,0.816] (0)
0031 -1.685678e+02 9.4e-03 9.5e-01 [0.1,3] [0.000,0.818] (0)
0032 -1.685282e+02 1.0e-02 9.5e-01 [0.1,3] [0.000,0.820] (1)
0033 -1.684927e+02 8.7e-03 9.5e-01 [0.1,3] [0.000,0.821] (1)
0034 -1.684607e+02 1.1e-02 9.5e-01 [0.1,3] [0.000,0.822] (1)
0035 -1.684311e+02 9.2e-03 9.5e-01 [0.1,3] [0.000,0.822] (1)
0036 -1.684038e+02 1.1e-02 9.5e-01 [0.1,3] [0.000,0.823] (1)
0037 -1.683780e+02 9.4e-03 9.5e-01 [0.1,3] [0.000,0.824] (1)
0038 -1.683537e+02 1.1e-02 9.5e-01 [0.1,3] [0.000,0.824] (1)
0039 -1.683303e+02 9.5e-03 9.5e-01 [0.1,3] [0.000,0.824] (1)
0040 -1.683080e+02 1.1e-02 9.5e-01 [0.1,3] [0.000,0.824] (1)
0041 -1.682862e+02 9.3e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0042 -1.682653e+02 1.1e-02 9.5e-01 [0.1,2] [0.000,0.825] (2)
0043 -1.682446e+02 9.0e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0044 -1.682246e+02 1.0e-02 9.5e-01 [0.1,2] [0.000,0.825] (2)
0045 -1.682048e+02 8.6e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0046 -1.681854e+02 9.5e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0047 -1.681662e+02 8.1e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0048 -1.681473e+02 9.2e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0049 -1.681285e+02 7.9e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0050 -1.681101e+02 9.4e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0051 -1.680917e+02 8.1e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0052 -1.680736e+02 9.5e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0053 -1.680555e+02 8.2e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0054 -1.680377e+02 9.6e-03 9.5e-01 [0.1,2] [0.000,0.825] (2)
0055 -1.680198e+02 8.3e-03 9.5e-01 [0.1,2] [0.000,0.825] (3)
0056 -1.680022e+02 9.7e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0057 -1.679846e+02 8.5e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0058 -1.679673e+02 9.7e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0059 -1.679499e+02 8.5e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0060 -1.679327e+02 9.7e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0061 -1.679154e+02 8.6e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0062 -1.678984e+02 9.7e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0063 -1.678813e+02 8.6e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0064 -1.678643e+02 9.6e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0065 -1.678473e+02 8.6e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0066 -1.678304e+02 9.5e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0067 -1.678133e+02 8.6e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0068 -1.677963e+02 9.9e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0069 -1.677789e+02 8.5e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0070 -1.677613e+02 1.0e-02 9.5e-01 [0.1,2] [0.000,0.824] (3)
0071 -1.677432e+02 8.6e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0072 -1.677245e+02 1.1e-02 9.5e-01 [0.1,2] [0.000,0.824] (3)
0073 -1.677045e+02 8.9e-03 9.5e-01 [0.1,2] [0.000,0.824] (3)
0074 -1.676831e+02 1.1e-02 9.5e-01 [0.1,2] [0.000,0.824] (3)
0075 -1.676587e+02 9.7e-03 9.5e-01 [0.1,1] [0.000,0.824] (4)
0076 -1.676301e+02 1.2e-02 9.5e-01 [0.1,1] [0.000,0.824] (4)
0077 -1.675931e+02 2.2e-02 9.6e-01 [0.1,1] [0.000,0.824] (4)
0078 -1.675408e+02 1.3e-02 9.6e-01 [0.1,1] [0.000,0.824] (4)
0079 -1.674499e+02 7.6e-02 9.6e-01 [0.1,1] [0.000,0.826] (4)
0080 -1.672504e+02 2.1e-02 9.6e-01 [0.1,1] [0.000,0.827] (4)
0081 -1.666602e+02 3.8e-01 9.8e-01 [0.1,1] [0.000,0.833] (4)
0082 -1.648648e+02 1.5e-01 9.6e-01 [0.1,1] [0.000,0.840] (4)
0083 -1.635905e+02 3.1e-01 9.5e-01 [0.1,1] [0.000,0.853] (4)
0084 -1.630913e+02 1.4e-02 9.5e-01 [0.1,1] [0.000,0.855] (4)
0085 -1.629968e+02 5.4e-02 9.5e-01 [0.1,1] [0.000,0.858] (4)
0086 -1.629772e+02 8.3e-03 9.5e-01 [0.1,1] [0.000,0.858] (4)
0087 -1.629724e+02 7.1e-03 9.5e-01 [0.1,1] [0.000,0.858] (5)
0088 -1.629700e+02 7.0e-03 9.5e-01 [0.1,1] [0.000,0.858] (5)
0089 -1.629680e+02 6.7e-03 9.5e-01 [0.1,1] [0.000,0.858] (5)
0090 -1.629662e+02 6.5e-03 9.5e-01 [0.1,1] [0.000,0.858] (6)
0091 -1.629645e+02 6.2e-03 9.5e-01 [0.1,1] [0.000,0.858] (7)
0092 -1.629628e+02 6.0e-03 9.5e-01 [0.1,1] [0.000,0.858] (7)
0093 -1.629612e+02 5.8e-03 9.5e-01 [0.1,1] [0.000,0.858] (7)
0094 -1.629597e+02 5.6e-03 9.5e-01 [0.1,1] [0.000,0.858] (7)
0095 -1.629583e+02 5.4e-03 9.5e-01 [0.1,1] [0.000,0.858] (7)
0096 -1.629569e+02 5.2e-03 9.5e-01 [0.1,1] [0.000,0.858] (8)
0097 -1.629556e+02 5.0e-03 9.5e-01 [0.1,1] [0.000,0.858] (8)
0098 -1.629543e+02 4.8e-03 9.5e-01 [0.1,1] [0.000,0.858] (8)
0099 -1.629531e+02 4.6e-03 9.5e-01 [0.1,1] [0.000,0.858] (8)
0100 -1.629520e+02 4.4e-03 9.5e-01 [0.1,1] [0.000,0.858] (8)
0101 -1.629509e+02 4.3e-03 9.5e-01 [0.1,1] [0.000,0.858] (8)
0102 -1.629498e+02 4.1e-03 9.5e-01 [0.1,1] [0.000,0.858] (9)
0103 -1.629488e+02 3.9e-03 9.5e-01 [0.1,1] [0.000,0.858] (9)
0104 -1.629479e+02 3.8e-03 9.5e-01 [0.1,1] [0.000,0.858] (9)
0105 -1.629470e+02 3.6e-03 9.5e-01 [0.1,1] [0.000,0.858] (9)
0106 -1.629461e+02 3.5e-03 9.5e-01 [0.3,1] [0.000,0.858] (10)
0107 -1.629453e+02 3.3e-03 9.5e-01 [0.3,1] [0.000,0.858] (11)
0108 -1.629445e+02 3.2e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0109 -1.629437e+02 3.0e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0110 -1.629430e+02 2.9e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0111 -1.629423e+02 2.8e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0112 -1.629416e+02 2.6e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0113 -1.629410e+02 2.5e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0114 -1.629404e+02 2.4e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0115 -1.629398e+02 2.3e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0116 -1.629393e+02 2.2e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0117 -1.629388e+02 2.1e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0118 -1.629383e+02 2.0e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0119 -1.629378e+02 1.9e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0120 -1.629374e+02 1.8e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0121 -1.629369e+02 1.7e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0122 -1.629365e+02 1.7e-03 9.5e-01 [0.4,1] [0.000,0.858] (12)
0123 -1.629361e+02 1.6e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0124 -1.629358e+02 1.5e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0125 -1.629354e+02 1.4e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0126 -1.629351e+02 1.4e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0127 -1.629348e+02 1.3e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0128 -1.629344e+02 1.3e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0129 -1.629342e+02 1.2e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0130 -1.629339e+02 1.1e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0131 -1.629336e+02 1.1e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0132 -1.629334e+02 1.0e-03 9.5e-01 [0.5,1] [0.000,0.858] (13)
0133 -1.629331e+02 9.8e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0134 -1.629329e+02 9.3e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0135 -1.629327e+02 8.9e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0136 -1.629325e+02 8.5e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0137 -1.629323e+02 8.1e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0138 -1.629321e+02 7.7e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0139 -1.629319e+02 7.3e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0140 -1.629318e+02 6.9e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0141 -1.629316e+02 6.6e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0142 -1.629315e+02 6.3e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0143 -1.629313e+02 6.0e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0144 -1.629312e+02 5.7e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0145 -1.629310e+02 5.4e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0146 -1.629309e+02 5.2e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0147 -1.629308e+02 4.9e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0148 -1.629307e+02 4.7e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0149 -1.629306e+02 4.4e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0150 -1.629305e+02 4.2e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0151 -1.629304e+02 4.0e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0152 -1.629303e+02 3.8e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0153 -1.629302e+02 3.6e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0154 -1.629301e+02 3.5e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0155 -1.629300e+02 3.3e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0156 -1.629300e+02 3.1e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0157 -1.629299e+02 3.0e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0158 -1.629298e+02 2.8e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0159 -1.629298e+02 2.7e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0160 -1.629297e+02 2.6e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0161 -1.629297e+02 2.4e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0162 -1.629296e+02 2.3e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0163 -1.629295e+02 2.2e-04 9.5e-01 [0.5,1] [0.000,0.858] (13)
0164 -1.629295e+02 2.1e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0165 -1.629294e+02 2.0e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0166 -1.629294e+02 1.9e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0167 -1.629294e+02 1.8e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0168 -1.629293e+02 1.7e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0169 -1.629293e+02 1.6e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0170 -1.629292e+02 1.5e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0171 -1.629292e+02 1.5e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0172 -1.629292e+02 1.4e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0173 -1.629291e+02 1.3e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0174 -1.629291e+02 1.3e-04 9.5e-01 [0.6,1] [0.000,0.858] (14)
0175 -1.629291e+02 1.2e-04 9.5e-01 [0.6,1] [0.000,0.858] (15)
0176 -1.629291e+02 1.1e-04 9.5e-01 [0.6,1] [0.000,0.858] (15)
0177 -1.629290e+02 1.1e-04 9.5e-01 [0.6,1] [0.000,0.858] (15)
0178 -1.629290e+02 1.0e-04 9.5e-01 [0.6,1] [0.000,0.858] (15)
0179 -1.629290e+02 9.8e-05 9.5e-01 [0.6,1] [0.000,0.858] (15)
plot(res2$lfsr)
Session information
sessionInfo()R version 3.3.2 (2016-10-31)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X El Capitan 10.11.6
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] varbvs_2.4-0
loaded via a namespace (and not attached):
[1] Rcpp_0.12.13 lattice_0.20-35 digest_0.6.12
[4] rprojroot_1.2 grid_3.3.2 backports_1.1.0
[7] git2r_0.19.0 magrittr_1.5 evaluate_0.10
[10] stringi_1.1.5 latticeExtra_0.6-28 Matrix_1.2-10
[13] rmarkdown_1.6 RColorBrewer_1.1-2 tools_3.3.2
[16] stringr_1.2.0 yaml_2.1.14 htmltools_0.3.6
[19] knitr_1.16 This R Markdown site was created with workflowr