compareKnownStructures: compareKnownStructures

In JLP-BioInf/comradesOO: Analysis of RNA Crosslinking Data

compareKnownStructures

Description

This method compares the predicted structures to a set of known interactions

Usage

compareKnownStructures(foldedCds, file)

Arguments

foldedCds	rnaCrosslinkDataSet after running foldrnaCrosslink
file	a two column file with column header i and j with numeric values showing nucleoide i binds to nucleotide j

Value

Returns a dataframe

a tables showing the number of predicted interactions and their agreement

Examples

## Not run: 
cds = makeExamplernaCrosslinkDataSet()
clusteredCds = clusterrnaCrosslink(cds = cds,
                               cores = 3,
                               stepCount = 2,
                               clusterCutoff = 1)


trimmedClusters = trimClusters(clusteredCds = clusteredCds,trimFactor = 1, clusterCutoff = 1)



fasta = paste(c(rep('A',25),
                rep('T',25),
                rep('A',10),
                rep('T',23)),collapse = "")

header = '>transcript1'


fastaFile = tempfile()
writeLines(paste(header,fasta,sep = "\n"),con = fastaFile)


rnaRefs = list()
rnaRefs[[rnas(cds)]] = read.fasta(fastaFile)
rnaRefs



foldedCds = foldrnaCrosslink(trimmedClusters,
                         rnaRefs = rnaRefs,
                         start = 1,
                         end = 83,
                         shape = 0,
                         ensembl = 5,
                         constraintNumber  = 1,
                         evCutoff = 1)


# make an example table of "know" interactions
file = data.frame(V1 = c(6), 
                  V2 = c(80))
compareKnownStructures(foldedCds,file)

## End(Not run)

JLP-BioInf/comradesOO documentation built on April 28, 2024, 4:22 a.m.