Hello, You may have more luck posting your question to the R-SIG-Geo mailing-list:
https://stat.ethz.ch/mailman/listinfo/R-SIG-Geo/ Be sure to use an appropriate "Subject" line, for example, the particular package/function that seems problematic. HTH, Bill. W. Michels, Ph.D. On Mon, Aug 30, 2021 at 1:36 AM SITI AISYAH ZAKARIA <aisyahzaka...@unimap.edu.my> wrote: > > Dear all, > > Can anyone help me? I'm using this coding to get the spatial gev model and > plotting the quantile plot to justify the model is fit but the plot look > like not fit.. How can I solve this? It is I need to change in the any > value at the coding.For example I want to change form.shape <-shape ~1 to > form.shape <- shape ~average shape but is not available for average coding. > > below is my coding > > #---------------------------------------------------------------------------------------------- > # fitspatgev > #---------------------------------------------------------------------------------------------- > # response surface model > Ozone<-data.matrix(Ozone_S) > Ozone > LotLatAlt<-data.matrix(OzoneLotLatAlt_S) > LotLatAlt > form.loc <- loc ~ Lon + Lat + Alt > form.scale <- scale ~ 1 > form.shape <- shape ~ 1 > dim(Ozone_S) > dim(OzoneLotLatAlt_S) > fit1 <- fitspatgev(Ozone, scale(LotLatAlt,scale=FALSE), form.loc, > form.scale, form.shape);fit1 > TIC(fit1) > fit1$param > #data(rain > #symbolplot(rain, coord, plot.border = swiss) > #check the fit of the model: compute QQplots for each station > par(mfrow=c(1,3)) > par(mar=c(3,2.5,1.5,0.5),mgp=c(1.5,0.5,0),font.main=1,cex=0.66,cex.main=1) > #calculation of confidance intervals > nc <- 10000 > M1 <- matrix(rfrechet(nc*nobs),nrow=nobs,ncol=nc) > M <- t(apply(M1,2,sort)) > E <- boot::envelope(mat=M) #compute 95% confidance bands > for (k in c(1:3,4,5,6)){ #choose some stations > park <- predict(fit1)[k,] > fk <- gev2frech(Ozone[,k],loc=park[4],scale=park[5],shape=park[6]) > qqplot(y=fk,x=qfrechet((1:nobs)/(nobs+1)),log='xy',main=k,ylab='Sample > Quantiles',xlab='Theoretical Quantiles',cex=0.7); > abline(0,1) > lines(y=E$overall[1,],x=qfrechet((1:nobs)/(nobs +1)),lty='dotted') > lines(y=E$overall[2,],x=qfrechet((1:nobs)/(nobs +1)),lty='dotted') > } > > the quantile plot is in attachment file. > please, can anyone help me? > > thank you > > -- > > > > > > "..Millions of trees are used to make papers, only to be thrown away > after a couple of minutes reading from them. Our planet is at stake. Please > be considerate. THINK TWICE BEFORE PRINTING THIS.." > > DISCLAIMER: This email > and any files transmitted with it are confidential and intended solely for > the use of the individual orentity to whom they are addressed. If you have > received this email in error please notify the UniMAP's Email > Administrator. Please note that any views or opinions presented in this > email are solely those of the author and do not necessarily represent those > of the university. Finally, the recipient should check this email and any > attachments for the presence of viruses.The university accepts no liability > for any damage caused by any virus transmitted by this email. > > Universiti > Malaysia Perlis (UniMAP) | Digital Management & Development Centre (DMDC), > Universiti Malaysia Perlis (UniMAP), Pauh Putra Campus, 02600 Arau, Perlis, > MALAYSIA | www.unimap.edu.my <http://www.unimap.edu.my/> > ______________________________________________ > R-help@r-project.org 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. ______________________________________________ R-help@r-project.org 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.