Hi,
I am trying to perform cross platform normalization for CytoScan and Omni
using MSCN but the program stops after smoothing. I imported the Affymetrix
raw intensity data and called the totalCn and Baf using CRMAv2 and CalMate.
LogR ratio (log2 scale) and BAF were extracted using GenomeStudio for Omni
and imported as 2^LRR*2 and BAF into .total.asb and .fracB.asb
respectively. Plotting the CN indicated similar level of mean intensities
for both platforms (CN~1) and thereby I proceeded to perform the mscn.
Smoothing looks normal but after that, the program seemed to stop while
backtransforming.
May I know how I can overcome the error of * smooth.spline(sKK, lambda,
...)*?? What's wrong in my data?
Thanks a lot!!!!!!
--Clara
*P.S. Log file is pasted as below,*
=====================================================
*dsNList <- process(mscn, verbose=log)*
20140427 17:06:46|Multi-source normalize all samples...
20140427 17:06:46| Number of unique samples in all sets: 8
chr [1:8] "9C" "9C-2" "2C" "2C-2" "4C" ...
20140427 17:06:46| Processing each array...
20140427 17:06:46| Sample #1 ('9C') of 8...
20140427 17:06:46| Identifying source data files...
20140427 17:06:46| Getting list tuple of data files for one sample...
20140427 17:06:46| Sample name: 9C
20140427 17:06:47| Number of arrays: 2
20140427 17:06:47| Getting list tuple of data files for one sample...done
[[1]]
AromaUnitTotalCnBinaryFile:
Name: 9C
Tags: total
Full name: 9C,total
Pathname:
totalAndFracBData/Cytoscan,ACC,ra,-XY,BPN,-XY,AVG,FLN,-XY,CMTN,v2/CytoScanHD_Array/9C,total.asb
File size: 10.77 MB (11288937 bytes)
RAM: 0.01 MB
Number of data rows: 2822125
File format: v1
Dimensions: 2822125x1
Column classes: double
Number of bytes per column: 4
Footer: <createdOn>20140421 18:35:53
HKT</createdOn><platform>Affymetrix</platform><chipType>CytoScanHD_Array</chipType><srcFile><srcDataSet>Cytoscan,ACC,ra,-XY,BPN,-XY,AVG,FLN,-XY</srcDataSet><srcChipType>CytoScanHD_Array,monocell</srcChipType><srcFullName>9C,chipEffects</srcFullName><srcChecksum>7d70e31c6aa60b7d8664373333ad5452</srcChecksum></srcFile>
Platform: Affymetrix
Chip type: CytoScanHD_Array
[[2]]
AromaUnitTotalCnBinaryFile:
Name: 9C
Tags:9C,ratio,total
Full name: 9C,9C,ratio,total
Pathname: totalAndFracBData/Omni/HumanOmni2.5/9C,9C,ratio,total.asb
File size: 9.08 MB (9519600 bytes)
RAM: 0.01 MB
Number of data rows: 2379855
File format: v1
Dimensions: 2379855x1
Column classes: double
Number of bytes per column: 4
Footer: <createdOn>20140427 16:14:31
HKT</createdOn><platform>Illumina</platform><chipType>HumanOmni2.5</chipType>
Platform: Illumina
Chip type: HumanOmni2.5
20140427 17:06:47| Identifying source data files...done
20140427 17:06:47| Check if all arrays are already normalized...
20140427 17:06:47| Is done: FALSE
20140427 17:06:47| Check if all arrays are already normalized...done
20140427 17:06:47| Fitting model...
20140427 17:06:47| Fitting one sample across multiple sources...
20140427 17:06:47| Number of arrays: 2
20140427 17:06:47| Sample name: 9C
List of 4
$ subsetToFit : int [1:28781] 1 2 3 4 5 6 7 8 9 10 ...
$ fitUgp :Classes 'AromaUgpFile',
'AromaUnitChromosomeTabularBinaryFile', 'AromaUnitTabularBinaryFile',
'UnitAnnotationDataFile', 'AromaMicroarrayTabularBinaryFile',
'AromaPlatformInterface', 'AromaTabularBinaryFile',
'FileCacheKeyInterface', 'CacheKeyInterface', 'GenericTabularFile',
'ColumnNamesInterface', 'GenericDataFile', 'FullNameInterface', 'Object'
atomic [1:1] NA
.. ..- attr(*, ".env")=<environment: 0x3c8f518>
$ align : chr "byChromosome"
$ targetDimension: int 1
20140427 17:06:47| Getting list tuple of data files for one sample...
20140427 17:06:47| Sample name: 9C
20140427 17:06:47| Number of arrays: 2
20140427 17:06:47| Getting list tuple of data files for one
sample...done
List of 2
$ :Classes 'AromaUnitTotalCnBinaryFile', 'CopyNumberDataFile',
'AromaUnitSignalBinaryFile', 'AromaPlatformInterface',
'AromaTabularBinaryFile', 'FileCacheKeyInterface', 'CacheKeyInterface',
'GenericTabularFile', 'ColumnNamesInterface', 'GenericDataFile',
'FullNameInterface', 'Object' atomic [1:1] NA
.. ..- attr(*, ".env")=<environment: 0x67be7f0>
$ :Classes 'AromaUnitTotalCnBinaryFile', 'CopyNumberDataFile',
'AromaUnitSignalBinaryFile', 'AromaPlatformInterface',
'AromaTabularBinaryFile', 'FileCacheKeyInterface', 'CacheKeyInterface',
'GenericTabularFile', 'ColumnNamesInterface', 'GenericDataFile',
'FullNameInterface', 'Object' atomic [1:1] NA
.. ..- attr(*, ".env")=<environment: 0x36d2f48>
20140427 17:06:47| Extracting data...
20140427 17:06:47| Subset of units used for fitting:
int [1:28781] 1 2 3 4 5 6 7 8 9 10 ...
num [1:28781, 1:2] 0.561 0.566 0.798 0.861 0.81 ...
V1 V2
Min. :-0.7285 Min. :-1.9090
1st Qu.: 0.9677 1st Qu.: 0.9403
Median : 0.9842 Median : 0.9838
Mean : 0.9805 Mean : 0.9898
3rd Qu.: 0.9995 3rd Qu.: 1.0340
Max. : 1.3343 Max. : 1.9075
NA's :1865 NA's :1934
20140427 17:06:48| Extracting data...done
20140427 17:06:48| Fitting across-source normalization function...
20140427 17:06:48| Estimator for principal curves:
function (Y, ...)
- attr(*, "smoother")=function (lambda, xj, ...)
- attr(*, "df")= num 5
20140427 17:07:52| Fitting time:
user system elapsed
63.884 0.008 64.062
20140427 17:07:52| Direction of fitted curve ('lambda') was flipped
such that it increases with the signal.
20140427 17:07:52| Processing time: 64.1 seconds
List of 7
$ s : num [1:28781, 1:2] 0.652 0.652 0.665 NA 0.836 ...
$ tag : num [1:28781] 15371 15372 7262 NA 15373 ...
$ lambda : num [1:28781] 1.99 1.99 1.93 NA 1.38 ...
$ dist : num 31.6
$ converged : logi TRUE
$ nbrOfIterations: int 6
$ call : language principal.curve(x = X, smoother = ..1,
trace = trace)
- attr(*, "class")= chr [1:2] "PrincipalCurve" "principal.curve"
- attr(*, "processingTime")=Class 'proc_time' Named num [1:5] 63.8
0.004 63.975 0 0
.. ..- attr(*, "names")= chr [1:5] "user.self" "sys.self" "elapsed"
"user.child" ...
20140427 17:07:52| Fitting across-source normalization function...done
20140427 17:07:52| Calculating shift for each chromosome...
20140427 17:07:52| align=byChromosome
20140427 17:07:52| Chromosomes: 1-25
20140427 17:07:52| Grouping units by chromosome...
List of 24
$ Chr01: int [1:2472] 1 2 3 4 5 6 7 8 9 10 ...
$ Chr02: int [1:2434] 2473 2474 2475 2476 2477 2478 2479 2480 2481
2482 ...
$ Chr03: int [1:1995] 4907 4908 4909 4910 4911 4912 4913 4914 4915
4916 ...
$ Chr04: int [1:1917] 6902 6903 6904 6905 6906 6907 6908 6909 6910
6911 ...
$ Chr05: int [1:1810] 8819 8820 8821 8822 8823 8824 8825 8826 8827
8828 ...
$ Chr06: int [1:1709] 10629 10630 10631 10632 10633 10634 10635
10636 10637 10638 ...
$ Chr07: int [1:1589] 12338 12339 12340 12341 12342 12343 12344
12345 12346 12347 ...
$ Chr08: int [1:1463] 13927 13928 13929 13930 13931 13932 13933
13934 13935 13936 ...
$ Chr09: int [1:1405] 15390 15391 15392 15393 15394 15395 15396
15397 15398 15399 ...
$ Chr10: int [1:1355] 16795 16796 16797 16798 16799 16800 16801
16802 16803 16804 ...
$ Chr11: int [1:1350] 18150 18151 18152 18153 18154 18155 18156
18157 18158 18159 ...
$ Chr12: int [1:1335] 19500 19501 19502 19503 19504 19505 19506
19507 19508 19509 ...
$ Chr13: int [1:1142] 20835 20836 20837 20838 20839 20840 20841
20842 20843 20844 ...
$ Chr14: int [1:1064] 21977 21978 21979 21980 21981 21982 21983
21984 21985 21986 ...
$ Chr15: int [1:1004] 23041 23042 23043 23044 23045 23046 23047
23048 23049 23050 ...
$ Chr16: int [1:900] 24045 24046 24047 24048 24049 24050 24051
24052 24053 24054 ...
$ Chr17: int [1:817] 24945 24946 24947 24948 24949 24950 24951
24952 24953 24954 ...
$ Chr18: int [1:778] 25762 25763 25764 25765 25766 25767 25768
25769 25770 25771 ...
$ Chr19: int [1:639] 26540 26541 26542 26543 26544 26545 26546
26547 26548 26549 ...
$ Chr20: int [1:637] 27179 27180 27181 27182 27183 27184 27185
27186 27187 27188 ...
$ Chr21: int [1:470] 27816 27817 27818 27819 27820 27821 27822
27823 27824 27825 ...
$ Chr22: int [1:496] 28286 28287 28288 28289 28290 28291 28292
28293 28294 28295 ...
$ Chr23: int [1:1550] 28782 28783 28784 28785 28786 28787 28788
28789 28790 28791 ...
$ Chr24: int [1:578] 30332 30333 30334 30335 30336 30337 30338
30339 30340 30341 ...
20140427 17:07:52| Grouping units by chromosome...done
20140427 17:07:52| Allocating matrix for smooth data...
20140427 17:07:52| RAM: 494760 bytes
20140427 17:07:52| Allocating matrix for smooth data...done
20140427 17:07:52| Loading and backtransforming *smoothed* data...
20140427 17:07:52| Source #1 ('9C,total') of 2...
20140427 17:07:52| Loading smoothed data...
num [1:30910] 0.561 0.566 0.798 0.861 0.81 ...
20140427 17:07:52| Loading smoothed data...done
20140427 17:07:52| Backtransforming smoothed data...
*Error in smooth.spline(sKK, lambda, ...) : *
missing or infinite values in inputs are not allowed
20140427 17:07:52| Backtransforming smoothed data...done
20140427 17:07:52| Source #1 ('9C,total') of 2...done
20140427 17:07:52| Loading and backtransforming *smoothed* data...done
20140427 17:07:52| Calculating shift for each chromosome...done
20140427 17:07:52| Fitting one sample across multiple sources...done
20140427 17:07:52| Fitting model...done
20140427 17:07:52| Sample #1 ('9C') of 8...done
20140427 17:07:52| Processing each array...done
20140427 17:07:52|Multi-source normalize all samples...done
--
--
When reporting problems on aroma.affymetrix, make sure 1) to run the latest
version of the package, 2) to report the output of sessionInfo() and
traceback(), and 3) to post a complete code example.
You received this message because you are subscribed to the Google Groups
"aroma.affymetrix" group with website http://www.aroma-project.org/.
To post to this group, send email to aroma-affymetrix@googlegroups.com
To unsubscribe and other options, go to http://www.aroma-project.org/forum/
---
You received this message because you are subscribed to the Google Groups
"aroma.affymetrix" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to aroma-affymetrix+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
