Re: [R] Lattices: Cloud: Background

On Friday 16 January 2004 00:51, Mulholland, Tom wrote: When I first started using lattice I found the colour schemes a bit confusing. So eventually I came up with the colours I wanted. The code below was one of those attempts. One thing that happened however was that I kept shutting down the

Re: [R] SIR

Vera Hofer wrote: Ich habe in R eine Version von SIR gefunden und ausprobiert. Leider kann diese multivariate Responses nicht verarbeiten. Gibt es in R eine ausgefeilte Version von SIR? Danke für die Hilfe. V.H. __ [EMAIL PROTECTED] mailing list

Re: [R] memory in R

Raphael Schoenle wrote: Hi, How can I assign memory in R? My simulations are very slow. What do you mean with assign memory in R? You might want to optimize your code... Uwe Ligges Thanks!, -R [[alternative HTML version deleted]] __ [EMAIL

[R] Naïve Bayes

I would like to know if there is a function to Naïve Bayes Classifier in R Almiro Moreira DRN/Núcleo de Estatísticas da Indústria Telf. 226072098 Fax. 226058204 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] SIR

I have found a version of SIR in R and I have tried it. But the problem with this file is the fact that it does not cope with multivariate response variables. Is there any version of SIR available that also works with multivariate responses? Thanks for help!

Re: [R] Naïve Bayes

Yes, look at library(e1071) christian Am Freitag, 16. Januar 2004 09:56 schrieb Almiro Moreira: I would like to know if there is a function to Naïve Bayes Classifier in R Almiro Moreira DRN/Núcleo de Estatísticas da Indústria Telf. 226072098 Fax. 226058204

[R] error message in plot(aov-object) -- repost

Hi all, I posted this question several days ago, but did not get any answer until now. Since I still have no clue about the source of this error message, I repost a description of the problem including some code: A student at our institute fitted an aov model, and got the following error

[R] Using dates for plotting

Dear R-Users! I have the following data: runs dates kms 1 2004-01-15 12:00:00 9 2 2004-01-16 12:00:00 10 3 2004-01-17 12:00:00 5 class(runs$dates) [1] POSIXt POSIXct Is it possible to _easily_ plot these data with real dates shown as tick marks on the x-axis? Perhaps even when the

Re: [R] Multiple comparisons in R; multicomp

Patrick I'm not familiar with the multcomp package from s-plus, but there is a package available through bioconductor (www.bioconductor.org) called multtest that has some functions for multiple testing procedures and adjusting p-values and computing false discovery rates. Sean -- Sean Davis,

Re: [R] Using dates for plotting

On Fri, 16 Jan 2004, Christoph Lange wrote: Dear R-Users! I have the following data: runs dates kms 1 2004-01-15 12:00:00 9 2 2004-01-16 12:00:00 10 3 2004-01-17 12:00:00 5 class(runs$dates) [1] POSIXt POSIXct Is it possible to _easily_ plot these data with real dates

[R] plot

Dear All, I'd like to 'highlight' say change font size/color 'some' point in a graph. I can do this with 'points', plot(x) points(x[1], col=red) Is there a more 'straightforward' way (e.g., option in plot) to do this? Thanks Marwan

Re: [R] plot

Marwan Khawaja wrote: Dear All, I'd like to 'highlight' say change font size/color 'some' point in a graph. I can do this with 'points', plot(x) points(x[1], col=red) Is there a more 'straightforward' way (e.g., option in plot) to do this? Thanks Marwan plot(x, col = colorvector) ??? Uwe

Re: [R] error message in plot(aov-object) -- repost

Pascal A. Niklaus wrote: Hi all, I posted this question several days ago, but did not get any answer until now. Since I still have no clue about the source of this error message, I repost a description of the problem including some code: A student at our institute fitted an aov model, and got

Re: [R] error message in plot(aov-object) -- repost

