search for: ci

...ines are statistically different when fitted with two different data sets; one corn, another soybean. I know I can do this in linear models by testing for interactions. See Introductory Statistics with R by Dallgaard p212-218 for an example. I have two different data sets I am comparing to lai. ci.re should have very little difference between corn and soybean, while ci.gr, there should be a difference. If I use the simplistic form described in Dallgaard (see example in code below) it does not work correctly. What do I need to do to test for this interaction? Thank you! My data is locate...

I am wanting to plot a 95% confidence band using segplot, yet I am wanting to have groups. For example if I have males and females, and then I have them in different races, I want the racial groups in different panels. I have this minor code, completely made up but gets at what I am wanting, 4 random samples and 4 samples of confidence, I know how to get A & B into one panel and C&D in to another but how do I get the x axis to label them properly and have it categorized as tw...

...fidence intervals.?I wanted to pass different sample sizes through the function, but for some reason it's not working.?n ? <- seq(from=40, to=300, by=2o) I was also wondering how I can return a plot for different sample sizes.?plot(m~d, main="n(i)")? Any help or suggestion is appreciated. ############################f<-function(n){? m = runif(n,50,200)? d = ?rnorm(n,0,1)?? ci.u<-mean(d)+1.96*sd(d)? ci.l<-mean(d)-1.96*sd(d)? ci.w<-ci.u-ci.l??? se=sd(d)/sqrt(n)? w.ci.uu<-ci.u+(qt(.975, df=n-1))*1.71*se? w.ci.ul<-ci.u-(qt(.975, df=n-1))*1.71*se? w.ci.upper<-...

