citrus.generateRegularizationThresholds | R Documentation |
Generate a range of regularization thresholds for model construction
citrus.generateRegularizationThresholds.classification(features, labels,
modelType, n = 100, ...)
citrus.generateRegularizationThresholds.continuous(features, labels, modelType,
n = 100, ...)
citrus.generateRegularizationThresholds(features, labels, modelType, family,
n = 100, ...)
features |
Features used to construct model |
labels |
Endpoint lables for samples and features |
modelType |
Method used to construct endpoint model. Valid options are: |
n |
Number of regularization thresholds to generate |
... |
Other arguments passed to model-fitting methods |
family |
Model family. Valid options are: |
A vector of regularization threshold values.
Robert Bruggner
# Where the data lives
dataDirectory = file.path(system.file(package = "citrus"),"extdata","example1")
# Create list of files to be analyzed
fileList = data.frame("unstim"=list.files(dataDirectory,pattern=".fcs"))
# Read the data
citrus.combinedFCSSet = citrus.readFCSSet(dataDirectory,fileList)
# List of columns to be used for clustering
clusteringColumns = c("Red","Blue")
# Cluster data
citrus.clustering = citrus.cluster(citrus.combinedFCSSet,clusteringColumns)
# Large enough clusters
largeEnoughClusters = citrus.selectClusters(citrus.clustering)
# Build features
abundanceFeatures = citrus.calculateFeatures(citrus.combinedFCSSet,clusterAssignments=citrus.clustering$clusterMembership,clusterIds=largeEnoughClusters)
# List disease group of each sample
labels = factor(rep(c("Healthy","Diseased"),each=10))
# Calculate regularization thresholds
regularizationThresholds = citrus.generateRegularizationThresholds(abundanceFeatures,labels,modelType="pamr",family="classification")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.