Last updated: 2018-10-09

workflowr checks: (Click a bullet for more information)
  • R Markdown file: up-to-date

    Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

  • Environment: empty

    Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.

  • Seed: set.seed(12345)

    The command set.seed(12345) was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.

  • Session information: recorded

    Great job! Recording the operating system, R version, and package versions is critical for reproducibility.

  • Repository version: 2756e44

    Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility. The version displayed above was the version of the Git repository at the time these results were generated.

    Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:
    
    Ignored files:
        Ignored:    .DS_Store
        Ignored:    .Rhistory
        Ignored:    .Rproj.user/
        Ignored:    analysis/.RData
        Ignored:    analysis/.Rhistory
        Ignored:    analysis/ALStruct_cache/
        Ignored:    data/.Rhistory
        Ignored:    docs/figure/.DS_Store
    
    Untracked files:
        Untracked:  .dropbox
        Untracked:  Icon
        Untracked:  _workflowr.yml
        Untracked:  analysis/GTEX-cogaps.Rmd
        Untracked:  analysis/SPCAvRP.rmd
        Untracked:  analysis/eQTL.perm.rand.pdf
        Untracked:  analysis/ieQTL.perm.rand.pdf
        Untracked:  analysis/mash_bhat_z.Rmd
        Untracked:  analysis/mash_ieqtl_permutations.Rmd
        Untracked:  analysis/pseudodata.Rmd
        Untracked:  analysis/sc_bimodal.Rmd
        Untracked:  analysis/susie_example.Rmd
        Untracked:  analysis/svd-timing.Rmd
        Untracked:  analysis/test_sparse.Rmd
        Untracked:  analysis/z.txt
        Untracked:  code/multivariate_testfuncs.R
        Untracked:  data/4matthew/
        Untracked:  data/4matthew2/
        Untracked:  data/E-MTAB-2805.processed.1/
        Untracked:  data/ENSG00000156738.Sim_Y2.RDS
        Untracked:  data/GDS5363_full.soft.gz
        Untracked:  data/GSE41265_allGenesTPM.txt
        Untracked:  data/Muscle_Skeletal.ACTN3.pm1Mb.RDS
        Untracked:  data/Thyroid.FMO2.pm1Mb.RDS
        Untracked:  data/bmass.HaemgenRBC2016.MAF01.Vs2.MergedDataSources.200kRanSubset.ChrBPMAFMarkerZScores.vs1.txt.gz
        Untracked:  data/bmass.HaemgenRBC2016.Vs2.NewSNPs.ZScores.hclust.vs1.txt
        Untracked:  data/bmass.HaemgenRBC2016.Vs2.PreviousSNPs.ZScores.hclust.vs1.txt
        Untracked:  data/finemap_data/fmo2.sim/b.txt
        Untracked:  data/finemap_data/fmo2.sim/dap_out.txt
        Untracked:  data/finemap_data/fmo2.sim/dap_out2.txt
        Untracked:  data/finemap_data/fmo2.sim/dap_out2_snp.txt
        Untracked:  data/finemap_data/fmo2.sim/dap_out_snp.txt
        Untracked:  data/finemap_data/fmo2.sim/data
        Untracked:  data/finemap_data/fmo2.sim/fmo2.sim.config
        Untracked:  data/finemap_data/fmo2.sim/fmo2.sim.k
        Untracked:  data/finemap_data/fmo2.sim/fmo2.sim.k4.config
        Untracked:  data/finemap_data/fmo2.sim/fmo2.sim.k4.snp
        Untracked:  data/finemap_data/fmo2.sim/fmo2.sim.ld
        Untracked:  data/finemap_data/fmo2.sim/fmo2.sim.snp
        Untracked:  data/finemap_data/fmo2.sim/fmo2.sim.z
        Untracked:  data/finemap_data/fmo2.sim/pos.txt
        Untracked:  data/logm.csv
        Untracked:  data/m.cd.RDS
        Untracked:  data/m.cdu.old.RDS
        Untracked:  data/m.new.cd.RDS
        Untracked:  data/m.old.cd.RDS
        Untracked:  data/mainbib.bib.old
        Untracked:  data/mat.csv
        Untracked:  data/mat.txt
        Untracked:  data/mat_new.csv
        Untracked:  data/paintor_data/
        Untracked:  data/temp.txt
        Untracked:  data/y.txt
        Untracked:  data/y_f.txt
        Untracked:  docs/figure/eigen.Rmd/
        Untracked:  docs/figure/fmo2.sim.Rmd/
        Untracked:  docs/figure/newVB.elbo.Rmd/
        Untracked:  docs/figure/rbc_zscore_mash2.Rmd/
        Untracked:  docs/figure/rbc_zscore_mash2_analysis.Rmd/
        Untracked:  docs/figure/rbc_zscores.Rmd/
        Untracked:  docs/trend_files/
        Untracked:  docs/z.txt
        Untracked:  explore_udi.R
        Untracked:  output/fit.varbvs.RDS
        Untracked:  output/glmnet.fit.RDS
        Untracked:  output/test.bv.txt
        Untracked:  output/test.gamma.txt
        Untracked:  output/test.hyp.txt
        Untracked:  output/test.log.txt
        Untracked:  output/test.param.txt
        Untracked:  output/test2.bv.txt
        Untracked:  output/test2.gamma.txt
        Untracked:  output/test2.hyp.txt
        Untracked:  output/test2.log.txt
        Untracked:  output/test2.param.txt
        Untracked:  output/test3.bv.txt
        Untracked:  output/test3.gamma.txt
        Untracked:  output/test3.hyp.txt
        Untracked:  output/test3.log.txt
        Untracked:  output/test3.param.txt
        Untracked:  output/test4.bv.txt
        Untracked:  output/test4.gamma.txt
        Untracked:  output/test4.hyp.txt
        Untracked:  output/test4.log.txt
        Untracked:  output/test4.param.txt
        Untracked:  output/test5.bv.txt
        Untracked:  output/test5.gamma.txt
        Untracked:  output/test5.hyp.txt
        Untracked:  output/test5.log.txt
        Untracked:  output/test5.param.txt
    
    Unstaged changes:
        Modified:   analysis/_site.yml
        Deleted:    analysis/chunks.R
        Modified:   analysis/eigen.Rmd
        Modified:   analysis/fmo2.sim.Rmd
        Modified:   analysis/newVB.Rmd
        Modified:   analysis/wSVD.Rmd
    
    
    Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.
Expand here to see past versions:
    File Version Author Date Message
    Rmd 2756e44 stephens999 2018-10-09 workflowr::wflow_publish(“plra.Rmd”)


Introduction

