Windows
sessionInfo() #R version 3.5.1 (2018-07-02)
Hi I am following along (using my own data) the processes described by
"glmulti: An R Package for Easy Automated Model Selection with (
Generalized ) Linear Models"
Journal of Statistical Software
May 2010, Volume 34, Issue 12. http://www.jstatsoft.org/
I will send the pdf if requested, however, I doubt it will get though attached
here.
I am having trouble understanding why I am getting these warnings:
ga1 <- glmulti(y ~ a+b+c+d, fitfunc = glm, crit = aic, method =
"h")
plot(ga1, type = "p", highlight = c("d:a","d:b"))
Warning messages:
1: In plot.window(...) : "highlight" is not a graphical parameter
2: In plot.xy(xy, type, ...) : "highlight" is not a graphical
parameter
3: In axis(side = side, at = at, labels = labels, ...) :
"highlight.3d" is not a graphical parameter
4: In axis(side = side, at = at, labels = labels, ...) :
"highlight.3d" is not a graphical parameter
5: In box(...) : "highlight" is not a graphical parameter
6: In title(...) : "highlight" is not a graphical parameter
The authors describe the plots on pg. 15
Also I am not getting the highlights in my plots as is shown in the author's
document?
I suppose I could reach out to the authors, assuming 8 years later the e-mail
address is the same, however before I did that I thought maybe someone out here
might know why this is happening for me.
Here is my structure
dim(r1)
[1] 23189 6
str(r1)
Classes 'data.table' and 'data.frame': 23189 obs. of 6
variables:
$ SavingsReversed: num 0 0 0 0 0 0 0 0 0 0 ...
$ productID : int 3 3 3 3 3 3 3 3 1 1 ...
$ ProviderID : int 113676 113676 113964 113964 114278 114278 114278 114278
114278 114278 ...
$ Editnumber : int 502 502 504 504 504 504 504 504 504 504 ...
$ PatientGender : Factor w/ 3 levels "F","M","U":
1 1 2 2 2 2 1 1 1 1 ...
$ Editnumber2 : num 0 0 1 1 1 1 1 1 1 1 ...
- attr(*, ".internal.selfref")=<externalptr
SavingsReversed productID ProviderID Editnumber PatientGender
Editnumber2
1: 0 3 113676 502 F 0
2: 0 3 113676 502 F 0
3: 0 3 113964 504 M 1
4: 0 3 113964 504 M 1
5: 0 3 114278 504 M 1
6: 0 3 114278 504 M 1
7: 0 3 114278 504 F 1
8: 0 3 114278 504 F 1
9: 0 1 114278 504 F 1
10: 0 1 114278 504 F 1
>
y <- r1$Editnumber2
a <- r1$productID
b <- r1$ProviderID
c <- r1$SavingsReversed
d <- r1$PatientGender
ga1 <- glmulti(y ~ a+b+c+d, fitfunc = glm, crit = aic, method =
"h")
I have looked at
https://cran.r-project.org/web/packages/scatterplot3d/scatterplot3d.pdf for more
info on highlight and it refers to highlight.3d which I also tried, no
difference.
Thank you for any insight.
Confidentiality Notice This message is sent from Zelis. ...{{dropped:15}}