Any thoughts on this error in robustSVD?
Thanks a lot!
Error in if (!all(tmp)) { : missing value where TRUE/FALSE needed
Enter a frame number, or 0 to exit
1: #73: pca(dTmp, method = "robustPca", nPcs = nNumFactors, center = FALSE)
2: robustPca(prepres$data, nPcs = nPcs, ...)
3: robustSvd(Matrix)
4: apply(x, 1, L1RegCoef, bk)
5: FUN(newX[, i], ...)
6: weightedMedian(x[keep]/a, abs(a), interpolate = FALSE)
7: weightedMedian.default(x[keep]/a, abs(a), interpolate = FALSE)
On Sun, Apr 22, 2012 at 6:43 PM, Michael <[email protected]> wrote:
> I actually tried "robustPca" in "pcaMethods" on bioconductor.
>
> It keeps giving me the warning "Input data is not complete"...
>
> Reading into the function:
>
> When there is no "NA"s, it will give this warning...
>
> It seems that there is a bug in this code...
>
> Is it reliable at all?
>
> ---------------------
>
>
> > robustPcafunction (Matrix, nPcs = 2, verbose = interactive(), ...)
> {
> nas <- is.na(Matrix)
> if (!any(nas) & verbose) {
> cat("Input data is not complete.\n")
> cat("Scores, R2 and R2cum may be inaccurate, handle with care\n")
> }
>
>
>
>
>
> On Fri, Apr 20, 2012 at 9:58 AM, Kevin Wright <[email protected]> wrote:
>
>> You can also have a look at the pcaMethods package on Bioconductor.
>>
>> Kevin
>>
>>
>> On Thu, Apr 19, 2012 at 11:20 PM, Michael <[email protected]> wrote:
>>
>>> Hi all,
>>>
>>> I found that the PCA gave chaotic results when there are big changes in a
>>> few data points.
>>>
>>> Are there "improved" versions of PCA in R that can help with this
>>> problem?
>>>
>>> Please give me some pointers...
>>>
>>> Thank you!
>>>
>>> [[alternative HTML version deleted]]
>>>
>>> ______________________________________________
>>> [email protected] mailing list
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide
>>> http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html>
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>
>>
>>
>> --
>> Kevin Wright
>>
>>
>
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
