Dennis, David, The code I'm running is from here: http://www.r-bloggers.com/nonparametric-copula-density-estimation/
s comes from here: s {mgcv} *Defining smooths in GAM formulae* *Package: * mgcv *Version: * 1.7-24 *Description* Function used in definition of smooth terms within *gam*<http://inside-r.org/r-doc/mgcv/gam> model formulae. The function does not evaluate a (spline) smooth - it exists purely to help set up a model using spline based smooths. *Usage* s(..., k=-1,fx=FALSE,bs="tp",m=NA,by=NA,xt=NULL,id=NULL,sp=NULL) The mgcv package is installed and loaded la function persp is in the package graphics persp {graphics} *perspective plots* *package: * graphics *version: * r 3.0.1 try installing graphics and tells me it can not be installed in version 3.0.1 of R 2013/8/3 marta Garcia <marta000gar...@gmail.com> > The code I'm running is from here: > > http://www.r-bloggers.com/nonparametric-copula-density-estimation/ > > s comes from here: > > > > s {mgcv} > > *Defining smooths in GAM formulae* > > *Package: * > > mgcv > > *Version: * > > 1.7-24 > > *Description* > > Function used in definition of smooth terms within > *gam*<http://inside-r.org/r-doc/mgcv/gam> model > formulae. The function does not evaluate a (spline) smooth - it exists > purely to help set up a model using spline based smooths. > > *Usage* > > s(..., k=-1,fx=FALSE,bs="tp",m=NA,by=NA,xt=NULL,id=NULL,sp=NULL) > > The mgcv package is installed and loaded > > > > > > > la function persp is in the package graphics persp {graphics} > > *perspective plots* > > *package: * > > graphics > > *version: * > > r 3.0.1 > try installing graphics and tells me it can not be installed in version > 3.0.1 of R > > > > > > > 2013/8/3 David Winsemius <dwinsem...@comcast.net> > >> >> >> Sent from my iPhone >> >> On Aug 3, 2013, at 8:27 AM, marta Garcia <marta000gar...@gmail.com> >> wrote: >> >> > ---------- Forwarded message ---------- >> > From: marta Garcia <marta000gar...@gmail.com> >> > Date: 2013/8/3 >> > Subject: (nonparametric) Copula density estimation >> > To: r-help@r-project.org >> > >> > >> > Hello community, sorry for the English, I speak Spanish ... >> > I try to run this code >> > >> > >> > install.packages("MASS") >> > >> > library(MASS) >> > install.packages ("evd") >> > library(evd) >> > >> > >> > X=lossalae >> > U=cbind(rank(X[,1])/(nrow(X)+1),rank(X[,2])/(nrow(X)+1)) >> > The standard tool to plot nonparametric estimators of densities is to >> use >> > multivariate kernels. We can look at the density using >> > >> > mat1=kde2d(U[,1],U[,2],n=35) >> > install.packages("mgcv") >> > library(mgcv) >> > >> > >> > persp(mat1$x,mat1$y,mat1$z,col="green" >> > , shade=TRUE,theta=s*5 >> > , xlab="",ylab="",zlab="",zlim=c(0,7)) >> > >> > >> > I get this error >> > >> > Error in s * 5 : non-numeric argument to binary operator >> > >> > Can anyone help? >> > >> The error message looks perfectly sensible: there is no numeric object >> named 's'. >> >> -- >> David >> > >> > I have R 3.0.1 >> > windows 8 >> > >> > Thanks in advance >> > >> > [[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. >> > > [[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.