##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function(filename,
con_group_name,tr_group_name,
con_gene_name, repeatnum, ... ){
rawdata<-read.csv(filename)
unknown<- subset(rawdata, Type != "NTC")
suppressWarnings(aggdata <-aggregate.data.frame(unknown, by=list(unknown$Fluor),FUN=mean, na.rm=TRUE))
#print(aggdata[1])
groups<-paste(aggdata[,1])
lengr <- length(groups)
temp<- data.frame()
i<- 1
for(i in 1: lengr)
{
if(groups[i] != con_gene_name ){
tr_gene_name <- groups[i] # assign group names
temp <- Calc.Expr(filename=filename, con_group_name, tr_group_name,
con_gene_name = con_gene_name,tr_gene_name = tr_gene_name, repeatnum = repeatnum)
plot.CT(temp)
}
}
}Run the code above in your browser using DataLab