search for: cint

...m inspecting the relevant bit of code in the function t.test) df <- length(x) - 1 stderr <- sqrt(vx/nx) tstat <- (mx - mu)/stderr pval <- 2 * pt(-abs(tstat), df) alpha <- 1 - conf.level cint <- c(tstat * stderr - qt((1 - alpha/2), df) * stderr, tstat * stderr + qt((1 - alpha/2), df) * stderr) that the confidence interval reported is for the expected value of (X-mu) rather than E(X). The help file says conf.int : a confidence int...

...rs, I'd like to make a graph of coefficients of the intercept, variable 1, and variable 2 (and possibly the interaction between variable 1 and variable 2). When I use the lmList function as attached below, it shows a nice coefficient graph. > PACRP.lis <- lmList(PAffect ~ CRPC + CRPT + CINT | ID, redinteract) > coef(PACRP.lis) > PACRPlis.coef <- coef(PACRP.lis) > plot(PACRPlis.coef) However, when I'd like to make a graph using the lmer function, it shows the error message below. > PACRP <- lmer(PAffect ~ CRPC + CRPT + CINT + (1 + CRPC + CRPT + CINT |ID), redi...

..., weights, conf.level = 0.95) { require(Hmisc) nx <- length(x) df <- nx - 1 vx <- wtd.var(x, weights, normwt = TRUE) ## From Hmisc mx <- weighted.mean(x, weights) stderr <- sqrt(vx/nx) tstat <- mx/stderr ## not mx - mu alpha <- 1 - conf.level cint <- qt(1 - alpha/2, df) cint <- tstat + c(-cint, cint) cint * stderr } However, in the below extreme case where the weights are 0 for 8 of the samples and 1 otherwise we get a much tighter confidence interval rather than a wider one that I would expect. So it seems that the above impl...

Hi. I'd like to do a t-test to compare the Delta values of items with Crit=1 with Delta values of items with Crit=0. What is the t.test syntax? It should produce a result like this below (I can't get in touch with the person who originally did this for me) Welch Two Sample t-test data: t1$Delta by Crit t = -3.4105, df = 8.674, p-value = 0.008173 alternative hypothesis: true

...conf.level < 0 || conf.level > 1)) stop(paste("conf.level must be a single number", "between 0 and 1")) conf.int <- TRUE z <- atanh(r) sigma <- 1 / sqrt(n - 3) cint <- switch(alternative, less = c(-Inf, z + sigma * qnorm(conf.level)), greater = c(z - sigma * qnorm(conf.level), Inf), two.sided = z + c(-1, 1) * sigma * qnorm((1 + conf.level) / 2))...

...n example. I'd like to revise the t.test function to return the stderr value to the user. We only need to change the "rval" in the third-from-the end line of stats:::t.test.default. Change this: rval <- list(statistic = tstat, parameter = df, p.value = pval, conf.int = cint, estimate = estimate, null.value = mu, alternative = alternative, method = method, data.name = dname) class(rval) <- "htest" return(rval) To this: rval <- list(statistic = tstat, parameter = df, stderr=stderr, p.value = pval, conf.int = cint, estimate = es...

...s A/B/C can be done with: javac -cp ABC.jar X.java This makes it easy to build java scripting engines (beanshell, jython, groovy) that need to be able to compile java classes that refs existing compiled classes. For C/C++ there are 2 tools that are close to fullfilling a similar goal: gcc+gdb and cint. gdb: If an .so/.exe is compiled with gcc -g there is enough type info in the debugging symbols that gdb can later print fields of structs. Theoretically you can use gdb to ptype all the relevant types and reproduce all the headers for another pass of compilation. cint: a C/C++ interpreter that in...

dear LLVMers, I am trying to use some of the LLVM alias analyses, and I would like to check two things with you: is scev-aa being maintained in LLVM 3.7? Second question: I run cfl-aa, and I got a very small number of pointer disambiguation (no alias) with it. My results for SPEC CINT 2006 follow below. Is this low number of no alias responses something to be excepted? Below the results that I got for SPEC with cfl-aa: TotalQueries | Name | NoAliasResponses 31944 | 470.lbm | 0 49133 | 429.mcf | 42 95098 | 473.astar | 0 146301 | 462.libquantum...

Full_Name: Wei Xu Version: 1.5.1 OS: WindowsME Submission from: (NULL) (63.215.238.92) The P-value for a two.sided test is not consistent with the confidence interval. For example, P-value=0.1372, but the 95% CI doesn't include the H0 value(1). > x [1] 0.80 0.83 1.89 1.04 1.45 1.38 1.91 1.64 0.73 1.46 > y [1] 1.15 0.88 0.90 0.74 1.21 >

...structions are scheduled? For example, I have this in our llvm IR: %0 = load i32 addrspace(4)* @answer, align 4 store i32 %0, i32* @xint, align 4 %1 = load i32 addrspace(2)* @seven, align 4 store i32 %1, i32* %bint, align 4 %2 = load i32 addrspace(3)* @two, align 4 store i32 %2, i32* %cint, align 4 Loading from addrspace(4) or addrspace(3) should have a much higher latency than loading from addrspace(2) or addrspace(1). But is that information available to the instruction scheduler in our overridden getOperandLatency method? Phil -------------- next part -------------- An HTML att...

