Nothing
`makeSplineData.fnc` <-
function(intr=0) {
set.seed(314)
if (!intr) {
stdev = 2.5
p = 2
}
else {
stdev = 0.5
p = 6
}
x = seq(0,p*pi,length=100)+2
if (intr==0) {
y = 30 + cos(x)
dfr = data.frame(y=y,X=x)
} else {
if (intr==1) {
z = seq(2,0.2,length=100)
y = 30 + cos(x)*z
dfr = data.frame(y=y,X=x,z=z)
} else {
z = 1
y = 30 + cos(x)*z
dfr = data.frame(y=y,X=x,z=1)
z = 2
y = 30 + cos(x)*z
dfr = rbind(dfr, data.frame(y=y,X=x,z=2))
z = 3
y = 30 + cos(x)*z
dfr = rbind(dfr, data.frame(y=y,X=x,z=3))
z = 4
y = 30 + cos(x)*z
dfr = rbind(dfr, data.frame(y=y,X=x,z=4))
}
}
nsubj = 10
ranefs = stats::rnorm(10, 0, 3)
dfr$Subject = rep("S1", nrow(dfr))
dfr$Ranef = rep(ranefs[1], nrow(dfr))
dfrCopy = dfr
for (i in 2:10){
tmp = dfrCopy
tmp$Subject = rep(paste("S", i, sep=""), nrow(tmp))
tmp$Ranef = ranefs[i]
dfr = rbind(dfr, tmp)
}
dfr$Subject = as.factor(dfr$Subject)
dfr$Error = stats::rnorm(nrow(dfr), 0, stdev)
dfr$Y = dfr$y+dfr$Ranef+dfr$Error
return(dfr)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.