Re: [R] partial least regression

From: Gabor Grothendieck <ggrothendieck_at_gmail.com>
Date: Wed, 16 May 2007 10:26:54 -0400

Your call to mvr does not correspond to the documentation. See ?mvr and try this:

Lines <- "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
"
library(pls)
quality.data <- read.table(textConnection(Lines), header = TRUE) pairs(quality.data[c("x4", "x5", "y1")]) quality.simpls <- mvr(y1 ~., ncomp = 6, quality.data, method = "simpls",   validation = "CV") # replace 6 with whatever you want

On 5/16/07, <zhaoy_at_ysu.edu.cn> wrote:
> 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
>
>
> ______________________________________________
> R-help_at_stat.math.ethz.ch 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-help_at_stat.math.ethz.ch 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 Wed 16 May 2007 - 14:45:13 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 Wed 16 May 2007 - 15:31:26 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.