...- pbinom(x - 1, n, p)) p2 <- max(pvec[pb - pvec >= - eps * pb], 0) min(1, pb + p2) } } ) p.L <- function(x, n, alpha) { if(x == 0) 0 else qbeta(alpha, x, n - x + 1) } p.U <- function(x, n, alpha) { if(x == n) 1 else qbeta(1 - alpha, x + 1, n - x) } CINT <- switch(alternative, less = c(0, p.U(x, n, 1 - conf.level)), greater = c(p.L(x, n, 1 - conf.level), 1), two.sided = { alpha <- (1 - conf.level)/2 c(p.L(x, n, alpha), p.U(x, n, alpha)) } ) attr(CINT, "conf.level") <- conf.level ESTIMATE <- x/n names(x) &l...

...l calculate age and display output first I do know the system date (month, day, year) I find with Sys.Date () I found a function. net to calculate the age Protected Sub calcul_age() Dim datetimenow = DateTime.Now Dim datedenaissance = CDate(TbDn.Text) Dim age As Integer age = CInt(Now.Year - datedenaissance.Year) If (datedenaissance.Month > Now.Month) Then age = age - 1 End If If ((datedenaissance.Month = Now.Month) And (datedenaissance.Day > Now.Day)) Then age = age - 1 End If Me.Label3.Text = " Age : " & age & " ans&...

Hi, what would be needed to make a C/C++ interpreter using the LLVM libraries. We have in our project (http://root.cern.ch) a C/C++ interpreter (http://root.cern.ch/twiki/bin/view/ROOT/CINT), but it has some limitations (the biggest being maintenance). I see there is a libLLVMInterpreter that can interpret the LLVM IR. Could this be used to interpret, or a starting point, to an real interpreter ? Cheers, Fons. -- Org: CERN, European Laboratory for Particle Physics. Mail: 12...

Dear R users, Here is a couple a quick questions, for which I was unable to not find any answer in the list archives and in the help: 1- Is there any R equivalents of the VB functions Cint, CStr, etc... (for non VB users, these functions transform the category of a specified variable and smartly adapt the value of this variable) ? I have tried to use the as.numeric, as.factor and as.vector commands but the result is not exactly what I want ([1] 1, 3, 5, 6) >a<-as.factor(cb...

...kinit and using C. This is where I need help. Here is the main function that needs to be translated into C: function FindRootDeviceByLabel(const RootFSLabel: string): boolean; const DEV_MAJOR: dev_t = 8; DEV_MODE: mode_t = S_IFBLK or &0660; var ActualFSLabel: string; FileDescriptor: cint; ReadBytes: TsSize; DeviceFound: boolean; dev_minor: dev_t; begin Result := False; SetLength(ActualFSLabel, Length(RootFSLabel)); //loop through all devices that might have our root filesystem sda2 to sdp2 dev_minor := 2; while dev_minor < 243 do begin // last one is 242 if...

...(2) global i32 7, align 4 @xint = common global i32 0, align 4 @y = global [1 x i32*] [i32* @xint], align 8 ; Function Attrs: nounwind uwtable define i32 @main(i32 %argc, i8** %argv) #0 { entry: %argc.addr = alloca i32, align 4 %argv.addr = alloca i8**, align 8 %bint = alloca i32, align 4 %cint = alloca i32, align 4 %a = alloca i32**, align 8 %b = alloca i32*, align 8 store i32 %argc, i32* %argc.addr, align 4 store i8** %argv, i8*** %argv.addr, align 8 %0 = load i32 addrspace(4)* @answer, align 4 store i32 %0, i32* @xint, align 4 ... Currently this produces the following asse...

Estimados Al dirigir la lectura de un folder tYA1.csv me da este error y no me percato del motivo, adjunto archivo. Espero su ayuda Saludos Jos? > setwd("D:/Public/Documents/R/bioimpedancia") > a<-read.csv("tYA1.csv",header=TRUE, sep=",", dec=".") Error in read.table(file = file, header = header, sep = sep, quote = quote, : more

Hello, I'm using var.test to do a simple F-test for equality of variances. I think I'm missing something small here: > m<-rnorm(10,sd=1) > n<-rnorm(5,sd=1) > var.test(m,n) F test to compare two variances data: m and n F = 13.7438, num df = 9, denom df = 4, p-value = 0.02256 alternative hypothesis: true ratio of variances is not equal to 1 95 percent confidence

Hi, we want to implement a C++ interpreter using LLVM and clang, to replace our existing one <http://root.cern.ch/drupal/content/cint> that is used as part of a data handling and analysis environment <http://root.cern.ch>. We plan to keep the set of features that our current interpreter offers, e.g. calls into and out of libraries, dynamic scoping, unloading of code, and a prompt. I will (hopefully :-) present the plans...

...s an issue. If the project stays with E.g. bash or Python it does not matter. Any decent OS will support it unless you want to do this in an embedded environment. I personally would love to see this done through C++ scriptlets. I love, love, love cling !!! Before cling I was dearly in love with CINT. But hey, that's just me :) Thanks again for this initiative. I think it is a good thing to give some one who wants to start out. Varol On 11/29/2012 04:45 PM, Marc J. Driftmeyer wrote: > So instead of being able to quickly extend scripting via Bash or C > Shell you'd prefer peo...