R: Predict oligomerization state of coiled-coil sequences.
LOGICOIL
R Documentation
Predict oligomerization state of coiled-coil sequences.
Description
Function for predicting the oligomeric state of a coiled-coil sequence. This function allows the user to separate between antiparallel dimer, parallel dimers, trimers and higher-order coiled coils.
Usage
LOGICOIL(id, seq, reg, plot.result = FALSE)
Arguments
id
A string identifier for the input coiled-coil sequence.
seq
The input coiled-coil sequence, which might a string or vector of amino-acids. The input must contain at least 14 characters. Valid characters should be uppercase letters and illegal terms include ‘B’, ‘J’, ‘O’, ‘U’, ‘Z’. Invalid characters will not be tolerated and their use will result in a failure of the program. If the amino-acid in a coiled-coil sequence is unknown, it is permitted to use the identifier "X" instead
reg
The register assignments for the amino-acids in the coiled-coil sequence. The input may be either a string or a vector and must always have the same length as the seq argument. If this is not the case then the program will fail. Valid characters are the lowercase letters ‘a’ to ‘g’. Register characters are not required to be in proper order
plot.result
A logical value indicating whether the results of the predictions should be plotted on a PNG file (default=FALSE).
Value
The function LOGICOIL is the parent function in the LOGICOIL package. It is used to call the LOGICOIL prediction algorithm in order to predict the oligomeric state of a new coiled-coil sequence. The LOGICOIL function returns a vector with 4 items, which gives the odd ratio that the input sequence belongs to an antiparallel dimer, parallel, trimer or higher-order coiled coiled. If a given sequence scores f(x)==1, then there is no evidence for a particular oligomeric state. If f(x) > 1 then there is strong evidence for a particular oligimeric state, and if f(x) < 1 then there is against a particular oligimeric state.
Thomas L. Vincent, Peter J. Green and Dek N. Woolfson. "LOGICOIL-multi-state prediction of coiled-coil oligomeric state", 29(1), pp69-76, Bioinformatics, (2013). DOI: 10.1093/bioinformatics/bts648
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(LOGICOIL)
Loading required package: nnet
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/LOGICOIL/LOGICOIL.Rd_%03d_medium.png", width=480, height=480)
> ### Name: LOGICOIL
> ### Title: Predict oligomerization state of coiled-coil sequences.
> ### Aliases: LOGICOIL
>
> ### ** Examples
>
> library(nnet)
> data(pssm)
> data(Model_Parameters)
> data(LOGICOILfit)
>
> score <- LOGICOIL("GCN4wt", "MKQLEDKVEELLSKNYHLENEVARLKKLV",
+ "abcdefgabcdefgabcdefgabcdefga", plot.result=FALSE)
Estimating oligomeric state of coiled-coil sequences
Analysis complete!
>
> # print output
> print(score)
1 2 3 4
0.6951839 1.5669293 0.5670830 2.0606170
>
>
>
>
>
> dev.off()
null device
1
>