scatterplotSexGenotypeBatchResult: Method "scatterplotSexGenotypeBatchResult"
In PhenStat: Statistical analysis of phenotypic data

Description Usage Arguments Value Author(s) References See Also Examples

Graph function for the phenotypic dataset. Creates a scatterplot split by sex, genotype and batch. refGenotype data points are shown in black and the testGenotype data points are shown in red. It is based on PhenTestResult object.

Note: the batches are not ordered with time but allow assessment of how the testGenotype data lie relative to the variation within the refGenotype data.

1
2
3

    scatterplotSexGenotypeBatchResult(phenTestResult,
                                  graphingName = NULL,
                                  outputMessages = TRUE)

`phenTestResult`	instance of the `PhenTestResult` class; mandatory argument
`graphingName`	a character string defining the name to be used on the graph for the dependent variable; mandatory argument
`outputMessages`	flag: "FALSE" value to suppress output messages; "TRUE" value to show output messages; default value TRUE

Empty.

Natalja Kurbatova, Natasha Karp, Jeremy Mason

Karp N, Melvin D, Sanger Mouse Genetics Project, Mott R (2012): Robust and Sensitive Analysis of Mouse Knockout Phenotypes. PLoS ONE 7(12): e52410. doi:10.1371/journal.pone.0052410

West B, Welch K, Galecki A (2007): Linear Mixed Models: A practical guide using statistical software New York: Chapman & Hall/CRC 353 p.

PhenList

    file <- system.file("extdata", "test1.csv", package="PhenStat")
    test <- PhenStat:::PhenList(dataset=read.csv(file,na.strings = '-'),
            testGenotype="Sparc/Sparc")
    result <- PhenStat:::testDataset(test,
            depVariable="Lean.Mass")
    # box plot for dataset with two sexes: males and females
    PhenStat:::scatterplotSexGenotypeBatchResult(result,
            graphingName="BMC" )

    file <- system.file("extdata", "test4.csv", package="PhenStat")
    test_1sex <- PhenStat:::PhenList(dataset=read.csv(file,na.strings = '-'),
            testGenotype="Mysm1/+")
    result_1sex <- PhenStat:::testDataset(test_1sex,
            depVariable="Lean.Mass")
    # box plot for females only dataset
    PhenStat:::scatterplotSexGenotypeBatchResult(result_1sex,
            graphingName="Lean Mass")