Uwe Ligges [EMAIL PROTECTED] writes: plot.lm() produces the error message, in particular the lines hii - lm.influence(x, do.coef = FALSE)$hat which appears to be 1 in some cases, and rs - r.w/(s * sqrt(1 - hii)) (you get a division by zero and the outcome is Inf, hence it's

[R] [R-pkgs] FW: Announce: RPy version 0.3.3

-Original Message- From: Warnes, Gregory R Sent: Friday, January 16, 2004 3:52 AM To: 'r-packages'; 'python-announce '; 'rsoap-talk'; 'Dan Nathan'; 'RPy List' Subject: Announce: RPy version 0.3.3 RPy Version 0.3.3 is now available from the RPy home page at

[R] reference to objects

Hi, is there a way to reference to a data object without copying it? For example I have a huge matrix called dist and I want two objects obj1 and obj2 to have a memeber dist that points to the matrix, but I don't want, for memory reasons, to copy the matrix twice. As far as I understand the

Re: [R] Exactness of ppois

Matthias == Matthias Kohl [EMAIL PROTECTED] on Thu, 15 Jan 2004 13:55:22 + writes: Matthias Hello, by checking the precision of a convolution Matthias algorithm, we found the following inexactness: Matthias We work with R Version 1.8.1 (2003-11-21) on Matthias Windows

[R] individual likelihoods

Dear all, is there a way to extract individual likelihoods from a glm/lrm object? By individual likelihoods, I mean the likelihoods whose product give the overall likelihood of the model. I guess the code in the base package, used to compute the Akaike Information Criterion may help me. However, I

RE: [R] SIR

SIR and other dimension reduction methods are available in the dr package on CRAN. (No idea if that's what Vera was referring to.) From a quick glance, I don't think it handles multivariate responses. Cheers, Andy From: Uwe Ligges Vera Hofer wrote: Ich habe in R eine Version von SIR

Re: [R] reference to objects

Hi, You will find some pieces of information about the way to handle such things at: http://www.maths.lth.se/help/R/ Programming with references Note that you will have to download and install the package developped by the author, Henrik Bengtsson HTH, Eric At 14:22 16/01/2004, Giampiero

Re: [R] reference to objects

Thanks for your answer, Yes, that will create two copies. What is it you want to do with the data? Do you want the capability of both of them changing the data? What type of processing are you going to do? The data should be read only (and all the objects share the same data values). One

Re: [R] memory in R

To expand a bit on Uwe's comments: R is vectorized and loops are notoriously inefficient. Consider the following: start.time - proc.time() set.seed(1) X - array(rnorm(1000), dim=c(100, 10)) X.1 - apply(X, 2, mean) print(var(X.1)) (elapsed.time - proc.time()-start.time) start.ttime -

Re: [R] reference to objects

Thank you, I'll read the documentation... Giampiero On Fri, 16 Jan 2004, Eric Lecoutre wrote: Hi, You will find some pieces of information about the way to handle such things at: http://www.maths.lth.se/help/R/ Programming with references Note that you will have to download and install

Re: [R] individual likelihoods

Dear Bruno, residuals(mod, type=deviance) gives you the deviance residuals for the model, so sum(residuals(mod, type=deviance)^2) is the deviance. I hope that this helps, John At 02:29 PM 1/16/2004 +0100, Bruno Giordano wrote: Dear all, is there a way to extract individual likelihoods from a

[R] analytic derivative of complex function

Dear R experts, I'd like to take analytic derivative of complex functions through D function. I succeded with simple functions like sin(x), but what if I want to work with function like A(x,y, ...) + i*B(x,y, ...)? Thank you. -- WBR, Timur. __ [EMAIL

[R] Saving jpg of plot with grid.arrows

I am having problems exporting a plot to a jpg file. I am first setting up a basic plot, and then opening up a viewport and overlaying grid.arrows onto the plot - I don't seem able to save this plot with the overlaid arrows to a file, can anyone offer any advice? I am using R-1.8.0 on Debian

[R] Nlm

Hi I use both R and Splus to fit miximum likelihood models. In general these are Gaussian mixture models of various kinds. The ms (in Splus) and nlm (in R) generally provide us with the same results However nlm is not nearly as stable and reliable with respect to finding the min of a surface.

Re: [R] SIR

This is strange; the sir for R I know (in package dr on CRAN, from S. Weisberg), last time I checked (about a year ago?) was able to handle multivariate responses. In fact, p. 6 of the documentation shows an example of SIR with a bivariate response, and I tried it, and it works. Best, R. On

Re: [R] invoking R scripts from a linux shell ?

Hi Enrico, As seen by other posts, there are many ways to accomplish this sort of thing. For simple tasks I've often found it easiest to call R from my script (in my case Python). I've not used BATCH mode and instead call R like this: R --slave --no-save script.R A common setup for me is

[R] anova repeated measure interpretation

Dear all, I have tried to use R for the repeated measures experiment design in a phonetic study, and in the resulting forms I met some problems which perturb the final interpretation of the results. I will explain the experiment design first: the same 7 subjects were answering a question about 25

re: [R] how to overlap plots

Try the distplot() function in package vcd. -Michael -- Michael Friendly Email: [EMAIL PROTECTED] Professor, Psychology Dept. York University Voice: 416 736-5115 x66249 Fax: 416 736-5814 4700 Keele Streethttp://www.math.yorku.ca/SCS/friendly.html Toronto, ONT M3J 1P3 CANADA

Re: [R] Saving jpg of plot with grid.arrows

On Friday 16 January 2004 10:02, Laura Quinn wrote: I am having problems exporting a plot to a jpg file. I am first setting up a basic plot, and then opening up a viewport and overlaying grid.arrows onto the plot - I don't seem able to save this plot with the overlaid arrows to a file, can

[R] Summary: Installing the Rcmdr and tcltk

Hello R-users, thanks to Brian Ripley, who guided me through this process, the tcltk package is now working and hence the Rcmdr package as well. The secrets were to first rebuild the Tcl and Tk packages from the sources files, making sure the same C compiler was used as the one used in

[R] Sparse matrix row/column access (SparseM)?

Hello list, I am creating sparse matrices using the SparseM package. These sparse structures are to be used as a basis for (a number of) distance matrices. But in order to create thoses, I must be able to read the elements/rows/columns of the sparse structures in order to operate on them. Thus,

Re: [R] Sparse matrix row/column access (SparseM)?

Given a matrix in any of the SparseM forms, you can access columns and rows with the usual R subsetting conventions, e.g. A[,k] gives the kth column, A[m,] gives the mth row... As a general principle, questions about packages should probably be directed first to package maintainers (in this case

[R] Weird problem with trying to change a variable

I have a dataframe called cvar, with two variables (among many others) called MSA and ACTUP. Both are numeric. This was working fine. Then I found out that for two MSAs, ACTUP should be 1, not 0. so I tried cvar$ACTUP[cvar$MSA == 6840] - 1 cvar$ACTUP[cvar$MSA == 5360] - 1 but when I try

Re: [R] Weird problem with trying to change a variable

On Fri, 16 Jan 2004, Peter Flom wrote: I have a dataframe called cvar, with two variables (among many others) called MSA and ACTUP. Both are numeric. This was working fine. Then I found out that for two MSAs, ACTUP should be 1, not 0. so I tried cvar$ACTUP[cvar$MSA == 6840] - 1

Re: [R] ecdf function

On Fri, 16 Jan 2004, Maya Sanders wrote: I am trying to use the ecdf function to find p-values (using a vector of numbers to represent my new distribution and a test specific t-statistic value). I am using : 1-ecdf(vector)(t-stat) vector-c(5.386, 3.701717, 3.8289, 3.602, 4.469, 5.2087,

Re: [R] random effects with lme() -- comparison with lm()

Jerome Asselin [EMAIL PROTECTED] writes: On Thu, 2004-01-15 at 16:30, Douglas Bates wrote: ...snip... (BTW, I wouldn't say that this is equivalent to a fixed effects model. It is still a random effects model with two variance components. It just doesn't have well-defined estimates for

[R] generate random number of any given distribution

Hello, Is there a function in R to generate random number of any given distribution (its pdf is given), besides uniform and gaussian distribution? Thanks, yan __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

RE: [R] generate random number of any given distribution

If you can compute the quantile function of the distribution (i.e., the inverse of the integral of the pdf), then you can use the probability integral transform: If U is a U(0,1) random variable and Q is the quantile function of the distribution F, then Q(U) is a random variable distributed as F.

Re: [R] generate random number of any given distribution

Hello, Yan: Are you aware that for all the standard distributions, R provides the probability density, the cumulative distribution function, the quantile function, and random number generation? When you said, besides uniform and gaussian, I wondered. The convention is that the first

Re: [R] Weird problem with trying to change a variable

Lists? Do you mean vectors? Columns of dataframes are not supposed to be lists. Ah. I've been using commands like x - list(c(1,2,3),c(2,1),c(6,6,1)) y - c(A,B,C) data.frame(I(x),y) x y 1 1, 2, 3 A 22, 1 B 3 6, 6, 1 C This sort of object behaves pretty much as I'd expect it