[R] Question about package "principal"

David Winsemius
2015-01-11T19:59:24Z

On Jan 11, 2015, at 3:55 AM, ??? wrote:

> Hi,
>
> I tried to run "principal" function in the 'psych'
package, but it failed
> to do..
>
> here is both my code and error message.
>
> I searched on the web, but couldn't find the exact answer I wanted.
>
>> data<-read.csv("
>
https://raw.githubusercontent.com/mylesmharrison/delta_PCA_kmeans/master/delta.csv
> ",row.names=1)
>
>> airpca<-principal(data,nfactors=33,rotate="none")
> Error in solve.default(model, r) :
>  system is computationally singular: reciprocal condition number >
7.05776e-17
> In addition: Warning messages:
> 1: In cor.smooth(model) :
>  Matrix was not positive definite, smoothing was done
> 2: In cor.smooth(r) : Matrix was not positive definite, smoothing was done

Looks like your data is highly correlated in some of those columns:

> dat[[1]]
 [1]  0.0 19.4  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
[25]  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0  0.0
0.0  0.0  0.0  0.0  0.0
> dat[[29]]
 [1] 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0
> dat[[2]]
 [1]  0 44  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
[42]  0  0  0
> dat[[3]]
 [1]  0 12  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0
[42]  0  0  0

I'd retry after removing the problem columns.

>
> Thank you for your help.
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.

David Winsemius
Alameda, CA, USA