search for: ldat

...I don't know much about the application, but apparently there are good empirical reasons to use a quadratic model. The goal is to find the X value which maximizes the response Y, and to find a confidence interval for this X value. Finding the maximizing X value is pretty straightforward: >Ldat <- data.frame("X"=c(3,7,14,24), "Y"=c(1,5,8,0)) >(LM<-lm(formula = Y ~ X + I(X^2), data = Ldat)) Call: lm(formula = Y ~ X + I(X^2), data = Ldat) Coefficients: (Intercept) X I(X^2) -3.86978 1.77762 -0.06729 > DZ<-function(B,C) { (-B)/...

..., { if (!is.null(userProb)) logp <- log(userProb(RAT, X)) else { y <- X %*% RAT if (!is.null(positive)) { dA.p <- dgamma(RAT[ind.pos], krat[ind.pos], lrat[ind.pos], log = TRUE) kdat <- y^2/sddat^2 ldat <- y/sddat^2 kdat[Dat == 0] <- 1 ldat[Dat == 0 & !is.na(Dat)] <- 1/y[Dat == 0 & !is.na(Dat)] dB.p <- dgamma(Dat[, positive], kdat[, positive], ldat[, positive], log = TRUE) } else { dA.p...

Hello All, I have a couple of stacked histograms which I need to compare/evaluate for similarity or difference. http://r.789695.n4.nabble.com/file/n4635668/Selection_011.png I believe rather than evaluating histograms is will be east to work with dataset used to plot these stacked histograms, which is in format: RED PURPLE BLUE GREY

...uot;R BATCH --vanilla <this>" > > if(paste(R.version$major, R.version$minor, sep=".") >= 1.7) > RNGversion("1.6") > set.seed(47) > > ## Here: Want usual "noisy" model; almost no printing > n <- 5000 > x <- rnorm(n) > ldat <- > data.frame(x1 = x, > x2 = sort(5*x - rnorm(n)), > f1 = factor(pmin(12, rpois(n, lam= 5))), > f2 = factor(pmin(20, rpois(n, lam= 9))), > f3 = factor(pmin(32, rpois(n, lam= 12)))) > with(ldat, > ldat$y...

...ally casts it and manipulates it using functions in pkg2.so. This works perfectly, because type1 is exactly the same in pkg1.so and pkg2.so. *But* UBSAN sanitizers give a runtime error, which arguably is a false positive. Real example on CRAN: https://www.stats.ox.ac.uk/pub/bdr/memtests/gcc-UBSAN/ldat/ldat-Ex.Rout A solution to this would be to dlopen pkg1.so with RTLD_GLOBAL, instead of RTLD_LOCAL, i.e., dyn.load(local=FALSE). So my proposal is to automatically set RTLD_GLOBAL for those packages that are listed at the same time in Depends/Imports/Suggests and LinkingTo, at least for those mach...

Dear All,I encountered a problem when I use RWeka for prediction. Specifically, I use the following: res=J48(X1~.,data=mydata); predict(res), #it worked fine but when I tried to use a different data set, i.e. predict(res,newdata=mynewdata); all the predictions I get is 0, which apparently is problematic. What is weird is, if I use the old data, but use the newdata option, i.e.

Good morning to you all, Sorry for taking your time from your research and teaching schedules.   If you have a non-stationary univariate time Series data that has the transformation: Say; l.dat<-log (series) d.ldat<-diff (l.dat, differences=1) and you fit say arima model. predit.arima<-predict (fit.series, n.ahead=10, xregnew= (n+1) :( n+10)) How could I re-transform "prediction$pred" to the level data since it has been differenced once? I know exp (prediction$pred) will bring the inverse of t...

Hi, I have 2 questions on list object:   1. Suppose I have a matrix like: dat <- matrix(1:9,3)   Now I want to replicate this entire matrix 3 times and put entire result in a list object. Means, if "res" is the resulting list then I should have:   res[[1]]=dat, res[[2]]=dat, res[[3]]=dat   How can I do that in the easilest manner?   2. Suppose I have 2 list objects: list1 <- list2

...sted that this thread would be a better fit for r-devel.) Thanks. In the example below I used rm() explicitly, but in general users wouldn't do that. One of the reasons for the large number of file handles is that sometimes unnamed temporary objects are created. For example: > library(ldat) > libraty(lvec) > > a <- lvec(10, "integer") OPENFILE '/tmp/RtmpVqkDsw/file32145169fb06/lvec3214753f2af0' > b <- as_rvec(a[1:3]) OPENFILE '/tmp/RtmpVqkDsw/file32145169fb06/lvec32146a50f383' OPENFILE '/tmp/RtmpVqkDsw/file32145169fb06/lvec3214484b6...

...>> Thanks. In the example below I used rm() explicitly, but in general >> users wouldn't do that. >> >> One of the reasons for the large number of file handles is that >> sometimes unnamed temporary objects are created. For example: >> >>> library(ldat) >>> libraty(lvec) >>> >>> a <- lvec(10, "integer") >> OPENFILE '/tmp/RtmpVqkDsw/file32145169fb06/lvec3214753f2af0' >>> b <- as_rvec(a[1:3]) >> OPENFILE '/tmp/RtmpVqkDsw/file32145169fb06/lvec32146a50f383' >> OPENF...

...etter fit for r-devel.) > > Thanks. In the example below I used rm() explicitly, but in general users > wouldn't do that. > > One of the reasons for the large number of file handles is that sometimes > unnamed temporary objects are created. For example: > >> library(ldat) >> libraty(lvec) >> >> a <- lvec(10, "integer") > OPENFILE '/tmp/RtmpVqkDsw/file32145169fb06/lvec3214753f2af0' >> b <- as_rvec(a[1:3]) > OPENFILE '/tmp/RtmpVqkDsw/file32145169fb06/lvec32146a50f383' > OPENFILE '/tmp/RtmpVqkDsw/fil...

The incoming (Class 5 switch to Asterisk PBX) side of a SIP trunk does not work unless I add "insecure=very" to my "Outgoing settings", but I don't want to do that. I do want to authenticate. Outgoing (Asterisk PBX to Class 5 switch) calls do authenticate and work. The Nortel CS 1500 I'm using as the PSTN-side of my SIP trunk has a username and password that it's