[R] Fw: Help with PLSR

Bjørn-Helge Mevik b.h.mevik at usit.uio.no
Fri May 13 09:57:19 CEST 2011


Amit Patel <amitrhelp at yahoo.co.uk> writes:

>>str(FullDataListTrans)
>  num [1:40, 1:94727] 42 40.9 65 56 61.7 ...
>  - attr(*, "dimnames")=List of 2
>   ..$ : chr [1:40] "X" "X.1" "X.12" "X.13" ...
>   ..$ : NULL
>
> I have also created a vector "GroupingList" which gives the groupnames for each 
> respective sample(row).
>
>> GroupingList
>  [1] 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3 3 3 4 4 4 4 4 4 4 4
> [39] 4 4
>> str(GroupingList)
>  int [1:40] 1 1 1 1 1 1 1 1 1 1 ...
>
> I am now stuck while conducting the plsr. I have tried various methods of 
> creating structured lists etc  and have got nowhere. I have also tried many 
> incarnations of 
>
>
> BHPLS1 <- plsr(GroupingList ~ PCIList, ncomp = FeaturePresenceExpected[1], data 
> = FullDataListTrans, validation = "LOO")
>
> Where am I going wrong.

You are not telling us what happens (or how you tried to make
"structured lists"), but from your description of the
data, FullDataListTrans is a matrix with only the the predictor
variables, and GroupingList is a vector with the response.  The data
argument of plsr() (as of most modelling functions in R) expects a
data.frame with both response and predictor variables.

Try this:

mydata <- data.frame(GroupingList = GroupingList, PCIList = I(FullDataListTrans))

(The I() is to prevent R from making the coloumns in FullDataListTrans
separate variables in the data frame.)

BHPLS1 <- plsr(GroupingList ~ PCIList, ncomp = FeaturePresenceExpected[1], data = FullDataListTrans, validation = "LOO")

-- 
Regards,
Bjørn-Helge Mevik



More information about the R-help mailing list