Untested guess -- are you reloading MASS in the new R session? I don't
believe loaded packages are saved in the session image... If MASS
isn't around, R won't know how to plot an lda object.
Best,
Michael
On Mon, Jun 11, 2012 at 4:55 PM, dga...@huskers.unl.edu
<dga...@huskers.unl.edu> wrote:
> Greetings R experts,
>
> I'm having some difficulty recovering lda objects that I've saved within
> sublists using the save.image() function. I am running a script that exports
> a variety of different information as a list, included within that list is an
> lda object. I then take that list and create a list of that with all the
> different replications I've run. Unfortunately I've been having difficulty
> accessing my lda data when I attempt to get it back after closing and
> reopening R. Here's an example that would give me the error I'm running into:
>
>> library(MASS)
>> ldaobject<-lda(Species~.,data=iris)
>> someotherresults<-c(1:5)
>> list1<-list(someotherresults,ldaobject)
>> list2<-list(list1,list1,list1)
>
>> plot(list2[[1]][[2]])#plots the ldaobject
>> save.image('ldalists.Rdata')
>
> ###Now if I close my R buffer and reopen it I get:
>
>> load('ldalists.Rdata')
>> plot(list2[[1]][[2]])
> Error in xy.coords(x, y, xlabel, ylabel, log) :
> 'x' is a list, but does not have components 'x' and 'y'
>
> ### And my lda obects appear to have changed and look like:
>> list2[[1]][[2]]
> $prior
> setosa versicolor virginica
> 0.3333333 0.3333333 0.3333333
>
> $counts
> setosa versicolor virginica
> 50 50 50
>
> $means
> Sepal.Length Sepal.Width Petal.Length Petal.Width
> setosa 5.006 3.428 1.462 0.246
> versicolor 5.936 2.770 4.260 1.326
> virginica 6.588 2.974 5.552 2.026
>
> $scaling
> LD1 LD2
> Sepal.Length 0.8293776 0.02410215
> Sepal.Width 1.5344731 2.16452123
> Petal.Length -2.2012117 -0.93192121
> Petal.Width -2.8104603 2.83918785
>
> $lev
> [1] "setosa" "versicolor" "virginica"
>
> $svd
> [1] 48.642644 4.579983
>
> $N
> [1] 150
>
> $call
> lda(formula = Species ~ ., data = iris)
>
> $terms
> Species ~ Sepal.Length + Sepal.Width + Petal.Length + Petal.Width
> attr(,"variables")
> list(Species, Sepal.Length, Sepal.Width, Petal.Length, Petal.Width)
> attr(,"factors")
> Sepal.Length Sepal.Width Petal.Length Petal.Width
> Species 0 0 0 0
> Sepal.Length 1 0 0 0
> Sepal.Width 0 1 0 0
> Petal.Length 0 0 1 0
> Petal.Width 0 0 0 1
> attr(,"term.labels")
> [1] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width"
> attr(,"order")
> [1] 1 1 1 1
> attr(,"intercept")
> [1] 1
> attr(,"response")
> [1] 1
> attr(,".Environment")
> <environment: R_GlobalEnv>
> attr(,"predvars")
> list(Species, Sepal.Length, Sepal.Width, Petal.Length, Petal.Width)
> attr(,"dataClasses")
> Species Sepal.Length Sepal.Width Petal.Length Petal.Width
> "factor" "numeric" "numeric" "numeric" "numeric"
>
> $xlevels
> named list()
>
> attr(,"class")
> [1] "lda"
>
>
> If anyone can help me out with a way to circumvent this problem or to recover
> my data that would be an immense help!
>
> Thanks again!
> -Dan
>
>
>
> [[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-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.
