[R-sig-Geo] help for show.vgms

[Saman Monfared]

Dear All,
I  try to plot some variogram models in a panel by show.vgms in gstat.
my program is:

show.vgms(min =0, max =200000, n = 50, sill =c(v.ols$psill[2],v.wls$psill[2],
v.gls$psill[2],v.rml$sigmasq,v.reml$sigmasq,v.opvme$psill[2],v.opv$psill[2],
v.opvab$psill[2],v.opvlin$psill[2]/3),range =c(v.ols$range[2],v.wls$range[2],
v.gls$range[2],v.rml$phi,v.reml$phi,v.opvme$range[2],v.opv$range[2]
,v.opvab$range[2],v.opvlin$range[2]),models ="Sph",nugget=c(v.ols$psill[1],
v.wls$psill[1],v.gls$psill[1],0,0,v.opvme$psill[1],v.opv$psill[1],
v.opvab$psill[1],v.opvlin$psill[1]/3),plot =T,as.groups =F)
I have some problem:
1-How can I replace my desired name
c("OLS","WLS","GLS","ML","REML",,.....)instead of model names?
2- How can I  create a legend with names:
c("OLS","WLS","GLS","ML","REML",,.....)?
3- How can I change lines to points and change pch?

Best.
-- 
Saman Monfared
Msc, Department of Statistics, Shiraz University,
Shiraz 71454, Iran
Email: Samanmonfared1 at gmail.com

Tel: +98 917 5305167