Description Usage Arguments Details Author(s) See Also Examples
Input data for this function should be assessed from any of the CG_KaCh() plots.
See argument details.
Very important: For this function to work you must have the file CG.Rdata, generated using patchwork.CG.plot(), in your working directory!
1 2 | patchwork.CG.copynumbers(cn2,delta,het,hom,maxCn=8,ceiling=1,name="copynumbers_",
CGfile=NULL,forcedelta=F)
|
cn2 |
The approximate position of copy number 2, diploid, on total intensity axis. |
delta |
The difference in total intensity between consecutive copy numbers. For example 1 and 2 or 2 and 3. If copy number 2 has total intensity ~0.6 and copy number 3 har total intensity ~0.8 then delta would be 0.2. |
het |
Allelic imbalance ratio of heterozygous copy number 2. |
hom |
Allelic imbalance ratio of Loss-of-Heterozygosity copy number 2. |
maxCn |
Highest copy number to calculate for. Default is 8. |
ceiling |
Default is 1. |
name |
Default is "copynumbers_". The name you want attached to generated plots. |
CGfile |
Default is NULL. If your CG.Rdata file is not in your working directory, and you dont wish to move it to your working directory, you can simply input the path here as CGfile = "path/to/file/CG.Rdata" so patchwork.CG.copynumbers() can find its data. |
forcedelta |
Default is FALSE. If TRUE the delta value will be absolute and not subject to adjustments. |
1. Run the patchwork.CG.plot correctly. This will generate CG.Rdata.
2. View your plots and decide the correct parameters for the other arguments of this function.
3. Make sure CG.Rdata is in your working directory, or set the CGfile parameter to point to CG.Rdata.
4. Run the function.
Upon completion two files will be generated, copynumbers.Rdata which contains the objects used to generate the plots and Copynumbers.csv. Other than this plots from CG_Ka_check and CG_KaChCN will be generated.
Markus Mayrhofer, markus.mayrhofer@medsci.uu.se
Sebastian DiLorenzo, sebastian.dilorenzo@medsci.uu.se
patchwork.CG.plot
subfunctions_copynumbers
CG_Ka_check
CG_KaChCN
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 | ## Not run:
>library(patchworkCG)
#With CG.Rdata in working directory. Values for parameters estimated from
#HCC2218 chromosome 1.
#See plot patchworkCG/example_plots/HCC2218_KaCh_chr1.png
>patchwork.CG.copynumbers(cn2=0.85,delta=0.35,het=0.25,hom=0.95,name="HCC2218_")
#Pointing to location of CG.Rdata file
>patchwork.CG.copynumbers(cn2=0.85,delta=0.35,het=0.25,hom=0.95,name="HCC2218_",
CGfile="path/to/CG.Rdata")
#Load the objects generated from patchwork.CG.copynumbers
#For comparison purposes only
>load("copynumbers.Rdata")
> head(ai)
$cn2m1
[1] 0.2460513
$cn2m0
[1] 0.95
$cn1m0
[1] 0.9690196
$cn0m0
[1] NA
$cn3m1
[1] 0.2677295
$cn3m0
[1] 0.9671712
> head(int)
$cn2
[1] 0.85
$cn1
[1] 0.5
$cn0
[1] 0.15
$cn3
[1] 1.2
$cn4
[1] 1.55
$cn5
[1] 1.9
> head(regions)
chr start end avgnormcov relcov ratio snvs ai np Cn mCn
1 chr1 10000 177417 79.0 1.05 0.9660335 24 0.5531915 17 2 1
2 chr1 227417 267719 79.0 1.05 0.9660335 2 0.3243243 4 2 1
3 chr1 317719 471368 79.0 1.05 0.9660335 2 0.7826087 15 2 0
4 chr1 521368 600000 170.8 2.26 2.0792721 2 0.1944444 8 6 3
5 chr1 600000 700000 74.6 0.99 0.9108316 2 0.5333333 10 2 1
6 chr1 700000 900000 51.9 0.69 0.6348220 109 0.4960325 20 1 0
fullCN
1 cn2m1
2 cn2m1
3 cn2m0
4 cn6m3
5 cn2m1
6 cn1m0
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.