lib/rnautil.c File Reference

#include "rnautil.h"
#include "common.h"

Include dependency graph for rnautil.c:

Go to the source code of this file.

Functions
void	dna2rna (char *s)
bool	rnaPair (char a, char b)
void	reverseFold (char *s)
void	fold2pairingList (char *fold, int len, int **p2pairList)
void	mkPairPartnerSymbols (int *pairList, char *pairSymbols, int size)
char *	projectString (char *s, char *ref, char refChar, char insertChar)
char *	gapAdjustFold (char *s, char *ref)
int *	projectIntArray (int *in, char *ref, char refChar, int insertInt)
int *	gapIntArrayAdjust (int *in, char *ref)
void	markCompensatoryMutations (char *s, char *ref, int *pairList, int *markList)
int	assignBin (double val, double minVal, double maxVal, int binCount)
Variables
const char *	RNA_PAIRS [] = {"AU","UA","GC","CG","GU","UG",0}

---

Function Documentation

int assignBin	(	double	val,
		double	minVal,
		double	maxVal,
		int	binCount
	)

Definition at line 203 of file rnautil.c.

{
double range = maxVal - minVal;
int maxBin   = binCount - 1;
int level    = (int) ( (val-minVal)*maxBin/range);
if (level <= 0) level = 0;
if (level > maxBin) level = maxBin;
return level;
}

void dna2rna

(

char *

s

)

Definition at line 8 of file rnautil.c.

Referenced by rnaPair().

{
for (;*s;s++)
    {
    if (*s == 't')
        *s = 'u';
    if (*s == 'T')
        *s = 'U';
    }
}

Here is the caller graph for this function:

void fold2pairingList	(	char *	fold,
		int	len,
		int **	p2pairList
	)

Definition at line 47 of file rnautil.c.

References needMem().

            : pairList[i] = j <=> i pair with j and pairList[i] = -1
   <=> i does not pair.*/
{
int i,j, stackSize = 0;
int *pairList      = needMem(len * sizeof(int));
*p2pairList        = pairList;

/* initialize array */
for (i = 0; i < len; i++)
    pairList[i] = -1;

/* fill out pairList */
for (i = 0; i < len; i++) 
    {
    if (fold[i] == '(')
        {
        stackSize = 1;
        for (j = i+1; j < len; j++) 
            {
            if (fold[j] == '(')
                stackSize += 1;
            else if (fold[j] == ')')
                stackSize -= 1;
            if (stackSize == 0)  /* found pair partner */
                {
                pairList[i] = j;
                pairList[j] = i;
                break;
                }
            }
        }
    }
}

Here is the call graph for this function:

char* gapAdjustFold	(	char *	s,
		char *	ref
	)

Definition at line 134 of file rnautil.c.

References projectString().

{
return projectString(s, ref, '-', ' ');
}

Here is the call graph for this function:

int* gapIntArrayAdjust	(	int *	in,
		char *	ref
	)

Definition at line 160 of file rnautil.c.

References projectIntArray().

{
return projectIntArray(in, ref, '-', 0);
}

Here is the call graph for this function:

void markCompensatoryMutations	(	char *	s,
		char *	ref,
		int *	pairList,
		int *	markList
	)

Definition at line 166 of file rnautil.c.

References rnaPair().

       : 
 * 0: not pairing, no substitution (default) 
 * 1: not pairing, single substitution
 * 2: pairing, no substitutions 
 * 3: pairing, single substitution (one of: CG<->TG, GC<->GT, TA<->TG, AT<->GT)
 * 4: pairing, double substitution (i.e. a compensatory change)
 * 5: annoated as pairing but dinucleotide cannot pair 
 */
{
int i, size = strlen(s);
for (i = 0; i < size; i++)
    {
    if (pairList[i] == -1)
        if (toupper(s[i]) != toupper(ref[i]) && s[i] != '.' && s[i] != '-' && ref[i] != '-')
            markList[i] = 1;
        else
            markList[i] = 0;
    else
        {
        if (s[i] == '.' || s[pairList[i]] == '.') /* treat missing data as possible pair partner */
            markList[i] = 2;
        else if (!rnaPair(s[i], s[pairList[i]]))
            markList[i] = 5;
        else if (toupper( s[i] ) != toupper( ref[i] ) && toupper( s[pairList[i]] ) != toupper( ref[pairList[i]] ) )
            markList[i] = 4;
        else if (toupper( s[i] ) != toupper( ref[i] ) || toupper( s[pairList[i]] ) != toupper( ref[pairList[i]] ) )
            markList[i] = 3;
        else
            markList[i] = 2;
        }
    }
}

