[R] error using method ls.ranking.capa.ident

From: Hazura Zulzalil <hazura.zulzalil_at_gmail.com>
Date: Tue, 08 Apr 2008 15:12:00 +0800

Good Day

I'm trying to run the generalized least square approach for my 2-additive problem,
unfortunately this error appeared. I have tried to figure out the error from the mailing list
but couldn't find the solution. Any help is highly appreciated. This is my source code:

 >a1 <- c(76.18, 61.84, 60.4, 69.09)

> a2 <- c(51.01, 50.39, 87.62, 52.03)
> a3 <- c(80.08, 48.49, 90.86, 76.11)
> a4 <- c(57.71, 38.99, 88.7, 53.12)
> a5 <- c(71.93, 82.04, 83.05, 85.99)
> a6 <- c(60.94, 71.12, 63.61, 69.47)
>
> ##preference threshold relative
> ## to the preorder of the alternatives
>
> C<- rbind (a1, a2, a3, a4, a5, a6)
> overall <- c (79, 78, 79.5, 77.5, 80, 78.5)
>
> rk.proto <- rbind (c(5,3), c(3,1), c(1,6), c(6,4), c(4,2))
>
>
> Asp <-rbind (c(1,2,-delta.S), c(2,1,-delta.S), c(3, 4,-delta.S), c(4, 3,
-delta.S))
>
> ## an interaction preorder constraint matrix
> delta.I <- 0.05
> Aii <- rbind (c(1,2,-1,-delta.I), c(1,3,delta.I, 1), c (1,4,-1,-delta.I),
+ c(2,3, delta.I,1), c(2, 4,-1,-delta.I), c(3,4,delta.I, 1))
>
>
> gls <-ls.ranking.capa.ident(4,2,C, rk.proto,0.5, A.Shapley.preorder = Asp,
A.interaction.interval = Aii)
Error in drop(.Call("La_dgesv", a, as.matrix(b), tol, PACKAGE = "base")) :   system is computationally singular: reciprocal condition number = 4.68861e-18

        [[alternative HTML version deleted]]



R-help_at_r-project.org 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. Received on Tue 08 Apr 2008 - 07:14:38 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 08 Apr 2008 - 07:30:27 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive