Skip to content
Snippets Groups Projects
Commit 56a04592 authored by Todor Kondić's avatar Todor Kondić
Browse files

Clean up some comments in `make_spec_ms2_plot'.

parent 4c3c6754
No related branches found
No related tags found
No related merge requests found
...@@ -457,25 +457,18 @@ make_spec_ms2_plot <- function(db,summ,kvals,labs,axis="linear",asp=1, colrdata= ...@@ -457,25 +457,18 @@ make_spec_ms2_plot <- function(db,summ,kvals,labs,axis="linear",asp=1, colrdata=
outcols=union(names(kvals), outcols=union(names(kvals),
colnames(summ)))[ms2_sel==T] colnames(summ)))[ms2_sel==T]
spectra = db$extr$spectra[db$extr$cgm$ms2[,unique(.SD), spectra = db$extr$spectra[mdata,on=.(precid,scan),.(catid,ce,mz,scan,rt=i.ms2_rt,intensity)]
.SDcols=c("ce","scan","rt")],
on=.(scan),c("ce","rt"):=.(i.ce,i.rt)]
subxdata = get_data_from_key(db=db,
tab=spectra,
kvals=kvals,
outcols=c("catid","ce","mz","scan","rt","intensity"))
if (NROW(mdata)==0L) return(NULL) if (NROW(mdata)==0L) return(NULL)
if (NROW(subxdata) == 0L) return(NULL) if (NROW(spectra) == 0L) return(NULL)
ans = data.table(scan=mdata[,unique(scan)],key="scan")
meta = db$cat[subxdata[,.(catid=unique(catid))], meta = db$cat[spectra[,.(catid=unique(catid))],
on=.(catid), on=.(catid),
.SD, .SD,
by=.EACHI, by=.EACHI,
.SDcols=labs] .SDcols=labs]
subxdata = meta[subxdata,on=.(catid),nomatch=NULL] subxdata = meta[spectra,on=.(catid),nomatch=NULL]
...@@ -485,7 +478,7 @@ make_spec_ms2_plot <- function(db,summ,kvals,labs,axis="linear",asp=1, colrdata= ...@@ -485,7 +478,7 @@ make_spec_ms2_plot <- function(db,summ,kvals,labs,axis="linear",asp=1, colrdata=
if (NROW(pdata)==0L) return(NULL) if (NROW(pdata)==0L) return(NULL)
# Aspect ratio. # Aspect ratio.
xrng = range(pdata$mz) xrng = range(pdata$mz,na.rm=T)
dx = abs(xrng[[2]]-xrng[[1]]) dx = abs(xrng[[2]]-xrng[[1]])
yrng = range(pdata$intensity) yrng = range(pdata$intensity)
dy = abs(yrng[[2]]-yrng[[1]]) dy = abs(yrng[[2]]-yrng[[1]])
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment