similar to: function ns()

I am fittling a spline to a variable in a regression model, I am then using the predict.glm funtion to make some predictions. When I use bs to fit the spline I don't have any problems using the predict.glm function however when I use ns I get the following error: Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : variable lengths differ (found for

I have been trying to find a command similiar to bs or ns to add a periodic/or cyclical splie term to a regression model in glm.nb I have been using ns to fit a spline to day of the year (1:365): glm.nb(CNT ~ WKDY + ns(DY,df=6) + HOLIDAY_FLAG + Trend) but I think a periodic spline might be more appropriate. Any suggestions would be appreciated. best, Spencer Jones, M.Stat. NLM Fellow Dept.

Hello: I have a question about spline: ns (or bs) in R (library: splines). What's the exact analytic form of the B splines generated by ns (or bs)? I need to calculate the penalty (square of second order derivative) of the basis functions. Is there any shortcut or I need to program myself? thanks, Shuangge Ma ********************* Department of Statistics Phone: 608-263-4782(O) University

In the following case, library(splines) tt <- c(55, 251, 380, 289, 210, 385, 361, 669) nn <- rep(0:7,tt) ns(nn,4) ## knots are located at (0.25,0.5,0.75); quantiles = (2,5,7) we get the error Error in qr.default(t(const)) : NA/NaN/Inf in foreign function call (arg 1) because the 75th quantile (the location of the last "interior" knot) ends up on the boundary. As a

I have an additive model of the following form : zmdlfit <- lm(z~ns(x,df=6)+ns(y,df=6)) I can get the fitted values and plot them against z easily enough, but I also want to both obtain and plot the two additive components (the estimates of the two additive terms on the RHS) I've been looking at manuals and searching on the internet and searching the archives, but I'm apparently

Hi, I never had a look on MS AD regarding that, to check if NS record is created for each DC running DNS service or not. Anyway Samba AD does not create them and that's not an issue as long as you don't plan to make your AD DNS zones public (available on the Internet). As far as I understood DNS, during DNS resolution, NS is used only when client's resolver does not know how to

On Tue, 18 Oct 2016 14:59:31 +0200 mathias dufresne via samba <samba at lists.samba.org> wrote: > Anyway NS records are used when DNS server speak to DNS server, not by > clients. So AD would work just fine without them. > > NS are used when a client ask something the configured resolver can't > resolve by himself and when the resolver is not configured to forward >

On Thu, 18 Aug 2016 16:59:51 +0200 mathias dufresne via samba <samba at lists.samba.org> wrote: > Hi, > > I never had a look on MS AD regarding that, to check if NS record is > created for each DC running DNS service or not. I couldn't find anything that explicitly says that each DC should have its own SOA in AD. What I could find was that each AD-integrated zone on a DC

Goodmorning, This is a documentation related question about the B-spline function in R. In the help file it is stated that: "df degrees of freedom; one can specify df rather than knots; bs() then chooses df-degree-1 knots at suitable quantiles of x (which will ignore missing values)." So if one were to specify a spline with 6 degrees of freedom (and no intercept) then a basis

2016-10-18 15:25 GMT+02:00 Rowland Penny via samba <samba at lists.samba.org>: > On Tue, 18 Oct 2016 14:59:31 +0200 > mathias dufresne via samba <samba at lists.samba.org> wrote: > > > Anyway NS records are used when DNS server speak to DNS server, not by > > clients. So AD would work just fine without them. > > > > NS are used when a client ask

Are there any general guidelines for writing functions that appear in model formulas for functions like lm()/glm()? The ones I can think of are ns()/bs() in `splines' and s() in `mgcv'. It seems that ns()/bs() basically create matrix with a number of attributes and a class, while s() returns a list of parameters. It seems that gam() in `mgcv' has its own way of dealing with

Hi - How can I 'manually' reproduce the results in 'pred1' below? My attempt is pred_manual, but is not correct. Any help is much appreciated. library(splines) set.seed(12345) y <- rgamma(1000, shape =0.5) age <- rnorm(1000, 45, 10) glm1 <- glm(y ~ ns(age, 4), family=Gamma(link=log)) dd <- data.frame(age = 16:80) mm <- model.matrix( ~ ns(dd$age, 4)) pred1 <-

Hi, I run a little router/firewall that also acts as a print server for the linux and win98 clients on the LAN. I have the harddisk configured so that if there is little activity, it spins down (hdparm -S 12 /dev/hda). I added smbd/nmbd today and with the two services running, the hard disk drive spins down but wakes up about every 12 minutes. Without the services, only a major event causes them

Hey everyone. I've had an AD domain running on a samba server for years now. I recently decided to add a second samba DC. (both 4.4.5) I ran through all the DNS checks mentioned here... https://wiki.samba.org/index.php/Check_and_fix_DNS_entries_on_DC_joins All the entries were there! However, I noticed that the domain only has one NS record present. "host mydomain.com" returns both

Anyway NS records are used when DNS server speak to DNS server, not by clients. So AD would work just fine without them. NS are used when a client ask something the configured resolver can't resolve by himself and when the resolver is not configured to forward request to relevant DNS server. IE: client search for toto.org and its resolver does not know anything about that zone. Resolver will

all.vars works fine, EXCEPT, it give a bit too much. I only want the regression variables, but in the following example I also get "k" the variable holding the chosen knots. Any machinery to find only "real" regression variables? cheers, Bendix library( splines ) y <- rnorm(100) x <- rnorm(100) k <- -1:1 ml <- lm( y ~ bs(x,knots=k) ) mg <- glm( y ~

Hello, I am attempting to use R as part of a time-series analysis investigating the influence of meteorological factors on health outcomes. The test csv dta file that i am working with contains a complete daily set of the variables ?DATE?, ?ADMINS?, ?NOX?, ?TEMP? across a 5 year period (1827 days). Within attempting to fit a GLM, I believe that I am having difficulty with the ns function,

I am fitting a cox PH model w/ 2 predictors, x1 = 0/1 treatment variable and x2=continuous variable. I am using natural splines (ns) to model the effect of x2. I would like to examine the estimated effect of x2 on the hazard. I have tried various approaches (below; let model.fit= fitted model using coxph in survival library): 1. The simplest method appears to be using termplot(model.fit).

Hi Akshay, (Please include r-help when replying) You have learned that PFC.NS and snl[[159]] are not identical. Now you have to figure out why they differ. This could also point to a bug or a logic error in your program. Figuring out how two objects differ can be a bit tricky, but with experience it becomes easier. (Some others may even have some suggestions for good ways to do it.) Basically

Hello, My domain was provisioned using a Windows 2008R2 server as the first DC, and I then subsequently joined a Samba 4.4.2 DC running on Ubuntu 14.04.4 and using BIND 9.9.5 as the DNS backend. Everything seems to be working properly after I added an NC replica for the zones as per the troubleshooting page. All zones resolve and replicate properly, but I noticed that there are no NS records