glmulti and plot highlight=

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

glmulti and plot highlight=

Bill Poling
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}}

______________________________________________
[hidden email] mailing list -- To UNSUBSCRIBE and more, see
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.