> des.nohome <- design(buli0506,var1="Heim", var2="Gast",
+ use.vars="Y3",home.advantage="no",
+ reference="GAMMA.MSV.Duisburg")
> head(des.nohome )
GAMMA.1.FC.Kaiserslautern GAMMA.1.FC.Koeln GAMMA.1.FC.Nuernberg
12795 0 0 0
12796 0 1 0
12797 0 0 0
12798 0 0 -1
12799 0 0 0
12800 0 0 0
GAMMA.1.FSV.Mainz.05 GAMMA.Arminia.Bielefeld GAMMA.Bayer.Leverkusen
12795 0 0 0
12796 -1 0 0
12797 0 0 0
12798 0 0 0
12799 0 0 0
12800 0 -1 0
GAMMA.Borussia.Dortmund GAMMA.Borussia.Moenchengladbach
12795 0 -1
12796 0 0
12797 0 0
12798 0 0
12799 -1 0
12800 0 0
GAMMA.Eintracht.Frankfurt GAMMA.FC.Bayern.Muenchen GAMMA.FC.Schalke.04
12795 0 1 0
12796 0 0 0
12797 0 0 0
12798 0 0 0
12799 0 0 0
12800 0 0 0
GAMMA.Hamburger.SV GAMMA.Hannover.96 GAMMA.Hertha.BSC.Berlin
12795 0 0 0
12796 0 0 0
12797 0 0 0
12798 1 0 0
12799 0 0 0
12800 0 0 0
GAMMA.SV.Werder.Bremen GAMMA.VfB.Stuttgart GAMMA.VfL.Wolfsburg Y3
12795 0 0 0 1
12796 0 0 0 1
12797 0 -1 0 2
12798 0 0 0 1
12799 0 0 1 2
12800 1 0 0 1
> utils:::menuInstallPkgs()
trying URL'https://mirrors.tuna.tsinghua.edu.cn/CRAN/bin/windows/contrib/3.2/gvcm.cat_1.9.zip'
Content type 'application/zip' length109702 bytes (107 KB)
downloaded 107 KB
package ‘gvcm.cat’ successfully unpackedand MD5 sums checked
The downloaded binary packages are in
C:\Users\Administrator\AppData\Local\Temp\RtmpMnDuOg\downloaded_packages
> library(gvcm.cat)
Loading required package: Matrix
Loading required package: MASS
Loading required package: mgcv
Loading required package: nlme
This is mgcv 1.8-12. For overview type'help("mgcv-package")'.
Attaching package: ‘mgcv’
The following object is masked from‘package:gtools’:
scat
The following object is masked from‘package:VGAM’:
s
Attaching package: ‘gvcm.cat’
The following object is masked from‘package
rdBTL’:
design
>f<-Y3~1+GAMMA.1.FC.Kaiserslautern+GAMMA.1.FC.Koeln+GAMMA.1.FC.Nuernberg+
+ GAMMA.1.FSV.Mainz.05+GAMMA.Arminia.Bielefeld+GAMMA.Bayer.Leverkusen+
+ GAMMA.Borussia.Dortmund+GAMMA.Borussia.Moenchengladbach+
+ GAMMA.Eintracht.Frankfurt+GAMMA.FC.Bayern.Muenchen+GAMMA.FC.Schalke.04+
+ GAMMA.Hamburger.SV +GAMMA.Hannover.96+GAMMA.Hertha.BSC.Berlin+
+ GAMMA.SV.Werder.Bremen+GAMMA.VfB.Stuttgart +GAMMA.VfL.Wolfsburg
>
> m1=gvcm.cat(f,des.nohome, family=gaussian ,plot=TRUE)
> summary(m1)
Call: gvcm.cat(formula = f, data = des.nohome, family = gaussian, plot = TRUE)
Deviance Residuals:
Min 1Q Median 3Q Max
-1.302 -0.58 -0.08 0.552 1.837
Coefficients:
coefficients.oml coefficients
Intercept 0.8301 0.83006536
GAMMA.1.FC.Kaiserslautern -0.0833 -0.08333333
GAMMA.1.FC.Koeln -0.0278 -0.02777778
GAMMA.1.FC.Nuernberg -0.2778 -0.27777778
GAMMA.1.FSV.Mainz.05 -0.1944 -0.19444444
GAMMA.Arminia.Bielefeld -0.1389 -0.13888889
GAMMA.Bayer.Leverkusen -0.4444 -0.44444444
GAMMA.Borussia.Dortmund -0.3611 -0.36111111
GAMMA.Borussia.Moenchengladbach -0.2778 -0.27777778
GAMMA.Eintracht.Frankfurt -0.1389 -0.13888889
GAMMA.FC.Bayern.Muenchen -0.8611 -0.86111111
GAMMA.FC.Schalke.04 -0.6389 -0.63888889
GAMMA.Hamburger.SV -0.6944 -0.69444444
GAMMA.Hannover.96 -0.2500 -0.25000000
GAMMA.Hertha.BSC.Berlin -0.3889 -0.38888889
GAMMA.SV.Werder.Bremen -0.7500 -0.75000000
GAMMA.VfB.Stuttgart -0.3333 -0.33333333
GAMMA.VfL.Wolfsburg -0.1389 -0.13888889
coefficients.reduced coefficients.refitted
Intercept 0.83006536 0.8301
GAMMA.1.FC.Kaiserslautern -0.08333333 -0.0833
GAMMA.1.FC.Koeln -0.02777778 -0.0278
GAMMA.1.FC.Nuernberg -0.27777778 -0.2778
GAMMA.1.FSV.Mainz.05 -0.19444444 -0.1944
GAMMA.Arminia.Bielefeld -0.13888889 -0.1389
GAMMA.Bayer.Leverkusen -0.44444444 -0.4444
GAMMA.Borussia.Dortmund -0.36111111 -0.3611
GAMMA.Borussia.Moenchengladbach -0.27777778 -0.2778
GAMMA.Eintracht.Frankfurt -0.13888889 -0.1389
GAMMA.FC.Bayern.Muenchen -0.86111111 -0.8611
GAMMA.FC.Schalke.04 -0.63888889 -0.6389
GAMMA.Hamburger.SV -0.69444444 -0.6944
GAMMA.Hannover.96 -0.25000000 -0.2500
GAMMA.Hertha.BSC.Berlin -0.38888889 -0.3889
GAMMA.SV.Werder.Bremen -0.75000000 -0.7500
GAMMA.VfB.Stuttgart -0.33333333 -0.3333
GAMMA.VfL.Wolfsburg -0.13888889 -0.1389
(Dispersion parameter for gaussian familytaken to be 0.5603667)
Null deviance: 412 on 306 degrees of freedom
Residual deviance: 161.3856 on 288 degreesof freedom
Removed parameters: 0 out of 0
Penalization parameter lambda = TRUE
Tuning: adapted.weights = FALSE,assured.intercept = TRUE
Number of iterations: 2
The model converged.
> plot(m1)
Error in plot.gvcm.cat(m1) :
Input argument 'plot' must be 'TRUE' for plotting coefficient paths.