Hello,
I tried on test data for both cnv.gr.rda and snap.gr.rda using
load(system.file("data", "cnv.gr.rda", package="BubbleTree"))
load(system.file("data", "snp.gr.rda", package="BubbleTree"))
It has the error message:
Error: bad restore file magic number (file may be corrupted) -- no data loaded
In addition: Warning message:
file ‘cnv.gr.rda’ has magic number '<'
Use of save versions prior to 2 is deprecated
I also tried my own dataset
vc.df = read.csv("F5.snp..inputBT.csv", header=T, sep=",")
vc.gr = GRanges(vc.df$chrom, IRanges(vc.df$start, vc.df$end),
freq=vc.df$freq, score=vc.df$freq)
cnv.df = read.csv("F5.cnv.inputBT.csv", header=T, sep=",")
cnv.gr = GRanges(cnv.df$chrom, IRanges(cnv.df$start, cnv.df$end),
num.mark=cnv.df$num.mark, seg.mean=cnv.df$seg.mean,
score=cnv.df$seg.mean)
When I initiated the object and create RBD, it met the errors again.
r <- new("RBD")
> new("RBD")
RBD object with 1 range and 0 metadata columns:
Error in `rownames<-`(`*tmp*`, value = "[1]") :
length of 'dimnames' [1] not equal to array extent
> rbd <- makeRBD(r, vc.gr, cnv.gr)
Error: could not find function "makeRBD"
Thanks for the help.
EL