Hi,
Maybe I'm not understanding you correctly, but wouldn't the following
be the fastest way to compute the sum of a column:
x - matrix(c(1:4462), nrow=4462, ncol=1) #setting up a column of
4462 elements to be able to make a comparison
sum(x[,1])
[1] 9956953
system.time(sum(x[,1]))
[1] 0 0 0
Robert,
Robert Keefe wrote:
Hello All,
Specifying 'type = response' when using predict() on a
model fit using glm(...,family=binomial) returns fitted
probabilities.
Is it possible to get the same from a model object
fit using glmmPQL() ?
glmmPQL returns an lmeObject which has no family
this is what my program read
dataa-data.frame(matrix(0, nrow=nrow(data), ncol=ncol(data))
i-1
j-1
while(i=nrow(data))
{
if(data$Index.Price==1)
data1[j,]-data[i,]
else
{
num-data$No.Primary[i]
data[j,]-apply(data[i:i+num-1,],2,sum)
i-i+num-1
}
j-j+1
On Sat, 2003-06-14 at 17:46, Peter Dalgaard BSA wrote:
Actually, the problem of making sure that N computers all have the
same set of RPMs installed is not terribly scalable either...
Centralized software maintenance is a pretty obvious thing to do, but
I often suspect RedHat (and several
On Sat, 14 Jun 2003 09:33:47 +0700
Philippe Glaziou [EMAIL PROTECTED] wrote:
Frank E Harrell Jr [EMAIL PROTECTED] wrote:
I tried this on the latest version of Hmisc (1.6-0):
library(Hmisc)
set.seed(1)
y - factor(sample(c('a','b','c'),100,T))
x - runif(100)
a - summary(y ~ x,
Thanks to Martin, Marc, Dennis, Rashid and Bill for comment.
The short answer is there is no factorial function within the base
package, so use the result that gamm(x + 1) = x!. This what the
factorial() function in the package gregmisc does:
factorial - function (x) gamma(1 + x)
Bill V.
Dear R-announce list members,
I've uploaded to CRAN a new version of my Rcmdr package, which provides a
cross-platform basic-statistics GUI for R. This version has a few new
features, but mostly addresses some problems that arose with the original
version. The relevant portion of the CHANGES
Frank E Harrell Jr [EMAIL PROTECTED] wrote:
Sorry to bother you again, I use the latex command with
Sweave, and I would like to see the reports come out well
without manual interventions on some tex files.
latex(object, file='') works with Sweave's results=tex= construct.
As in my