Some time ago I was allowed to discuss "Diamond Graphs", and whether
they would be useful in R, in this mailing list.
The August 2003 issue of The American Statistician has finally arrived
here and I have been able to read the article. A number of points of
interest arise.
1. The article is
Oh, disregard the second argument to cv.tree, it's one of my
gazillion attempts to make the .lost work by saying silly
things a la local=list(X=X) or some such, all to no avail...
Cheers,
Alexy
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.eth
Greetings -- I'm programming a data mining system
in R for protein structural data. As a seasoned
Perl and Python and Ada and ML, et al., programmer,
I am severely befuddled by the environment problem,
where data is not found in a 3rd level nested
function. I did peruse the parent frame not on t
Zitan Broth wrote:
That sounds cool, it also avoids the file IO. I'm an R newbie, but could I
achieve something similar with pure PHP.
I don't know of a way to achieve it with pure PHP short of writing your
own C extension to PHP. I had considered that at one point, but have
found it unnecessary
Hi all,
I have an experiment where 5 raters assessed the quality of 24 web sites. (each rater
rated each site once). I want to come up with a measure of reliability of the ratings
for the web sites ie to what extent does each rater give the same (or similar) rating
to each web site. My idea was
Hi Joe,
Thanks for your message,
> I'm doing something similar using PL/R (an R procedural language handler
> extension to Postgres that I wrote) with Postgres, R, and PHP. In
> Postgres 7.4 (currently at beta3) or with a back-patched copy of 7.3,
> you can preload the R interpreter when the Post
Opps I meant "incur a couple of seconds overhead for starting R" rather than
file IO.
It would be nice to be able to preload the R interpreter like Joe is doing
with Postgres thereby running it as a daemon .
Z.
- Original Message -
> Yeah I agree with what you are saying here for sur
On Wed, 2003-09-24 at 18:36, Liaw, Andy wrote:
> There was a related thread on R-help, probably last year. The question was
> getting density() to numerically integrate to 1. The answer is, "yes". If
> you do fine enough partitions, you will see that it integrates to one. And
> yes, a kernel de
Hi Paul,
Thanks greatly for your message.
- Original Message -
> Below is the test I ran awhile back on invoking R as a system call. It
> might be faster if you had a c-extension to R but before I went that route
I
> would want to know 1) roughly how fast Python and Perl are in returning
There was a related thread on R-help, probably last year. The question was
getting density() to numerically integrate to 1. The answer is, "yes". If
you do fine enough partitions, you will see that it integrates to one. And
yes, a kernel density estimate is theoretically a true density (assumin
Andy -
help("assign") says:
"Value: This function is invoked for its side effect, which is
assigning `value' to the variable `x'. ..."
Gosh. The help page isn't very specific about what the return
value of assign() IS, but it's not a named object. Your basic
strategy of assigning, writin
On Thu, 2003-09-25 at 11:20, Fernando Henrique Ferraz wrote:
>Hi, I'm having a little issue with \Sexpr{bla} relating to the number of digits
> it is using to print its output.
Try \Sexpr{format(bla,digits=3)}
Cheers
Jason
--
Indigo Industrial Controls Ltd.
http://www.indigoindustrial.co.n
Hi, I'm having a little issue with \Sexpr{bla} relating to the number of digits it
is using to print its output. I understand that what \Sexpr{bla} does is parsing 'bla'
inside R and then returns it coerced into a character string. The problem I'm having
is that I'm losing control over the n
I have an object which I want to save into an external file, using save, on a time
step in loop. In the loop I have an if statment that will save when the appropraite
time, t, is reached.
I want to change the object and the name of the saved object to reflect which time
step is being written. I
Visual inspection of the plot of a density() function vs a normal with
the same mean and variance suggests the area under the density curve is
bigger than under the normal curve. The two curves are very close over
most of the domain. Assuming the normal curve does integrate to 1, this
implies the
Hi,
Has anyone seen this problem before, and does anyone have a solution?
When I query a database using DBI and ROracle to retrieve a date in a format
other than the system default, I get a core dump.
System: Redhat Linux 7.3
uname: Linux 2.4.20-20.7smp #1 SMP Mon Aug 18 14:46:14 EDT 2003 i686
Greetings, Fellow R-ians:
I'm working with a few different quasi-likelihood formulations for some data
I'm analyzing. I'd like to implement the Information Matrix Test (see,
e.g., White, 1982, or Lancaster, 1984) for each of them to determine which
of the models is more likely. Since the null di
Jason Turner <[EMAIL PROTECTED]> writes:
> On Thu, 2003-09-25 at 05:34, Christoph Lehmann wrote:
> > Hi
> >
> > I have calculated lots (>1000) of linear models and would like to store
> > each single result as an element of a 3D matrix or a similar storage:
> > something like
> >
> > glm[i][j][k
Hi,
That would be great! I will wait.
Thanks.
Weiming
On Wed, 2003-09-24 at 13:59, Duncan Temple Lang wrote:
>
> Hi.
>
> Over the past few weeks, I have spent some time enhancing the
> SJava package and ironing out some issues (syncrhonization
> and locking of the R evaluator, etc.). Along th
Hi.
Over the past few weeks, I have spent some time enhancing the
SJava package and ironing out some issues (syncrhonization
and locking of the R evaluator, etc.). Along the way,
I did encounter the problem you mention about java.lang.Hashtable.
I can't quite remember what the problem was, but
Dear all,
I have a question on the dose-response estimation with clustered/ correlated binary
data.
I would like to estimate the hit rate for a certain test at various concentration
levels. The test is used
on 5 subjects, and each subject is tested 20 times. If we assume that the 100 samples
ar
Hi, Henrik,
Thank you very much for replying.
I actually have the JAVA_HOME environment variable set in a linux shell
script .bash_profile which set up all of my development environment
before I install this SJava package and it points to the current java
cdk I am using. So it is weird that R cou
Hi, I do not have Linux myself but looking at your error message
> > .JavaInit()
> [1] error initializing manager class can't find class
> java/lang/Hashtable Error in .JavaInit() : Couldn't start
> Java Virtual Machine: can't find class java/lang/Hashtable
> >
it looks like you have not set t
On Wed, 2003-09-24 at 19:31, Vito Muggeo wrote:
> Dear all,
> I'm interested in fitting time-series linear models with I(1) errors. Namely
> given
> y_t=a+b*t+u_t
> the random term u_t are such that
> u_t-u_{t-1}=e_t~iid N(0,\sigma)
>
library(nlme)
help(lme) #note the optional correlation argume
As you noted, the output of, e.g., "lm" is a list. When I've done many
fits like this, I've usely been interested in only a few components of
the list. In that case, I first construct a matrix or higher
dimensional array and then store what I want in that array. Example:
> DF0 <- data.frame
On Thu, 2003-09-25 at 05:34, Christoph Lehmann wrote:
> Hi
>
> I have calculated lots (>1000) of linear models and would like to store
> each single result as an element of a 3D matrix or a similar storage:
> something like
>
> glm[i][j][k] = lm(...)
>
> Since I read that results are lists: Is i
Hi
I have calculated lots (>1000) of linear models and would like to store
each single result as an element of a 3D matrix or a similar storage:
something like
glm[i][j][k] = lm(...)
Since I read that results are lists: Is it possible to define a matrix
of type list?
Or what do you recommend?
Hi,
I installed SJava 0.66-1 on my linux RH 7.2 machine. I am using R-1.7.1.
I set LD_LIBRARY_PATH up then called library(SJava) and it was fine. But
when I called .JavaInit(), I got weird error that states it could not
find java/lang/Hashtable. The session is below and shows my java
configuratio
Hi All:
I am clustering 500 genes using hclust of R. Visualizing cluster
membership becomes difficult with so many genes in each cluster...Is there
a way of printing the dendrogram in multiple pages so that I can clearly
see what is in each cluster?
Thanks in advance.
Karthi.
_
You need to declare the drand48 function within your C source file.
Either
#include
or
extern double drand48(void);
[EMAIL PROTECTED] writes:
> Hello
>
> I'm trying to use C random generator function drand48() which return
> floating-point values, uniformly distributed over [0,1],
> (http:
Hi,
running gc() helped me a lot
Regards,
Vladimir
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of
> ZABALZA-MEZGHANI Isabelle
> Sent: Wednesday, September 24, 2003 7:27 PM
> To: help R (E-mail)
> Subject: [R] Problem with memory for large datasets
>
What have you done to remove "temporary objects"? Does this include
"rm(list='ls()')" or "remove(list=objects())"?
hope this helps. spencer graves
ZABALZA-MEZGHANI Isabelle wrote:
Hello,
I would like to know if there is a possibility to "clean" the R memory
during a R session. In fact, I r
ZABALZA-MEZGHANI Isabelle wrote:
Hello,
I would like to know if there is a possibility to "clean" the R memory
during a R session. In fact, I realize a lot of instruction with large
objects (matrix of 500*5000), and I can not manage to achieve the end of my
script due to memory lack. Of course, I
Hello,
I would like to know if there is a possibility to "clean" the R memory
during a R session. In fact, I realize a lot of instruction with large
objects (matrix of 500*5000), and I can not manage to achieve the end of my
script due to memory lack. Of course, I've tried to remove all "temporary
Hi,
I'm ploting a multivariate time series but I want to change the y
labels. I've tried with ylab but I get an error
> plot(hke.ts,ylab=letters[1:5])
Error in plot.ts(x[, i], axes = FALSE, xlab = "", ylab = "", log = log,
:
formal argument "ylab" matched by multiple actual arguments
I'
Hi,
this is clearly not a problem of ROracle but very likely one of the
client configuration/connection stuff.
Is TNS on the client configured correctly (looks like the client cannot
map the server and, hence, does not find the listener; very common
Oracle mistake)?
Regards
Michael
Joerg Schabe
ZABALZA-MEZGHANI Isabelle wrote:
Hello,
I would like to know where I can find a Windows binary version of R 1.6.2,
since I want to test the last SJava library (which has been built under
1.6.2) which does not work on R 1.7.1
I've already search on the CRAN web site but on the r1070 page, if I cli
Hi,
In lattice, is it possible to get qqplot for t
distribution in different panels with different dfs.
For example,
# in panel 1
qqmath(~x, distribution=function(p) qt(p, df=3))
# in panel 2
qqmath(~x, distribution=function(p) qt(p, df=5))
etc...
Thanks,
Mahbub.
Hi,
I am running the 64bit version of R on a AIX RS6000 cluster on which is
also running the ' SQL*Plus: Release 9.2.0.2.0' in the
'ORACLE-environment Version 64 bit AIX 5.1'.
I have problems setting up an Oralce database connection.
I can instantiate a driver using drv <- Oracle(). But when I
This is sort of like calling up an auto mechanic and ask "My car just won't
start. Can you fix it?" The answer is almost surely "yes", but not without
looking at the car (or at least get answers on quite a few questions)!
Showing us the code would be a start.
Since you can use gdb to debug the c
Maybe like this:
X.ts <- matrix(rnorm(4000),1000,4)
par(mfrow=c(2,2))
apply(X.ts, 2, plot.ts, plot.type=c("single"))
or
apply(X.ts, 2, ts.plot)
HTH
Thomas
>Hi,
>How can I force a plot.ts to draw a 2x2 plot matrix instead of 4x1 ?
>Regards
>EJ
_
Hello,
I would like to know where I can find a Windows binary version of R 1.6.2,
since I want to test the last SJava library (which has been built under
1.6.2) which does not work on R 1.7.1
I've already search on the CRAN web site but on the r1070 page, if I click
on 'old', I get an error.
Than
[EMAIL PROTECTED] writes:
> This is a relatively recent article that is somewhat accessible.
> Jensen, D. R., and Solomon, Herbert (1994), "Approximations to joint
> distributions of definite quadratic forms", Journal of the American
> Statistical Association, 89 , 480-486
> It has references to p
Hi,
How can I force a plot.ts to draw a 2x2 plot matrix instead of 4x1 ?
Regards
EJ
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[EMAIL PROTECTED] writes:
> Dear Peter,
>
> I don't know if this is the proper way to ask for help installing R, but if
> not, I presume you can pass this on to the appropriate place.
>
> I'm trying to install the latest binary on Redhat 9, and keep getting the same
> error messag
> Is there an alternative way of "dispatching" the printing, such that
> the usual print method for loglm is used after doing what is special
> for hllm?
You might want to have a look at `NextMethod()'
Best,
David
__
[EMAIL PROTECTED] mailing list
https:
Hello !
I try to plot some data with lattice function xyplot. These data are
separated in two groups by a factor. If I let the default color and symbol
settings, I get the right legent with auto.key. But if I change them by col
and pch arguments, the plot is right but the legend still reflects
Obviously, when you call the function "trellis.device(postscript, ..)"
you have set current output device to postscript, look at the file *.ps !
Try the same call with "windows" at the place of "postscript"!
Best!
A.S.
Alessandro Semeria
Models and Simulations
Hi all,
I have currently been using the book Modern Applied Statistics with S from
Venables and Ripley. At chapter 6 on Linear Statistical models I wanted to
produce the plot as shown by "Figure 6.1" using Whiteside's data. xyplot
command seems not to work on my version of "R"(version 1.7.0) runni
On Wed, Sep 24, 2003 at 02:47:21PM +0200, Martin Olivier wrote:
> Hi all,
>
> The function heatmap uses the functions dist and hclust with default
> parameters.
> How to change these parameters? For example, i want to use the ward
> criterion for hierarchical
> clustering with binary distance.
Søren:
> I take a an object of class loglm and specialize it into class
> c("hllm", "loglm"). I would like to define print.hllm such that it
> prints some special information for the hllm object but ALSO prints
> the loglm-information using the print method for loglm. There is no
> print.loglm met
Hi all,
The function heatmap uses the functions dist and hclust with default
parameters.
How to change these parameters? For example, i want to use the ward
criterion for hierarchical
clustering with binary distance.
Best regards,
Olivier.
__
[EMAIL
Hello
I'm trying to use C random generator function drand48() which return
floating-point values, uniformly distributed over [0,1],
(http://www.opengroup.org/onlinepubs/007908799/xsh/drand48.html). When
values are returned to R, they are not in [0,1]. A simple C program using
drand48() gives v
Dear all,
I take a an object of class loglm and specialize it into class c("hllm", "loglm"). I
would like to define print.hllm such that it prints some special information for the
hllm object but ALSO prints the loglm-information using the print method for loglm.
There is no print.loglm method (
Hi!
> is there a exstisting function (..i found nothing until now.)
> what makes it possible transfrom a dataset:
>
> ID AGE V.MAI V.JUNE
> 11 20 100 120
> 12 30 200 90
>
> into
>
> IDAGEV
> 1120 100
> 1120 120
> 1230200
> 1230
Dear all,
I'm interested in fitting time-series linear models with I(1) errors. Namely
given
y_t=a+b*t+u_t
the random term u_t are such that
u_t-u_{t-1}=e_t~iid N(0,\sigma)
Please, could anyone suggest me any reference (book, article, R functions)
dealing with such models?
Many thanks in advance,
56 matches
Mail list logo