Re: [R] package mgcv - predict with bam: Error in X[ind, ] : subscript out of bounds

the random effects for site set to zero > then one trick is to use terms like s(site,bs="re",by=dum) in fitting with > dum set to 1. Then in prediction you can set 'site' to any existing level, > and dum to zero, in order to get a prediction for the missing lev

Re: [R] package mgcv - predict with bam: Error in X[ind, ] : subscript out of bounds

Hi Simon, thank you for your reply, I really appreciate any help to understand the problem here... Unluckily the package upgrade didn't help with this issue. An example reproducing the error, and a current sessionInfo() Output can be found below. Many thanks once again, Katharina R Code

[R] package mgcv - predict with bam: Error in X[ind, ] : subscript out of bounds

Dear R-Community, I`m trying to apply the mgcv package to fill gaps in sensor data from different sites (9 sites, 2 sensors per site) and do the filling on a site-wise level. Based on http://r.789695.n4.nabble.com/mgcv-gamm-predict-to-reflect-random-s-effects-td3622738.html my model looks like th

Re: [R] Package zoo - na.rm ignored?

many thanks for the enlightenment - I guess I completely misunderstood the na.rm option here... On 18 June 2013 00:50, Gabor Grothendieck wrote: > On Mon, Jun 17, 2013 at 6:23 PM, Katharina May > wrote: > > Dear R Users, > > > > I've got a strange problem, whi

[R] Package zoo - na.rm ignored?

Dear R Users, I've got a strange problem, which I do not really understand: when I use na.spline from the zoo package, the option na.rm is just being ignored, see this adjusted example from the na.spline help page: d0 <- as.Date("2000-01-01") z <- zoo(c(NA, 11, 13, NA, 15

[R] round Date object to 10 minutes intervals?

Hello *, I'm trying to round time stamps, consisting of a date and time part, to 10 minutes intervals The timestamps have the following format: %d.%m.%Y %H:%M *timeStamp <- c("31.03.2011 09:30", "31.03.2011 09:41", "31.03.2011 10:04", "31.03.2011 10:10", "31.03.2011 10:28", "31.03.2011 10:35", "3

Re: [R] linear model: Test difference between coefficients and given values (t.test?)

ubling you: I struggling with the data analysis of my bachelor thesis and just want to compare the coefficients of my linear regression to published ones... a method often used I guess, but still I cannot find any appropriate documentations. Thanks, Katharina > On Aug 8, 2009, Kathari

[R] linear model: Test difference between coefficients and given values (t.test?)

Hi there, I've got a question which is really trivial for sure but still I have to ask as I'm not making any progress solving it by myself (please be patient with an undergraduate student): I've got a linear model (lm and lmer fitted with method="ML"). Now I want to compare the coefficients (slop

Re: [R] package lmodel2: p-value RMA fitting?

2009/7/20 Katharina May > Hi *, > > is there a way to obtain some kind of p-value for a model fitted with RMA > using the lmodel2 package? > I know that p-values are discussed and criticized a lot and as you can > image from my question I'm not > very much of a statisti

[R] package lmodel2: p-value RMA fitting?

Hi *, is there a way to obtain some kind of p-value for a model fitted with RMA using the lmodel2 package? I know that p-values are discussed and criticized a lot and as you can image from my question I'm not very much of a statistican (only writing my bachelor thesis). As fare as I understood th

Re: [R] distinguish regression lines in grouped, black and white lattice xyplot

think more creatively about how to make a > more meaningful display to provide viewers interpretable information. > > Bert Gunter > Genentech Nonclinical Biostatistics > > > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] O

[R] distinguish regression lines in grouped, black and white lattice xyplot

Hi, I've got the following problem which I cannot think of a solution right now: if got a lattice xyplot in black and white and a grouping variable with many (more than 8 values) and I plot it as regression lines (type="r"), just like this one (not reproducable but that's I guess not the point he

Re: [R] lattice logaritmic scale (basis "e" ), rewriting labels using xscale.component

thanks Deepayan, that works great! 2009/6/19 Deepayan Sarkar : > On 6/18/09, Katharina May wrote: >> Hi there, >> >>  sorry for troubling everybody once again, I've got a problem rewriting >>  Sarkar's function for >>  rewriting the tick locations

Re: [R] lattice: axis ticks, axis alignment and remove axis from plot

, Katharina 2009/6/18 Jim Lemon : > Katharina May wrote: >> >> What I sort of looking for is a xyplot with both the x axis at the bottom >> and y axis at the >> left going through 0, but continuing  in the positive and negative area >> (forming a kind

[R] lattice logaritmic scale (basis "e" ), rewriting labels using xscale.component

Hi there, sorry for troubling everybody once again, I've got a problem rewriting Sarkar's function for rewriting the tick locations in a logaritmic way (s. http://lmdvr.r-forge.r-project.org/code/Chapter08.R): His example works for log 2 but I need log e (natural logarithm). My problem is that if

Re: [R] lattice: axis ticks, axis alignment and remove axis from plot

s. > axs > character, "r" or "i". In the latter case, the axis limits are calculated > as the exact data range, instead of being padded on either side. (May not > always work as expected.) > > > > ---

[R] lattice: axis ticks, axis alignment and remove axis from plot

Hi there, I'm a bit confused concerning the axis tck setting in the lattice package as the ticks on left sided axis aren't drawn at all with the following setting: dados <- data.frame(varsep = factor(rep(1:2,10)), i = runif(20)) library(lattice) my.theme <- list( axi

[R] color in grouped lattice xyplot plot

Hi, I've got a data.frame object which I would like to plot in a lattice xyplot grouped by a variable and using a value for the color of each point which is stored in a variable in the same data.frame. I tried the following code: my.panel.xy <- function(x, y, ...){ panel.xyplot(x, y, ...)

[R] wrong labels and colors of points in graph/plot

Hi there, I trying to solve this problem for the whole day not going anywhere, so I really hope maybe somebody can help me in this community... I've got an object coefficient2 which I want to plot in differerent ways, with colors and labels added to the points, but somehow there seems to be a pro

[R] write values of points beside points in plot

Hi everyone, is it possible to write a certain value beside a point in a plot? I'm plotting the following: plot(coefficient$intercept ~ coefficient$average_BM_leaves_needles, main="intercepts ::: BM_leaves_needles", ylab="intercepts", xlab="average BM_leaves_needles per site [kg]") and I wa

[R] create string of comma-separated content of vector

Hi, how do I create a string of the comma-separated content of a vector? I've got the vector i with several numeric values as content: >str(i) num 99 and want to create a SQL statement to look like the following where the part '(2, 4, 6, 7)' should be the content of the vector i: select * from

Re: [R] Remove objects names like character String

thanks to all your solutions, works out perfectly! 2009/5/19 Henrique Dallazuanna : > Try this: > > rm(list=ls(patt="site[0-9]$")) > > On Tue, May 19, 2009 at 7:47 AM, Katharina May > wrote: >> >> Hi, >> >> how can I use rm() on objects

Re: [R] remove empty objects from workspace

t))[1]) || dim((get(object)))[1] == 0) rm(list=object) Thanks a lot! 2009/5/19 Jim Lemon : > Katharina May wrote: >> >> Hi, >> >> how can I remove all empty objects (which are NA or have zero rows) >> from my workspace? >> >> > > Hi Katharina, > T

[R] Remove objects names like character String

Hi, how can I use rm() on objects named like: paste("site",i,"_data",sep="") while looping through i? I tried rm(paste("site",i,"_data",sep="")) but I get the error that rm() must contain names or text strings which is confusing me as I thought paste() would create something like that...? Thanks,

[R] remove empty objects from workspace

Hi, how can I remove all empty objects (which are NA or have zero rows) from my workspace? Thanks, Katharina __ 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.o

Re: [R] create objects in a loop and adding sqlQuery content to them

nnel, paste("select * from [biomass_data$] where site_no = ",i,sep=""))) } Thanks for all your solution proposals! Katharina On 21.04.2009, at 02:43, David Winsemius wrote: On Apr 20, 2009, at 8:04 PM, Katharina May wrote: Hi there, I've got a database or rath

[R] create objects in a loop and adding sqlQuery content to them

Hi there, I've got a database or rather spreadsheet with several columns and rows. For one column named sites I want to loop through all possible values and retrieve all data out of the database where site = x and write it into an object named 'sitex_data'. Somehow I'm really missing someth