Entering edit mode
francesco.brundu@gmail.com
▴
40
@francescobrundugmailcom-5985
Last seen 6.7 years ago
Hi all,
I am trying to run ComBat on a dataset. I installed sva bioconductor
library and I am calling ComBat in this way:
ComBat(dat='dataset_13294/dChip_signal_13294.csv',batch='sif.tsv',
mod=NULL)
The output of the script is the following:
Found 1 batches
Found 1 categorical covariate(s)
Error in matrix(0, length(vec), nlevels(vec) - start + 1) :
invalid 'ncol' value (< 0)
I have searched a lot but I did not find any clue to resolve this
issue
and, since I am an R newbie, I am unable to inspect the code to know
what I
am doing wrong.
Here are the first two lines of dat and batch files passed as
parameters.
Thanks
dat:
probe set GSM335547 GSM335547 call GSM335558 GSM335558 call GSM335552
GSM335552
call GSM335597 GSM335597 call GSM335549 GSM335549 call GSM335608
GSM335608
call GSM335529 GSM335529 call GSM335661 GSM335661 call GSM335560
GSM335560
call GSM335655 GSM335655 call GSM335653 GSM335653 call GSM335629
GSM335629
call GSM335562 GSM335562 call GSM335524 GSM335524 call GSM335635
GSM335635
call GSM335663 GSM335663 call GSM335517 GSM335517 call GSM335639
GSM335639
call GSM335591 GSM335591 call GSM335563 GSM335563 call GSM335616
GSM335616
call GSM335652 GSM335652 call GSM335599 GSM335599 call GSM335581
GSM335581
call GSM335515 GSM335515 call GSM335630 GSM335630 call GSM335565
GSM335565
call GSM335659 GSM335659 call GSM335641 GSM335641 call GSM335569
GSM335569
call GSM335542 GSM335542 call GSM335574 GSM335574 call GSM335514
GSM335514
call GSM335516 GSM335516 call GSM335617 GSM335617 call GSM335525
GSM335525
call GSM335541 GSM335541 call GSM335527 GSM335527 call GSM335627
GSM335627
call GSM335657 GSM335657 call GSM335575 GSM335575 call GSM335618
GSM335618
call GSM335571 GSM335571 call GSM335650 GSM335650 call GSM335556
GSM335556
call GSM335590 GSM335590 call GSM335621 GSM335621 call GSM335662
GSM335662
call GSM335619 GSM335619 call GSM335534 GSM335534 call GSM335656
GSM335656
call GSM335593 GSM335593 call GSM335622 GSM335622 call GSM335628
GSM335628
call GSM335626 GSM335626 call GSM335606 GSM335606 call GSM335645
GSM335645
call GSM335568 GSM335568 call GSM335636 GSM335636 call GSM335559
GSM335559
call GSM335526 GSM335526 call GSM335640 GSM335640 call GSM335561
GSM335561
call GSM335522 GSM335522 call GSM335642 GSM335642 call GSM335611
GSM335611
call GSM335637 GSM335637 call GSM335658 GSM335658 call GSM335664
GSM335664
call GSM335532 GSM335532 call GSM335570 GSM335570 call GSM335647
GSM335647
call GSM335537 GSM335537 call GSM335555 GSM335555 call GSM335614
GSM335614
call GSM335588 GSM335588 call GSM335651 GSM335651 call GSM335521
GSM335521
call GSM335564 GSM335564 call GSM335530 GSM335530 call GSM335543
GSM335543
call GSM335572 GSM335572 call GSM335646 GSM335646 call GSM335598
GSM335598
call GSM335553 GSM335553 call GSM335551 GSM335551 call GSM335615
GSM335615
call GSM335578 GSM335578 call GSM335531 GSM335531 call GSM335654
GSM335654
call GSM335583 GSM335583 call GSM335595 GSM335595 call GSM335567
GSM335567
call GSM335536 GSM335536 call GSM335592 GSM335592 call GSM335585
GSM335585
call GSM335643 GSM335643 call GSM335631 GSM335631 call GSM335605
GSM335605
call GSM335576 GSM335576 call GSM335510 GSM335510 call GSM335582
GSM335582
call GSM335604 GSM335604 call GSM335573 GSM335573 call GSM335638
GSM335638
call GSM335612 GSM335612 call GSM335644 GSM335644 call GSM335603
GSM335603
call GSM335577 GSM335577 call GSM335540 GSM335540 call GSM335589
GSM335589
call GSM335535 GSM335535 call GSM335566 GSM335566 call GSM335607
GSM335607
call GSM335610 GSM335610 call GSM335609 GSM335609 call GSM335586
GSM335586
call GSM335648 GSM335648 call GSM335660 GSM335660 call GSM335554
GSM335554
call GSM335546 GSM335546 call GSM335545 GSM335545 call GSM335512
GSM335512
call GSM335579 GSM335579 call GSM335557 GSM335557 call GSM335633
GSM335633
call GSM335600 GSM335600 call GSM335620 GSM335620 call GSM335587
GSM335587
call GSM335519 GSM335519 call GSM335548 GSM335548 call GSM335584
GSM335584
call GSM335528 GSM335528 call GSM335625 GSM335625 call GSM335623
GSM335623
call GSM335596 GSM335596 call GSM335602 GSM335602 call GSM335513
GSM335513
call GSM335544 GSM335544 call GSM335613 GSM335613 call GSM335511
GSM335511
call GSM335550 GSM335550 call GSM335580 GSM335580 call GSM335649
GSM335649
call GSM335533 GSM335533 call GSM335624 GSM335624 call GSM335601
GSM335601
call GSM335594 GSM335594 call GSM335523 GSM335523 call GSM335538
GSM335538
call GSM335539 GSM335539 call GSM335632 GSM335632 call GSM335518
GSM335518
call GSM335520 GSM335520 call GSM335634 GSM335634 call
AFFX-BioB-5_at 9.17 P 8.96 P 9.09 P 8.76 P 8.74 P 8.89 P 6.76 P 9.31 P
8.99
P 7.86 P 7.98 P 8.24 P 7.84 P 6.78 P 8.12 P 8.58 P 7.12 P 7.92 P 8.33
P 8.94
P 8.45 P 7.83 P 8.70 P 9.94 P 6.92 P 8.54 P 8.11 P 9.41 P 8.08 P 8.96
P 7.20
P 7.98 P 6.98 P 6.72 P 8.53 P 7.08 P 7.44 P 7.23 P 8.06 P 9.01 P 9.06
P 8.54
P 8.97 P 8.66 P 8.77 P 8.88 P 7.60 P 7.65 P 8.36 P 7.07 P 8.13 P 8.80
P 8.34
P 8.16 P 9.11 P 8.96 P 7.71 P 8.43 P 8.35 P 8.55 P 6.97 P 8.06 P 8.77
P 8.68
P 8.38 P 9.67 P 9.54 P 8.36 P 8.80 P 6.75 P 9.23 P 9.98 P 7.26 P 8.46
P 8.80
P 8.58 P 8.19 P 7.38 P 8.96 P 6.85 P 7.49 P 9.64 P 8.98 P 8.33 P 9.24
P 8.95
P 8.42 P 8.88 P 7.59 P 8.77 P 8.50 P 9.44 P 8.94 P 6.77 P 8.89 P 9.07
P 8.15
P 7.48 P 8.31 P 8.89 P 6.94 P 9.33 P 8.94 P 8.96 P 9.82 P 8.81 P 8.11
P 8.89
P 8.84 P 7.20 P 9.05 P 7.42 P 8.95 P 8.50 P 9.27 P 8.17 P 8.69 P 10.19
P
8.13 P 8.59 P 7.20 P 7.65 P 6.95 P 8.78 P 8.93 P 8.79 P 8.64 P 8.93 P
8.58 P
7.04 P 8.74 P 8.94 P 8.38 P 8.85 P 8.85 P 8.62 P 8.62 P 7.25 P 7.44 P
8.54 P
7.08 P 8.86 P 9.03 P 8.21 P 7.05 P 8.72 P 8.40 P 8.63 P 7.16 P 7.08 P
7.08 P
7.64 P 6.83 P 8.15 P 7.98 P
batch:
Array name Sample name Batch
GSM335547 GSM335547 1
--
~Francesco Brundu
[[alternative HTML version deleted]]