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

2008-04-08 Thread Hazura Zulzalil
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@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.


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

2008-04-07 Thread Hazura Zulzalil
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@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.