Displaying 20 results from an estimated 400 matches similar to: "bit set or bit test"
2003 Feb 15
2
(no subject)
Hi,
Are there some packages which can generate multi-normal, multi-t, etc
multivariate sampling? thanks!
Best wishes,
Peng
*******************************
Peng Zhang
Department of Biostatistics
Harvard School of Public Health
655 Huntington Avenue
Boston, Massachusetts 02115
*******************************
I believe I can fly
I believe I can touch the sky
2003 Feb 19
1
How to use Cox PH model to select genes from DNA gene expression profiles?
I'm doing prediction of the survival cases using gene expression
profiles(Affymetrix chips). Can somebody tell me how to use the Cox PH
model to select genes and make a prediction of survival?
Thanks.
Guangchun
2003 Feb 27
2
epoch time conversion in R
I have a data file where each entry is indexed by the time in seconds since
epoch (e.g. 1046315697). Is there an easy way to convert this time value into
a more friendly time (such as Month-Year) when plotting it?
I searched through the manual, mailing lists, and functions like as.POSIXct
and strptime, but didn't find what I need.
Thanks,
Sharad.
2003 Feb 24
2
"trace" argument in legend() (PR#2578)
Full_Name: Jerome Asselin
Version: 1.6.2
OS: RedHat Linux 7.2
Submission from: (NULL) (142.103.173.179)
Should be an easy fix...
Consider the examble below:
plot(0,0)
legend(0,0,c("Hello!","Hi!"),pch=1:2,lty=1:2,trace=T)
It gives the following trace:
> plot(0,0)
> legend(0,0,c("Hello!","Hi!"),pch=1:2,lty=1:2,trace=T)
xchar= 0.05178 ;
2003 Feb 24
2
fill prob. in legend
Hi,
I'm trying to construct a legend which has four lines of text and associated
symbols. The first two symbols need to be normal lines which vary only in
colour. The second two symbols should have filled boxes. How do I suppress
the fill boxes in the first two lines?
J.
2003 Apr 14
7
functions in a package
Hello,
I have a beginner's Q:
if i want to know all the functions provided by a package,
what is command for that?
in another word,
Is there a command to list all the commands available in a packege?
thanks a lot,
yan
2003 May 16
2
Efficient subsetting
Hi, I'm facing this problem quite a lot, so it seems worthwhile
to check to see what the most efficient solution is.
I've two vectors x (values ordered) and y. I've ranges
x < x0, x0 <= x < x1, x1 <= x < x2, x2 <= x < x3, x > xn
and want to construct a subvector yprime of y which consists
of the first/last value of y whose x values are in the range.
For
2003 Aug 07
2
segmentation fault: formula() with long variable names (PR#3680)
R version: 1.7.1
OS: Red Hat Linux 7.2
In this example, I would expect an error for the overly long variable
name. This is always reproducable for me.
> formula(paste("y~",paste(rep("x",50000),collapse="")))
Segmentation fault
Sincerely,
Jerome Asselin
--
Jerome Asselin (Jérôme), Statistical Analyst
British Columbia Centre for Excellence in HIV/AIDS
St.
2003 Aug 08
1
covmat argument in princomp() (PR#3682)
R version: 1.7.1
OS: Red Hat Linux 7.2
When "covmat" is supplied in princomp(), the output value "center" is all
NA's, even though the input matrix was indeed centered. I haven't read
anything about this in the help file for princomp(). See code below for an
example: pc2$center is all NA's.
Jerome Asselin
x <- rnorm(6)
y <- rnorm(6)
X <- cbind(x,y)
2005 Jun 14
1
RGui crashes on wle call
Hi all --
I'm seeing the following commands reliably produce a crash in RGui,
version 2.0.1, for both my home and office machine:
> rm(list = ls(all = TRUE));
> load("dataset.R");
> library("wle");
> data.wle = wle.lm(abortion ~ year * lib.con + age + gender +
+ urbanism + census + income + church.att + children + educ +
+ religion.imp, data =
2003 Aug 30
3
fisher.test() gives wrong confidence interval (PR#4019)
The problem occurs when the sample odds ratio is Inf, such as in the
following example. Given the fact that both upper bounds of the two 95%
confidence intervals are Inf, I would have expected that the two lower
bounds be equal, but they aren't.
x <- matrix(c(9,4,0,2),2,2)
x
# [,1] [,2]
#[1,] 9 0
#[2,] 4 2
rbind("two.sided.95CI"=fisher.test(x)$conf.int,
2000 Dec 12
1
[Fwd: R code and robust regression]
-------------- next part --------------
An embedded message was scrubbed...
From: Hella Heikki <Heikki.Hella at bof.fi>
Subject: R code and robust regression
Date: Tue, 12 Dec 2000 12:12:14 +0200
Size: 2187
Url: https://stat.ethz.ch/pipermail/r-help/attachments/20001212/c3361c7d/attachment-0001.mht
2003 Aug 07
2
model.frame() call from inside a function (PR#3671)
R version: 1.7.1
OS: Red Hat Linux 7.2
Hi all,
The formula object in model.frame() is not retrieved properly when
model.frame() is called from within a function and the "subset" argument
is supplied.
foo <- function(formula,data,subset=NULL)
{
cat("\n*****Does formula[-3] == ~y ?**** TRUE *****\n")
print(formula[-3] == ~y)
cat("\n*****Result of model.frame()
2003 Apr 22
7
Subject: Eliminate repeated components from a vector
X-Mailer: VM 7.00 under 21.4 (patch 6) "Common Lisp" XEmacs Lucid
Reply-To: fjmolina at lbl.gov
FCC: /home/f/.xemacs/mail/sent
Does anyone know how I can eliminate repeated elements from a vector?
2010 Mar 12
2
find "CP"
R users
Thanks in advance:
I need to use the function "mle.cp" in R. When I run this function I recive
result <- mle.cp(y.hald~x.hald)
Error: could not find function "mle.cp"
could some one give me a help?
thanks
[[alternative HTML version deleted]]
2003 Nov 24
1
mle in the gamma model
Dear [R]-list,
I'm looking for a classic equivalent of the wle.gamma function (library
wle) that estimate robustly the shape and the scale parameters of gamma
data.
I have a vector of iid gamma rv :
>data=rgamma(100,shape=10,scale=3)
and a vector of their weights:
>weights=c(rep(.5/70,70),rep(.25/20,20),rep(.25/10,10))
and want to estimate the scale and shape of the gamma
2002 Apr 06
2
packages in OS X
=======================================================================
Simple CRAN packages which do not compile without modifications (all
others do)
=======================================================================
-- akima
/usr/bin/ld: multiple definitions of symbol _idlc_
-- fracdiff
/usr/bin/ld: multiple definitions of symbol _gammfd_
(and others)
-- odesolve
--
2005 Dec 08
1
mle.stepwise versus step/stepAIC
Hello,
I have a question pertaining to the stepwise regression which I am trying to
perform. I have a data set in which I have 14 predictor variables
accompanying my response variable. I am not sure what the difference is
between the function "mle.stepwise" found in the wle package and the
functions "step" or "stepAIC"? When would one use
2003 May 20
1
legend() with option adj=1
Hi there,
I want to justify to right the text of my legend. Consider this short
reproducable example.
x <- 1:5
y1 <- 1/x
y2 <- 2/x
plot(rep(x,2),c(y1,y2),type="n",xlab="x",ylab="y")
lines(x,y1)
lines(x,y2,lty=2)
legend(5,2,c("1,000","1,000,000"),lty=1:2,xjust=1,yjust=1)
2003 May 07
0
frailty models in survreg() -- survival package (PR#2933)
I am confused on how the log-likelihood is calculated in a parametric
survival problem with frailty. I see a contradiction in the frailty() help
file vs. the source code of frailty.gamma(), frailty.gaussian() and
frailty.t().
The function frailty.gaussian() appears to calculate the penalty as the
negative log-density of independent Gaussian variables, as one would
expect:
>