-
Notifications
You must be signed in to change notification settings - Fork 15
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
unable to plot using ggplot2, object 'count' not found #44
Comments
Hi, Did you do a nearBait analysis? - nb_results=nearBaitAnalysis(my_obj,k=5) |
yes! i'm sorry. i did, nb_results=nearBaitAnalysis(my_obj,k=10). I get
that error.
… Hi,
Did you do a nearBait analysis? - nb_results=nearBaitAnalysis(my_obj,k=5)
--
You are receiving this because you authored the thread.
Reply to this email directly or view it on GitHub:
#44 (comment)
|
maybe it doesnt work with no replicates? is that the case?? |
Yes sorry the plotting doesnt work with one replicate - i would recommend visualizing the bedGraph file on a browser. |
Hi, In that case, what if i give the same file copy as a replicate? the results wouldn't be different from there being just one replicate would it? or will i give incorrect results? |
Hello,
I have successfully loaded the 4Cker object as:
my_obj = createR4CkerObjectFromFiles(files = "hs75l6_13_DN-4C_9p21_DD1_CTCF_view_GACCAG_L006_R1_001_aligned.bedGraph",bait_chr="chr9",bait_coord= 22231303,bait_name = "CTCF-9p21locus", primary_enz = "AAGCTT", samples = c("CTCF_DD1"), conditions = "DD1", replicates = 1,species = "hg",output_dir = "/home/dnotani/ssanaa/software/R.4Cker-master/folder_for_analysis",enz_file=enz_file)
however, upon plotting the object i am getting the following error:
Don't know how to automatically pick scale for object of type Coord/ggproto. Defaulting to continuous.
Error in eval(expr, envir, enclos) : object 'Count' not found
the command for plotting used is the default command given in readme:
ggplot(nb_results$norm_counts_avg, aes(x=Coord, y=Count, colour=Condition))+
theme_bw()+
geom_line()+xlab(paste("Chromosome coordinates (", my_obj@bait_chr, ")", sep =""))+
ylab("Normalized counts")+
ggtitle(paste("Near bait analysis (", my_obj@bait_name, " bait)", sep = ""))
Please help!
Thank you.
The text was updated successfully, but these errors were encountered: