Hi Julie,
I am using offTargetAnalysis to find gRNAs with CFDscores. I followed the manual and included the subPAM.activity option as follows:
results=offTargetAnalysis("exon.fasta", findgRNAsWithREcutOnly = FALSE, annotateExon = TRUE,findPairedgRNAOnly = FALSE, txdb= txdb, max.mismatch = 3, BSgenomeName = ss, outputDir = ".", overwrite = TRUE,scoring.method="CFDscore", subPAM.activity = hash(AA =0,AC = 0,AG = 0.259259259,AT = 0,CA = 0,CC = 0,CG = 0.107142857,CT = 0,GA = 0.069444444,GC = 0.022222222,GG = 1,GT = 0.016129032,TA = 0,TC = 0,TG = 0.038961039,TT = 0),subPAM.position = c(22, 23),PAM.location = "3prime",mismatch.activity.file = system.file("extdata","NatureBiot2016SuppTable19DoenchRoot.csv",package = "CRISPRseek"))
I encountered the following error message at the score-calculation step:
Error in hash(AA = 0, AC = 0, AG = 0.259259259, AT = 0, CA = 0, CC = 0, :
could not find function "hash"
Could you please tell me how to fix this?
Thanks!