1 | computeProteinPI(keyAA)
|
keyAA |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (keyAA)
{
AspNumber <- keyAA[1]
GluNumber <- keyAA[2]
CysNumber <- keyAA[3]
TyrNumber <- keyAA[4]
HisNumber <- keyAA[5]
LysNumber <- keyAA[6]
ArgNumber <- keyAA[7]
NQ <- 0
QN1 <- 0
QN2 <- 0
QN3 <- 0
QN4 <- 0
QN5 <- 0
QP1 <- 0
QP2 <- 0
QP3 <- 0
QP4 <- 0
pH <- 6.5
pHprev <- 0
pHnext <- 14
X <- 0
E <- 0.01
temp <- 0
tabl <- 650
while (TRUE) {
QN1 <- tab[tabl, 1]
QN2 <- -AspNumber/tab[tabl, 2]
QN3 <- -GluNumber/tab[tabl, 3]
QN4 <- -CysNumber/tab[tabl, 4]
QN5 <- -TyrNumber/tab[tabl, 5]
QP1 <- HisNumber/tab[tabl, 6]
QP2 <- tab[tabl, 7]
QP3 <- LysNumber/tab[tabl, 8]
QP4 <- ArgNumber/tab[tabl, 9]
NQ <- QN1 + QN2 + QN3 + QN4 + QN5 + QP1 + QP2 + QP3 +
QP4
if (NQ < 0) {
temp <- pH
pH <- pH - ((pH - pHprev)/2)
pHnext <- temp
}
else {
temp <- pH
pH <- pH + ((pHnext - pH)/2)
pHprev <- temp
}
if ((pH - pHprev < E) & (pHnext - pH < E)) {
break
}
n <- pH * 100
obc <- n
ulamek <- n - obc
if (ulamek >= 0.5) {
obc <- obc + 1
}
tabl <- obc
pH <- obc
pH <- pH/100
}
return(pH)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.