Displaying 7 results from an estimated 7 matches for "prob2".
Did you mean:
probe
2010 Mar 06
1
Plotting Comparisons with Missing Data
...file containing some simple information about the performance of 4
algorithms. The columns are the name of the algorithm, the problem instance
and the resulting score on that problem (if it wasn't solved I mark that
with NA).
solver instance result
A prob1 40
B prob1 NA
C prob1 39
D prob1 35
A prob2 100
B prob2 50
C prob2 NA
D prob2 NA
A prob3 75
B prob3 80
C prob3 60
D prob3 70
A prob4 80
B prob4 NA
C prob4 85
D prob4 75
I've managed to read in the data as follows:
data <- table.read("./test.txt", header = TRUE, colClasses =
c("factor","factor","nume...
2010 Sep 06
1
calculating area between plot lines
Hi everyone. I have these data:
probClass<-seq(0,0.9,0.1)
prob1<-c(0.0070,0.0911,0.1973,0.2949,0.3936,0.5030,0.5985,0.6869,0.7820,0.8822)
prob2<-c(0.0066,0.0791,0.2358,0.3478,0.3714,0.3860,0.6667,0.6400,0.7000,1.0000)
# which I'm plotting as follows:
plot(probClass,prob1,xlim=c(0,1),ylim=c(0,1),xaxs='i',yaxs='i',type="n")
lines(probClass,prob1)
lines(probClass,prob2)
polygon(c(probClass,rev(probClass)),c(...
2011 Nov 29
0
[SOLVED]looking for beta parameters
...> flag = 1
> }
> return(m0)
> }
> p1 = quantile1$p
> x1 = quantile1$x
> p2 = quantile2$p
> x2 = quantile2$x
> logK = seq(-3, 8, length = 100)
> K = exp(logK)
> m = sapply(K, betaprior1, x1, p1)
> prob2 = pbeta(x2, K * m, K * (1 - m))
> ind = ((prob2 > 0) & (prob2 < 1))
> app = approx(prob2[ind], logK[ind], p2)
> K0 = exp(app$y)
> m0 = betaprior1(K0, x1, p1)
> return(round(K0 * c(m0, (1 - m0)), 2))
> }
>
> I assume one could change this code to...
2009 Apr 21
4
My surprising experience in trying out REvolution's R
...chol(var.root) );
u <- matrix( rnorm(n.random*n.simu), nrow=n.random, ncol=n.simu );
log.prob1 <- -colSums(u*u)/2;
u <- u.mean.initial + var.root %*% u;
ada <- exp( ada.part + z %*% u );
ada <- ada/(1+ada); #it is probability now
log.prob2 <- colSums( y*log(ada) + (1-y)*log(1-ada) );
log.prob2 <- -quadratic.form2(D.u.inv, u)/2 + log.prob2;
weight <- exp(log.prob2 - log.prob1);
weight <- weight/sum(weight);
ada <- t(ada);
pi <- colSums(ada*weight);
product <- colSums(...
2012 Apr 04
0
multivariate ordered probit regression---use standard bivariate normal distribution?
...verse = TRUE), data = refdata)
Coefficients:
(Intercept):1 (Intercept):2 age:1 age:2
-1.65895567 -2.14755951 0.06688242 0.04055919
Degrees of Freedom: 1492 Total; 1488 Residual
Residual Deviance: 1188.909
Log-likelihood: -594.4543
##########################################################
unvar.prob2<-vglm(auric4~age,data=refdata,cumulative(link="probit",parallel=FALSE,reverse=TRUE))
> unvar.prob2
Call:
vglm(formula = auric4 ~ age, family = cumulative(link = "probit",
parallel = FALSE, reverse = TRUE), data = refdata)
Coefficients:
(Intercept):1 (Intercept):2 (Interce...
2002 Aug 06
3
hard to believe speed difference
...public normal_univariate_truncated(double mean, double var,
double lower, double upper)
{
double sd = Math.sqrt(var);
lower -= .5;
upper += .5;
double left = (lower - mean)/sd;
double right = (upper - mean)/sd;
double prob2 = Statistics.normalCdf(right) -
Statistics.normalCdf(left);
while(true)
{
harvest = mean + ran_obj.nextGaussian()*sd;
if(harvest>lower && harvest<upper) break;
}
harvest = Sfun.near...
2019 Jun 24
1
Calculation of e^{z^2/2} for a normal deviate z
>>>>> jing hua zhao
>>>>> on Mon, 24 Jun 2019 08:51:43 +0000 writes:
> Hi All,
> Thanks for all your comments which allows me to appreciate more of these in Python and R.
> I just came across the matrixStats package,
> ## EXAMPLE #1
> lx <- c(1000.01, 1000.02)
> y0 <- log(sum(exp(lx)))
> print(y0) ## Inf