Hi
[I'm not sure if this is intelligible to non-programmers or not]
R-exts (section 4.2) gives an example of the .C() function whose
third argument is
as.integer(length(a)),
and urges the user to coerce all the arguments to the correct form
(on pain of hard-to-catch errors
which I now know
On Mon, 8 Nov 2004, Robin Hankin wrote:
[I'm not sure if this is intelligible to non-programmers or not]
Not, I suspect: R-exts is not intended to be.
R-exts (section 4.2) gives an example of the .C() function whose
third argument is
as.integer(length(a)),
and urges the user to
On Sun, 2004-11-07 at 02:54, Saiwing Yeung wrote:
It seems like a number of people on this list can install rgl but have
problem loading it. I found myself in the same situation too.
I have tried the workaround of removing /usr/X11R6/lib from
DYLD_LIBRARY_PATH, but it doesn't seem to work
Hi,
I posted some days ago a question concerning the computation of lambda in the
smooth.spline function (which I repreat at the bottom of the mail) given _df_ .
Unfortunately the documentation is not clear to me. Maybee someone can help to
answer in my view the basic question:
If the
Hello R-users,
When I didn't know about the internal 'choose'
function, I made such function, 'my.choose' below. But
when I used them instead of choose(6000,20), they
didn't give me any answer.
What is the difference between 'choose', 'my.choose1',
and 'my.choose2' below? That is, what is behind
I am trying to use the Windows COM interface under R 2.0, and have
encountered the following difficulties:
- the package RDCOMClient installs, loads and works under R 1.9.1,
installs under R2.0, but does not load or work under R2.0
- the package SWinTypeLibs does not install or load under
I'm doing some work in C with the R_ExternalPointer
interface, and having some seg fault problems. I expect the
crash is my fault, bad pointer in my code causing a fault
later etc, but I'm curious about the point of failure.
R almost always falls over in a call to unzReadCurrentFile
following a
Lindveld, Charles k.lindveld at imperial.ac.uk writes:
:
: I am trying to use the Windows COM interface under R 2.0, and have
: encountered the following difficulties:
: - the package RDCOMClient installs, loads and works under R 1.9.1,
: installs under R2.0, but does not load or work under
John wrote:
Hello R-users,
When I didn't know about the internal 'choose'
function, I made such function, 'my.choose' below. But
when I used them instead of choose(6000,20), they
didn't give me any answer.
What is the difference between 'choose', 'my.choose1',
and 'my.choose2' below? That is, what
Lindveld, Charles wrote:
I am trying to use the Windows COM interface under R 2.0, and have
encountered the following difficulties:
- the package RDCOMClient installs, loads and works under R 1.9.1,
installs under R2.0, but does not load or work under R2.0
- the package SWinTypeLibs does not
Try with less ambitious numbers such as my.choose1(60,20). I think it
works fine.
I think the problem is that gamma(6001) and prod(1:6000) are so large
that it gives Inf as the answer. Hence the numerator and denominator
approaches Inf and division of two Inf gives NaN.
You could use the natural
On Mon, 8 Nov 2004, Lindveld, Charles wrote:
I am trying to use the Windows COM interface under R 2.0, and have
encountered the following difficulties:
- the package RDCOMClient installs, loads and works under R 1.9.1,
installs under R2.0, but does not load or work under R2.0
- the
On Mon, 8 Nov 2004, Mark White wrote:
I'm doing some work in C with the R_ExternalPointer
interface, and having some seg fault problems. I expect the
crash is my fault, bad pointer in my code causing a fault
later etc, but I'm curious about the point of failure.
R almost always falls over
Eryk == Eryk Wolski [EMAIL PROTECTED]
on Mon, 08 Nov 2004 10:39:26 +0100 writes:
Eryk Hi,
Eryk I posted some days ago a question concerning the
Eryk computation of lambda in the smooth.spline function
Eryk (which I repreat at the bottom of the mail) given _df_.
I had
On Mon, 8 Nov 2004 11:58:04 + (GMT), John [EMAIL PROTECTED]
wrote :
What is the difference between 'choose', 'my.choose1',
and 'my.choose2' below? That is, what is behind
'choose' function and what's the problem using 'prod'
or 'gamma' function?
Thanks a lot.
John
##
I am just realized that survival has the facility to do survival time
splitting survSplit
after read some postings about time dependency in the list.
Is it survSplit only for the survival data input (time,status) and not
for the 'counting process' input (start,stop,status)?
I take one
Prof Brian Ripley [EMAIL PROTECTED] writes:
On Mon, 8 Nov 2004, Mark White wrote:
I'm doing some work in C with the R_ExternalPointer
interface, and having some seg fault problems. I expect the
crash is my fault, bad pointer in my code causing a fault
later etc, but I'm curious about
On 08-Nov-04 John wrote:
Hello R-users,
When I didn't know about the internal 'choose'
function, I made such function, 'my.choose' below. But
when I used them instead of choose(6000,20), they
didn't give me any answer.
What is the difference between 'choose', 'my.choose1',
and
Hi there,
I'm trying to fit a growth curve to some data and need to use a weighted least
squares estimator to account for heteroscedasticity in the data. A weights
argument is available in nls that would appear to be appropriate for this
purpose, but it is listed as 'not yet implemented'. Is
On 08-Nov-04 Ted Harding wrote:
[...]
choose(600,31)
[1] 1.612899e+300
[...]
my.choose3-function(x,y){
if((x==y)||(y==0)) return(1);
m - min(y,x-y)
prod(seq(x,(x-m+1),by=-1)/(seq(m,1,by=-1)))
}
when
my.choose3(600,31)
[1] 1.613121e+300
Hi Robert
You can try gnls() in the nlme package or if that is not suitable, try
to
follow the example in the nls() help page for weighted regression that has
an example of weighted regression from MASS.
from the nls help page:
## weighted nonlinear regression
Treated -
On Mon, 8 Nov 2004, Danardono wrote:
I am just realized that survival has the facility to do survival time
splitting survSplit
after read some postings about time dependency in the list.
Is it survSplit only for the survival data input (time,status) and not for
the 'counting process' input
This question has been asked before on the list, but I'm not sure if the
answer were posted. Basically, the trick is to write the formula a bit
differently in nls() so that it does weighted least squares. nls() tries to
minimize the sum of squared differences between the two sides of ~. If you
On 08-Nov-04 Ted Harding wrote:
On 08-Nov-04 Ted Harding wrote:
[...]
choose(600,31)
[1] 1.612899e+300
[...]
my.choose3-function(x,y){
if((x==y)||(y==0)) return(1);
m - min(y,x-y)
prod(seq(x,(x-m+1),by=-1)/(seq(m,1,by=-1)))
}
when
I am trying to write a function that will run a linear
model and plot the regression coeficients with their
corresponding means. I am having two problems. I can
get the plot with the function below, but I am having
trouble labeling the points.
function(y,x1,x2,x3,x4){
outlm-lm(y~x1+x2+x3+x4)
In addition to John Fox's recommendations, some further hints and resources
are:
The tools of TCL/TK allow one to develop a comprehensive GUI. Note that the
Bioconductor project, http://www.bioconductor.org/ , has some packages
(tkWidgets, widgetTools are at least two that I know of) that claim
Thanks to both Duncan Murdoch and Brian Ripley. LEM indeed
is the program that Brian referred to in his response, a
program for latent class/trait loglinear and event history
analysis compiled with Borland Pascal. There's a dos
version and a windows version. Not sure the program will
Dear R users,
I have a data frame which I create with read.csv and then order by
date:
d - na.omit(read.csv(...))
d - d[order(as.Date(as.character(d$Date), format=%d-%b-%y),
decreasing=F, na.last=F),]
My problem is that even though the data frame is ordered as
requested, the old row
Dear R users:
I'm generating the following survival data:
set.seed(123)
n=200 #sample size
x=rbinom(n,size=1,prob=.5) #binomial treatment
v=rgamma(n,shape=1,scale=1) #gamma frailty
w=rweibull(n,shape=1,scale=1) #Weibull deviates
b=-log(2)
On Mon, 8 Nov 2004, bogdan romocea wrote:
My problem is that even though the data frame is ordered as
requested, the old row numbers are preserved. For example:
* Before sorting:
d[1:3,]
Date Amt
1 5-Nov-04 87.07
2 4-Nov-04 85.80
3 3-Nov-04 82.90
* After sorting:
d[1:3,]
Date Amt
Dear All,
I am trying to install Rmpi on AMD 64 with Fedora Core 2 for x86_64, but I
receive the following error:
/usr/bin/ld:
/usr/lib/gcc-lib/x86_64-redhat-linux/3.3.3/../../../../lib64/libmpi.a(abort.o):
relocation R_X86_64_32 can not be used when making a shared object;
recompile with -fPIC
It's your misinterpretion that is misleading, not the output.
Data frames have row *names* and not *numbers*.
On Mon, 8 Nov 2004, bogdan romocea wrote:
Dear R users,
I have a data frame which I create with read.csv and then order by
date:
d - na.omit(read.csv(...))
d -
You wrote:
Is there a way to update the row numbers as well? It's not that
important, but I find it a bit confusing.
They're not actually row numbers, they're row ***names***.
These default to row numbers.
If they were real-live names you'd want them to be
On Mon, 8 Nov 2004, Krzysztof Kolanek wrote:
Dear All,
I am trying to install Rmpi on AMD 64 with Fedora Core 2 for x86_64, but I
receive the following error:
/usr/bin/ld:
/usr/lib/gcc-lib/x86_64-redhat-linux/3.3.3/../../../../lib64/libmpi.a(abort.o):
relocation R_X86_64_32 can not be
Two groups with a difference in mean of 20.6 have a between-group variance
of 212.18 (as Minitab was given), not 424.36.
var(c(0, 20.6))
[1] 212.18
(and the example on the help page shows this *is* what is meant by the
variance).
On Mon, 8 Nov 2004, Peter Levy wrote:
I think there is a bug
Hi,
I want to insert a key into each panel of a trellis plot, which I can do
with a custom panel function that calles draw.key. The problem arises
because I want the top right hand corner of the key to start in the top
right hand corner of the panel. If you run my code below, you can see that
Hi,
In my previous posting I forgot my system information. Sorry. It is listed
below.
I want to insert a key into each panel of a trellis plot, which I can do
with a custom panel function that calles draw.key. The problem arises
because I want the top right hand corner of the key to start
How about something like this?
my.func - function(y, x1, x2, x3, x4 = NULL){
my.formula - as.formula(y ~ x1 + x2 + x3 + x4)
if(is.null(x4)) { my.formula - as.formula(y ~ x1 + x2 + x3) }
outlm - lm(my.formula)
meanvec-c(mean(x1),mean(x2),mean(x3))
if(is.null(x4) == F) {
Hi,
I've searched the archives for a discussion on the use of R for developing
small world or scale-free network models but have been unable to dig
anything up. Is anyone working on these types of models using R, and if so,
have you found this language amenable to their development?
thanks,
Hello,
I have the following problem:
test is a data frame with 9 fields. The field test$Date is factorized with
dates. The format is dd-mm- (using Oracle notation). I want to convert this
to Date in '%Y-%m-%d format.
What I am doing is:
for (i in 1:nrow(test))
{
On Mon, 08 Nov 2004 15:56:49 - (GMT), (Ted Harding)
[EMAIL PROTECTED] wrote :
Further investigation indicates that there is an integer truncation
problem with 'choose':
print(choose(600,1),digits=20)
[1] 60001679906.0
print(my.choose3(600,1),digits=20)
[1] 6e+10
Hi there. OSX/R2.0
We are trying to implement a bootstrap of the coeffecients of a mixed
effect model. In particular, we are interested in the intercept and
slope of the random effects.
Following from the basics for a linear model, we construct our lme
models and a boot function:
On 08-Nov-04 Duncan Murdoch wrote:
On Mon, 08 Nov 2004 15:56:49 - (GMT), (Ted Harding)
[EMAIL PROTECTED] wrote :
Further investigation indicates that there is an integer truncation
problem with 'choose':
print(choose(600,1),digits=20)
[1] 60001679906.0
Hi
Here's a modification of your panel function that I think does what you
want (comments embedded):
tmp.xyplot - function(x,y,subscripts=subscripts,cdata=cdata){
# plot data points
panel.xyplot(x,y)
# extract parameter values
right - as.character(cdata[subscripts,][1,c(4,5,6)])
Dear Colleagues,
Does anyone know if the data sets used in the book Graphical Methods
for Data Analysis by Chambers, Cleveland, Kleiner and Tukey are
available on the web?
[[alternative HTML version deleted]]
__
[EMAIL PROTECTED] mailing
Hi
Steven Lacey wrote:
In my previous posting I forgot my system information. Sorry. It is listed
below.
snip
platform i386-pc-mingw32
arch i386
os mingw32
system i386, mingw32
status
major1
minor8.1
year
You want as.Date, not strptime that gives you a *list*.
Please do read the help pages before posting, as the the posting guide
asks you to.
On Mon, 8 Nov 2004, João Mendes Moreira wrote:
Hello,
I have the following problem:
test is a data frame with 9 fields. The field test$Date is
Where would I find a list of valid options for \docType{}?
I found the following types using grep:
- class
- data
- genericFunction
- methods
Are there others?
--
SIGSIG -- signature too long (core dumped)
On Mon, 8 Nov 2004, Roberto Perrelli wrote:
Dear R users:
I'm generating the following survival data:
set.seed(123)
n=200 #sample size
x=rbinom(n,size=1,prob=.5) #binomial treatment
v=rgamma(n,shape=1,scale=1) #gamma frailty
w=rweibull(n,shape=1,scale=1)
I'd like to insert (for example) the current datetime into a comment so
it goes into the history.
I can of course cut and paste the results of date() into a comment line,
but it would be easier and more powerful to be able to type something
like
hstamp()
and have it go into the history.
Ted == Ted Harding [EMAIL PROTECTED]
on Mon, 08 Nov 2004 19:10:55 - (GMT) writes:
Ted On 08-Nov-04 Duncan Murdoch wrote:
On Mon, 08 Nov 2004 15:56:49 - (GMT), (Ted Harding)
[EMAIL PROTECTED] wrote :
Further investigation indicates that there is an integer
Have you perchance looked at
?paste
For example:
paste(The current time is , date())
[1] The current time is Mon Nov 08 17:06:25 2004
[EMAIL PROTECTED]
Sent by: [EMAIL PROTECTED]
11/08/2004 04:47 PM
To: [EMAIL PROTECTED]
cc:
Subject:[R] can one
But that doesn't put the result into the history buffer, to be written
to a file only later when I savehistory(filename).
Bert Gunter also suggested ?capture.output and ?textConnection,
but I cannot see how to get text into the history buffer as comments,
but with evaluated expressions (values).
Hello,
I have a large data frame and am aiming to create a summary data frame in
order to plot quarterly means by age. However, it seems as though I am
taking a very long winded approach to this - can anybody point me in the
direction of something neater?
#original table is dat2,
dat5- tapply
On Mon, 8 Nov 2004 [EMAIL PROTECTED] wrote:
But that doesn't put the result into the history buffer, to be written
to a file only later when I savehistory(filename).
Bert Gunter also suggested ?capture.output and ?textConnection,
but I cannot see how to get text into the history buffer as
(Ted Harding) [EMAIL PROTECTED] writes:
On 08-Nov-04 Duncan Murdoch wrote:
Are you sure you're using the standard R choose()? I get different
results from you:
print(choose(600,1),digits=20)
[1] 6e+10
...
Well, that's good to learn! As usual, I'm out of date, but this
time
--- mattf [EMAIL PROTECTED] wrote:
Hi,
I've searched the archives for a discussion on the use of R for
developing
small world or scale-free network models but have been unable to
dig
anything up. Is anyone working on these types of models using R, and
if so,
have you found this
On Mon, 8 Nov 2004 [EMAIL PROTECTED] wrote:
Is there an add-on package in R for QTL interval mapping for outbred
population, eg. Haley-Knott regression method ?
Stella
I believe Karl Broman's R/QTL
(http://biosun01.biostat.jhsph.edu/~kbroman/software/)
and Brian Yandell's bim
Hi,
I have a data.frame with several variables and 50,000 observations.
i.e.
data[1:2,1:7]
Iteration Day Production.Type tsUSusc tsASusc tsULat tsALat
1 0 Generic 17965 8833053 0 0
1 1 Generic 17965 8833053 0 0
.
.
Hello sir:
I am a data analysist of a bio-company.
Here's a question about microarray for identifying the differentially expressed
genes.Thanks for your help.
First,I wanna make sure whether I've grasped the difference between
theunadjusted p-value and adjusted p-value.
1.Unadjusted
F Z gerifalte28 at hotmail.com writes:
:
: Hi,
:
: I have a data.frame with several variables and 50,000 observations.
: i.e.
: data[1:2,1:7]
: Iteration Day Production.Type tsUSusc tsASusc tsULat tsALat
: 1 0 Generic 17965 8833053 0 0
: 1 1
Yes, I am a newbie at R, but it is not the complex commands in R that
have me baffled, but simple data commands. For example, why does
something like:
plot(Girth ~ Height)
*not* work after a command that allegedly loads the data:
data(trees)
with the error message:
Error in eval(expr, envir,
On Mon, 2004-11-08 at 21:56, Dr. John R. Vokey wrote:
Yes, I am a newbie at R, but it is not the complex commands in R that
have me baffled, but simple data commands. For example, why does
something like:
plot(Girth ~ Height)
*not* work after a command that allegedly loads the data:
Hi, I have a number of SPLUS scripts that I wish to run in R. Can you point
me to the right FAQ or help area to do this?
Cheers
Andy
Andy Donaldson
Head Of Market Risk
Institutional Markets Services
Bank of New Zealand
Phone: +64 4 474 6937
Mob No: +64 29 222 0050
CAUTION - This message
Jari, thanks so much for the advice. I followed your instruction and could
install rgl too. It's really fun playing with it.
I have written up a more detailed instruction, hopefully this will help less
experienced R users.
1) go to apple.com, download Developer Tools and X11. You need to sign
hi all
does anyone know of a package that does neural networks?
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
On Tue, 9 Nov 2004, UCT Staff Member - Allan wrote:
does anyone know of a package that does neural networks?
Yes, very many people know.
Do please read the posting guide and do your homework. Try
help.search(neural net) for example. Or, there is a list of packages in
the R FAQ so why not
[EMAIL PROTECTED] wrote:
Hi, I have a number of SPLUS scripts that I wish to run in R. Can you point
me to the right FAQ or help area to do this?
What about reading the FAQ and the manual An Introduction to R
yourself? If you don't want to read this minimal amount of documentation
before asking
68 matches
Mail list logo