User Tools


Differences

This shows you the differences between two versions of the page.

Link to this comparison view

output [2016/04/13 18:02]
output [2016/04/13 18:02] (current)
Line 1: Line 1:
 +===== SEQPower Output =====
 +==== Output Contents ====
 +Two types of output are generated in power calculation process: association group statistics as well as summary site specific attributes.
  
 + *  **Association group statistics**:​ summary statistic for the group of variants on which power analysis was carried out, including power or required sample size, cumulative MAF (by case control group, by functionality categories),​ data properties (missing / error genotype call rates), etc.
 + *  **Site specific attributes summary**: attributes of each variant sites within the group, including MAF, function score, etc, from input data, and model specific attributes such as odds ratio, population attributable risk, MAF in cases/ctrls for each variant site, etc.
 +
 +==== Output Format ====
 +Output are written to two separate ''​.csv''​ files and to a [[http://​sqlite.org|SQLite database]]. The ''​.csv''​ output will be overwritten each time the program is executed, yet the results in database will be updated with the latest calculation while keeping all existing information it has. The ''​.csv''​ can be considered as an intermediate output to readily view current calculation results. The database contains all calculations ever performed, which serves as useful resource for future references.
 +
 +=== Extract information from output files ===
 +''​spower show''​ command provides features to extract information from output text or database files.
 +
 +== Text file ==
 +For example,
 +
 +<code bash>
 +spower show test.csv
 +</​code>​\\
 +shows all column names in the output file ''​test.csv''​. To view values of particular columns, for example "​method"​ and "​power":​
 +
 +<code bash>
 +spower show test.csv method power
 +</​code>​\\
 +<hidden initialState="​visible"​ -noprint>​
 +  +--------+-------+
 +  | method | power |
 +  +--------+-------+
 +  |  BQt1  |  0.82 |
 +  |  BQt5  |  0.78 |
 +  +--------+-------+
 +</​hidden>​\\
 +Wildcard character "​*"​ is supported, for example:
 +
 +<code bash>
 +spower show test.csv *deleterious*
 +</​code>​\\
 +<hidden initialState="​visible"​ -noprint>​
 +  +------------------+-------------------------+----------------------+
 +  | cmaf_deleterious | cmaf_deleterious_median | cmaf_deleterious_std |
 +  +------------------+-------------------------+----------------------+
 +  | 0.00550143284332 |     ​0.00133893662132 ​   |   ​0.00416249625512 ​  |
 +  | 0.00179535530657 |     ​0.00157300161663 ​   |  0.000222353711192 ​  |
 +  +------------------+-------------------------+----------------------+
 +</​hidden>​\\
 +== Database file ==
 +For example
 +
 +<code bash>
 +spower show test.SEQPower.DB
 +</​code>​\\
 +shows all table names in the database, e.g.,
 +
 +<hidden initialState="​visible"​ -noprint>​
 +  +--------+
 +  | TABLES |
 +  +--------+
 +  |  BLNR  |
 +  |  ELNR  |
 +  |  LNR   |
 +  | LOGIT  |
 +  |  PAR   |
 +  +--------+
 +</​hidden>​\\
 +To show column names in a particular table:
 +
 +<code bash>
 +spower show test.SEQPowerDB LOGIT
 +</​code>​\\
 +To show values for a particular column for a table:
 +
 +<code bash>
 +spower show test.SEQPowerDB LOGIT power
 +</​code>​\\
 +Note that wildcard character is not supported for database output. However an additional ''​%%--%%condition''​ command is supported for database. For example to select only the tests having power greater than 0.8:
 +
 +<code bash>
 +spower show test.SEQPowerDB LOGIT title method power --condition "where power > 0.8"
 +</​code>​\\
 +<hidden initialState="​visible"​ -noprint>​
 +  +-------+---------+----------------+
 +  | title |  method |     ​power ​     |
 +  +-------+---------+----------------+
 +  |   ​EA ​ | default | 0.999938410481 |
 +  +-------+---------+----------------+
 +</​hidden>​\\
 +=== Graphical interface for SEQPower databases ===
 +We recommend the use of an SQLite GUI software such as [[http://​sqlitestudio.pl/​|SQLite Studio]] to browse the result SQLite database. These programs display the data as spreadsheets,​ and provides simple "click & run" features to filter or query the current database.
 +
 +Below is a screen shot of SQLite Studio displaying the power analysis result of a whole exome scan (left) and the tests with power greater than 80% (right), extracted via running SQLite query from the GUI.
 +
 +{{:​sqlite-manager-demo.png?​900.0}}