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.