...rkeley Breathed in his "Bloom County" comic strip ) On Mon, Sep 25, 2017 at 3:28 PM, Farnoosh Sheikhi <farnoosh_81 at yahoo.com> wrote: > Hi Bert, > > Here is the code: > f<-function(n){ > m=runif(n, 50, 100) > d=rnorm(n, 0, 1) > se=sd(d)/sqrt (n) > > ci.u<-mean(d)+1.96*sd(d) > ci.l<-mean(d)-1.96*sd(d) > ci.w<-ci.u-ci.l > > w.ci.uu<-ci.u+(qt (0.975, df=n-1)*1.71(se) > w.ci.ul<-ci.u-(qt (0.975, df=n-1)*1.71(se) > w.ci.upper<-w.ci.uu-w.ci.ul > > w.ci.lu<-ci.l+(qt (0.975, df=n-1)*1.71(se) > w.ci.ll<-c...

I took the easy route of crippling what I couldn't get working, on the grounds that partial coverage is better than none now that we have Cirrus CI checking commits on additional platforms. This series got me to a green checkmark: https://gitlab.com/ebblake/nbdkit/-/pipelines/793156983 but depends on an as-yet uncommitted patch in libvirt-ci: https://gitlab.com/libvirt/libvirt-ci/-/merge_requests/360 Eric Blake (5): ci: Expose more...

hi List and Manuel, I have encounter the following problem with the function "lineplot.CI".? I'm running R 2.10.1, sciplot 1.0-7 on Win XP.? It seems like it's a scoping issue, but I couldn't figure it out. Thanks! ...Tao > lineplot.CI(x.factor = dose, response = len, data = ToothGrowth)??? ## fine > lineplot.CI(x.factor = dose, response = len, data = ToothGr...

On 3/1/23 17:54, Eric Blake wrote: > I took the easy route of crippling what I couldn't get working, on the > grounds that partial coverage is better than none now that we have > Cirrus CI checking commits on additional platforms. > > This series got me to a green checkmark: > https://gitlab.com/ebblake/nbdkit/-/pipelines/793156983 > > but depends on an as-yet uncommitted patch in libvirt-ci: > https://gitlab.com/libvirt/libvirt-ci/-/merge_requests/360 >...

Hello! I am trying to duplicate the studentized bootstrap confidence interval calculations in boot.ci() without using boot.ci(). My w/o boot.ci() calculations are close to boot.ci() using the same object from implementation of boot(), but not quite the same (I think the differences are due to how quantiles are computed for z*). Through examining the boot.ci() code, I can see there is a call out t...

On Wed, Mar 01, 2023 at 06:43:11PM +0100, Laszlo Ersek wrote: > On 3/1/23 17:54, Eric Blake wrote: > > I took the easy route of crippling what I couldn't get working, on the > > grounds that partial coverage is better than none now that we have > > Cirrus CI checking commits on additional platforms. > > > > This series got me to a green checkmark: > > https://gitlab.com/ebblake/nbdkit/-/pipelines/793156983 > > > > but depends on an as-yet uncommitted patch in libvirt-ci: > > https://gitlab.com/libvirt/libvir...

Hello, When using bargraph.CI in package sciplot can the bars for each group be different colors? How do I select the color for each group? When I use this instead of the default (SD vs SE): bargraph.CI(x.factor = dose, response = len, data = ToothGrowth, ci.fun= function(x) c(mean(x)-sd(x), mean(x) + sd(x))...

Hi, I am trying to calculate confidence intervals using ci.numeric from epicalc package. If I generate a normal set of data and find the 99% and 95% CI, they seem too narrow to me. Am I doing something wrong?? The IQR goes from -0.62 to 0.62, so I thought the CI limits should be more extreme than these values. x<- rnorm(200,0,1) ci.numeric(x=mean(x),n...

...ls::runOnFunction(). If you browse the code a bit and backtrack the call stack to see what happens with the return value of PrintFOpt::OptimizeFixedFormatString(), you will stumble upon this segment in SimplifyLibCalls.cpp:1703ff.: // Try to optimize this call. Value *Result = LCO->OptimizeCall(CI, TD, TLI, Builder); if (Result == 0) continue; DEBUG(dbgs() << "SimplifyLibCalls simplified: " << *CI; dbgs() << " into: " << *Result << "\n"); // Something changed! Changed = true; ++NumSimplified; // Inspect the instruction aft...

Hi there, I am a beginner. I would like to change the error bars in the bargraph.CI function from the default (se) to (sd). The help file says ci.fun= function(x) c(fun(x)-se(x), fun(x)+se(x)) Is there a simple way of telling the function what (x) precisely is - I already define in in the of the bargraph.CI function and assume that is should be able to use that information. chee...

Dear All I am trying to graph a proportion and CI95% by a factor with ooplot (any other better solution ?) It works well until I try to add the confidence interval. this is the error message and and a description of the data: > dat1 PointEst TT1 1 3.6 TT2 2 5.0 TT3 3 5.8 TT4 4 11.5 TT5 5 7.5 TT5...

How can I condition any error bar function that use the arrows() function, such as 'CI.plot' (see example below) or 'error.bars', to draw only upper error bar (upper CI) if the bar value (mean) is positive and the lower error bar (lower CI) if bar value is negative? CI.plot <- function(mean, se,length, ylim=c(-5, max(CI.H)), ...) { CI.H <- mean+se CI.L <...

On Wed, Oct 14, 2020 at 10:37:26AM +0100, Richard W.M. Jones wrote: > On Tue, Oct 13, 2020 at 03:21:35PM -0600, alan somers wrote: > > Would libguestfs be willing to enable CI for the nbdkit project?  It's very > > easy to set up, at least for the Rust portion.  I'm comfortable with both > > Cirrus CI and Github's native CI.  I can get it started if you agree to do it. > > Sure. > > Be aware we may move to gitlab (same as libvirt) at...

I have a file in the following format: Scenario1 Scenario1CIL Scenario1CIU Scenario2 Scenario2CIL Scenario2CIU 60 57 62 45 48 50 110 101 111 51 50 52 120 117 122 64 62 66 192 190 194 79 75 79 where: First column = Scenario1 mean value Second column = Scenario1 Low Confidence Interval Third column = Scenario1 Upper Confidence Interval Fourth column =...

...ot;greater",conf.int=T) gives Wilcoxon rank sum test data: x and y W = 94, p-value = 0.0001624 alternative hypothesis: true mu is greater than 0 95 percent confidence interval: -3.344416 NA which is apparently contradictory (p-value is much less than 0.05 while exact CI contains 0). I examined the code and I found that while the exact wilcoxon test is applied on ranks of x (thus testing that mu = x – y is greater than 0), the exact CI is computed based on differences y - x, hence finding a CI for mu = y – x). The same applies, in my opinion, to approximative CI, w...

I'm confused about name spaces. This morning I installed the boot package because I wanted to look at bca.ci. So I did library(boot), but then I had, > bca.ci Error: Object "bca.ci" not found I had a look in the boot R directory and bca.ci was there as expected. So then I took a look at the NAMESPACE file for the boot package and saw that bca.ci wasn't exported. I tried example(boot....

There's something peculiar that I do not understand here. However, did you realize that the thing you are assigning into parts of `a' is NULL? Check you're my.test.boot.ci.1: It's NULL. Be that as it may, I get: > a <- data.frame(matrix(1:4, nrow=2), X3=NA, X4=NA) > a X1 X2 X3 X4 1 1 3 NA NA 2 2 4 NA NA > a[a$X1 == 1,]$X3 <- NULL > a X1 X2 X3 X4 1 1 3 NA 1 2 2 4 NA NA > a[a$X1 == 1,]$X4 <- NULL > a X1 X2 X3 X4 1 1 3...