I’m going to try to do a version of Poisson low rank approximation in R so I can play with it, based on Abhishek Sarkar’s python code (here)[https://github.com/aksarkar/wlra/blob/master/wlra/wlra.py].

Poisson low rank approximation

# return truncated svd (truncated to rank principal components)
lra = function(x, rank){
  x.svd = svd(x,nu=rank,nv=rank)
  
  return(x.svd$u %*% diag(x.svd$d[1:rank],nrow=rank) %*% t(x.svd$v))
}
  
# Return the weighted low rank approximation of x
# Minimize the weighted Frobenius norm between x and the approximation z using EM [SJ03].
#
#' @param x input data (n, p)
#' @param w input weights (n, p)
#' @param rank rank of the approximation (non-negative)
#' @param max_iters - maximum number of EM iterations
#' @param atol - minimum absolute difference in objective function for convergence
#' @param verbose - print objective function updates
#' @return an n by p matrix
wlra = function(x, w, rank, max_iters=1000, atol=1e-3,verbose=FALSE){
  n = nrow(x)
  p = ncol(x)

  # Important: WLRA requires weights 0 <= w <= 1
  w = w/max(w)
  
    # Important: the procedure is deterministic, so initialization
  # matters.
  #
  # Srebro and Jaakkola suggest the best strategy is to initialize
  # from zero, but go from a full rank down to a rank k approximation in
  # the first iterations
  #
  # For now, take the simpler strategy of just initializing to zero. Srebro and
  # Jaakkola suggest this can underfit.
  z = matrix(0,nrow=n,ncol=p)
  obj = mean(w * x^2)
  if(verbose){
    print(paste0("wsvd [0]=",obj))
  }

  for(i in 0:max_iters){
    z1 = lra(w * x + (1 - w) * z, rank)
    update = mean(w * (x - z1)^2)
    if(verbose){print(paste0("wsvd [",i+1,"]=",update))}
    if(update > obj){
      stop("objective increased")
    } else if(max(abs(update-obj))<atol){
      return(z1)
    } else {
      z = z1
      obj = update
    }
  }
  stop("failed to converge")
}

# return log(p(Y|lambda=exp(eta))) for Y \sim Poi(lambda)
pois_llik= function(y, eta){
  sum(dpois(y,exp(eta),log=TRUE))
}

#' @details Assume x_ij ~ Poisson(exp(eta_ij)), eta_ij = L_ik F_kj
#'  Maximize the log likelihood by using Taylor approximation to rewrite the problem as WLRA.
#' @param x: input data (n, p)
#' @param rank: rank of the approximation
#' @param max_outer_iters: maximum number of calls to WLRA
#' @param max_iters: maximum number of EM iterations in WLRA
#' @param verbose: print objective function updates
#' @return low rank approximation (n, p) matrix
pois_lra= function(x, rank, max_outer_iters=50, max_iters=1000, atol=1e-3, verbose=FALSE){
  n = nrow(x)
  p = ncol(x)
  nmf = NNLM::nnmf(x, rank)
  eta = log(nmf$W %*% nmf$H)
  obj = mean(pois_llik(x, eta))
  if(verbose)
    print(paste0("pois_lra [0]:",obj))
  for(i in 0:max_outer_iters){
    lam = exp(eta)
    w = lam
    target = eta + x / lam - 1
    w[is.na(x)]=0 # Mark missing data with weight 0
    
      # Now we can go ahead and fill in the missing values with something
      # computationally convenient, because the WLRA EM update will ignore the
      # value for weight zero.
    target[is.na(x)] = 0
    eta1 = wlra(target, w, rank, max_iters=max_iters, atol=atol, verbose=verbose)
    update = mean(pois_llik(x, eta1))
    if(verbose){
      print(paste0("pois_lra [",i + 1,"]:",update))
    }
    if(max(abs(update-obj))<atol){
      return(list(fit=eta1,w=w,target=target))
    } else {
      eta = eta1
      obj = update
    }
  }
  stop("failed to converge")
}

# this just does a simple TSE about lam and runs wSVD once without any iteration
pois_lra1 = function(x, rank, lam = ifelse(x>0,x,0.5), max_iters=1000, atol=1e-3, verbose=FALSE){
  target = log(lam) + x / lam - 1
  w = lam
  w[is.na(x)]=0 # Mark missing data with weight 0
    
      # Now we can go ahead and fill in the missing values with something
      # computationally convenient, because the WLRA EM update will ignore the
      # value for weight zero.
    target[is.na(x)] = 0
    eta1 = wlra(target, w, rank, max_iters=max_iters, atol=atol, verbose=verbose)
    return(list(fit=eta1,w=w,target=target))
}

First simulate some data:

set.seed(1)
l = rnorm(100)
f = rnorm(100)
eta = outer(l,f)
lambda = exp(eta)
x = matrix(rpois(length(lambda),lambda),nrow=nrow(lambda))

Now fit various models: the plra, plra1 with the “naive” expansion (around x, or 0.5 for x=0) and around the true value of lambda:

x.plra=pois_lra(x,rank=1,verbose=TRUE)
[1] "pois_lra [0]:-56393.8383721883"
[1] "wsvd [0]=9.11187664142563"
[1] "wsvd [1]=9.09767330084035"
[1] "wsvd [2]=9.0942294282996"
[1] "wsvd [3]=9.09252076154906"
[1] "wsvd [4]=9.09146781601123"
[1] "wsvd [5]=9.09073663864168"
[1] "pois_lra [1]:-20730.5265718132"
[1] "wsvd [0]=0.0254032600526617"
[1] "wsvd [1]=0.018730302153258"
[1] "wsvd [2]=0.0173472213350827"
[1] "wsvd [3]=0.0166205649043371"
[1] "pois_lra [2]:-24476.2310334721"
[1] "wsvd [0]=0.0342958650072393"
[1] "wsvd [1]=0.0302492341697698"
[1] "wsvd [2]=0.0295946545829222"
[1] "pois_lra [3]:-26631.5708456328"
[1] "wsvd [0]=0.0594742391599983"
[1] "wsvd [1]=0.0558114666026009"
[1] "wsvd [2]=0.0551636349735505"
[1] "pois_lra [4]:-26775.5968189102"
[1] "wsvd [0]=0.0654476543467692"
[1] "wsvd [1]=0.0617080003798833"
[1] "wsvd [2]=0.0609634116407461"
[1] "pois_lra [5]:-26663.0720352895"
[1] "wsvd [0]=0.0639419959608682"
[1] "wsvd [1]=0.0601941802731442"
[1] "wsvd [2]=0.0594430040900086"
[1] "pois_lra [6]:-26653.5190878185"
[1] "wsvd [0]=0.0637737656022713"
[1] "wsvd [1]=0.0600262104698666"
[1] "wsvd [2]=0.0592754713652051"
[1] "pois_lra [7]:-26653.9093042342"
[1] "wsvd [0]=0.063778632660418"
[1] "wsvd [1]=0.0600311044827268"
[1] "wsvd [2]=0.0592803862918574"
[1] "pois_lra [8]:-26653.9404182592"
[1] "wsvd [0]=0.0637791754208205"
[1] "wsvd [1]=0.0600316465363189"
[1] "wsvd [2]=0.0592809270276872"
[1] "pois_lra [9]:-26653.9393089239"
[1] "wsvd [0]=0.0637791618254631"
[1] "wsvd [1]=0.0600316328646063"
[1] "wsvd [2]=0.0592809132945037"
[1] "pois_lra [10]:-26653.9392199259"
x.plra1.naive=pois_lra1(x,rank=1,verbose=TRUE)
[1] "wsvd [0]=0.0357370934894206"
[1] "wsvd [1]=0.0216973393395743"
[1] "wsvd [2]=0.0183559836838193"
[1] "wsvd [3]=0.0166671902679346"
[1] "wsvd [4]=0.0156058206162763"
[1] "wsvd [5]=0.0148556454379835"
x.plra1.true = pois_lra1(x,1,lam=lambda) 
pois_llik(x,eta) #log likelihood at true value of eta
[1] -13269.62
pois_llik(x,x.plra$fit)
[1] -26653.94
pois_llik(x,x.plra1.naive$fit)
[1] -20770.27
pois_llik(x,x.plra1.true$fit)
[1] -20830.33

This was weird that expansion around the truth was so bad. So I tried more stringent convergence:

x.plra1.true2 = pois_lra1(x,1,lam=lambda,max_iters = 10000,atol=1e-9,verbose=TRUE) 
[1] "wsvd [0]=0.0345547809715475"
[1] "wsvd [1]=0.0207263959799866"
[1] "wsvd [2]=0.0174407325265061"
[1] "wsvd [3]=0.0157745000299274"
[1] "wsvd [4]=0.0147294810989384"
[1] "wsvd [5]=0.0139947047339968"
[1] "wsvd [6]=0.0134405646816002"
[1] "wsvd [7]=0.0130030394359756"
[1] "wsvd [8]=0.0126464909227861"
[1] "wsvd [9]=0.012349219468963"
[1] "wsvd [10]=0.0120970680644167"
[1] "wsvd [11]=0.0118802897240533"
[1] "wsvd [12]=0.0116918827807944"
[1] "wsvd [13]=0.0115266446027147"
[1] "wsvd [14]=0.0113806024851364"
[1] "wsvd [15]=0.0112506546653568"
[1] "wsvd [16]=0.0111343344795743"
[1] "wsvd [17]=0.0110296499093777"
[1] "wsvd [18]=0.0109349710709728"
[1] "wsvd [19]=0.0108489492265877"
[1] "wsvd [20]=0.0107704571407363"
[1] "wsvd [21]=0.0106985442718116"
[1] "wsvd [22]=0.0106324025160935"
[1] "wsvd [23]=0.010571339613379"
[1] "wsvd [24]=0.0105147582173843"
[1] "wsvd [25]=0.0104621392222806"
[1] "wsvd [26]=0.0104130283325444"
[1] "wsvd [27]=0.0103670251352177"
[1] "wsvd [28]=0.010323774124084"
[1] "wsvd [29]=0.0102829572609871"
[1] "wsvd [30]=0.010244287757877"
[1] "wsvd [31]=0.0102075048355588"
[1] "wsvd [32]=0.0101723692691997"
[1] "wsvd [33]=0.0101386595716251"
[1] "wsvd [34]=0.0101061686969389"
[1] "wsvd [35]=0.0100747011715784"
[1] "wsvd [36]=0.0100440705794027"
[1] "wsvd [37]=0.010014097343138"
[1] "wsvd [38]=0.00998460675745698"
[1] "wsvd [39]=0.0099554272398796"
[1] "wsvd [40]=0.00992638877511808"
[1] "wsvd [41]=0.00989732153685665"
[1] "wsvd [42]=0.00986805467858832"
[1] "wsvd [43]=0.00983841529224545"
[1] "wsvd [44]=0.00980822754011714"
[1] "wsvd [45]=0.00977731197201945"
[1] "wsvd [46]=0.00974548504587779"
[1] "wsvd [47]=0.00971255887571109"
[1] "wsvd [48]=0.00967834123630306"
[1] "wsvd [49]=0.00964263585833216"
[1] "wsvd [50]=0.00960524305102644"
[1] "wsvd [51]=0.00956596069101885"
[1] "wsvd [52]=0.00952458561540583"
[1] "wsvd [53]=0.0094809154533755"
[1] "wsvd [54]=0.00943475092345007"
[1] "wsvd [55]=0.00938589861167933"
[1] "wsvd [56]=0.00933417422944722"
[1] "wsvd [57]=0.00927940632756793"
[1] "wsvd [58]=0.00922144041609716"
[1] "wsvd [59]=0.00916014340735611"
[1] "wsvd [60]=0.00909540826434175"
[1] "wsvd [61]=0.00902715870005596"
[1] "wsvd [62]=0.00895535373823257"
[1] "wsvd [63]=0.00887999191611026"
[1] "wsvd [64]=0.00880111488941361"
[1] "wsvd [65]=0.0087188101927575"
[1] "wsvd [66]=0.00863321291898988"
[1] "wsvd [67]=0.00854450611109585"
[1] "wsvd [68]=0.00845291971097561"
[1] "wsvd [69]=0.00835872797914897"
[1] "wsvd [70]=0.00826224538415828"
[1] "wsvd [71]=0.00816382105367115"
[1] "wsvd [72]=0.00806383197274062"
[1] "wsvd [73]=0.00796267519922755"
[1] "wsvd [74]=0.00786075943327471"
[1] "wsvd [75]=0.00775849631983691"
[1] "wsvd [76]=0.00765629187623029"
[1] "wsvd [77]=0.00755453841945491"
[1] "wsvd [78]=0.00745360732314718"
[1] "wsvd [79]=0.00735384286694553"
[1] "wsvd [80]=0.00725555735945015"
[1] "wsvd [81]=0.00715902762842662"
[1] "wsvd [82]=0.00706449288680204"
[1] "wsvd [83]=0.00697215390741367"
[1] "wsvd [84]=0.00688217337851253"
[1] "wsvd [85]=0.00679467726853902"
[1] "wsvd [86]=0.00670975700330715"
[1] "wsvd [87]=0.00662747225023184"
[1] "wsvd [88]=0.00654785411004353"
[1] "wsvd [89]=0.00647090853322582"
[1] "wsvd [90]=0.0063966198026211"
[1] "wsvd [91]=0.00632495395190205"
[1] "wsvd [92]=0.00625586201901082"
[1] "wsvd [93]=0.00618928306196224"
[1] "wsvd [94]=0.00612514689001128"
[1] "wsvd [95]=0.00606337648514052"
[1] "wsvd [96]=0.00600389010671461"
[1] "wsvd [97]=0.00594660308596605"
[1] "wsvd [98]=0.00589142932700378"
[1] "wsvd [99]=0.00583828253773391"
[1] "wsvd [100]=0.00578707721799839"
[1] "wsvd [101]=0.00573772943393736"
[1] "wsvd [102]=0.00569015740760073"
[1] "wsvd [103]=0.00564428194964958"
[1] "wsvd [104]=0.00560002676100251"
[1] "wsvd [105]=0.00555731862682069"
[1] "wsvd [106]=0.00551608752354144"
[1] "wsvd [107]=0.0054762666569499"
[1] "wsvd [108]=0.0054377924466528"
[1] "wsvd [109]=0.00540060446987247"
[1] "wsvd [110]=0.00536464537526318"
[1] "wsvd [111]=0.00532986077548909"
[1] "wsvd [112]=0.00529619912559655"
[1] "wsvd [113]=0.005263611592755"
[1] "wsvd [114]=0.0052320519217122"
[1] "wsvd [115]=0.0052014762992891"
[1] "wsvd [116]=0.00517184322040291"
[1] "wsvd [117]=0.00514311335742961"
[1] "wsvd [118]=0.00511524943417583"
[1] "wsvd [119]=0.00508821610530338"
[1] "wsvd [120]=0.00506197984171756"
[1] "wsvd [121]=0.00503650882217724"
[1] "wsvd [122]=0.00501177283119364"
[1] "wsvd [123]=0.00498774316314558"
[1] "wsvd [124]=0.00496439253243934"
[1] "wsvd [125]=0.00494169498947336"
[1] "wsvd [126]=0.0049196258421246"
[1] "wsvd [127]=0.00489816158244848"
[1] "wsvd [128]=0.00487727981827311"
[1] "wsvd [129]=0.00485695920936781"
[1] "wsvd [130]=0.00483717940787239"
[1] "wsvd [131]=0.00481792100268496"
[1] "wsvd [132]=0.00479916546752084"
[1] "wsvd [133]=0.00478089511237186"
[1] "wsvd [134]=0.0047630930381128"
[1] "wsvd [135]=0.00474574309401977"
[1] "wsvd [136]=0.00472882983798274"
[1] "wsvd [137]=0.00471233849921164"
[1] "wsvd [138]=0.00469625494325128"
[1] "wsvd [139]=0.00468056563913587"
[1] "wsvd [140]=0.00466525762852754"
[1] "wsvd [141]=0.0046503184966968"
[1] "wsvd [142]=0.00463573634521435"
[1] "wsvd [143]=0.00462149976623485"
[1] "wsvd [144]=0.00460759781826321"
[1] "wsvd [145]=0.00459402000330287"
[1] "wsvd [146]=0.00458075624529387"
[1] "wsvd [147]=0.004567796869756"
[1] "wsvd [148]=0.00455513258455877"
[1] "wsvd [149]=0.00454275446174662"
[1] "wsvd [150]=0.00453065392035257"
[1] "wsvd [151]=0.00451882271013928"
[1] "wsvd [152]=0.00450725289621047"
[1] "wsvd [153]=0.00449593684444013"
[1] "wsvd [154]=0.00448486720767052"
[1] "wsvd [155]=0.00447403691263356"
[1] "wsvd [156]=0.00446343914755314"
[1] "wsvd [157]=0.00445306735038899"
[1] "wsvd [158]=0.00444291519768516"
[1] "wsvd [159]=0.0044329765939886"
[1] "wsvd [160]=0.00442324566180582"
[1] "wsvd [161]=0.00441371673206721"
[1] "wsvd [162]=0.00440438433507084"
[1] "wsvd [163]=0.0043952431918792"
[1] "wsvd [164]=0.00438628820614386"
[1] "wsvd [165]=0.00437751445633477"
[1] "wsvd [166]=0.00436891718835201"
[1] "wsvd [167]=0.00436049180849943"
[1] "wsvd [168]=0.00435223387680052"
[1] "wsvd [169]=0.00434413910063832"
[1] "wsvd [170]=0.00433620332870181"
[1] "wsvd [171]=0.00432842254522281"
[1] "wsvd [172]=0.00432079286448764"
[1] "wsvd [173]=0.00431331052560928"
[1] "wsvd [174]=0.00430597188754628"
[1] "wsvd [175]=0.00429877342435543"
[1] "wsvd [176]=0.00429171172066608"
[1] "wsvd [177]=0.00428478346736455"
[1] "wsvd [178]=0.00427798545747773"
[1] "wsvd [179]=0.00427131458224563"
[1] "wsvd [180]=0.00426476782737319"
[1] "wsvd [181]=0.00425834226945212"
[1] "wsvd [182]=0.00425203507254413"
[1] "wsvd [183]=0.00424584348491733"
[1] "wsvd [184]=0.00423976483592813"
[1] "wsvd [185]=0.00423379653304117"
[1] "wsvd [186]=0.00422793605898048"
[1] "wsvd [187]=0.00422218096900524"
[1] "wsvd [188]=0.00421652888830404"
[1] "wsvd [189]=0.00421097750950155"
[1] "wsvd [190]=0.00420552459027236"
[1] "wsvd [191]=0.00420016795105644"
[1] "wsvd [192]=0.00419490547287146"
[1] "wsvd [193]=0.00418973509521716"
[1] "wsvd [194]=0.00418465481406731"
[1] "wsvd [195]=0.00417966267994508"
[1] "wsvd [196]=0.00417475679607782"
[1] "wsvd [197]=0.00416993531662737"
[1] "wsvd [198]=0.0041651964449925"
[1] "wsvd [199]=0.00416053843217985"
[1] "wsvd [200]=0.00415595957524034"
[1] "wsvd [201]=0.00415145821576786"
[1] "wsvd [202]=0.00414703273845744"
[1] "wsvd [203]=0.00414268156972007"
[1] "wsvd [204]=0.00413840317635165"
[1] "wsvd [205]=0.00413419606425348"
[1] "wsvd [206]=0.0041300587772021"
[1] "wsvd [207]=0.00412598989566609"
[1] "wsvd [208]=0.00412198803566783"
[1] "wsvd [209]=0.00411805184768819"
[1] "wsvd [210]=0.00411418001561219"
[1] "wsvd [211]=0.0041103712557139"
[1] "wsvd [212]=0.00410662431567878"
[1] "wsvd [213]=0.00410293797366188"
[1] "wsvd [214]=0.00409931103738037"
[1] "wsvd [215]=0.00409574234323881"
[1] "wsvd [216]=0.00409223075548585"
[1] "wsvd [217]=0.00408877516540107"
[1] "wsvd [218]=0.0040853744905105"
[1] "wsvd [219]=0.00408202767382982"
[1] "wsvd [220]=0.00407873368313396"
[1] "wsvd [221]=0.00407549151025205"
[1] "wsvd [222]=0.00407230017038665"
[1] "wsvd [223]=0.00406915870145627"
[1] "wsvd [224]=0.00406606616346029"
[1] "wsvd [225]=0.00406302163786529"
[1] "wsvd [226]=0.00406002422701191"
[1] "wsvd [227]=0.00405707305354163"
[1] "wsvd [228]=0.00405416725984232"
[1] "wsvd [229]=0.00405130600751221"
[1] "wsvd [230]=0.00404848847684126"
[1] "wsvd [231]=0.00404571386630938"
[1] "wsvd [232]=0.0040429813921009"
[1] "wsvd [233]=0.00404029028763453"
[1] "wsvd [234]=0.00403763980310828"
[1] "wsvd [235]=0.0040350292050589"
[1] "wsvd [236]=0.00403245777593502"
[1] "wsvd [237]=0.00402992481368371"
[1] "wsvd [238]=0.00402742963134994"
[1] "wsvd [239]=0.00402497155668828"
[1] "wsvd [240]=0.0040225499317866"
[1] "wsvd [241]=0.00402016411270127"
[1] "wsvd [242]=0.00401781346910336"
[1] "wsvd [243]=0.00401549738393553"
[1] "wsvd [244]=0.00401321525307922"
[1] "wsvd [245]=0.00401096648503178"
[1] "wsvd [246]=0.00400875050059306"
[1] "wsvd [247]=0.00400656673256141"
[1] "wsvd [248]=0.00400441462543841"
[1] "wsvd [249]=0.00400229363514234"
[1] "wsvd [250]=0.00400020322872984"
[1] "wsvd [251]=0.00399814288412569"
[1] "wsvd [252]=0.00399611208986031"
[1] "wsvd [253]=0.00399411034481466"
[1] "wsvd [254]=0.0039921371579725"
[1] "wsvd [255]=0.0039901920481796"
[1] "wsvd [256]=0.0039882745439096"
[1] "wsvd [257]=0.00398638418303661"
[1] "wsvd [258]=0.00398452051261395"
[1] "wsvd [259]=0.00398268308865908"
[1] "wsvd [260]=0.00398087147594446"
[1] "wsvd [261]=0.00397908524779412"
[1] "wsvd [262]=0.0039773239858857"
[1] "wsvd [263]=0.00397558728005802"
[1] "wsvd [264]=0.00397387472812364"
[1] "wsvd [265]=0.00397218593568663"
[1] "wsvd [266]=0.00397052051596509"
[1] "wsvd [267]=0.00396887808961845"
[1] "wsvd [268]=0.00396725828457931"
[1] "wsvd [269]=0.00396566073588976"
[1] "wsvd [270]=0.00396408508554192"
[1] "wsvd [271]=0.00396253098232272"
[1] "wsvd [272]=0.00396099808166269"
[1] "wsvd [273]=0.00395948604548862"
[1] "wsvd [274]=0.00395799454208011"
[1] "wsvd [275]=0.00395652324592973"
[1] "wsvd [276]=0.00395507183760675"
[1] "wsvd [277]=0.00395364000362443"
[1] "wsvd [278]=0.00395222743631054"
[1] "wsvd [279]=0.00395083383368128"
[1] "wsvd [280]=0.00394945889931825"
[1] "wsvd [281]=0.00394810234224864"
[1] "wsvd [282]=0.00394676387682824"
[1] "wsvd [283]=0.0039454432226275"
[1] "wsvd [284]=0.00394414010432037"
[1] "wsvd [285]=0.00394285425157581"
[1] "wsvd [286]=0.00394158539895207"
[1] "wsvd [287]=0.00394033328579347"
[1] "wsvd [288]=0.00393909765612977"
[1] "wsvd [289]=0.00393787825857792"
[1] "wsvd [290]=0.00393667484624622"
[1] "wsvd [291]=0.00393548717664077"
[1] "wsvd [292]=0.00393431501157418"
[1] "wsvd [293]=0.00393315811707651"
[1] "wsvd [294]=0.00393201626330821"
[1] "wsvd [295]=0.00393088922447527"
[1] "wsvd [296]=0.00392977677874625"
[1] "wsvd [297]=0.0039286787081714"
[1] "wsvd [298]=0.00392759479860356"
[1] "wsvd [299]=0.00392652483962099"
[1] "wsvd [300]=0.00392546862445197"
[1] "wsvd [301]=0.00392442594990114"
[1] "wsvd [302]=0.0039233966162776"
[1] "wsvd [303]=0.00392238042732465"
[1] "wsvd [304]=0.0039213771901511"
[1] "wsvd [305]=0.0039203867151642"
[1] "wsvd [306]=0.00391940881600416"
[1] "wsvd [307]=0.00391844330948003"
[1] "wsvd [308]=0.00391749001550716"
[1] "wsvd [309]=0.00391654875704602"
[1] "wsvd [310]=0.00391561936004241"
[1] "wsvd [311]=0.00391470165336903"
[1] "wsvd [312]=0.00391379546876832"
[1] "wsvd [313]=0.00391290064079668"
[1] "wsvd [314]=0.00391201700676978"
[1] "wsvd [315]=0.00391114440670923"
[1] "wsvd [316]=0.00391028268329038"
[1] "wsvd [317]=0.00390943168179122"
[1] "wsvd [318]=0.00390859125004251"
[1] "wsvd [319]=0.00390776123837892"
[1] "wsvd [320]=0.00390694149959127"
[1] "wsvd [321]=0.00390613188887986"
[1] "wsvd [322]=0.0039053322638087"
[1] "wsvd [323]=0.00390454248426088"
[1] "wsvd [324]=0.00390376241239479"
[1] "wsvd [325]=0.00390299191260134"
[1] "wsvd [326]=0.00390223085146209"
[1] "wsvd [327]=0.00390147909770826"
[1] "wsvd [328]=0.00390073652218069"
[1] "wsvd [329]=0.00390000299779055"
[1] "wsvd [330]=0.00389927839948093"
[1] "wsvd [331]=0.0038985626041893"
[1] "wsvd [332]=0.0038978554908107"
[1] "wsvd [333]=0.00389715694016168"
[1] "wsvd [334]=0.00389646683494512"
[1] "wsvd [335]=0.00389578505971563"
[1] "wsvd [336]=0.00389511150084584"
[1] "wsvd [337]=0.00389444604649324"
[1] "wsvd [338]=0.0038937885865678"
[1] "wsvd [339]=0.00389313901270026"
[1] "wsvd [340]=0.00389249721821106"
[1] "wsvd [341]=0.00389186309807986"
[1] "wsvd [342]=0.00389123654891582"
[1] "wsvd [343]=0.00389061746892833"
[1] "wsvd [344]=0.00389000575789846"
[1] "wsvd [345]=0.00388940131715096"
[1] "wsvd [346]=0.00388880404952678"
[1] "wsvd [347]=0.00388821385935623"
[1] "wsvd [348]=0.00388763065243264"
[1] "wsvd [349]=0.00388705433598651"
[1] "wsvd [350]=0.0038864848186603"
[1] "wsvd [351]=0.00388592201048362"
[1] "wsvd [352]=0.00388536582284892"
[1] "wsvd [353]=0.00388481616848775"
[1] "wsvd [354]=0.00388427296144744"
[1] "wsvd [355]=0.00388373611706821"
[1] "wsvd [356]=0.00388320555196078"
[1] "wsvd [357]=0.00388268118398446"
[1] "wsvd [358]=0.00388216293222557"
[1] "wsvd [359]=0.00388165071697637"
[1] "wsvd [360]=0.00388114445971437"
[1] "wsvd [361]=0.00388064408308207"
[1] "wsvd [362]=0.00388014951086704"
[1] "wsvd [363]=0.00387966066798245"
[1] "wsvd [364]=0.00387917748044798"
[1] "wsvd [365]=0.00387869987537103"
[1] "wsvd [366]=0.00387822778092835"
[1] "wsvd [367]=0.00387776112634804"
[1] "wsvd [368]=0.00387729984189186"
[1] "wsvd [369]=0.00387684385883786"
[1] "wsvd [370]=0.00387639310946342"
[1] "wsvd [371]=0.00387594752702854"
[1] "wsvd [372]=0.00387550704575948"
[1] "wsvd [373]=0.00387507160083272"
[1] "wsvd [374]=0.00387464112835917"
[1] "wsvd [375]=0.0038742155653688"
[1] "wsvd [376]=0.0038737948497954"
[1] "wsvd [377]=0.00387337892046176"
[1] "wsvd [378]=0.00387296771706505"
[1] "wsvd [379]=0.00387256118016251"
[1] "wsvd [380]=0.00387215925115741"
[1] "wsvd [381]=0.00387176187228523"
[1] "wsvd [382]=0.00387136898660018"
[1] "wsvd [383]=0.00387098053796184"
[1] "wsvd [384]=0.00387059647102221"
[1] "wsvd [385]=0.00387021673121285"
[1] "wsvd [386]=0.00386984126473236"
[1] "wsvd [387]=0.00386947001853403"
[1] "wsvd [388]=0.00386910294031376"
[1] "wsvd [389]=0.00386873997849814"
[1] "wsvd [390]=0.00386838108223284"
[1] "wsvd [391]=0.00386802620137113"
[1] "wsvd [392]=0.00386767528646263"
[1] "wsvd [393]=0.00386732828874229"
[1] "wsvd [394]=0.00386698516011953"
[1] "wsvd [395]=0.00386664585316763"
[1] "wsvd [396]=0.00386631032111324"
[1] "wsvd [397]=0.00386597851782614"
[1] "wsvd [398]=0.00386565039780914"
[1] "wsvd [399]=0.00386532591618822"
[1] "wsvd [400]=0.00386500502870276"
[1] "wsvd [401]=0.00386468769169599"
[1] "wsvd [402]=0.00386437386210568"
[1] "wsvd [403]=0.00386406349745482"
[1] "wsvd [404]=0.00386375655584265"
[1] "wsvd [405]=0.00386345299593573"
[1] "wsvd [406]=0.00386315277695923"
[1] "wsvd [407]=0.00386285585868829"
[1] "wsvd [408]=0.00386256220143967"
[1] "wsvd [409]=0.00386227176606341"
[1] "wsvd [410]=0.00386198451393469"
[1] "wsvd [411]=0.00386170040694587"
[1] "wsvd [412]=0.00386141940749856"
[1] "wsvd [413]=0.003861141478496"
[1] "wsvd [414]=0.00386086658333534"
[1] "wsvd [415]=0.00386059468590031"
[1] "wsvd [416]=0.0038603257505538"
[1] "wsvd [417]=0.00386005974213067"
[1] "wsvd [418]=0.0038597966259307"
[1] "wsvd [419]=0.0038595363677116"
[1] "wsvd [420]=0.00385927893368217"
[1] "wsvd [421]=0.00385902429049559"
[1] "wsvd [422]=0.00385877240524277"
[1] "wsvd [423]=0.00385852324544591"
[1] "wsvd [424]=0.00385827677905205"
[1] "wsvd [425]=0.00385803297442681"
[1] "wsvd [426]=0.00385779180034824"
[1] "wsvd [427]=0.00385755322600069"
[1] "wsvd [428]=0.00385731722096889"
[1] "wsvd [429]=0.00385708375523204"
[1] "wsvd [430]=0.00385685279915805"
[1] "wsvd [431]=0.0038566243234979"
[1] "wsvd [432]=0.00385639829937996"
[1] "wsvd [433]=0.00385617469830461"
[1] "wsvd [434]=0.00385595349213877"
[1] "wsvd [435]=0.00385573465311058"
[1] "wsvd [436]=0.00385551815380424"
[1] "wsvd [437]=0.0038553039671548"
[1] "wsvd [438]=0.00385509206644316"
[1] "wsvd [439]=0.00385488242529107"
[1] "wsvd [440]=0.00385467501765623"
[1] "wsvd [441]=0.00385446981782751"
[1] "wsvd [442]=0.00385426680042019"
[1] "wsvd [443]=0.00385406594037132"
[1] "wsvd [444]=0.00385386721293514"
[1] "wsvd [445]=0.00385367059367855"
[1] "wsvd [446]=0.00385347605847669"
[1] "wsvd [447]=0.00385328358350861"
[1] "wsvd [448]=0.00385309314525289"
[1] "wsvd [449]=0.00385290472048353"
[1] "wsvd [450]=0.00385271828626569"
[1] "wsvd [451]=0.00385253381995165"
[1] "wsvd [452]=0.00385235129917679"
[1] "wsvd [453]=0.00385217070185561"
[1] "wsvd [454]=0.00385199200617782"
[1] "wsvd [455]=0.00385181519060454"
[1] "wsvd [456]=0.00385164023386448"
[1] "wsvd [457]=0.00385146711495025"
[1] "wsvd [458]=0.00385129581311469"
[1] "wsvd [459]=0.00385112630786728"
[1] "wsvd [460]=0.00385095857897059"
[1] "wsvd [461]=0.00385079260643675"
[1] "wsvd [462]=0.0038506283705241"
[1] "wsvd [463]=0.00385046585173371"
[1] "wsvd [464]=0.00385030503080614"
[1] "wsvd [465]=0.00385014588871808"
[1] "wsvd [466]=0.00384998840667919"
[1] "wsvd [467]=0.00384983256612886"
[1] "wsvd [468]=0.00384967834873316"
[1] "wsvd [469]=0.00384952573638167"
[1] "wsvd [470]=0.0038493747111845"
[1] "wsvd [471]=0.0038492252554693"
[1] "wsvd [472]=0.00384907735177831"
[1] "wsvd [473]=0.00384893098286545"
[1] "wsvd [474]=0.00384878613169349"
[1] "wsvd [475]=0.00384864278143123"
[1] "wsvd [476]=0.00384850091545071"
[1] "wsvd [477]=0.00384836051732454"
[1] "wsvd [478]=0.00384822157082316"
[1] "wsvd [479]=0.00384808405991223"
[1] "wsvd [480]=0.00384794796875001"
[1] "wsvd [481]=0.0038478132816848"
[1] "wsvd [482]=0.00384767998325242"
[1] "wsvd [483]=0.00384754805817371"
[1] "wsvd [484]=0.00384741749135208"
[1] "wsvd [485]=0.0038472882678711"
[1] "wsvd [486]=0.00384716037299212"
[1] "wsvd [487]=0.00384703379215192"
[1] "wsvd [488]=0.00384690851096041"
[1] "wsvd [489]=0.00384678451519835"
[1] "wsvd [490]=0.00384666179081511"
[1] "wsvd [491]=0.00384654032392644"
[1] "wsvd [492]=0.00384642010081234"
[1] "wsvd [493]=0.00384630110791487"
[1] "wsvd [494]=0.00384618333183606"
[1] "wsvd [495]=0.00384606675933582"
[1] "wsvd [496]=0.0038459513773299"
[1] "wsvd [497]=0.00384583717288785"
[1] "wsvd [498]=0.00384572413323105"
[1] "wsvd [499]=0.00384561224573071"
[1] "wsvd [500]=0.00384550149790597"
[1] "wsvd [501]=0.00384539187742199"
[1] "wsvd [502]=0.00384528337208805"
[1] "wsvd [503]=0.00384517596985571"
[1] "wsvd [504]=0.00384506965881697"
[1] "wsvd [505]=0.0038449644272025"
[1] "wsvd [506]=0.00384486026337983"
[1] "wsvd [507]=0.00384475715585163"
[1] "wsvd [508]=0.00384465509325396"
[1] "wsvd [509]=0.00384455406435461"
[1] "wsvd [510]=0.00384445405805138"
[1] "wsvd [511]=0.00384435506337044"
[1] "wsvd [512]=0.00384425706946473"
[1] "wsvd [513]=0.00384416006561232"
[1] "wsvd [514]=0.00384406404121485"
[1] "wsvd [515]=0.00384396898579595"
[1] "wsvd [516]=0.00384387488899974"
[1] "wsvd [517]=0.00384378174058926"
[1] "wsvd [518]=0.003843689530445"
[1] "wsvd [519]=0.00384359824856346"
[1] "wsvd [520]=0.00384350788505562"
[1] "wsvd [521]=0.00384341843014557"
[1] "wsvd [522]=0.00384332987416909"
[1] "wsvd [523]=0.00384324220757222"
[1] "wsvd [524]=0.00384315542090991"
[1] "wsvd [525]=0.00384306950484467"
[1] "wsvd [526]=0.00384298445014522"
[1] "wsvd [527]=0.00384290024768516"
[1] "wsvd [528]=0.00384281688844172"
[1] "wsvd [529]=0.00384273436349442"
[1] "wsvd [530]=0.00384265266402382"
[1] "wsvd [531]=0.0038425717813103"
[1] "wsvd [532]=0.0038424917067328"
[1] "wsvd [533]=0.00384241243176761"
[1] "wsvd [534]=0.00384233394798719"
[1] "wsvd [535]=0.00384225624705899"
[1] "wsvd [536]=0.00384217932074424"
[1] "wsvd [537]=0.00384210316089686"
[1] "wsvd [538]=0.0038420277594623"
[1] "wsvd [539]=0.0038419531084764"
[1] "wsvd [540]=0.00384187920006435"
[1] "wsvd [541]=0.00384180602643954"
[1] "wsvd [542]=0.00384173357990249"
[1] "wsvd [543]=0.00384166185283984"
[1] "wsvd [544]=0.00384159083772326"
[1] "wsvd [545]=0.00384152052710845"
[1] "wsvd [546]=0.00384145091363408"
[1] "wsvd [547]=0.00384138199002083"
[1] "wsvd [548]=0.00384131374907039"
[1] "wsvd [549]=0.00384124618366447"
[1] "wsvd [550]=0.00384117928676383"
[1] "wsvd [551]=0.00384111305140736"
[1] "wsvd [552]=0.00384104747071112"
[1] "wsvd [553]=0.00384098253786741"
[1] "wsvd [554]=0.00384091824614386"
[1] "wsvd [555]=0.00384085458888256"
[1] "wsvd [556]=0.0038407915594991"
[1] "wsvd [557]=0.00384072915148176"
[1] "wsvd [558]=0.00384066735839061"
[1] "wsvd [559]=0.00384060617385666"
[1] "wsvd [560]=0.00384054559158102"
[1] "wsvd [561]=0.00384048560533406"
[1] "wsvd [562]=0.00384042620895461"
[1] "wsvd [563]=0.00384036739634909"
[1] "wsvd [564]=0.0038403091614908"
[1] "wsvd [565]=0.00384025149841905"
[1] "wsvd [566]=0.00384019440123839"
[1] "wsvd [567]=0.00384013786411791"
[1] "wsvd [568]=0.00384008188129039"
[1] "wsvd [569]=0.00384002644705159"
[1] "wsvd [570]=0.00383997155575952"
[1] "wsvd [571]=0.0038399172018337"
[1] "wsvd [572]=0.00383986337975441"
[1] "wsvd [573]=0.00383981008406203"
[1] "wsvd [574]=0.00383975730935629"
[1] "wsvd [575]=0.00383970505029561"
[1] "wsvd [576]=0.00383965330159638"
[1] "wsvd [577]=0.00383960205803231"
[1] "wsvd [578]=0.00383955131443377"
[1] "wsvd [579]=0.00383950106568708"
[1] "wsvd [580]=0.00383945130673393"
[1] "wsvd [581]=0.00383940203257067"
[1] "wsvd [582]=0.00383935323824773"
[1] "wsvd [583]=0.00383930491886896"
[1] "wsvd [584]=0.00383925706959101"
[1] "wsvd [585]=0.00383920968562275"
[1] "wsvd [586]=0.00383916276222462"
[1] "wsvd [587]=0.00383911629470808"
[1] "wsvd [588]=0.003839070278435"
[1] "wsvd [589]=0.00383902470881707"
[1] "wsvd [590]=0.00383897958131526"
[1] "wsvd [591]=0.0038389348914392"
[1] "wsvd [592]=0.00383889063474669"
[1] "wsvd [593]=0.00383884680684307"
[1] "wsvd [594]=0.00383880340338075"
[1] "wsvd [595]=0.00383876042005861"
[1] "wsvd [596]=0.00383871785262152"
[1] "wsvd [597]=0.00383867569685976"
[1] "wsvd [598]=0.00383863394860856"
[1] "wsvd [599]=0.00383859260374753"
[1] "wsvd [600]=0.00383855165820021"
[1] "wsvd [601]=0.00383851110793353"
[1] "wsvd [602]=0.00383847094895732"
[1] "wsvd [603]=0.00383843117732384"
[1] "wsvd [604]=0.00383839178912729"
[1] "wsvd [605]=0.00383835278050335"
[1] "wsvd [606]=0.00383831414762867"
[1] "wsvd [607]=0.00383827588672045"
[1] "wsvd [608]=0.00383823799403596"
[1] "wsvd [609]=0.00383820046587207"
[1] "wsvd [610]=0.00383816329856486"
[1] "wsvd [611]=0.00383812648848911"
[1] "wsvd [612]=0.00383809003205791"
[1] "wsvd [613]=0.00383805392572223"
[1] "wsvd [614]=0.00383801816597045"
[1] "wsvd [615]=0.00383798274932799"
[1] "wsvd [616]=0.00383794767235687"
[1] "wsvd [617]=0.00383791293165528"
[1] "wsvd [618]=0.00383787852385725"
[1] "wsvd [619]=0.00383784444563214"
[1] "wsvd [620]=0.00383781069368434"
[1] "wsvd [621]=0.00383777726475283"
[1] "wsvd [622]=0.00383774415561081"
[1] "wsvd [623]=0.00383771136306529"
[1] "wsvd [624]=0.00383767888395677"
[1] "wsvd [625]=0.00383764671515882"
[1] "wsvd [626]=0.00383761485357772"
[1] "wsvd [627]=0.00383758329615212"
[1] "wsvd [628]=0.00383755203985264"
[1] "wsvd [629]=0.00383752108168156"
[1] "wsvd [630]=0.00383749041867244"
[1] "wsvd [631]=0.00383746004788978"
[1] "wsvd [632]=0.00383742996642869"
[1] "wsvd [633]=0.00383740017141453"
[1] "wsvd [634]=0.0038373706600026"
[1] "wsvd [635]=0.00383734142937778"
[1] "wsvd [636]=0.00383731247675424"
[1] "wsvd [637]=0.00383728379937508"
[1] "wsvd [638]=0.00383725539451205"
[1] "wsvd [639]=0.00383722725946521"
[1] "wsvd [640]=0.0038371993915626"
[1] "wsvd [641]=0.00383717178815999"
[1] "wsvd [642]=0.00383714444664052"
[1] "wsvd [643]=0.00383711736441444"
[1] "wsvd [644]=0.00383709053891879"
[1] "wsvd [645]=0.00383706396761709"
[1] "wsvd [646]=0.00383703764799912"
[1] "wsvd [647]=0.00383701157758054"
[1] "wsvd [648]=0.00383698575390267"
[1] "wsvd [649]=0.00383696017453221"
[1] "wsvd [650]=0.00383693483706093"
[1] "wsvd [651]=0.00383690973910542"
[1] "wsvd [652]=0.0038368848783068"
[1] "wsvd [653]=0.00383686025233048"
[1] "wsvd [654]=0.0038368358588659"
[1] "wsvd [655]=0.00383681169562622"
[1] "wsvd [656]=0.00383678776034811"
[1] "wsvd [657]=0.00383676405079149"
[1] "wsvd [658]=0.00383674056473925"
[1] "wsvd [659]=0.00383671729999702"
[1] "wsvd [660]=0.00383669425439295"
[1] "wsvd [661]=0.00383667142577742"
[1] "wsvd [662]=0.00383664881202281"
[1] "wsvd [663]=0.0038366264110233"
[1] "wsvd [664]=0.00383660422069459"
[1] "wsvd [665]=0.00383658223897368"
[1] "wsvd [666]=0.00383656046381867"
[1] "wsvd [667]=0.00383653889320849"
[1] "wsvd [668]=0.0038365175251427"
[1] "wsvd [669]=0.00383649635764127"
[1] "wsvd [670]=0.00383647538874435"
[1] "wsvd [671]=0.00383645461651206"
[1] "wsvd [672]=0.00383643403902427"
[1] "wsvd [673]=0.0038364136543804"
[1] "wsvd [674]=0.0038363934606992"
[1] "wsvd [675]=0.00383637345611853"
[1] "wsvd [676]=0.00383635363879521"
[1] "wsvd [677]=0.00383633400690475"
[1] "wsvd [678]=0.00383631455864119"
[1] "wsvd [679]=0.00383629529221689"
[1] "wsvd [680]=0.00383627620586233"
[1] "wsvd [681]=0.00383625729782595"
[1] "wsvd [682]=0.0038362385663739"
[1] "wsvd [683]=0.0038362200097899"
[1] "wsvd [684]=0.00383620162637504"
[1] "wsvd [685]=0.0038361834144476"
[1] "wsvd [686]=0.00383616537234284"
[1] "wsvd [687]=0.00383614749841286"
[1] "wsvd [688]=0.00383612979102639"
[1] "wsvd [689]=0.00383611224856863"
[1] "wsvd [690]=0.00383609486944107"
[1] "wsvd [691]=0.00383607765206133"
[1] "wsvd [692]=0.00383606059486295"
[1] "wsvd [693]=0.00383604369629528"
[1] "wsvd [694]=0.00383602695482327"
[1] "wsvd [695]=0.00383601036892733"
[1] "wsvd [696]=0.00383599393710313"
[1] "wsvd [697]=0.00383597765786149"
[1] "wsvd [698]=0.00383596152972818"
[1] "wsvd [699]=0.00383594555124378"
[1] "wsvd [700]=0.00383592972096354"
[1] "wsvd [701]=0.00383591403745718"
[1] "wsvd [702]=0.00383589849930878"
[1] "wsvd [703]=0.00383588310511661"
[1] "wsvd [704]=0.003835867853493"
[1] "wsvd [705]=0.00383585274306415"
[1] "wsvd [706]=0.00383583777247004"
[1] "wsvd [707]=0.00383582294036425"
[1] "wsvd [708]=0.00383580824541382"
[1] "wsvd [709]=0.00383579368629912"
[1] "wsvd [710]=0.00383577926171371"
[1] "wsvd [711]=0.00383576497036418"
[1] "wsvd [712]=0.00383575081097007"
[1] "wsvd [713]=0.00383573678226365"
[1] "wsvd [714]=0.00383572288298988"
[1] "wsvd [715]=0.00383570911190619"
[1] "wsvd [716]=0.00383569546778241"
[1] "wsvd [717]=0.00383568194940065"
[1] "wsvd [718]=0.00383566855555509"
[1] "wsvd [719]=0.00383565528505195"
[1] "wsvd [720]=0.00383564213670932"
[1] "wsvd [721]=0.00383562910935702"
[1] "wsvd [722]=0.00383561620183653"
[1] "wsvd [723]=0.00383560341300081"
[1] "wsvd [724]=0.00383559074171422"
[1] "wsvd [725]=0.0038355781868524"
[1] "wsvd [726]=0.00383556574730214"
[1] "wsvd [727]=0.00383555342196125"
[1] "wsvd [728]=0.0038355412097385"
[1] "wsvd [729]=0.00383552910955344"
[1] "wsvd [730]=0.00383551712033633"
[1] "wsvd [731]=0.00383550524102802"
[1] "wsvd [732]=0.00383549347057985"
[1] "wsvd [733]=0.00383548180795352"
[1] "wsvd [734]=0.00383547025212101"
[1] "wsvd [735]=0.00383545880206443"
[1] "wsvd [736]=0.00383544745677597"
[1] "wsvd [737]=0.00383543621525776"
[1] "wsvd [738]=0.00383542507652179"
[1] "wsvd [739]=0.00383541403958978"
[1] "wsvd [740]=0.00383540310349312"
[1] "wsvd [741]=0.00383539226727272"
[1] "wsvd [742]=0.00383538152997897"
[1] "wsvd [743]=0.0038353708906716"
[1] "wsvd [744]=0.00383536034841958"
[1] "wsvd [745]=0.00383534990230106"
[1] "wsvd [746]=0.00383533955140327"
[1] "wsvd [747]=0.00383532929482239"
[1] "wsvd [748]=0.0038353191316635"
[1] "wsvd [749]=0.00383530906104048"
[1] "wsvd [750]=0.00383529908207588"
[1] "wsvd [751]=0.00383528919390089"
[1] "wsvd [752]=0.00383527939565524"
[1] "wsvd [753]=0.00383526968648706"
[1] "wsvd [754]=0.00383526006555287"
[1] "wsvd [755]=0.00383525053201742"
[1] "wsvd [756]=0.00383524108505368"
[1] "wsvd [757]=0.00383523172384271"
[1] "wsvd [758]=0.00383522244757356"
[1] "wsvd [759]=0.00383521325544326"
[1] "wsvd [760]=0.00383520414665666"
[1] "wsvd [761]=0.0038351951204264"
[1] "wsvd [762]=0.00383518617597282"
[1] "wsvd [763]=0.00383517731252387"
[1] "wsvd [764]=0.00383516852931505"
[1] "wsvd [765]=0.00383515982558931"
[1] "wsvd [766]=0.00383515120059699"
[1] "wsvd [767]=0.00383514265359577"
[1] "wsvd [768]=0.00383513418385054"
[1] "wsvd [769]=0.00383512579063336"
[1] "wsvd [770]=0.0038351174732234"
[1] "wsvd [771]=0.00383510923090683"
[1] "wsvd [772]=0.0038351010629768"
[1] "wsvd [773]=0.00383509296873331"
[1] "wsvd [774]=0.00383508494748319"
[1] "wsvd [775]=0.00383507699854"
[1] "wsvd [776]=0.003835069121224"
[1] "wsvd [777]=0.00383506131486201"
[1] "wsvd [778]=0.00383505357878744"
[1] "wsvd [779]=0.00383504591234015"
[1] "wsvd [780]=0.0038350383148664"
[1] "wsvd [781]=0.00383503078571882"
[1] "wsvd [782]=0.0038350233242563"
[1] "wsvd [783]=0.00383501592984395"
[1] "wsvd [784]=0.00383500860185305"
[1] "wsvd [785]=0.00383500133966096"
[1] "wsvd [786]=0.00383499414265108"
[1] "wsvd [787]=0.00383498701021277"
[1] "wsvd [788]=0.0038349799417413"
[1] "wsvd [789]=0.00383497293663781"
[1] "wsvd [790]=0.00383496599430923"
[1] "wsvd [791]=0.0038349591141682"
[1] "wsvd [792]=0.00383495229563306"
[1] "wsvd [793]=0.00383494553812778"
[1] "wsvd [794]=0.00383493884108187"
[1] "wsvd [795]=0.00383493220393037"
[1] "wsvd [796]=0.00383492562611375"
[1] "wsvd [797]=0.0038349191070779"
[1] "wsvd [798]=0.00383491264627406"
[1] "wsvd [799]=0.00383490624315874"
[1] "wsvd [800]=0.0038348998971937"
[1] "wsvd [801]=0.00383489360784589"
[1] "wsvd [802]=0.0038348873745874"
[1] "wsvd [803]=0.00383488119689538"
[1] "wsvd [804]=0.00383487507425203"
[1] "wsvd [805]=0.00383486900614455"
[1] "wsvd [806]=0.00383486299206503"
[1] "wsvd [807]=0.00383485703151049"
[1] "wsvd [808]=0.00383485112398276"
[1] "wsvd [809]=0.00383484526898847"
[1] "wsvd [810]=0.00383483946603896"
[1] "wsvd [811]=0.0038348337146503"
[1] "wsvd [812]=0.0038348280143432"
[1] "wsvd [813]=0.00383482236464295"
[1] "wsvd [814]=0.0038348167650794"
[1] "wsvd [815]=0.00383481121518693"
[1] "wsvd [816]=0.00383480571450436"
[1] "wsvd [817]=0.00383480026257495"
[1] "wsvd [818]=0.00383479485894632"
[1] "wsvd [819]=0.00383478950317042"
[1] "wsvd [820]=0.00383478419480351"
[1] "wsvd [821]=0.00383477893340609"
[1] "wsvd [822]=0.00383477371854285"
[1] "wsvd [823]=0.00383476854978267"
[1] "wsvd [824]=0.00383476342669853"
[1] "wsvd [825]=0.0038347583488675"
[1] "wsvd [826]=0.0038347533158707"
[1] "wsvd [827]=0.00383474832729324"
[1] "wsvd [828]=0.0038347433827242"
[1] "wsvd [829]=0.00383473848175658"
[1] "wsvd [830]=0.00383473362398727"
[1] "wsvd [831]=0.00383472880901698"
[1] "wsvd [832]=0.00383472403645028"
[1] "wsvd [833]=0.00383471930589545"
[1] "wsvd [834]=0.00383471461696456"
[1] "wsvd [835]=0.00383470996927333"
[1] "wsvd [836]=0.00383470536244117"
[1] "wsvd [837]=0.00383470079609111"
[1] "wsvd [838]=0.00383469626984975"
[1] "wsvd [839]=0.00383469178334728"
[1] "wsvd [840]=0.00383468733621736"
[1] "wsvd [841]=0.00383468292809718"
[1] "wsvd [842]=0.00383467855862736"
[1] "wsvd [843]=0.00383467422745193"
[1] "wsvd [844]=0.0038346699342183"
[1] "wsvd [845]=0.00383466567857725"
[1] "wsvd [846]=0.00383466146018287"
[1] "wsvd [847]=0.00383465727869251"
[1] "wsvd [848]=0.00383465313376681"
[1] "wsvd [849]=0.0038346490250696"
[1] "wsvd [850]=0.00383464495226791"
[1] "wsvd [851]=0.00383464091503192"
[1] "wsvd [852]=0.00383463691303495"
[1] "wsvd [853]=0.00383463294595341"
[1] "wsvd [854]=0.00383462901346676"
[1] "wsvd [855]=0.00383462511525753"
[1] "wsvd [856]=0.00383462125101121"
[1] "wsvd [857]=0.00383461742041631"
[1] "wsvd [858]=0.00383461362316426"
[1] "wsvd [859]=0.00383460985894942"
[1] "wsvd [860]=0.00383460612746904"
[1] "wsvd [861]=0.00383460242842323"
[1] "wsvd [862]=0.00383459876151492"
[1] "wsvd [863]=0.00383459512644987"
[1] "wsvd [864]=0.00383459152293662"
[1] "wsvd [865]=0.00383458795068643"
[1] "wsvd [866]=0.00383458440941333"
[1] "wsvd [867]=0.003834580898834"
[1] "wsvd [868]=0.00383457741866784"
[1] "wsvd [869]=0.00383457396863686"
[1] "wsvd [870]=0.00383457054846572"
[1] "wsvd [871]=0.00383456715788165"
[1] "wsvd [872]=0.00383456379661446"
[1] "wsvd [873]=0.00383456046439651"
[1] "wsvd [874]=0.00383455716096268"
[1] "wsvd [875]=0.00383455388605033"
[1] "wsvd [876]=0.00383455063939932"
[1] "wsvd [877]=0.00383454742075194"
[1] "wsvd [878]=0.00383454422985289"
[1] "wsvd [879]=0.00383454106644931"
[1] "wsvd [880]=0.00383453793029068"
[1] "wsvd [881]=0.00383453482112885"
[1] "wsvd [882]=0.003834531738718"
[1] "wsvd [883]=0.00383452868281462"
[1] "wsvd [884]=0.00383452565317749"
[1] "wsvd [885]=0.00383452264956765"
[1] "wsvd [886]=0.00383451967174838"
[1] "wsvd [887]=0.00383451671948518"
[1] "wsvd [888]=0.00383451379254576"
[1] "wsvd [889]=0.0038345108907"
[1] "wsvd [890]=0.00383450801371994"
[1] "wsvd [891]=0.00383450516137975"
[1] "wsvd [892]=0.00383450233345572"
[1] "wsvd [893]=0.00383449952972625"
[1] "wsvd [894]=0.00383449674997179"
[1] "wsvd [895]=0.00383449399397485"
[1] "wsvd [896]=0.00383449126152"
[1] "wsvd [897]=0.00383448855239379"
[1] "wsvd [898]=0.00383448586638479"
[1] "wsvd [899]=0.00383448320328355"
[1] "wsvd [900]=0.00383448056288256"
[1] "wsvd [901]=0.00383447794497626"
[1] "wsvd [902]=0.00383447534936102"
[1] "wsvd [903]=0.00383447277583508"
[1] "wsvd [904]=0.00383447022419859"
[1] "wsvd [905]=0.00383446769425358"
[1] "wsvd [906]=0.00383446518580389"
[1] "wsvd [907]=0.00383446269865521"
[1] "wsvd [908]=0.00383446023261504"
[1] "wsvd [909]=0.00383445778749268"
[1] "wsvd [910]=0.00383445536309921"
[1] "wsvd [911]=0.00383445295924746"
[1] "wsvd [912]=0.003834450575752"
[1] "wsvd [913]=0.00383444821242915"
[1] "wsvd [914]=0.00383444586909691"
[1] "wsvd [915]=0.00383444354557501"
[1] "wsvd [916]=0.00383444124168482"
[1] "wsvd [917]=0.00383443895724941"
[1] "wsvd [918]=0.00383443669209346"
[1] "wsvd [919]=0.0038344344460433"
[1] "wsvd [920]=0.00383443221892688"
[1] "wsvd [921]=0.00383443001057374"
[1] "wsvd [922]=0.003834427820815"
[1] "wsvd [923]=0.00383442564948335"
[1] "wsvd [924]=0.00383442349641305"
[1] "wsvd [925]=0.00383442136143988"
[1] "wsvd [926]=0.00383441924440115"
[1] "wsvd [927]=0.00383441714513568"
[1] "wsvd [928]=0.00383441506348379"
[1] "wsvd [929]=0.00383441299928728"
[1] "wsvd [930]=0.00383441095238941"
[1] "wsvd [931]=0.0038344089226349"
[1] "wsvd [932]=0.0038344069098699"
[1] "wsvd [933]=0.00383440491394201"
[1] "wsvd [934]=0.0038344029347002"
[1] "wsvd [935]=0.00383440097199488"
[1] "wsvd [936]=0.00383439902567782"
[1] "wsvd [937]=0.00383439709560217"
[1] "wsvd [938]=0.00383439518162244"
[1] "wsvd [939]=0.00383439328359449"
[1] "wsvd [940]=0.00383439140137549"
[1] "wsvd [941]=0.00383438953482397"
[1] "wsvd [942]=0.00383438768379973"
[1] "wsvd [943]=0.00383438584816388"
[1] "wsvd [944]=0.00383438402777883"
[1] "wsvd [945]=0.00383438222250823"
[1] "wsvd [946]=0.00383438043221702"
[1] "wsvd [947]=0.00383437865677135"
[1] "wsvd [948]=0.00383437689603865"
[1] "wsvd [949]=0.00383437514988753"
[1] "wsvd [950]=0.00383437341818785"
[1] "wsvd [951]=0.00383437170081063"
[1] "wsvd [952]=0.00383436999762813"
[1] "wsvd [953]=0.00383436830851373"
[1] "wsvd [954]=0.00383436663334203"
[1] "wsvd [955]=0.00383436497198874"
[1] "wsvd [956]=0.00383436332433075"
[1] "wsvd [957]=0.00383436169024607"
[1] "wsvd [958]=0.00383436006961383"
[1] "wsvd [959]=0.00383435846231427"
[1] "wsvd [960]=0.00383435686822873"
[1] "wsvd [961]=0.00383435528723967"
[1] "wsvd [962]=0.0038343537192306"
[1] "wsvd [963]=0.0038343521640861"
[1] "wsvd [964]=0.00383435062169183"
[1] "wsvd [965]=0.0038343490919345"
[1] "wsvd [966]=0.00383434757470184"
[1] "wsvd [967]=0.00383434606988263"
[1] "wsvd [968]=0.00383434457736666"
[1] "wsvd [969]=0.00383434309704475"
[1] "wsvd [970]=0.00383434162880869"
[1] "wsvd [971]=0.0038343401725513"
[1] "wsvd [972]=0.00383433872816636"
[1] "wsvd [973]=0.00383433729554862"
[1] "wsvd [974]=0.00383433587459382"
[1] "wsvd [975]=0.00383433446519863"
[1] "wsvd [976]=0.00383433306726069"
[1] "wsvd [977]=0.00383433168067854"
[1] "wsvd [978]=0.0038343303053517"
[1] "wsvd [979]=0.00383432894118057"
[1] "wsvd [980]=0.00383432758806648"
[1] "wsvd [981]=0.00383432624591166"
[1] "wsvd [982]=0.00383432491461924"
[1] "wsvd [983]=0.00383432359409322"
[1] "wsvd [984]=0.0038343222842385"
[1] "wsvd [985]=0.00383432098496082"
[1] "wsvd [986]=0.00383431969616682"
[1] "wsvd [987]=0.00383431841776396"
[1] "wsvd [988]=0.00383431714966057"
[1] "wsvd [989]=0.00383431589176579"
[1] "wsvd [990]=0.00383431464398962"
[1] "wsvd [991]=0.00383431340624286"
[1] "wsvd [992]=0.00383431217843713"
[1] "wsvd [993]=0.00383431096048486"
[1] "wsvd [994]=0.00383430975229929"
[1] "wsvd [995]=0.00383430855379441"
[1] "wsvd [996]=0.00383430736488505"
[1] "wsvd [997]=0.00383430618548678"
[1] "wsvd [998]=0.00383430501551595"
[1] "wsvd [999]=0.00383430385488966"
[1] "wsvd [1000]=0.0038343027035258"
[1] "wsvd [1001]=0.00383430156134296"
[1] "wsvd [1002]=0.0038343004282605"
[1] "wsvd [1003]=0.00383429930419853"
[1] "wsvd [1004]=0.00383429818907784"
[1] "wsvd [1005]=0.00383429708281998"
[1] "wsvd [1006]=0.00383429598534721"
[1] "wsvd [1007]=0.00383429489658247"
[1] "wsvd [1008]=0.00383429381644943"
[1] "wsvd [1009]=0.00383429274487244"
[1] "wsvd [1010]=0.00383429168177654"
[1] "wsvd [1011]=0.00383429062708745"
[1] "wsvd [1012]=0.00383428958073157"
[1] "wsvd [1013]=0.00383428854263597"
[1] "wsvd [1014]=0.00383428751272836"
[1] "wsvd [1015]=0.00383428649093713"
[1] "wsvd [1016]=0.00383428547719132"
[1] "wsvd [1017]=0.0038342844714206"
[1] "wsvd [1018]=0.00383428347355528"
x.plra1.naive2 = pois_lra1(x,1,max_iters = 10000,atol=1e-9,verbose=TRUE) 
[1] "wsvd [0]=0.0357370934894206"
[1] "wsvd [1]=0.0216973393395743"
[1] "wsvd [2]=0.0183559836838193"
[1] "wsvd [3]=0.0166671902679346"
[1] "wsvd [4]=0.0156058206162763"
[1] "wsvd [5]=0.0148556454379835"
[1] "wsvd [6]=0.0142861298681759"
[1] "wsvd [7]=0.013833147538475"
[1] "wsvd [8]=0.0134611235799792"
[1] "wsvd [9]=0.0131484370771326"
[1] "wsvd [10]=0.012880983456594"
[1] "wsvd [11]=0.0126490337851493"
[1] "wsvd [12]=0.0124455735801941"
[1] "wsvd [13]=0.0122653630701558"
[1] "wsvd [14]=0.0121043743099799"
[1] "wsvd [15]=0.0119594370517658"
[1] "wsvd [16]=0.0118280063032763"
[1] "wsvd [17]=0.011708004069002"
[1] "wsvd [18]=0.0115977081385263"
[1] "wsvd [19]=0.0114956717781811"
[1] "wsvd [20]=0.011400664365524"
[1] "wsvd [21]=0.0113116266182692"
[1] "wsvd [22]=0.011227636252621"
[1] "wsvd [23]=0.0111478812675515"
[1] "wsvd [24]=0.0110716389255055"
[1] "wsvd [25]=0.0109982590761071"
[1] "wsvd [26]=0.0109271508586666"
[1] "wsvd [27]=0.0108577720882793"
[1] "wsvd [28]=0.010789620819991"
[1] "wsvd [29]=0.0107222287215383"
[1] "wsvd [30]=0.0106551559838807"
[1] "wsvd [31]=0.0105879875706331"
[1] "wsvd [32]=0.0105203306593684"
[1] "wsvd [33]=0.0104518131639629"
[1] "wsvd [34]=0.0103820832504919"
[1] "wsvd [35]=0.0103108097714377"
[1] "wsvd [36]=0.0102376835453458"
[1] "wsvd [37]=0.0101624194024742"
[1] "wsvd [38]=0.0100847589023808"
[1] "wsvd [39]=0.0100044736079557"
[1] "wsvd [40]=0.00992136877377786"
[1] "wsvd [41]=0.00983528727707797"
[1] "wsvd [42]=0.00974611358989482"
[1] "wsvd [43]=0.00965377756471187"
[1] "wsvd [44]=0.00955825778683369"
[1] "wsvd [45]=0.00945958423899753"
[1] "wsvd [46]=0.00935784003081019"
[1] "wsvd [47]=0.0092531619702242"
[1] "wsvd [48]=0.00914573979760668"
[1] "wsvd [49]=0.00903581396424226"
[1] "wsvd [50]=0.00892367191338431"
[1] "wsvd [51]=0.00880964290810675"
[1] "wsvd [52]=0.00869409153934431"
[1] "wsvd [53]=0.00857741013175579"
[1] "wsvd [54]=0.00846001033647807"
[1] "wsvd [55]=0.00834231425155498"
[1] "wsvd [56]=0.00822474543790227"
[1] "wsvd [57]=0.00810772019884672"
[1] "wsvd [58]=0.00799163946522883"
[1] "wsvd [59]=0.00787688157921606"
[1] "wsvd [60]=0.00776379620395342"
[1] "wsvd [61]=0.00765269950992438"
[1] "wsvd [62]=0.00754387070971126"
[1] "wsvd [63]=0.00743754993749412"
[1] "wsvd [64]=0.00733393740361531"
[1] "wsvd [65]=0.00723319370168147"
[1] "wsvd [66]=0.00713544110796183"
[1] "wsvd [67]=0.00704076569054782"
[1] "wsvd [68]=0.00694922003773357"
[1] "wsvd [69]=0.00686082641923442"
[1] "wsvd [70]=0.0067755802074884"
[1] "wsvd [71]=0.00669345340651867"
[1] "wsvd [72]=0.00661439815994234"
[1] "wsvd [73]=0.0065383501353162"
[1] "wsvd [74]=0.00646523170720401"
[1] "wsvd [75]=0.00639495488473385"
[1] "wsvd [76]=0.00632742395007832"
[1] "wsvd [77]=0.00626253779174328"
[1] "wsvd [78]=0.00620019193064968"
[1] "wsvd [79]=0.0061402802478447"
[1] "wsvd [80]=0.00608269643056267"
[1] "wsvd [81]=0.0060273351586537"
[1] "wsvd [82]=0.00597409305653323"
[1] "wsvd [83]=0.0059228694372047"
[1] "wsvd [84]=0.00587356686497122"
[1] "wsvd [85]=0.00582609156253562"
[1] "wsvd [86]=0.00578035368659729"
[1] "wsvd [87]=0.00573626749403727"
[1] "wsvd [88]=0.00569375141853934"
[1] "wsvd [89]=0.00565272807517424"
[1] "wsvd [90]=0.005613124208189"
[1] "wsvd [91]=0.00557487059506764"
[1] "wsvd [92]=0.00553790191791545"
[1] "wsvd [93]=0.0055021566113923"
[1] "wsvd [94]=0.00546757669479487"
[1] "wsvd [95]=0.00543410759446247"
[1] "wsvd [96]=0.0054016979614494"
[1] "wsvd [97]=0.00537029948835545"
[1] "wsvd [98]=0.00533986672831934"
[1] "wsvd [99]=0.00531035691844026"
[1] "wsvd [100]=0.00528172980928286"
[1] "wsvd [101]=0.00525394750162384"
[1] "wsvd [102]=0.00522697429119773"
[1] "wsvd [103]=0.00520077652188098"
[1] "wsvd [104]=0.00517532244750396"
[1] "wsvd [105]=0.00515058210228814"
[1] "wsvd [106]=0.00512652717976084"
[1] "wsvd [107]=0.00510313091989329"
[1] "wsvd [108]=0.00508036800413224"
[1] "wsvd [109]=0.00505821445794513"
[1] "wsvd [110]=0.00503664756046735"
[1] "wsvd [111]=0.00501564576082491"
[1] "wsvd [112]=0.00499518860070152"
[1] "wsvd [113]=0.00497525664272427"
[1] "wsvd [114]=0.00495583140425306"
[1] "wsvd [115]=0.0049368952961752"
[1] "wsvd [116]=0.00491843156632551"
[1] "wsvd [117]=0.00490042424717311"
[1] "wsvd [118]=0.00488285810743825"
[1] "wsvd [119]=0.00486571860732497"
[1] "wsvd [120]=0.00484899185707717"
[1] "wsvd [121]=0.00483266457858789"
[1] "wsvd [122]=0.00481672406981185"
[1] "wsvd [123]=0.00480115817175139"
[1] "wsvd [124]=0.00478595523780479"
[1] "wsvd [125]=0.00477110410528297"
[1] "wsvd [126]=0.0047565940689176"
[1] "wsvd [127]=0.00474241485619808"
[1] "wsvd [128]=0.00472855660438936"
[1] "wsvd [129]=0.00471500983909516"
[1] "wsvd [130]=0.00470176545424309"
[1] "wsvd [131]=0.00468881469337887"
[1] "wsvd [132]=0.00467614913216681"
[1] "wsvd [133]=0.00466376066200272"
[1] "wsvd [134]=0.0046516414746536"
[1] "wsvd [135]=0.00463978404784583"
[1] "wsvd [136]=0.00462818113173063"
[1] "wsvd [137]=0.00461682573616145"
[1] "wsvd [138]=0.00460571111872367"
[1] "wsvd [139]=0.00459483077346211"
[1] "wsvd [140]=0.00458417842025635"
[1] "wsvd [141]=0.00457374799479819"
[1] "wsvd [142]=0.00456353363912928"
[1] "wsvd [143]=0.00455352969270034"
[1] "wsvd [144]=0.00454373068391679"
[1] "wsvd [145]=0.00453413132213805"
[1] "wsvd [146]=0.0045247264901008"
[1] "wsvd [147]=0.00451551123673849"
[1] "wsvd [148]=0.00450648077037172"
[1] "wsvd [149]=0.00449763045224599"
[1] "wsvd [150]=0.00448895579039524"
[1] "wsvd [151]=0.00448045243381087"
[1] "wsvd [152]=0.00447211616689796"
[1] "wsvd [153]=0.00446394290420117"
[1] "wsvd [154]=0.00445592868538455"
[1] "wsvd [155]=0.00444806967045027"
[1] "wsvd [156]=0.00444036213518251"
[1] "wsvd [157]=0.00443280246680363"
[1] "wsvd [158]=0.00442538715983074"
[1] "wsvd [159]=0.00441811281212126"
[1] "wsvd [160]=0.00441097612109737"
[1] "wsvd [161]=0.00440397388013931"
[1] "wsvd [162]=0.0043971029751386"
[1] "wsvd [163]=0.00439036038120258"
[1] "wsvd [164]=0.00438374315950236"
[1] "wsvd [165]=0.00437724845425651"
[1] "wsvd [166]=0.00437087348984374"
[1] "wsvd [167]=0.00436461556803774"
[1] "wsvd [168]=0.00435847206535808"
[1] "wsvd [169]=0.00435244043053136"
[1] "wsvd [170]=0.00434651818205714"
[1] "wsvd [171]=0.00434070290587335"
[1] "wsvd [172]=0.00433499225311651"
[1] "wsvd [173]=0.00432938393797209"
[1] "wsvd [174]=0.00432387573561062"
[1] "wsvd [175]=0.00431846548020562"
[1] "wsvd [176]=0.00431315106302933"
[1] "wsvd [177]=0.00430793043062276"
[1] "wsvd [178]=0.00430280158303644"
[1] "wsvd [179]=0.00429776257213882"
[1] "wsvd [180]=0.00429281149998904"
[1] "wsvd [181]=0.00428794651727129"
[1] "wsvd [182]=0.00428316582178792"
[1] "wsvd [183]=0.00427846765700871"
[1] "wsvd [184]=0.00427385031067386"
[1] "wsvd [185]=0.00426931211344821"
[1] "wsvd [186]=0.00426485143762463"
[1] "wsvd [187]=0.0042604666958743"
[1] "wsvd [188]=0.00425615634004198"
[1] "wsvd [189]=0.00425191885998428"
[1] "wsvd [190]=0.00424775278244916"
[1] "wsvd [191]=0.00424365666999489"
[1] "wsvd [192]=0.00423962911994684"
[1] "wsvd [193]=0.00423566876339054"
[1] "wsvd [194]=0.00423177426419953"
[1] "wsvd [195]=0.00422794431809653"
[1] "wsvd [196]=0.00422417765174661"
[1] "wsvd [197]=0.00422047302188111"
[1] "wsvd [198]=0.00421682921445101"
[1] "wsvd [199]=0.00421324504380861"
[1] "wsvd [200]=0.00420971935191642"
[1] "wsvd [201]=0.00420625100758218"
[1] "wsvd [202]=0.00420283890571894"
[1] "wsvd [203]=0.00419948196662938"
[1] "wsvd [204]=0.00419617913531321"
[1] "wsvd [205]=0.004192929380797"
[1] "wsvd [206]=0.00418973169548547"
[1] "wsvd [207]=0.00418658509453339"
[1] "wsvd [208]=0.0041834886152374"
[1] "wsvd [209]=0.00418044131644707"
[1] "wsvd [210]=0.00417744227799424"
[1] "wsvd [211]=0.00417449060014029"
[1] "wsvd [212]=0.00417158540304041"
[1] "wsvd [213]=0.00416872582622442"
[1] "wsvd [214]=0.00416591102809353"
[1] "wsvd [215]=0.00416314018543232"
[1] "wsvd [216]=0.00416041249293562"
[1] "wsvd [217]=0.00415772716274957"
[1] "wsvd [218]=0.00415508342402653"
[1] "wsvd [219]=0.00415248052249314"
[1] "wsvd [220]=0.00414991772003131"
[1] "wsvd [221]=0.00414739429427156"
[1] "wsvd [222]=0.00414490953819821"
[1] "wsvd [223]=0.00414246275976627"
[1] "wsvd [224]=0.00414005328152934"
[1] "wsvd [225]=0.00413768044027837"
[1] "wsvd [226]=0.00413534358669071"
[1] "wsvd [227]=0.00413304208498941"
[1] "wsvd [228]=0.00413077531261203"
[1] "wsvd [229]=0.00412854265988899"
[1] "wsvd [230]=0.004126343529731"
[1] "wsvd [231]=0.00412417733732524"
[1] "wsvd [232]=0.00412204350984005"
[1] "wsvd [233]=0.00411994148613793"
[1] "wsvd [234]=0.00411787071649641"
[1] "wsvd [235]=0.00411583066233672"
[1] "wsvd [236]=0.00411382079595983"
[1] "wsvd [237]=0.00411184060028985"
[1] "wsvd [238]=0.00410988956862431"
[1] "wsvd [239]=0.00410796720439128"
[1] "wsvd [240]=0.0041060730209131"
[1] "wsvd [241]=0.00410420654117636"
[1] "wsvd [242]=0.00410236729760818"
[1] "wsvd [243]=0.00410055483185839"
[1] "wsvd [244]=0.00409876869458754"
[1] "wsvd [245]=0.0040970084452605"
[1] "wsvd [246]=0.00409527365194553"
[1] "wsvd [247]=0.00409356389111864"
[1] "wsvd [248]=0.00409187874747308"
[1] "wsvd [249]=0.00409021781373375"
[1] "wsvd [250]=0.0040885806904765"
[1] "wsvd [251]=0.00408696698595208"
[1] "wsvd [252]=0.00408537631591461"
[1] "wsvd [253]=0.00408380830345451"
[1] "wsvd [254]=0.00408226257883568"
[1] "wsvd [255]=0.00408073877933682"
[1] "wsvd [256]=0.00407923654909679"
[1] "wsvd [257]=0.00407775553896392"
[1] "wsvd [258]=0.00407629540634908"
[1] "wsvd [259]=0.00407485581508248"
[1] "wsvd [260]=0.00407343643527402"
[1] "wsvd [261]=0.00407203694317716"
[1] "wsvd [262]=0.00407065702105619"
[1] "wsvd [263]=0.00406929635705671"
[1] "wsvd [264]=0.00406795464507942"
[1] "wsvd [265]=0.00406663158465696"
[1] "wsvd [266]=0.00406532688083376"
[1] "wsvd [267]=0.00406404024404885"
[1] "wsvd [268]=0.00406277139002157"
[1] "wsvd [269]=0.00406152003963994"
[1] "wsvd [270]=0.00406028591885187"
[1] "wsvd [271]=0.00405906875855892"
[1] "wsvd [272]=0.00405786829451261"
[1] "wsvd [273]=0.00405668426721323"
[1] "wsvd [274]=0.00405551642181113"
[1] "wsvd [275]=0.00405436450801025"
[1] "wsvd [276]=0.00405322827997402"
[1] "wsvd [277]=0.00405210749623348"
[1] "wsvd [278]=0.00405100191959753"
[1] "wsvd [279]=0.00404991131706534"
[1] "wsvd [280]=0.00404883545974079"
[1] "wsvd [281]=0.00404777412274889"
[1] "wsvd [282]=0.0040467270851542"
[1] "wsvd [283]=0.00404569412988111"
[1] "wsvd [284]=0.00404467504363594"
[1] "wsvd [285]=0.00404366961683094"
[1] "wsvd [286]=0.0040426776435099"
[1] "wsvd [287]=0.00404169892127559"
[1] "wsvd [288]=0.00404073325121878"
[1] "wsvd [289]=0.00403978043784893"
[1] "wsvd [290]=0.00403884028902645"
[1] "wsvd [291]=0.00403791261589643"
[1] "wsvd [292]=0.00403699723282399"
[1] "wsvd [293]=0.00403609395733096"
[1] "wsvd [294]=0.00403520261003409"
[1] "wsvd [295]=0.00403432301458457"
[1] "wsvd [296]=0.00403345499760895"
[1] "wsvd [297]=0.00403259838865132"
[1] "wsvd [298]=0.00403175302011687"
[1] "wsvd [299]=0.00403091872721659"
[1] "wsvd [300]=0.00403009534791327"
[1] "wsvd [301]=0.00402928272286864"
[1] "wsvd [302]=0.00402848069539174"
[1] "wsvd [303]=0.00402768911138832"
[1] "wsvd [304]=0.00402690781931143"
[1] "wsvd [305]=0.00402613667011305"
[1] "wsvd [306]=0.00402537551719679"
[1] "wsvd [307]=0.00402462421637155"
[1] "wsvd [308]=0.00402388262580631"
[1] "wsvd [309]=0.00402315060598575"
[1] "wsvd [310]=0.0040224280196669"
[1] "wsvd [311]=0.00402171473183673"
[1] "wsvd [312]=0.00402101060967061"
[1] "wsvd [313]=0.00402031552249164"
[1] "wsvd [314]=0.0040196293417309"
[1] "wsvd [315]=0.0040189519408885"
[1] "wsvd [316]=0.00401828319549545"
[1] "wsvd [317]=0.00401762298307638"
[1] "wsvd [318]=0.00401697118311297"
[1] "wsvd [319]=0.0040163276770082"
[1] "wsvd [320]=0.00401569234805135"
[1] "wsvd [321]=0.00401506508138367"
[1] "wsvd [322]=0.00401444576396482"
[1] "wsvd [323]=0.00401383428454003"
[1] "wsvd [324]=0.00401323053360781"
[1] "wsvd [325]=0.00401263440338851"
[1] "wsvd [326]=0.00401204578779333"
[1] "wsvd [327]=0.00401146458239415"
[1] "wsvd [328]=0.00401089068439383"
[1] "wsvd [329]=0.00401032399259719"
[1] "wsvd [330]=0.00400976440738259"
[1] "wsvd [331]=0.004009211830674"
[1] "wsvd [332]=0.00400866616591379"
[1] "wsvd [333]=0.00400812731803585"
[1] "wsvd [334]=0.0040075951934395"
[1] "wsvd [335]=0.00400706969996371"
[1] "wsvd [336]=0.00400655074686194"
[1] "wsvd [337]=0.0040060382447775"
[1] "wsvd [338]=0.00400553210571933"
[1] "wsvd [339]=0.0040050322430383"
[1] "wsvd [340]=0.004004538571404"
[1] "wsvd [341]=0.00400405100678192"
[1] "wsvd [342]=0.00400356946641117"
[1] "wsvd [343]=0.00400309386878256"
[1] "wsvd [344]=0.00400262413361713"
[1] "wsvd [345]=0.00400216018184514"
[1] "wsvd [346]=0.00400170193558542"
[1] "wsvd [347]=0.00400124931812515"
[1] "wsvd [348]=0.00400080225389999"
[1] "wsvd [349]=0.00400036066847468"
[1] "wsvd [350]=0.00399992448852389"
[1] "wsvd [351]=0.00399949364181358"
[1] "wsvd [352]=0.0039990680571826"
[1] "wsvd [353]=0.00399864766452469"
[1] "wsvd [354]=0.00399823239477083"
[1] "wsvd [355]=0.00399782217987195"
[1] "wsvd [356]=0.00399741695278185"
[1] "wsvd [357]=0.00399701664744063"
[1] "wsvd [358]=0.00399662119875826"
[1] "wsvd [359]=0.00399623054259856"
[1] "wsvd [360]=0.00399584461576346"
[1] "wsvd [361]=0.00399546335597752"
[1] "wsvd [362]=0.0039950867018728"
[1] "wsvd [363]=0.00399471459297397"
[1] "wsvd [364]=0.0039943469696837"
[1] "wsvd [365]=0.00399398377326832"
[1] "wsvd [366]=0.00399362494584378"
[1] "wsvd [367]=0.00399327043036183"
[1] "wsvd [368]=0.00399292017059647"
[1] "wsvd [369]=0.00399257411113063"
[1] "wsvd [370]=0.00399223219734315"
[1] "wsvd [371]=0.00399189437539592"
[1] "wsvd [372]=0.00399156059222133"
[1] "wsvd [373]=0.00399123079550988"
[1] "wsvd [374]=0.00399090493369808"
[1] "wsvd [375]=0.0039905829559565"
[1] "wsvd [376]=0.00399026481217811"
[1] "wsvd [377]=0.00398995045296678"
[1] "wsvd [378]=0.00398963982962598"
[1] "wsvd [379]=0.00398933289414775"
[1] "wsvd [380]=0.00398902959920176"
[1] "wsvd [381]=0.00398872989812469"
[1] "wsvd [382]=0.00398843374490968"
[1] "wsvd [383]=0.00398814109419606"
[1] "wsvd [384]=0.0039878519012592"
[1] "wsvd [385]=0.00398756612200058"
[1] "wsvd [386]=0.00398728371293803"
[1] "wsvd [387]=0.00398700463119609"
[1] "wsvd [388]=0.00398672883449663"
[1] "wsvd [389]=0.00398645628114956"
[1] "wsvd [390]=0.0039861869300437"
[1] "wsvd [391]=0.00398592074063792"
[1] "wsvd [392]=0.00398565767295225"
[1] "wsvd [393]=0.00398539768755933"
[1] "wsvd [394]=0.00398514074557586"
[1] "wsvd [395]=0.00398488680865431"
[1] "wsvd [396]=0.00398463583897469"
[1] "wsvd [397]=0.00398438779923652"
[1] "wsvd [398]=0.0039841426526509"
[1] "wsvd [399]=0.00398390036293273"
[1] "wsvd [400]=0.00398366089429306"
[1] "wsvd [401]=0.00398342421143159"
[1] "wsvd [402]=0.00398319027952927"
[1] "wsvd [403]=0.00398295906424104"
[1] "wsvd [404]=0.00398273053168867"
[1] "wsvd [405]=0.00398250464845376"
[1] "wsvd [406]=0.00398228138157085"
[1] "wsvd [407]=0.00398206069852061"
[1] "wsvd [408]=0.00398184256722316"
[1] "wsvd [409]=0.00398162695603156"
[1] "wsvd [410]=0.00398141383372531"
[1] "wsvd [411]=0.00398120316950401"
[1] "wsvd [412]=0.00398099493298118"
[1] "wsvd [413]=0.00398078909417805"
[1] "wsvd [414]=0.00398058562351759"
[1] "wsvd [415]=0.00398038449181852"
[1] "wsvd [416]=0.00398018567028953"
[1] "wsvd [417]=0.00397998913052354"
[1] "wsvd [418]=0.00397979484449199"
[1] "wsvd [419]=0.00397960278453939"
[1] "wsvd [420]=0.00397941292337775"
[1] "wsvd [421]=0.00397922523408132"
[1] "wsvd [422]=0.00397903969008121"
[1] "wsvd [423]=0.00397885626516026"
[1] "wsvd [424]=0.00397867493344789"
[1] "wsvd [425]=0.0039784956694151"
[1] "wsvd [426]=0.00397831844786947"
[1] "wsvd [427]=0.00397814324395035"
[1] "wsvd [428]=0.00397797003312404"
[1] "wsvd [429]=0.00397779879117907"
[1] "wsvd [430]=0.00397762949422157"
[1] "wsvd [431]=0.00397746211867073"
[1] "wsvd [432]=0.00397729664125426"
[1] "wsvd [433]=0.00397713303900402"
[1] "wsvd [434]=0.00397697128925162"
[1] "wsvd [435]=0.0039768113696242"
[1] "wsvd [436]=0.00397665325804016"
[1] "wsvd [437]=0.00397649693270505"
[1] "wsvd [438]=0.0039763423721075"
[1] "wsvd [439]=0.00397618955501514"
[1] "wsvd [440]=0.00397603846047074"
[1] "wsvd [441]=0.00397588906778824"
[1] "wsvd [442]=0.00397574135654896"
[1] "wsvd [443]=0.00397559530659783"
[1] "wsvd [444]=0.00397545089803966"
[1] "wsvd [445]=0.00397530811123549"
[1] "wsvd [446]=0.00397516692679903"
[1] "wsvd [447]=0.00397502732559306"
[1] "wsvd [448]=0.003974889288726"
[1] "wsvd [449]=0.00397475279754843"
[1] "wsvd [450]=0.00397461783364974"
[1] "wsvd [451]=0.00397448437885481"
[1] "wsvd [452]=0.00397435241522071"
[1] "wsvd [453]=0.00397422192503349"
[1] "wsvd [454]=0.003974092890805"
[1] "wsvd [455]=0.00397396529526976"
[1] "wsvd [456]=0.00397383912138189"
[1] "wsvd [457]=0.00397371435231205"
[1] "wsvd [458]=0.0039735909714445"
[1] "wsvd [459]=0.00397346896237411"
[1] "wsvd [460]=0.00397334830890351"
[1] "wsvd [461]=0.00397322899504018"
[1] "wsvd [462]=0.0039731110049937"
[1] "wsvd [463]=0.00397299432317294"
[1] "wsvd [464]=0.00397287893418336"
[1] "wsvd [465]=0.00397276482282432"
[1] "wsvd [466]=0.00397265197408639"
[1] "wsvd [467]=0.00397254037314885"
[1] "wsvd [468]=0.00397243000537699"
[1] "wsvd [469]=0.00397232085631971"
[1] "wsvd [470]=0.00397221291170694"
[1] "wsvd [471]=0.00397210615744722"
[1] "wsvd [472]=0.00397200057962529"
[1] "wsvd [473]=0.0039718961644997"
[1] "wsvd [474]=0.00397179289850046"
[1] "wsvd [475]=0.00397169076822673"
[1] "wsvd [476]=0.00397158976044457"
[1] "wsvd [477]=0.00397148986208466"
[1] "wsvd [478]=0.00397139106024009"
[1] "wsvd [479]=0.00397129334216422"
[1] "wsvd [480]=0.0039711966952685"
[1] "wsvd [481]=0.00397110110712037"
[1] "wsvd [482]=0.00397100656544117"
[1] "wsvd [483]=0.00397091305810408"
[1] "wsvd [484]=0.00397082057313213"
[1] "wsvd [485]=0.00397072909869616"
[1] "wsvd [486]=0.00397063862311288"
[1] "wsvd [487]=0.00397054913484295"
[1] "wsvd [488]=0.00397046062248902"
[1] "wsvd [489]=0.00397037307479392"
[1] "wsvd [490]=0.00397028648063873"
[1] "wsvd [491]=0.00397020082904101"
[1] "wsvd [492]=0.00397011610915299"
[1] "wsvd [493]=0.00397003231025977"
[1] "wsvd [494]=0.00396994942177758"
[1] "wsvd [495]=0.00396986743325209"
[1] "wsvd [496]=0.00396978633435664"
[1] "wsvd [497]=0.00396970611489065"
[1] "wsvd [498]=0.0039696267647779"
[1] "wsvd [499]=0.00396954827406492"
[1] "wsvd [500]=0.00396947063291943"
[1] "wsvd [501]=0.00396939383162867"
[1] "wsvd [502]=0.00396931786059793"
[1] "wsvd [503]=0.00396924271034896"
[1] "wsvd [504]=0.00396916837151846"
[1] "wsvd [505]=0.0039690948348566"
[1] "wsvd [506]=0.00396902209122554"
[1] "wsvd [507]=0.00396895013159798"
[1] "wsvd [508]=0.00396887894705571"
[1] "wsvd [509]=0.00396880852878823"
[1] "wsvd [510]=0.00396873886809134"
[1] "wsvd [511]=0.00396866995636575"
[1] "wsvd [512]=0.00396860178511576"
[1] "wsvd [513]=0.0039685343459479"
[1] "wsvd [514]=0.00396846763056964"
[1] "wsvd [515]=0.00396840163078804"
[1] "wsvd [516]=0.00396833633850855"
[1] "wsvd [517]=0.00396827174573368"
[1] "wsvd [518]=0.00396820784456176"
[1] "wsvd [519]=0.00396814462718578"
[1] "wsvd [520]=0.00396808208589209"
[1] "wsvd [521]=0.00396802021305927"
[1] "wsvd [522]=0.00396795900115691"
[1] "wsvd [523]=0.00396789844274451"
[1] "wsvd [524]=0.00396783853047024"
[1] "wsvd [525]=0.00396777925706992"
[1] "wsvd [526]=0.00396772061536581"
[1] "wsvd [527]=0.00396766259826559"
[1] "wsvd [528]=0.00396760519876121"
[1] "wsvd [529]=0.00396754840992787"
[1] "wsvd [530]=0.00396749222492294"
[1] "wsvd [531]=0.00396743663698493"
[1] "wsvd [532]=0.00396738163943247"
[1] "wsvd [533]=0.00396732722566329"
[1] "wsvd [534]=0.00396727338915322"
[1] "wsvd [535]=0.0039672201234552"
[1] "wsvd [536]=0.00396716742219834"
[1] "wsvd [537]=0.00396711527908692"
[1] "wsvd [538]=0.00396706368789949"
[1] "wsvd [539]=0.00396701264248789"
[1] "wsvd [540]=0.00396696213677637"
[1] "wsvd [541]=0.00396691216476067"
[1] "wsvd [542]=0.0039668627205071"
[1] "wsvd [543]=0.00396681379815172"
[1] "wsvd [544]=0.00396676539189938"
[1] "wsvd [545]=0.00396671749602296"
[1] "wsvd [546]=0.00396667010486243"
[1] "wsvd [547]=0.00396662321282409"
[1] "wsvd [548]=0.00396657681437968"
[1] "wsvd [549]=0.00396653090406562"
[1] "wsvd [550]=0.00396648547648216"
[1] "wsvd [551]=0.00396644052629265"
[1] "wsvd [552]=0.00396639604822266"
[1] "wsvd [553]=0.00396635203705932"
[1] "wsvd [554]=0.00396630848765047"
[1] "wsvd [555]=0.00396626539490396"
[1] "wsvd [556]=0.00396622275378688"
[1] "wsvd [557]=0.00396618055932486"
[1] "wsvd [558]=0.0039661388066013"
[1] "wsvd [559]=0.00396609749075671"
[1] "wsvd [560]=0.00396605660698798"
[1] "wsvd [561]=0.00396601615054767"
[1] "wsvd [562]=0.00396597611674336"
[1] "wsvd [563]=0.00396593650093695"
[1] "wsvd [564]=0.00396589729854402"
[1] "wsvd [565]=0.00396585850503312"
[1] "wsvd [566]=0.0039658201159252"
[1] "wsvd [567]=0.00396578212679287"
[1] "wsvd [568]=0.00396574453325986"
[1] "wsvd [569]=0.00396570733100034"
[1] "wsvd [570]=0.00396567051573833"
[1] "wsvd [571]=0.00396563408324707"
[1] "wsvd [572]=0.00396559802934844"
[1] "wsvd [573]=0.00396556234991235"
[1] "wsvd [574]=0.00396552704085618"
[1] "wsvd [575]=0.00396549209814418"
[1] "wsvd [576]=0.00396545751778689"
[1] "wsvd [577]=0.00396542329584064"
[1] "wsvd [578]=0.00396538942840689"
[1] "wsvd [579]=0.0039653559116318"
[1] "wsvd [580]=0.0039653227417056"
[1] "wsvd [581]=0.0039652899148621"
[1] "wsvd [582]=0.00396525742737814"
[1] "wsvd [583]=0.00396522527557309"
[1] "wsvd [584]=0.00396519345580832"
[1] "wsvd [585]=0.00396516196448673"
[1] "wsvd [586]=0.00396513079805217"
[1] "wsvd [587]=0.00396509995298905"
[1] "wsvd [588]=0.00396506942582175"
[1] "wsvd [589]=0.00396503921311423"
[1] "wsvd [590]=0.0039650093114695"
[1] "wsvd [591]=0.00396497971752914"
[1] "wsvd [592]=0.00396495042797289"
[1] "wsvd [593]=0.00396492143951816"
[1] "wsvd [594]=0.00396489274891956"
[1] "wsvd [595]=0.0039648643529685"
[1] "wsvd [596]=0.00396483624849273"
[1] "wsvd [597]=0.00396480843235587"
[1] "wsvd [598]=0.00396478090145705"
[1] "wsvd [599]=0.00396475365273042"
[1] "wsvd [600]=0.00396472668314479"
[1] "wsvd [601]=0.00396469998970317"
[1] "wsvd [602]=0.00396467356944238"
[1] "wsvd [603]=0.00396464741943266"
[1] "wsvd [604]=0.00396462153677727"
[1] "wsvd [605]=0.00396459591861206"
[1] "wsvd [606]=0.00396457056210513"
[1] "wsvd [607]=0.00396454546445644"
[1] "wsvd [608]=0.00396452062289741"
[1] "wsvd [609]=0.00396449603469056"
[1] "wsvd [610]=0.00396447169712914"
[1] "wsvd [611]=0.00396444760753677"
[1] "wsvd [612]=0.00396442376326707"
[1] "wsvd [613]=0.00396440016170334"
[1] "wsvd [614]=0.00396437680025815"
[1] "wsvd [615]=0.00396435367637305"
[1] "wsvd [616]=0.00396433078751821"
[1] "wsvd [617]=0.00396430813119205"
[1] "wsvd [618]=0.00396428570492096"
[1] "wsvd [619]=0.00396426350625896"
[1] "wsvd [620]=0.00396424153278733"
[1] "wsvd [621]=0.00396421978211435"
[1] "wsvd [622]=0.00396419825187494"
[1] "wsvd [623]=0.00396417693973036"
[1] "wsvd [624]=0.00396415584336793"
[1] "wsvd [625]=0.00396413496050067"
[1] "wsvd [626]=0.00396411428886705"
[1] "wsvd [627]=0.00396409382623064"
[1] "wsvd [628]=0.00396407357037988"
[1] "wsvd [629]=0.00396405351912774"
[1] "wsvd [630]=0.00396403367031145"
[1] "wsvd [631]=0.00396401402179223"
[1] "wsvd [632]=0.00396399457145496"
[1] "wsvd [633]=0.00396397531720799"
[1] "wsvd [634]=0.00396395625698279"
[1] "wsvd [635]=0.0039639373887337"
[1] "wsvd [636]=0.0039639187104377"
[1] "wsvd [637]=0.0039639002200941"
[1] "wsvd [638]=0.00396388191572431"
[1] "wsvd [639]=0.00396386379537157"
[1] "wsvd [640]=0.00396384585710069"
[1] "wsvd [641]=0.00396382809899783"
[1] "wsvd [642]=0.00396381051917022"
[1] "wsvd [643]=0.00396379311574594"
[1] "wsvd [644]=0.00396377588687366"
[1] "wsvd [645]=0.0039637588307224"
[1] "wsvd [646]=0.00396374194548133"
[1] "wsvd [647]=0.00396372522935949"
[1] "wsvd [648]=0.00396370868058559"
[1] "wsvd [649]=0.00396369229740779"
[1] "wsvd [650]=0.00396367607809343"
[1] "wsvd [651]=0.00396366002092888"
[1] "wsvd [652]=0.00396364412421925"
[1] "wsvd [653]=0.00396362838628823"
[1] "wsvd [654]=0.00396361280547783"
[1] "wsvd [655]=0.0039635973801482"
[1] "wsvd [656]=0.00396358210867744"
[1] "wsvd [657]=0.00396356698946132"
[1] "wsvd [658]=0.00396355202091317"
[1] "wsvd [659]=0.0039635372014636"
[1] "wsvd [660]=0.00396352252956037"
[1] "wsvd [661]=0.00396350800366813"
[1] "wsvd [662]=0.00396349362226827"
[1] "wsvd [663]=0.00396347938385872"
[1] "wsvd [664]=0.00396346528695374"
[1] "wsvd [665]=0.00396345133008377"
[1] "wsvd [666]=0.00396343751179521"
[1] "wsvd [667]=0.00396342383065027"
[1] "wsvd [668]=0.00396341028522675"
[1] "wsvd [669]=0.0039633968741179"
[1] "wsvd [670]=0.00396338359593223"
[1] "wsvd [671]=0.00396337044929334"
[1] "wsvd [672]=0.00396335743283974"
[1] "wsvd [673]=0.00396334454522467"
[1] "wsvd [674]=0.00396333178511597"
[1] "wsvd [675]=0.00396331915119589"
[1] "wsvd [676]=0.00396330664216091"
[1] "wsvd [677]=0.00396329425672162"
[1] "wsvd [678]=0.00396328199360253"
[1] "wsvd [679]=0.00396326985154191"
[1] "wsvd [680]=0.00396325782929165"
[1] "wsvd [681]=0.00396324592561711"
[1] "wsvd [682]=0.00396323413929694"
[1] "wsvd [683]=0.00396322246912297"
[1] "wsvd [684]=0.00396321091390001"
[1] "wsvd [685]=0.00396319947244577"
[1] "wsvd [686]=0.00396318814359065"
[1] "wsvd [687]=0.00396317692617763"
[1] "wsvd [688]=0.00396316581906211"
[1] "wsvd [689]=0.00396315482111182"
[1] "wsvd [690]=0.0039631439312066"
[1] "wsvd [691]=0.00396313314823834"
[1] "wsvd [692]=0.00396312247111079"
[1] "wsvd [693]=0.00396311189873947"
[1] "wsvd [694]=0.0039631014300515"
[1] "wsvd [695]=0.00396309106398549"
[1] "wsvd [696]=0.00396308079949141"
[1] "wsvd [697]=0.00396307063553047"
[1] "wsvd [698]=0.00396306057107497"
[1] "wsvd [699]=0.00396305060510821"
[1] "wsvd [700]=0.00396304073662432"
[1] "wsvd [701]=0.00396303096462818"
[1] "wsvd [702]=0.00396302128813531"
[1] "wsvd [703]=0.00396301170617168"
[1] "wsvd [704]=0.00396300221777368"
[1] "wsvd [705]=0.00396299282198795"
[1] "wsvd [706]=0.00396298351787129"
[1] "wsvd [707]=0.00396297430449052"
[1] "wsvd [708]=0.00396296518092239"
[1] "wsvd [709]=0.00396295614625348"
[1] "wsvd [710]=0.00396294719958007"
[1] "wsvd [711]=0.00396293834000803"
[1] "wsvd [712]=0.00396292956665272"
[1] "wsvd [713]=0.00396292087863889"
[1] "wsvd [714]=0.00396291227510059"
[1] "wsvd [715]=0.00396290375518101"
[1] "wsvd [716]=0.00396289531803246"
[1] "wsvd [717]=0.00396288696281619"
[1] "wsvd [718]=0.00396287868870235"
[1] "wsvd [719]=0.00396287049486986"
[1] "wsvd [720]=0.00396286238050633"
[1] "wsvd [721]=0.00396285434480793"
[1] "wsvd [722]=0.00396284638697935"
[1] "wsvd [723]=0.00396283850623367"
[1] "wsvd [724]=0.00396283070179227"
[1] "wsvd [725]=0.00396282297288473"
[1] "wsvd [726]=0.00396281531874878"
[1] "wsvd [727]=0.00396280773863015"
[1] "wsvd [728]=0.00396280023178255"
[1] "wsvd [729]=0.00396279279746751"
[1] "wsvd [730]=0.00396278543495436"
[1] "wsvd [731]=0.0039627781435201"
[1] "wsvd [732]=0.00396277092244932"
[1] "wsvd [733]=0.00396276377103416"
[1] "wsvd [734]=0.00396275668857415"
[1] "wsvd [735]=0.00396274967437621"
[1] "wsvd [736]=0.00396274272775451"
[1] "wsvd [737]=0.00396273584803042"
[1] "wsvd [738]=0.00396272903453242"
[1] "wsvd [739]=0.00396272228659603"
[1] "wsvd [740]=0.00396271560356372"
[1] "wsvd [741]=0.00396270898478485"
[1] "wsvd [742]=0.00396270242961559"
[1] "wsvd [743]=0.00396269593741883"
[1] "wsvd [744]=0.00396268950756413"
[1] "wsvd [745]=0.00396268313942762"
[1] "wsvd [746]=0.00396267683239197"
[1] "wsvd [747]=0.00396267058584627"
[1] "wsvd [748]=0.003962664399186"
[1] "wsvd [749]=0.00396265827181292"
[1] "wsvd [750]=0.00396265220313506"
[1] "wsvd [751]=0.00396264619256658"
[1] "wsvd [752]=0.00396264023952777"
[1] "wsvd [753]=0.00396263434344492"
[1] "wsvd [754]=0.00396262850375034"
[1] "wsvd [755]=0.00396262271988218"
[1] "wsvd [756]=0.00396261699128448"
[1] "wsvd [757]=0.00396261131740703"
[1] "wsvd [758]=0.00396260569770534"
[1] "wsvd [759]=0.00396260013164056"
[1] "wsvd [760]=0.00396259461867946"
[1] "wsvd [761]=0.0039625891582943"
[1] "wsvd [762]=0.00396258374996284"
[1] "wsvd [763]=0.00396257839316824"
[1] "wsvd [764]=0.003962573087399"
[1] "wsvd [765]=0.00396256783214893"
[1] "wsvd [766]=0.00396256262691706"
[1] "wsvd [767]=0.00396255747120763"
[1] "wsvd [768]=0.00396255236452997"
[1] "wsvd [769]=0.0039625473063985"
[1] "wsvd [770]=0.00396254229633265"
[1] "wsvd [771]=0.0039625373338568"
[1] "wsvd [772]=0.00396253241850025"
[1] "wsvd [773]=0.00396252754979714"
[1] "wsvd [774]=0.00396252272728644"
[1] "wsvd [775]=0.00396251795051183"
[1] "wsvd [776]=0.00396251321902172"
[1] "wsvd [777]=0.00396250853236915"
[1] "wsvd [778]=0.00396250389011178"
[1] "wsvd [779]=0.00396249929181179"
[1] "wsvd [780]=0.00396249473703588"
[1] "wsvd [781]=0.0039624902253552"
[1] "wsvd [782]=0.0039624857563453"
[1] "wsvd [783]=0.0039624813295861"
[1] "wsvd [784]=0.0039624769446618"
[1] "wsvd [785]=0.00396247260116089"
[1] "wsvd [786]=0.00396246829867608"
[1] "wsvd [787]=0.00396246403680424"
[1] "wsvd [788]=0.00396245981514638"
[1] "wsvd [789]=0.00396245563330758"
[1] "wsvd [790]=0.00396245149089698"
[1] "wsvd [791]=0.0039624473875277"
[1] "wsvd [792]=0.00396244332281684"
[1] "wsvd [793]=0.00396243929638538"
[1] "wsvd [794]=0.0039624353078582"
[1] "wsvd [795]=0.003962431356864"
[1] "wsvd [796]=0.00396242744303526"
[1] "wsvd [797]=0.00396242356600823"
[1] "wsvd [798]=0.00396241972542285"
[1] "wsvd [799]=0.00396241592092273"
[1] "wsvd [800]=0.00396241215215513"
[1] "wsvd [801]=0.00396240841877089"
[1] "wsvd [802]=0.0039624047204244"
[1] "wsvd [803]=0.00396240105677357"
[1] "wsvd [804]=0.00396239742747979"
[1] "wsvd [805]=0.0039623938322079"
[1] "wsvd [806]=0.00396239027062612"
[1] "wsvd [807]=0.00396238674240608"
[1] "wsvd [808]=0.0039623832472227"
[1] "wsvd [809]=0.00396237978475424"
[1] "wsvd [810]=0.00396237635468218"
[1] "wsvd [811]=0.00396237295669127"
[1] "wsvd [812]=0.00396236959046944"
[1] "wsvd [813]=0.00396236625570777"
[1] "wsvd [814]=0.00396236295210048"
[1] "wsvd [815]=0.00396235967934488"
[1] "wsvd [816]=0.00396235643714135"
[1] "wsvd [817]=0.00396235322519329"
[1] "wsvd [818]=0.0039623500432071"
[1] "wsvd [819]=0.00396234689089216"
[1] "wsvd [820]=0.00396234376796077"
[1] "wsvd [821]=0.00396234067412813"
[1] "wsvd [822]=0.00396233760911234"
[1] "wsvd [823]=0.0039623345726343"
[1] "wsvd [824]=0.00396233156441778"
[1] "wsvd [825]=0.00396232858418928"
[1] "wsvd [826]=0.00396232563167809"
[1] "wsvd [827]=0.00396232270661621"
[1] "wsvd [828]=0.00396231980873835"
[1] "wsvd [829]=0.00396231693778187"
[1] "wsvd [830]=0.0039623140934868"
[1] "wsvd [831]=0.00396231127559574"
[1] "wsvd [832]=0.00396230848385393"
[1] "wsvd [833]=0.00396230571800912"
[1] "wsvd [834]=0.00396230297781162"
[1] "wsvd [835]=0.00396230026301424"
[1] "wsvd [836]=0.00396229757337227"
[1] "wsvd [837]=0.00396229490864344"
[1] "wsvd [838]=0.00396229226858794"
[1] "wsvd [839]=0.00396228965296833"
[1] "wsvd [840]=0.00396228706154956"
[1] "wsvd [841]=0.00396228449409894"
[1] "wsvd [842]=0.00396228195038609"
[1] "wsvd [843]=0.00396227943018296"
[1] "wsvd [844]=0.00396227693326375"
[1] "wsvd [845]=0.00396227445940495"
[1] "wsvd [846]=0.00396227200838524"
[1] "wsvd [847]=0.00396226957998556"
[1] "wsvd [848]=0.00396226717398899"
[1] "wsvd [849]=0.0039622647901808"
[1] "wsvd [850]=0.0039622624283484"
[1] "wsvd [851]=0.00396226008828131"
[1] "wsvd [852]=0.00396225776977115"
[1] "wsvd [853]=0.00396225547261163"
[1] "wsvd [854]=0.00396225319659848"
[1] "wsvd [855]=0.00396225094152951"
[1] "wsvd [856]=0.00396224870720449"
[1] "wsvd [857]=0.00396224649342522"
[1] "wsvd [858]=0.00396224429999545"
[1] "wsvd [859]=0.00396224212672089"
[1] "wsvd [860]=0.00396223997340917"
[1] "wsvd [861]=0.00396223783986984"
[1] "wsvd [862]=0.00396223572591433"
[1] "wsvd [863]=0.00396223363135594"
[1] "wsvd [864]=0.00396223155600984"
[1] "wsvd [865]=0.003962229499693"
[1] "wsvd [866]=0.00396222746222424"
[1] "wsvd [867]=0.00396222544342414"
[1] "wsvd [868]=0.00396222344311507"
[1] "wsvd [869]=0.00396222146112117"
[1] "wsvd [870]=0.0039622194972683"
[1] "wsvd [871]=0.00396221755138404"
[1] "wsvd [872]=0.0039622156232977"
[1] "wsvd [873]=0.00396221371284025"
[1] "wsvd [874]=0.00396221181984434"
[1] "wsvd [875]=0.00396220994414427"
[1] "wsvd [876]=0.00396220808557597"
[1] "wsvd [877]=0.003962206243977"
[1] "wsvd [878]=0.0039622044191865"
[1] "wsvd [879]=0.00396220261104522"
[1] "wsvd [880]=0.00396220081939548"
[1] "wsvd [881]=0.00396219904408111"
[1] "wsvd [882]=0.00396219728494753"
[1] "wsvd [883]=0.00396219554184165"
[1] "wsvd [884]=0.0039621938146119"
[1] "wsvd [885]=0.00396219210310818"
[1] "wsvd [886]=0.00396219040718189"
[1] "wsvd [887]=0.00396218872668587"
[1] "wsvd [888]=0.00396218706147443"
[1] "wsvd [889]=0.00396218541140327"
[1] "wsvd [890]=0.00396218377632954"
[1] "wsvd [891]=0.00396218215611178"
[1] "wsvd [892]=0.00396218055060991"
[1] "wsvd [893]=0.00396217895968524"
[1] "wsvd [894]=0.00396217738320041"
[1] "wsvd [895]=0.00396217582101944"
[1] "wsvd [896]=0.00396217427300765"
[1] "wsvd [897]=0.00396217273903169"
[1] "wsvd [898]=0.00396217121895952"
[1] "wsvd [899]=0.00396216971266039"
[1] "wsvd [900]=0.00396216822000481"
[1] "wsvd [901]=0.00396216674086458"
[1] "wsvd [902]=0.00396216527511273"
[1] "wsvd [903]=0.00396216382262355"
[1] "wsvd [904]=0.00396216238327253"
[1] "wsvd [905]=0.0039621609569364"
[1] "wsvd [906]=0.00396215954349308"
[1] "wsvd [907]=0.00396215814282167"
[1] "wsvd [908]=0.00396215675480247"
[1] "wsvd [909]=0.00396215537931694"
[1] "wsvd [910]=0.00396215401624767"
[1] "wsvd [911]=0.00396215266547842"
[1] "wsvd [912]=0.00396215132689406"
[1] "wsvd [913]=0.00396215000038061"
[1] "wsvd [914]=0.00396214868582516"
[1] "wsvd [915]=0.00396214738311592"
[1] "wsvd [916]=0.00396214609214218"
[1] "wsvd [917]=0.0039621448127943"
[1] "wsvd [918]=0.00396214354496372"
[1] "wsvd [919]=0.0039621422885429"
[1] "wsvd [920]=0.00396214104342539"
[1] "wsvd [921]=0.00396213980950572"
[1] "wsvd [922]=0.00396213858667948"
[1] "wsvd [923]=0.00396213737484326"
[1] "wsvd [924]=0.00396213617389465"
[1] "wsvd [925]=0.00396213498373222"
[1] "wsvd [926]=0.00396213380425554"
[1] "wsvd [927]=0.00396213263536513"
[1] "wsvd [928]=0.00396213147696251"
[1] "wsvd [929]=0.0039621303289501"
[1] "wsvd [930]=0.00396212919123129"
[1] "wsvd [931]=0.00396212806371042"
[1] "wsvd [932]=0.00396212694629271"
[1] "wsvd [933]=0.00396212583888433"
[1] "wsvd [934]=0.00396212474139235"
[1] "wsvd [935]=0.00396212365372471"
[1] "wsvd [936]=0.00396212257579027"
[1] "wsvd [937]=0.00396212150749875"
[1] "wsvd [938]=0.00396212044876074"
[1] "wsvd [939]=0.00396211939948769"
[1] "wsvd [940]=0.00396211835959191"
[1] "wsvd [941]=0.00396211732898655"
[1] "wsvd [942]=0.00396211630758559"
[1] "wsvd [943]=0.00396211529530385"
[1] "wsvd [944]=0.00396211429205694"
[1] "wsvd [945]=0.00396211329776132"
pois_llik(x,x.plra1.naive2$fit)
[1] -13163.59
pois_llik(x,x.plra1.true2$fit)
[1] -13145.19

Session information

sessionInfo()
R version 3.5.1 (2018-07-02)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: OS X El Capitan 10.11.6

Matrix products: default
BLAS: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRlapack.dylib

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     

loaded via a namespace (and not attached):
 [1] workflowr_1.1.1   Rcpp_0.12.19      NNLM_0.4.2       
 [4] digest_0.6.17     rprojroot_1.3-2   R.methodsS3_1.7.1
 [7] backports_1.1.2   git2r_0.23.0      magrittr_1.5     
[10] evaluate_0.11     stringi_1.2.4     whisker_0.3-2    
[13] R.oo_1.22.0       R.utils_2.7.0     rmarkdown_1.10   
[16] tools_3.5.1       stringr_1.3.1     yaml_2.2.0       
[19] compiler_3.5.1    htmltools_0.3.6   knitr_1.20       



This reproducible R Markdown analysis was created with workflowr 1.1.1