similar to: NA's in segmented

Am I trying to perform multiple linear regressions on each 'VARIABLE2'. I figured out that there are different ways, using the following code: (data is given at the end of this message) reg <- lapply(split(TRY, VARIABLE2), function(X){lm(X2 ~ X3, data=X)}) lapply(reg, summary) Which produces the following: $`1` Call: lm(formula = X2 ~ X3, data = X) Residuals: Min

Hello all, I have 3 time series (tt) that I've fitted segmented regression models to, with 3 breakpoints that are common to all, using code below (requires segmented package). However I wish to specifiy a zero coefficient, a priori, for the last segment of the KW series (green) only. Is this possible to do with segmented? If not, could someone point in a direction? The final goal is to

Hi everyone, while trying to use 'segmented' (R i386 2.15.0 for Windows 32bit OS) to determine the breakpoint I got stuck with an error message and I can't find solution. It is connected with psi value, and the error says: Error in seg.glm.fit(y, XREG, Z, PSI, weights, offs, opz) : (Some) estimated psi out of its range This is the code I am using:

Hi everyone, I have encountered this problem while using 'segmented' plugin for R i386 2.15.0 (for Windows 32bit OS) and I just cannot find neither explanation nor solution for it. I am trying to run this data gpp temp 1.661 5 5.028 10 9.772 15 8.692 20 5.693 25 6.293 30 7.757 5 4.604 10 8.763 15 8.134 20 4.616 25 8.417 30 3.483 5 5.046 10 8.306 15 9.142 20 4.686 25 7.301 30 and with

Hi Folks, I've been playing around with SYSLINUX for a little bit; here's what I'm trying to do... My disk will have three partitions: [BOOT: (FAT16 w/ SYSLINUX)] [OS1] [OS2] OS1 and OS2 will both contain an XP Embedded system which is read-only and which will always boot from a hibernate file. What I'd like to be able to do if flip-flop between them at runtime (that is, run

Hi all, Suppose I want to set the values in a column to the log of the values of another column like so: object$LogDist <- log10(object$Distance) How do I loop through the objects if I have object1, object2, etc to perform this function? object1$LogDist <- log10(object1$Distance) object2$LogDist <- log10(object2$Distance) object3$LogDist <- log10(object3$Distance) I was trying to

Hi Peter, On Mon, May 14, 2018 at 8:14 AM Peter Smith via llvm-dev < llvm-dev at lists.llvm.org> wrote: > My understanding from the RFC is: > - All global objects in the bitcode file will be assigned a section name. > ... which is equal to the section name that they would have been emitted to if this was a regular compilation. In addition to allowing the linker to read section

I would expect that several math operations should always return values with a class of numeric. If the input is defined with multiple classes, however, the class attribute is preserved. I would think this may have some unintended side-effects. Here's an example: > sessionInfo()$R.version$version.string [1] "R version 3.4.0 (2017-04-21)" > x <- seq.int(5) >

Hello all, This is the error message that I get. > hyp.res <- nls(log(y)~log(pdf.hyperb(theta,X)), data=dataModel, + start=list(theta=thetaE0), + trace=TRUE) 45.54325 : 0.1000000 1.3862944 -4.5577142 0.0005503 3.728302 : 0.0583857346 0.4757772859 -4.9156128701 0.0005563154 1.584317 : 0.0194149477 0.3444648833 -4.9365149150 0.0004105426 1.569333 :

quite some speed improvement indeed. attached the updated patch to apply to svn/trunk. j -------------- next part -------------- A non-text attachment was scrubbed... Name: theora-mmx.patch.gz Type: application/x-gzip Size: 8648 bytes Desc: not available Url : http://lists.xiph.org/pipermail/theora-dev/attachments/20040824/5a5f2731/theora-mmx.patch-0001.bin

I think you may be confusing (S3) class and ?mode. > x <- seq.int(1:3) > class(x) [1] "integer" > mode(x) [1] "numeric" > class(x+.5) ## coercion [1] "numeric" > mode(x+.5) [1] "numeric" But note: > y <- as.integer(1) > class(y) [1] "integer" > class(y) <- "foo" > mode(y) [1] "numeric" >

On Friday 17 October 2014 11:25:03 Olaf Hering wrote: > On Thu, Oct 16, Olaf Hering wrote: > > I can probably add a check somewhere to catch the whatever=="@". A > > better fix would be to check if a given subvolume is for the entire > > partition. > > So after talking to David Sterba there is no way to tell if a given > subvolume is for the entire

I read a dataset with times in them, e.g., "09:31:29.18761". I then parse them: > all$X.Time <- strptime(all$X.Time, format = "%H:%M:%OS6"); and get a vector of NAs (how do I check that except for a visual inspection?) then I do > options("digits.secs"=6); > all$X.Time <- strptime(all$X.Time, format = "%H:%M:%OS"); and it, apparently, works:

In the vcdExtra package on R-Forge, I have functions and generic methods for calculating log odds ratios for R x C x strata tables. I'd like to define methods for fitting weighted lm()s to the resulting loddsratio objects, but I'm having problems figuring out how to do this generally. # install.packages("vcdExtra", repos="http://R-Forge.R-Project.org")

Thanks Bert, I think we agree on the current behaviour, but I'm still not sure if it's desirable. The mode isn't used for method dispatch. In the following example, I have to write `log.foo` in order for the correct method to be called. > x <- seq.int(5) > class(x) <- c("integer", "foo") > half <- function(x) UseMethod("half") >

Listmates, I have a fresh openSuSE 10.3 install with all updates installed, the video card is a MSI nVidia 8600GT (pci-e 16). After following the openSuSE nVidia install and the Compiz Fusion install (and after log outs and reboots) compiz refused to start. As shown on the openSuSE Compiz page he file .config/compiz/compiz-manager was created with: COMPIZ_OPTIONS="--no-libgl-fallback

Hello, first of all, thanks to LT for \pkg{codeutils}. I agree that it is indeed very useful to identify errors and also to encourage re-thinking past solutions. My problem: I want to compare different sets of related sub-functions which should be used alternatively by the same top-level function. Sets of related functions should be bound together (as lists) and the workspace should be as clean

I have a matrix of frequency counts from 0-160. x<-as.matrix(c(0,1,0,0,1,0,0,0,1,0,0,0,0,1)) I would like to apply a function creating a new column (x[,2])containing values equal to: a) log(x[m,1]) if x[m,1] > 0; and b) for all x[m,1]= 0, log(next x[m,1] > 0 / count of preceding zero values +1) for example, x[1,2] should equal log(x[2,1]/2) = log(1/2) = -0.6931472 whereas x[3,2] should

This is a very simple question but I couldn't form a site search quesry that would return a reasonable result set. Say I have a vector: x <- c(0,2,3,4,5,-1,-2) I want to replace all of the values in 'x' with the log of x. Naturally this runs into problems since some of the values are negative or zero. So how can I replace all of the positive elements of x with the log(x) and the

Hi- I'm gettting the following error message when trying to use the segmented function to look for breakpoints in my data. Error in segmented.glm(glm, seg.Z = ~segmentdist, psi = 2, control = seg.control(display = F), : (Some) estimated psi out of its range Here are some real data and the models I'm calling which gives the error above. > segmentdist [1] 0.000000 8.547576