GroupPlot: ~~PLOT all graph in the group ... ~~

In qPCR.CT: qPCR data analysis and plot package

Description

~~ Plot all in barplot with Error Bar(se or sd) ~~

Usage

GroupPlot(filename, con_group_name, tr_group_name, con_gene_name, repeatnum, ...)

Arguments

filename	~~specific filename such as c: pcr.csv~~
con_group_name	~~set con_group_name ~~
tr_group_name	~~set tr_group_name here~~
con_gene_name	~~set con_gene_name here~~
repeatnum	~~set repeatnum here~~
...

Author(s)

~~yz pan, xy yan, jx li~~

Examples

##---- 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)
   }
  }
  }

Example output

Loading required package: RColorBrewer
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))
    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]
            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)
        }
    }
}

qPCR.CT documentation built on May 30, 2017, 3:51 a.m.