similar to: Piecewise continuous logistic regression with one knot

Windows Vista R 2.10 - I know it is old, I will update later today. How might I perform a piece-wise linear regression where two linear segments are separated by a single knot? In addition to estimating the slopes of the two segments (or the slope in one segment and the difference between the slope of the first and second segment), I would like the analysis to select the optimum knot. My first

windows 7, R 2.12 I am trying to run a piecewise linear regression with a single knot, i.e. a regression composed of two straight lines where the two lines intersect at an x value given by the variable knot. I wish to estimate the slope of both lines, the value of knot, the x value where the two lines intersect, and an intercept. I am using the nls code below, and get the following error

>>>>> Anupam Tyagi >>>>> on Tue, 9 Jul 2024 16:16:43 +0530 writes: > How can I do automatic knot selection while fitting piecewise linear > splines to two variables x and y? Which package to use to do it simply? I > also want to visualize the splines (and the scatter plot) with a graph. > Anupam NB: linear splines, i.e. piecewise

dear all, I apologize for my delay in replying you. Here my contribution, maybe just for completeness: Similar to "earth", "segmented" also fits piecewise linear relationships with the number of breakpoints being selected by the AIC or BIC (recommended). #code (example and code from Martin Maechler previous email) library(segmented) o<-selgmented(y, ~x, Kmax=20,

How can I do automatic knot selection while fitting piecewise linear splines to two variables x and y? Which package to use to do it simply? I also want to visualize the splines (and the scatter plot) with a graph. Anupam [[alternative HTML version deleted]]

??I would like to apply a function, fract, to the columns of a dataframe. I tried the following apply(data5NonEventEpochs,2,fract) but, no surprise it did not work as apply works on matrices not data frames. How can I apply a fuction to the columns of a data frame? (I can't covert data5NonEventsEpochs to a matrix as it contains character data). Thank you, John John David Sorkin M.D., Ph.D.

David, predict.glm and se.fit were exactly what I was looking for. Many thanks! John John David Sorkin M.D., Ph.D. Professor of Medicine Chief, Biostatistics and Informatics University of Maryland School of Medicine Division of Gerontology and Geriatric Medicine Baltimore VA Medical Center 10 North Greene Street GRECC (BT/18/GR) Baltimore, MD 21201-1524 (Phone) 410-605-7119 (Fax)

Dear All, I am trying to simulate data for a spline/piecewise regression model. I am missing something fundamental in my simulation procedure because when I try to fit my simulated data using the Gauss-Newton method in SAS, I am getting some wacky parameter estimates. Can anyone please check my simulation code and tell me what mistake I am making in generating data for spline model? Thank you

David, Again you have my thanks!. You are correct. What I want is not technically a contrast. What I want is the estimate for "regional" and its SE. I don't mind if I get these on the log scale; I can get the anti-log. Can you suggest how I can get the point estimate and its SE for "regional"? The predict function will give the point estimate, but not (to my knowledge)

Windows 7 (local computer) Windows server (server I am trying to reach) I need to read a file whose windows path is of the form \\Theserver\mydirectory\data.csv You will note that as per windows standards the server name is preceded by two backslashes. I am not sure how to specify this in R. One usually needs to specify an escape characters in a path new which would suggest my path should

> On Oct 22, 2017, at 5:01 PM, Sorkin, John <jsorkin at som.umaryland.edu> wrote: > > David, > Again you have my thanks!. > You are correct. What I want is not technically a contrast. What I want is the estimate for "regional" and its SE. There needs to be a reference value for the contrast. Contrasts are differences. I gave you the choice of two references

Windows 10 64-bit, R-Studio, R version 3.4.3 Several questions relating to groupedData: (1) I am trying to create a groupedData object that can be used to run an analysis that I have been able to urn using lmer. When I include the interaction terms in the groupedData opbject I get an error message stating that + is not meaningful for factors. How do I include factors in my model? See code and

I am trying to get the mean of a row of a data frame. My code follows: roop <- data.frame(x=c(1,2,3),y=c(4,5,2),z=c(0,9,4)) roop mean(roop[1,]) mean(roop[1,c("x","y","z")]) I get the following output: > roop x y z 1 1 4 0 2 2 5 9 3 3 2 4 > mean(roop[1,]) [1] NA Warning message: In mean.default(roop[1, ]) : argument is not numeric or logical:

Jeff, Perhaps I was sloppy with my notation: I want groups >=0 <10 >=10 <20 >=20<30 ...... >=90 <100 In any event, my question remains, why did the four different versions of cut give me the same results? I hope someone can explain to me the function of include.lowest and right in the call to cut. As demonstrated in my example below, the parameters do not seem to alter

I have installed rJava into my Windows 10 (64-bit) R instillation using the Tools > Install Packages command of my RStudion IDE. When I issued the R command in my R code library(rJava) I received the following error: library(rJava) Error: package or namespace load failed for ?rJava? in get(Info[i, 1], envir = env): lazy-load database 'C:/Users/John

How can I get the p values from a glm ? I want to get the p values so I can add them to a custom report fitwean<- glm(data[,"JWean"]~data[,"Group"],data=data,family=binomial(link ="logit")) summary(fitwean) # This lists the coefficeints, SEs, z and p values, but I can't isolate the pvalues. names(summary(fitwean)) # I see the coefficients,

Colleagues, I am trying to learn to use uniroot to solve a simple linear equation. I define the function, prove the function and a call to the function works. When I try to use uniroot to solve the equation I get an error message, Error in yfu n(x,10,20) : object 'x' not found. I hope someone can tell we how I can fix the problem

Dear John You are only allowed to have the primary covariate on the left hand side of the vertical bar. Other covariates go in inner or outer. Michael On 22/03/2018 12:59, Sorkin, John wrote: > Windows 10 64-bit, R-Studio, R version 3.4.3 > > > Several questions relating to groupedData: > > (1) I am trying to create a groupedData object that can be used to run an analysis

> On Oct 22, 2017, at 3:56 PM, Sorkin, John <jsorkin at som.umaryland.edu> wrote: > > David, > Thank you for responding to my post. > > Please consider the following output (typeregional is a factor having two levels, "regional" vs. "general"): > Call: > glm(formula = events ~ type, family = poisson(link = log), data = data, > offset =

David, Thank you for responding to my post. Please consider the following output (typeregional is a factor having two levels, "regional" vs. "general"): Call: glm(formula = events ~ type, family = poisson(link = log), data = data, offset = log(SS)) Deviance Residuals: Min 1Q Median 3Q Max -43.606 -17.295 -4.651 4.204 38.421 Coefficients: