fExpr = expression(x1+x2)
x.mu = c(1,1)
x.u = c(0.1,1.0)
x.pdf = c('norm','norm')
x.df = c(Inf, Inf)
names(x.u) = names(x.mu) = names(x.pdf) = names(x.df) = c('x1','x2')
test_that('vgSA results fit Combination of Variances for linear model',
{
nrun = 1000
Sref = gumCV(fExpr,x.mu,x.u,silent=TRUE)
S = gumS1(fExpr,x.mu,x.u,x.pdf,x.df,nrun=nrun,silent=TRUE)
SA = vgSA(fExpr,x.mu,x.u,S$X,S$Y,silent=TRUE)
expect_equal(SA$vg[['x1']],Sref$anova[['x1']],
tolerance = 3/sqrt(nrun))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.