Displaying 2 results from an estimated 2 matches for "nonscalar".
2007 Sep 10
2
[Fwd: buglet (?) in de.restore()]
I'm resending this after a decent interval of 20 days -- any
opinions? Should I file it as a bug report? Is it my mistake?
cheers
Ben Bolker
-------- Original Message --------
Subject: buglet (?) in de.restore()
Date: Tue, 21 Aug 2007 13:29:33 -0400
From: Ben Bolker <bolker at zoo.ufl.edu>
To: r-devel at r-project.org
If one calls data.entry() with a matrix:
A =
2009 Dec 08
4
lower.tail option in pnorm
Hi,
I would have thought that these two constructions would
produce the same result but they do not.
Resp <- rbinom(10, 1, 0.5)
Stim <- rep(0:1, 5)
mm <- model.matrix(~ Stim)
Xb <- mm %*% c(0, 1)
ifelse(Resp, log(pnorm(Xb)), log(1 - pnorm(Xb)))
pnorm(as.vector(Xb), lower.tail = Resp, log.p = TRUE)
> ifelse(Resp, log(pnorm(Xb)), log(1 - pnorm(Xb)))
[1] -0.6931472 -1.8410216