similar to: GAM Poisson

I am fitting a model in which the response variable y is a function of two independent, quantitative variables x1 and x2; thus: y = f(x1, x2). For reasons I do not believe to be important for the purpose of this post, I find it desirable to find f by means of GAM; also, I require principal effects and interactions to be specified separately, so I am using using te and ti tensors. Thus, I am using

There probably is a way, but it involves some programming. You would need to clone a smooth constructor (e.g. for the "cr" class), and then modify it to add a linear constraint matrix C to the returned smooth object. If b are the smooth coefficients then C should? be the matrix such that s(0) = Cb (you can get this from the Predict.matrix method for the class). Then the constraint

Thanks Simon, by cloning a smooth construct do you mean copying and modifying the smooth constructor code? Could you pleas elaborate on your answer? Which is the Predict.matrix method? 2018-01-12 3:20 GMT-06:00 Simon Wood <simon.wood at bath.edu>: > There probably is a way, but it involves some programming. You would need to > clone a smooth constructor (e.g. for the "cr"

Hello R users, I have a quick question I was hoping to get your input on. I am new to R and the smooth statistical regression world, and am trying to wrap my mind around the issues concerning using splines for mixed effect modeling. My question is the following: in the ?gamm? function, generalized additive mixed models can be estimated by including random components. These can be explicitly

Hello Jim, Thanks again. I am getting the two-headed arrow but I cannot seem to get the coordinates right for the arrow to appear beneath the map. These coordinates puts the arrow on the left hand side. Thanks again! Sincerely, Milu On Tue, Apr 12, 2016 at 1:15 PM, Jim Lemon <drjimlemon at gmail.com> wrote: > Hi Milu, > There is a two-headed arrow on the image you sent, and it

Hi Milu, My fault here. As I don't have the data to make the map and try out my suggestions I mixed up the x and y coordinates. Try this: par(xpd=TRUE) arrows(-19.75966,53,33.60000,53,code=3) par(xpd=FALSE) Jim On Tue, Apr 12, 2016 at 10:11 PM, Miluji Sb <milujisb at gmail.com> wrote: > Hello Jim, > > Thanks again. I am getting the two-headed arrow but I cannot seem to get

Dear David, The device was the issue. The quartz() device works fine but pdf() does not. Now I just need to figure out the limits for map for Europe. Thanks for all your help and patience. Sincerely, Milu On Sun, Apr 10, 2016 at 7:10 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > > On Apr 10, 2016, at 4:12 AM, Miluji Sb <milujisb at gmail.com> wrote: > >

> On Apr 10, 2016, at 1:45 PM, Miluji Sb <milujisb at gmail.com> wrote: > > Dear David, > > The device was the issue. The quartz() device works fine but pdf() does not. Now I just need to figure out the limits for map for Europe. Thanks for all your help and patience. After plotting a map of Europe with base graphics the coordinates of the lower-left and upper-right

Hello David, This is exactly what I want but I still can't get the arrows. R and R studio is updated. Thanks again! Sincerely, Milu On Sat, Apr 9, 2016 at 10:29 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > > On Apr 9, 2016, at 1:27 PM, David Winsemius <dwinsemius at comcast.net> > wrote: > > > > > >> On Apr 9, 2016, at 11:18 AM,

> On Apr 10, 2016, at 4:12 AM, Miluji Sb <milujisb at gmail.com> wrote: > > Hello David, > > This is exactly what I want but I still can't get the arrows. R and R studio is updated. Thanks again! I didn't try it in Rstudio until just now (and I don't remember that you ever mentioned RStudio as a possible issue.) The plotting I see in the default graphics Rstudio

This is what I was looking for. Thank you everyone! Sincerely, Milu <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Mail priva di virus. www.avast.com <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>

Ista, et. al: efficiency? (Note: I needed to correct my previous post: do.call() is required for pmax() over the data frame) > x <- data.frame(matrix(runif(12e6), ncol=12)) > system.time(r1 <- do.call(pmax,x)) user system elapsed 0.049 0.000 0.049 > identical(r1,r2) [1] FALSE > system.time(r2 <- apply(x,1,max)) user system elapsed 2.162 0.045 2.207 ##

Dear David, Thank you very much for your replies! I didn't know about par('usr'). I get different coordinates though: [1] -19.75966 54.75966 33.60000 71.40000 But the arrow is not at the bottom of the map. I will keep playing with this. Thanks again! Sincerely, Milu On Mon, Apr 11, 2016 at 12:00 AM, David Winsemius <dwinsemius at comcast.net> wrote: > > > >

Forgot to copy the list Dear Jim, Thank you for your reply. I must be doing something wrong, If this is my command to plot a map of Europe: eps_europe <- mapCountryData(n, nameColumnToPlot="eps_score", mapTitle="EPS Score - Europe",colourPalette=colourPalette, catMethod="fixedWidth", missingCountryCol = "white",

Hi Milu, There is a two-headed arrow on the image you sent, and it seems to be where you specified. Did you want it beneath the map, as: par(xpd=TRUE) arrows(-22,54.75,-22,74,code=3) par(xpd=FALSE) Jim On Tue, Apr 12, 2016 at 7:58 PM, Miluji Sb <milujisb at gmail.com> wrote: > Dear Jim, > > Thanks again! I do want the arrows at the bottom (beneath the map). This is > what I am

Hi Milu, I just realized that by "the bottom of the map" you may mean "beneath the map", in which case you should use: par(xpd=TRUE) arrows(...) par(xpd=FALSE) Jim On Mon, Apr 11, 2016 at 11:50 PM, Miluji Sb <milujisb at gmail.com> wrote: > Dear David, > > Thank you very much for your replies! I didn't know about par('usr'). > > I get

Dear Jim, Thanks again! I do want the arrows at the bottom (beneath the map). This is what I am doing: # Draw the map eps_europe <- mapCountryData(n, nameColumnToPlot="eps_score", mapTitle="EPS Score - Europe",colourPalette=colourPalette, catMethod="fixedWidth", missingCountryCol = "white", mapRegion="Europe", addLegend=FALSE) # ISO3 codes on

> On Apr 9, 2016, at 11:18 AM, David Winsemius <dwinsemius at comcast.net> wrote: > > >> On Apr 9, 2016, at 10:46 AM, Miluji Sb <milujisb at gmail.com> wrote: >> >> Dear David, >> >> Thank you for your answer. Sorry for the embarrassing mistake. However, even with when I generate a map for the whole world using: >> >> eps <-

Dear Bert and Huzefa, Apologies for the late reply, my account got hacked and I have just managed to recover it. Thank you very much for your replies and the solutions. Both work well. I was wondering if there was any way to ensure (force) that all possible combinations show up in the output. The full dataset has 25 cities but of course people have not moved from Boston to all the other 24

Dear David, Thank you for your answer. Sorry for the embarrassing mistake. However, even with when I generate a map for the whole world using: eps <- mapCountryData(n, nameColumnToPlot="eps_score", mapTitle="EPS Score",colourPalette=colourPalette, catMethod="fixedWidth", missingCountryCol = "white", addLegend=FALSE) And then use: