Description Usage Arguments Details Value Author(s) References See Also Examples
This procedure quantifies the accuracy with which one can predict a given genotypes (SNPs or SAAPs) from the corresponding phenotypes using linear support vector machines (LSVM).
1  runGenphenSvm(genotype, phenotype, cv.fold, cv.steps, hdi.level)

genotype 
Character matrix or data frame, containing SNPs/SAAPs as columns or alternatively as DNAMultipleAlignment or AAMultipleAlignment Biostrings object. 
phenotype 
Numerical vector, where each element is a measured phenotype corresponding to the observations of the genotype data. 
cv.fold 
The crossvalidation fraction (0, 1) of the data which is used to train the classifier (recommended = 0.66). The remaining fraction (1cv.fold) of the data is used to test the classifier. 
cv.steps 
Number of steps in the crossvalidation to be performed to estimate the classification accuracy and the corresponding highest density intervals(recommended >= 100). 
hdi.level 
Highest density interval (default = 0.99). 
This procedure takes two types of data as input: first a genotype data composed of a set of single nucleotide polymorphisms (SNPs) or alternatively single amino acid polymorphisms (SAAPs), each of which is represented by a column of character amino acids; second a numerical phenotype vector, where the elements sorted to correspond to the rows of the genotype data. This method quantifies the association between the polymorphic site (SNP or SAAP) and the phenotype via a classification analysis using linear support vector machines. The analysis results in a classification accuracy score between 0 and 1, where 1 indicates a perfect association between the genotype and the phenotype. To validate the classification accuracy, the tool also computes the Cohen's kappa statistics (Cohen 1960) which compares the observed classification accuracy with the expected classification accuracy. If the expected and observed classification accuracies are in concordance, the computed association can be taken seriously, otherwise it can be discarded as noise.
The function runGenphenSvm also computes statistics such as Cohen's d (effect size) and the Pvalue resulting from a twosample Ttest, allowing the user to compare the linear support vector based results with those computed with simpler techniques which are frequently used for genetic association studies.
Five classes of results are computed for each SAAP with respect to the phenotype, resulting in a 18 element vector which is stored as a row in the final data frame:
effect.size, effect.CI.low, effect.CI.high 
Cohen's effect size and CI. 
ca, ca.hdi.low, ca.hdi.high, ca.hdi.length 
Mean classification accuracy and its HDI. 
kappa, kappa.hdi.low, kappa.hdi.high, kappa.hdi.length 
Cohen's kappa statistics and its HDI. 
site, g.1, g.2, count.1, count.2 
General information about the genotype. 
t.test.pvalue 
Pvalue score from an twosample Ttest. 
Simo Kitanovski <[email protected]>
Cortes, Corinna, and Vladimir Vapnik. Supportvector networks. Machine learning 20.3 (1995): 273297.
Cohen, Jacob. Statistical power analysis for the behavior science. Lawrance Eribaum Association (1988).
Cohen, Jacob. A coefficient of agreement for nominal scales (1960).
runGenpenRf, runGenpenBayes, plotGenphenRfSvm, plotGenphenBayes, plotSpecificGenotype, plotManhattan
1 2 3 4 5 6 7 8  data(genotype.saap)
#or data(genotype.saap.msa) in this case you cannot subset genotype.saap[, 1:5]
data(phenotype.saap)
genphen.svm < runGenphenSvm(genotype = genotype.saap[, 1:5],
phenotype = phenotype.saap,
cv.fold = 0.66,
cv.steps = 100,
hdi.level = 0.99)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.