Nothing
familyDesign_c <-
function(n=1000, affectnum=0, m.carrier= 0, dominant.m = TRUE, dominant.s = TRUE,
interaction = TRUE, depend=1,
base.dist="Weibull", frailty.dist="gamma",
vbeta= c(-1.126, 2.55, 1.6), parms=c(0.016, 3), variation="none",
allelefreq=c(0.02, 0.2), mrate=0.1, probandage=c(45,1.5),
agemin=20, agemax=100)
{
data<-numeric()
cumind<-0
i<- 1
j<- 0
while (i <= n) {
j <- j + 1
dat<-try(familyStructure_c(i,cumind=cumind, m.carrier=m.carrier,
base.dist=base.dist, frailty.dist=frailty.dist, interaction=interaction,
depend=depend, parms=parms, vbeta=vbeta, dominant.m=dominant.m, dominant.s=dominant.s,
variation=variation, allelefreq=allelefreq, mrate=mrate, probandage=probandage,
agemin=agemin, agemax=agemax))
if(is.null(attr(dat, "class"))){
# At least one parent in first gen and two sibs in the second gen should be affected
#[,7]=generation, [,6]=proband, [,13]=status
# affectnum = 3 "cli+"
# affectnum = 2 "twostage"
# affectnum = 1 "pop+"
# affectnum = 0 "pop"
if(affectnum==3) until <- ifelse(sum(dat[dat[,"generation"]==1,"status"] %in% c(1,2))>0 & sum(dat[dat[,"generation"]==2,"status"] %in% c(1,2)) > 1, TRUE, FALSE)
else if(affectnum==2) until <- ifelse(sum(dat[, "status"] %in% c(1,2)) >= affectnum, TRUE, FALSE)
else if(affectnum==1) until <- ifelse( dat[dat[,"proband"]==1,"status"] %in% c(1,2), TRUE, FALSE)
else until <- TRUE
if(!is.null(dim(dat))){
if(nrow(dat)>0 ){
if(until){
data<-rbind(data, dat)
cumind<-cumind+nrow(dat)
i<-i+1
}
}
}
} # close "is.null(attr(dat, "class"))"
} # close while
data
}
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.