[R] partial least regression

=?gb2312?Q?=D5=D4=D3=F1=D6=D2?= zhaoy at ysu.edu.cn
Wed May 16 15:47:15 CEST 2007


hello r-helpers:
there is a .txt file:
x1	x2	x3	x4	x5	x6	x7	x8	x9	x10	x11	y1
17	5	77	18	19	24	7	24	24	72	52	100
2	6	72	18	17	15	4	12	18	35	42	97.2
17	2	58	10	5	3	4	3	3	40	28	98
17	2	69	14	13	12	4	6	6	50	37	93
2	3	75	20	38	18	6	12	18	73	67	99
14	4	59	16	18	9	4	3	15	47	40	99.95
17	4	87	18	17	12	4	15	12	69	46	100
14	3	74	15	9	12	1	15	12	44	35	98
17	6	76	15	33	21	15	9	18	46	41	100
17	5	76	17	22	18	1	15	9	54	38	100
20	16	98	23	43	27	27	24	27	87	67	99.9
20	6	110	25	36	18	7	18	12	99	88	90
20	6	100	20	19	18	4	15	21	84	60	99.2
17	6	72	13	31	18	5	21	12	52	43	95.6
20	15	104	25	53	30	31	30	30	100	90	99.93
17	13	68	21	42	24	14	18	12	75	69	99.1
17	8	63	14	18	9	10	15	25	39	38	95
20	7	88	19	28	18	7	15	24	61	49	100
20	14	61	18	48	15	11	21	18	48	48	99.91
20	9	94	22	38	21	13	31	12	84	80	99.43
17	7	70	16	28	12	18	12	12	52	38	100
17	6	80	17	46	18	18	24	18	70	69	99.2 

and i write this program:
library(pls)
quality.data<-read.table("d:\\work\\12.txt",col.names=c("x1","x2","x3","x4",
"x5","x6","x7","x8","x9","x10","x11","y1"))
attach(quality.data)
quality.data              #
pairs(cbind(x4,x5,y1))
X<-data.frame(x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,y1)
Y<-data.frame(y1)
quality.simpls<-mvr(X,Y,1:11,method="simpls",validation="CV")  #####error
happens in this row ,and the error information is about: error in
:eval(expr, envir, enclos) :the length of number'envir' is not 1. 

i don't know how to solve this problem.
who can help me make this program correct,please. thank you very much!

                                      kenneth



More information about the R-help mailing list