Last data update: 2014.03.03

R: Score test of association
ASP.Score R Documentation

Score test of association

Description

Calculate score statistics and the associated P-value for each SNPs

Usage

 ASP.Score(Tem_Gen, Index_Gen, IBD) 

Arguments

Tem_Gen

Genotypes of controls (denoted by the number of alternative allele)

Index_Gen

Genotypes of Index cases

IBD

IBD states for each affected sib pair

Details

Give the values of statistic and p-value of the association score test.

Value

List of 2 vectors of length the number of SNPs:

Value

Statistic values for each SNPs

Pvalue

P-values of the score test for each SNPs

Author(s)

Claire Dandine-Roulland

References

Perdry, Herve, Muller-Myhsok, Bertram, et Clerget-Darpoux, Francoise. Using affected sib-pairs to uncover rare disease variants. Hum Hered, 2013.

See Also

ASP.Selection

Examples

data(ASPData)
ASP.Score(ASPData$Control, ASPData$Index, ASPData$IBD )

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(ASPBay)
Loading required package: hexbin
Loading required package: Rcpp
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ASPBay/ASP.Score.Rd_%03d_medium.png", width=480, height=480)
> ### Name:  ASP.Score 
> ### Title: Score test of association
> ### Aliases: ' ASP.Score '
> ### Keywords: Affected sib pairs Score test
> 
> ### ** Examples
> 
> data(ASPData)
> ASP.Score(ASPData$Control, ASPData$Index, ASPData$IBD )
$Value
      SNP1       SNP2       SNP3       SNP4       SNP5       SNP6       SNP7 
10.5697498 -2.6353276 -0.7334053 -0.7617181  6.3666141 -2.4239793 -0.3127143 
      SNP8       SNP9      SNP10      SNP11      SNP12      SNP13      SNP14 
-2.5229966  9.2774756 -3.7236812 -0.9274910  9.4499055 11.2370354  9.4499055 
     SNP15      SNP16      SNP17      SNP18      SNP19      SNP20      SNP22 
11.2370354 -0.9820782 -0.9820782 -1.5493967 -0.4868593 -0.1776845 -0.6855409 

$Pvalue
        SNP1         SNP2         SNP3         SNP4         SNP5         SNP6 
0.000000e+00 8.405611e-03 4.633113e-01 4.462283e-01 1.932465e-10 1.535148e-02 
        SNP7         SNP8         SNP9        SNP10        SNP11        SNP12 
7.544978e-01 1.163595e-02 0.000000e+00 1.963389e-04 3.536717e-01 0.000000e+00 
       SNP13        SNP14        SNP15        SNP16        SNP17        SNP18 
0.000000e+00 0.000000e+00 0.000000e+00 3.260613e-01 3.260613e-01 1.212864e-01 
       SNP19        SNP20        SNP22 
6.263581e-01 8.589708e-01 4.930027e-01 

> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>