Here is the call graph for this function:

void mkPairPartnerSymbols	(	int *	pairList,
		char *	pairSymbols,
		int	size
	)

Definition at line 83 of file rnautil.c.

{
/* Make a symbol string indicating pairing partner */
int i;
char symbols[] = "abcdefghiklmnopqrstuvwxyzABCDEFGHIKLMNOPQRSTUVWXYZ1234567890!@#$%^&*=+{}|[]\\;'<>"; /* length 80 */
int  symbolMax = strlen(symbols);
int index;
for (i = 0, index = 0; i < size; i++)
    {
    pairSymbols[i] = ' ';
    if (pairList[i] >= 0)
        {
        if (pairList[i] < i)
            {
            --index;
            if (index<0)
                index = symbolMax-1;
            }
        pairSymbols[i] = symbols[index];
        if (pairList[i] > i)
            {
            ++index;
            if (index>=symbolMax)
                index=0;
            }
        }
    }
}

int* projectIntArray	(	int *	in,
		char *	ref,
		char	refChar,
		int	insertInt
	)

Definition at line 141 of file rnautil.c.

References needMem().

Referenced by gapIntArrayAdjust().

{
int i,j,size = strlen(ref);
int   *copy = (int *) needMem(size *sizeof(int) );

for (i = 0, j = 0; i < size; i++)
    {
    if (ref[i] == refChar)
        copy[i] = insertInt;
    else
        {       
        copy[i] = in[j];
        j++;
        }
    }
return copy;
}

Here is the call graph for this function:

Here is the caller graph for this function:

char* projectString	(	char *	s,
		char *	ref,
		char	refChar,
		char	insertChar
	)

Definition at line 112 of file rnautil.c.

References countChars(), errAbort(), and needMem().

Referenced by gapAdjustFold().

{
int i,j,size = strlen(ref);
char *copy = (char *) needMem(size + 1);

if (strlen(s) != strlen(ref) - countChars(ref, refChar))
  errAbort("ERROR from rnautil::projectString: Input string 's' has wrong length.\n"); 

for (i = 0, j = 0; i < size; i++)
    {
    if (ref[i] == refChar)
        copy[i] = insertChar;
    else
        {       
        copy[i] = s[j];
        j++;
        }
    }
return copy;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void reverseFold

(

char *

s

)

Definition at line 34 of file rnautil.c.

References reverseBytes().

{
reverseBytes(s, strlen(s));
for (;*s;s++)
    {
    if (*s == '(')
        *s = ')';
    else if (*s == ')')
        *s = '(';
    }
}

Here is the call graph for this function:

bool rnaPair	(	char	a,
		char	b
	)

Definition at line 20 of file rnautil.c.

References dna2rna(), FALSE, RNA_PAIRS, touppers(), and TRUE.

Referenced by markCompensatoryMutations().

{
char pair[] = {a,b,'\0'};
int i;
dna2rna(pair);
touppers(pair);

for (i=0;RNA_PAIRS[i] != 0; i++)
    if (pair[0] == RNA_PAIRS[i][0] && pair[1] == RNA_PAIRS[i][1] )
        return TRUE;
return FALSE;
}

Here is the call graph for this function:

Here is the caller graph for this function:

---

Variable Documentation

const char* RNA_PAIRS[] = {"AU","UA","GC","CG","GU","UG",0}

Definition at line 5 of file rnautil.c.

Referenced by rnaPair().

---