qpcrREPEATED {rtpcr}

Fold change (\Delta \Delta C_T method) analysis of repeated measure qPCR data

Description

qpcrREPEATED function performs fold change (\Delta \Delta C_T method) analysis of observations repeatedly taken over different time courses, Data may be obtained over time from a uni- or multi-factorial experiment. The bar plot of the fold changes (FC) values along with the standard error (se) or confidence interval (ci) is also returned by the qpcrREPEATED function.

Usage

qpcrREPEATED(
  x,
  numberOfrefGenes,
  factor,
  block = NULL,
  width = 0.5,
  fill = "#BFEFFF",
  y.axis.adjust = 1,
  y.axis.by = 1,
  ylab = "Fold Change",
  xlab = "none",
  fontsize = 12,
  fontsizePvalue = 7,
  axis.text.x.angle = 0,
  axis.text.x.hjust = 0.5,
  x.axis.labels.rename = "none",
  letter.position.adjust = 0,
  p.adj = "none",
  errorbar = "se"
)

Arguments

x

input data frame in which the first column is id, followed by the factor(s) which include at least time factor. The first level of time factor in data set is used as calibrator or reference level. Additional factor(s) may also be present. Other columns are efficiency and Ct values of target and reference genes. In the "id" column, a unique number is assigned to each individual from which samples have been takes over time, for example in the data_repeated_measure_1, all the three number 1 indicate one individual which has been sampled over three different time courses. To prepare a data frame from a repeated measure analysis, please refer to the vignette.

numberOfrefGenes

number of reference genes which is 1 or 2 (Up to two reference genes can be handled). as reference or calibrator which is the reference level or sample that all others are compared to. Examples are untreated of time 0. The FC value of the reference or calibrator level is 1 because it is not changed compared to itself. If NULL, the first level of the main factor column is used as calibrator.

factor

the factor for which the FC values is analysed. The first level of the specified factor in the input data frame is used as calibrator.

block

column name of the block if there is a blocking factor (for correct column arrangement see example data.). Block effect is usually considered as random and its interaction with any main effect is not considered.

width

a positive number determining bar width in the output bar plot.

fill

specify the fill color for the columns in the bar plot. If a vector of two colors is specified, the reference level is differentialy colored.

y.axis.adjust

a negative or positive value for reducing or increasing the length of the y axis.

y.axis.by

determines y axis step length

ylab

the title of the y axis

xlab

the title of the x axis

fontsize

font size of the plot

fontsizePvalue

font size of the pvalue labels

axis.text.x.angle

angle of x axis text

axis.text.x.hjust

horizontal justification of x axis text

x.axis.labels.rename

a vector replacing the x axis labels in the bar plot

letter.position.adjust

adjust the distance between the signs and the error bars.

p.adj

Method for adjusting p values

errorbar

Type of error bar, can be se or ci.

Details

The qpcrREPEATED function performs fold change (FC) analysis of observations repeatedly taken over time. The intended factor (could be time or any other factor) is defined for the analysis by the factor argument, then the function performs FC analysis on its levels so that the first levels (as appears in the input data frame) is used as reference or calibrator. the function returns FC values along with confidence interval and standard error for the FC values.

Value

A list with 5 elements:

Final_data

Input data frame plus the weighted Delat Ct values (wDCt)

lm

lm of factorial analysis-tyle

ANOVA_table

ANOVA table

FC Table

Table of FC values, significance, confidence interval and standard error for the selected factor levels.

Bar plot of FC values

Bar plot of the fold change values for the main factor levels.

Author(s)

Ghader Mirzaghaderi

Examples

qpcrREPEATED(data_repeated_measure_1,
            numberOfrefGenes = 1,
            factor = "time")

qpcrREPEATED(data_repeated_measure_2,
             numberOfrefGenes = 1,
             factor = "time")