similar to: use of glm and cubic splines

Here is an answer to a 1999 post. I didn't find a direct answer anywhere on the Web, perhaps because it is "obvious" once one sees it. Suppose you have data from a longitudinal study, where each subject was measured *up to* four times, with missing measurements, so that the data look like this: > MAT<- structure(c(23, 24, 6, 19, 16, 20, 13, 11, NA, 8, NA, 21, 19, 15, 11,

Hi - I am fitting various Cox PH models with spline predictors. After fitting the model, I would like to use termplot() to examine the functional form of the fitted model (e.g., to obtain a plot of the relative risk (or log r.r.) versus the predictors). When there is only 1 predictor in the model, termplot returns a "?". In this case, I have not been able to figure out how to create

Hello everybody, I would like to use the SpDep-package (especially the Local Moran index analysis and the Getis-Ord statistics) in R for analysing my data. However, I don't have x-y coordinates, but my data is in a distance matrix format. Is it possible to use the SpDep package with predefined distances as well instead of letting the program determining the distance itself by the function

Hi, the first thing that comes to mind, did you set ip_forward to enable in /etc/sysctl.conf ? net.ipv4.ip_forward = 1 Should explain why you IP on the bridge works but not on the vms. Regards, Michel On Wed, 2020-06-17 at 09:43 -0400, Lamar Owen wrote: > As part of my initial KVM host on C8 deployment, I decided to set up > some HA features on the new host, specifically NIC

Dear forum,   I am working in R 2.9.1 and I am trying to sample locations from a network file. Reading in nor plotting is a problem, however when I am trying to sample from the file I get the following message:   nwlim<-readShapeLines("C:/Limburg_nwshape", proj4string=CRS("+init=epsg:31300")) plot(nwlim)   randacc<-spsample(nwlim,n=1000,"random")

Hello, I am using R and the smooth.Pspline function in the pspline package to smooth some data by using natural cubic splines. After fitting a sufficiently smooth spline using the following call: (ps=smooth.Pspline(x,y,norder=2,spar=0.8,method=1) [the values of x are age in years from 1 to 100] I tried to check that R in fact had fitted a natural cubic spline by checking that the resulting

re: Cubic splines in package "mgcv" I don't have access to Gu (2002) but clearly the function R(x,z) defined on p126 of Simon Wood's book is piecewise quartic, not piecewise cubic. Like Kunio Takezawa (below) I was puzzled by the word "cubic" on p126. As Simon Wood writes, this basis is not actually used by mgcv when specifying bs="cr". Maybe the point is

R-users E-mail: r-help@r-project.org My understanding is that package "mgcv" is based on "Generalized Additive Models: An Introduction with R (by Simon N. Wood)". On the page 126 of this book, eq(3.4) looks a quartic equation with respect to "x", not a cubic equation. I am wondering if all routines which uses cubic splines in mgcv are based on this quartic

Dear R-helpers, I need to fit natural cubic spline with specified number of knots. I expected 'splines' package will be helpful, but I am confused by its help. Is more detailed documentation available for it or could you recommend another R function? Best regards Ondrej Mikula

You wrote: "It may be best to either write to the package maintainer (me, as you did) or post to the group but not both." This is just a note that I disagree wrt my own packages: I go on vacation or trips, or have other projects so won't always answer Other folks on the list often have good ideas that I'd miss My preferred standard is "ask the list, with a copy

In the book, "Regression Modeling Strategies", Frank Harrell uses the function rcs() often. The current version or R and Hmisc library don't appear to have this function. What is an appropriate substitute? Thank You -- Noah Silverman UCLA Department of Statistics 8208 Math Sciences Building Los Angeles, CA 90095 [[alternative HTML version deleted]]

Hi, On http://www.speex.org/docs/manual/speex-manual/node11.html#SECTION001140000000000000000 the bitrate for Speex Wideband mode 9 and 10 is 9 34,400 10 42,400 On http://tools.ietf.org/html/draft-ietf-avt-rtp-speex-03 the bitrate for Speex Wideband mode 9 and 10 is 9 34.2 kbit/s 10 42.2 kbit/s I believe the draft contains the right values: E.g. mode 10 Narrowband mode 7 -> 492

Dear R-users, Is there an easy way to determine all possible vectors of length 21 with each entry having permutations from 0 to 4, instead of doing it like this? It really takes up too much time, and I am convinced that there exists something easier. Can you help me? Thank you in advance! Kind regards, Elke for (a in 0:4){ for (b in 0:4){ for (c in 0:4){ for (d in 0:4){

Hi all I was playing with termplot(), and came across what appears to be an inconsistency. It would appreciate if someone could enlighten me: > # First, generate some data: > y <- rnorm(100) > x <- runif(length(y),1,2) > # Now find the log of x: > logx <- log(x) > > # Now fit two models that are exactly the same, but specified differently: > m1 <-

Thanks for your prompt reply Duncan. I had indeed assumed they were what the help file says until observation raised doubts, which is why I queried it. >From reading the code for termplot(), it seems that either the predict() function doesn't return the 1x standard error, or the curves plotted by the termplot() function are not 1x standard errors. If they're not 1x standard errors,

Hi, we now have an old Windows NT4.0 domain served by Samba 4.2.x (using Samba and LDAP) and want to move to Windows AD. The reason we need to do that is because of the clients (Windows 10 and MacOS) and because of a third party device which does not want to talk to Samba AD (Isilon OneFS). I did a 'classicupgrade' to Samba AD from our Samba/LDAP config and then I use this guide :

The bug exists on R-1.9.0-alpha compiled the 10/3. Termplot has a problem if either the model only contains a single term or if asked to plot a single term. In addition there are problems with the option se = TRUE. Analysis: termplot starts with terms <- if (is.null(terms)) predict(model, type = "terms", se = se) else predict(model, type = "terms", se = se,

Hi, I think it would be nice to have par(ask=T) set by default in termplot(), like it is in plot.lm(). Best, Emmanuel Paradis -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !)

Hi together, I always get an error message with using ylim in termplot(), like this: > x<-(1:10) > y<-(10:1) > l<-lm(y~x) > termplot(l,ylim=c(1,2)) Is this a bug, or is there another possibility to do that? Especially, I would like to use term.plot() for gamlss objects. Thanks for your help! Andreas -- Echte DSL-Flatrate dauerhaft f?r 0,- Euro*!

Hello, A reviewer queried what the intervals were on the termplot I provided in a report. The help file for termplot() suggests they're standard errors (se=T), but in the code the se.fit values from predict() are multiplied by 2, suggesting it's a rough 95% confidence interval, is that right? Many thanks, Eric Goodwin Scientific data analyst | Coastal and Freshwater Group Cawthron