similar to: BCa Bootstrap confidence intervals

Björn Pettersson A <bjorn.a.pettersson at ericsson.com> writes: > In the new monorepo UC1 may or may not be a parent to UL1. > We could actually have something like this: > > UL4->UC2->UL3->UL2->UL1->UL0->UC1 > > Our DL1 commit should preferably have UL1 as parent after > conversion > > UL4->UC2->UL3->UL2->UL1->UL0->UC1 >

In linear regression, regression weights of x1 on Y given x2 and x3 should be mathematically identical to the semipartial correlations between x1 and Y, given x2 and x3. However, I do not obtain identical results, so apparently I'm doing something wrong in R. Data preparation: data<-read.csv("file.csv", head=T) data <- subset(data, select=c(age,s1,s2,s3,s5,s10,WSAS01)) data

Björn Pettersson A <bjorn.a.pettersson at ericsson.com> writes: > In llvm (split) we have: > > UL4->UL3->UL2->UL1->UL0 > \ > ...->DL2->DL1 > > In clang (split) we have: > > UC4->UC3->UC2->UC1->UC0 > \ > ...->DC2->DC1 > > > DL1 is a commit that updates the

Dear helpers, I am using R version 2.5.1 to estimate a multinomial logit model using my own maximum likelihood function (I work with share data and the default function of R cannot deal with that). However, the computer (I have an Athlon XP 3200+ with 512 GB ram) takes quite a while to estimate the model. With 3 categories, 5 explanatory variables and roughly 5000 observations it takes 2-3 min.

I like to bootstrap regression models, saving the entire set of bootstrapped regression coefficients for later use so that I can get confidence limits for a whole set of contrasts derived from the coefficients. I'm finding that ordinary bootstrap percentile confidence limits can provide poor coverage for odds ratios for binary logistic models with small N. So I'm exploring BCa confidence

Dear All, I am trying to use the bca.ci function on a boot.out object which consists a few NaN's and I want to ignore those NaN's, and get a ci only for the "normal" values. boot.out$t has R number of values for 3000 different statistics, so when I use boot.ci(boot.out, index=i) and i happens to be a column in boot.out$t with some NaN's in there I get an error message. I

Hello everybody, when I create an object of class boot with the function tsboot() from the package boot and try to compute several types of confidence intervals with boot.ci("object of class boot created with tsboot") I obtain the warning message, that "BCa-intervals are not defined for time-series bootstraps". Does that hold in general? Or is it just not defined in

i was searching for the source of bca.ci, a function of the package boot. I tried require(boot, keep.source=TRUE) but again the source was not viewable. How should i do? Best regards Meinhard Ploner ---------------- PS > version _ platform i386-apple-darwin8.8.1 arch i386 os darwin8.8.1 system i386, darwin8.8.1 status major 2

Dear List, Couple of issues while using functions from ?BCA? library: 1. I am trying to use ?lift.chart? function from ?BCA? library, but facing issues while using model where model formula is passed as formula object in glm. When model formula is written as text, then it works fine. In my case input variables and target variables are going to change dynamically, so have to used formula as

Hello List! I asked this before (with no solution), but maybe this time... I'm trying to project a surface to the XY under a 3d cloud using lattice. I can project contour lines following the code for fig 13.7 in Deepayan Sarkar's "Lattice, Multivariate Data Visualization with R", but it fails when I try to "color them in" using panel.levelplot. ?utilities.3d says there

I'm using R 2.15.1 on a 64-bit machine with Windows 7 Home Premium and package 'boot'. I've found that using a number of bootstrap resamples in boot() that is less than the number of data results in a fatal error. Once the number of resamples meets or exceeds the number of data, the error disappears. Sample problem (screwy subscripted syntax is a relic of edited down a more

Sometimes one might like to obtain pointwise bootstrap bias-corrected, accelerated (BCA) confidence intervals for a large number of statistics computed from a single dataset. For instance, one might like to get (so as to plot graphically) bootstrap confidence bands for the fitted values in a regression model. (Example: Chiu S et al., Early Acceleration of Head Circumference in Children with

I'm using R 2.15.1 on a 64-bit machine with Windows 7 Home Premium. Sample problem (screwy subscripted syntax is a relic of edited down a more complex script): > N <- 25 > s <- rlnorm(N, 0, 1) > require("boot") Loading required package: boot > v <- NULL # hold sample variance estimates > i <- 1 > v[i] <- var(s) # get sample variance >

Hi everyone, I use the boot package within the MBESS package to automate much of the hard part for folks interested in performing a (simple) mediation model. The function mediation() works well, except for a (probably) unrealistic artificial data set. When I apply the bootstrap I sometimes get: Error in t.star[r, ] <- res[[r]] which is (I think) due to Inf and -Inf and NaN values produced for

Hi, I'm trying to make a ts object that has both NA values and a frequency other than 1 (so I can use stl). I've tried all permutations I can think of, but cannot get the desired (expected?) results. The values live in x and the corresponding semi-regular time stamps are in t: > library('zoo') > z = zoo(x, order.by=t, frequency=24) > zzr = as.zooreg(z, start=0) > zr

HLVM is a garbage collected virtual machine built upon LLVM: http://forge.ocamlcore.org/projects/hlvm/ The development of HLVM was described in a series of three OCaml Journal articles that turned out to be among our most popular. Consequently, we have decided to run another series of related articles that build upon this foundation in order to develop complete compilers. The first article

Hello We are trying to implement SSO with mod_auth_ntlm_winbind. We followed the instructions on [1], but have the issue that users can not authenticate with web browsers. In addition to that document we did the following extra steps: - chown root:winbind /var/lib/samba/winbindd_privileged/ - apache user (vagrant) is in group winbind - net setauthuser -U vagrant - smb.conf has: winbind use

With the exception of "L-BFGS-B", all of the other optim() methods return the value of the function when they are given a trivial function (i.e., one with no variable arguments) to optimize. I don't think this is a "bug" in L-BFGS-B (more like a response to an undefined condition), but it leads to a bug in stats4::mle -- a spurious error saying that a better fit has been

Hi, I'd like to simply add column-wise using Matrix objects (Csparse). It looks like one can apply mosty any base function to these objects (i.e., apply, colSums), but there is a nasty conversion to traditional matrix objects if one does that. Is there any workaround? I can see colSum listed in the help for Class 'CsparseMatrix' , but I wonder whether I'm using the default

Hello, Can anyone explain the following behaviour? To me it seems a bug, but maybe it is intentional. It seems that a diff on a zooreg class that is not _strictly_ regular only considers those entries that are 'deltat' apart. In the following, diff on the zooreg class only returns values where the index was one second apart. The example replicates by dev code, but I've also tested