[R] A problem with "copy()" in the svIO package

2006-12-20 Thread Kenneth Cabrera
Hi R-users: When I type the following code in R, using the svIO package I got a problem only with the "latex" type option.(With "raw", "ascii" and "html" options, it works fine). x<-1:50 copy('x', type='latex', objname='x') Erro en file.info(fn <- c(...)) : argumento del nombre de archivo inválid

Re: [R] Upgrading

2006-12-20 Thread Ivailo Stoyanov
On Wed, 2006-12-20 at 11:37 -0500, DEEPANKAR BASU wrote: > Thanks for the comments. I had got Thomas' message wrong the first time. Now, > I have downloaded and re-installed the latest version of R. After reading the > "administration and maintenance" manual, I added > > deb http://cran.R-projec

Re: [R] call by reference

2006-12-20 Thread Henrik Bengtsson
There is no support for 'call by reference' in the S language, and this is intentionally, but you can use environments to imitate it, cf. ?environment. See also the R.oo package. /Henrik On 12/20/06, biter bilen <[EMAIL PROTECTED]> wrote: > Can anyone help me about pass by reference of arguments

Re: [R] Problems with Optimization

2006-12-20 Thread Tobias
My apologies if this came across as unfair. This was by no means intended. I find R to be a fantastic software and at least for my needs faster and more comfortable to use than commercial ones such as Matlab. I am just wondering, because I implemented the same code on GAUSS and it seemed to hand

Re: [R] Problems with Optimization

2006-12-20 Thread Prof Brian Ripley
On Wed, 20 Dec 2006, Tobias wrote: > > Dear R-helpers, > > I am having following problem: > > Let P be an observed quantity, F(...) a function describing P, and e = P - > F(...) the error. > > F(...) is essentially a truncated mean whose value is obtained via > integrating from some value X to inf

[R] Problems with Optimization

2006-12-20 Thread Tobias
Dear R-helpers, I am having following problem: Let P be an observed quantity, F(...) a function describing P, and e = P - F(...) the error. F(...) is essentially a truncated mean whose value is obtained via integrating from some value X to inf over a probability density with six parameters. Th

Re: [R] call by reference

2006-12-20 Thread biter bilen
a = 9 b = 8 print (a) print (b) fa<-function(bS, aV) { aV <<- 3 assign(bS,3,env = .GlobalEnv) } fa("b", a) print (a) #9 print (b) #3 #only this call is ok - Original Message From: Prof Brian Ripley <[EMAIL PROTECTED]> To: biter b

[R] FW: R Version Problem in using write.foreign+SAS

2006-12-20 Thread Shubha Vishwanath Karanth
Could I find a solution for this? Please please please... How can I change the built in code of "foreign:::writeForeignSAS" for the format? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Shubha Vishwanath Karanth Sent: Wednesday, December 20, 2006 6:06 PM

Re: [R] \Sexpr in MikTex

2006-12-20 Thread Duncan Murdoch
On 12/20/2006 6:26 PM, Blanchard, Suzette wrote: > Greetings, \Sexpr{} has worked on MikTex with earlier versions of R, I can > not seem to get it to work with > R-2.4.0. I run Sweave with MikTex using the following statement in Run under > Accessories. > latex -include-directory="C:\Program

Re: [R] How to use strings from a data.frame as the argument of an expression() for plot

2006-12-20 Thread Gabor Grothendieck
Just store them as strings in the data frame and then perform a conversion to expressions prior to the loop: DF <- data.frame(s = c("O[3]", "NO", "NO[2]"), stringsAsFactors = FALSE) e <- parse(text = DF$s) ... continue with prior solution ... On 12/20/06, MrJ Man <[EMAIL PROTECTED]> wrote: > Tha

[R] Spline models in sspir

2006-12-20 Thread Simon Frost
Dear R-Help, I'm trying to learn the sspir package for state space modeling. Has anyone coded a cubic spline smoother (continuous time) in state space format in sspir? The syntax for setting up the various matrices would be really helpful. Best Simon -- Simon D.W. Frost, D.Phil. Assistant Adjunc

[R] lda plotting: labeling x axis and changing y-axis scale

2006-12-20 Thread Wade Wall
Hi all, I have performed an lda on two groups and have plotted using plot(x.lda), with x.lda being my lda results. I have forgotten how to change the labels of the of the x-axes (they are currently listed as Group1 and Group 13), and to rescale the y-axis to reflect frequency. If anyone knows how

Re: [R] how to change the date (Year, Month, Day) to Julian date in R language?

2006-12-20 Thread jim holtman
Check out date.ddmm in the 'date' package. On 12/20/06, Zhang Jian <[EMAIL PROTECTED]> wrote: > > How to do the change? > And how to change Julian date to the date (Year,Month,Day)? > Thanks > >[[alternative HTML version deleted]] > > __ > R-

[R] how to change the date (Year, Month, Day) to Julian date in R language?

2006-12-20 Thread Zhang Jian
How to do the change? And how to change Julian date to the date (Year,Month,Day)? Thanks [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] How to use strings from a data.frame as the argument of an expression() for plot

2006-12-20 Thread MrJ Man
Thanks for your response; one more thing: Is it possible to use a data.frame for the same effect, as these values are associated with some others (I could of course use them separately, but merging them in a data.frame should be possible in R)? > Try this: > > > e <- expression(O[3], NO, NO[2])

Re: [R] How to use strings from a data.frame as the argument of an expression() for plot

2006-12-20 Thread Gabor Grothendieck
Try this: e <- expression(O[3], NO, NO[2]) opar <- par(mfrow = c(2,2)) for(i in 1:3) plot(1, 1, type = "b", main = bquote(.(e[[i]]) ~ Year ~ 2005)) par(opar) Also please read the last line to every post to r-help and particularly note the part about reproducible examples. x and y.were undefine

Re: [R] \Sexpr in MikTex

2006-12-20 Thread Doran, Harold
Suzette I have not experienced any problems with \Sexpr{} and new R versions. It might be helpful if you could provied a minimal example of your .Rnw and how \Sexpr{} is used within. It is obvious what OS you're using given the path dirctories below, but normally it might be useful to be expli

[R] \Sexpr in MikTex

2006-12-20 Thread Blanchard, Suzette
Greetings, \Sexpr{} has worked on MikTex with earlier versions of R, I can not seem to get it to work with R-2.4.0. I run Sweave with MikTex using the following statement in Run under Accessories. latex -include-directory="C:\Program Files\R\R-2.4.0\share\texmf" "C:\Documents and Settings\Su

[R] How to use strings from a data.frame as the argument of an expression() for plot

2006-12-20 Thread MrJ Man
Greetings, I would like to use a data.frame with strings to feed the expression() in the title of a plot. The way I did this is: molecules <-data.frame(name=c("o3","no","no2"),expression=c("quote(O[3])","quote(NO)","quote(NO[2])")) for (mol in c(5,7,9)) { plot(x, y, type="b", main=eval(substitu

Re: [R] writing R extension

2006-12-20 Thread ahmad ajakh
Thanks Barry, I tested this solution and it works. Thanks also to Sarah Goslee for bringing up alternative ideas. I guess I need to get into building proper packages now. AA. - Original Message From: Barry Rowlingson <[EMAIL PROTECTED]> Cc: r-help@stat.math.ethz.ch Sent: Wednesday, Decemb

[R] mclust priors

2006-12-20 Thread Finny Kuruvilla
I am currently exploring mclust to do clustering by Gaussian mixture modeling. There are typically 1-5 clusters in my datasets. In the latest version of mclust (version 3), there is now the ability to specify priors. However, I am confused about the implementation. I've read through the document

Re: [R] Reformat meteorological data

2006-12-20 Thread hadley wickham
Hi Toby, Have a look at the reshape package, http://had.co.nz. Code something like the following should do what you want: library(reshape) dfm <- melt(data, id=1:3) dfm <- rename(dfm, c(variable = "month")) # to make it more obvious cast(dfm, ... ~ var) Hadley On 12/20/06, Toby Gass <[EMAIL P

[R] daisy() and agnes() output

2006-12-20 Thread Marion Echard
Dear All, I used the function daisy() and agnes() in the package cluster to, respectively, generate a dissimilarity matrix and a dendrogram (via plot function). But I have 405 species plotted on the dendrogram. This means that the result is not easily readable. Do anyone know how arrange the

[R] lme4 mcmcsamp matrix not PD

2006-12-20 Thread Martin Henry H. Stevens
Hi folks, I am not providing a small replicable example, because I assume the problem is related to my quirky data (~650 obs.). I am using the latest lme4, matrix and coda, and R 2.4.0. I frequently get the following error message for this particular lmer model. I do not get this message for t

[R] Reformat meteorological data

2006-12-20 Thread Toby Gass
Dear HelpeRs: I have a data set in the following format, which will be familiar to those of you working with NCDC climate data. Example: Year <- rep(1:3, each = 3) Year <- c(Year,Year) ID <- rep(1:2, each = 9) Jan <- runif(18, min = 0, max = 20) Feb <- runif(18, min = 0, max = 20) Mar <- runif(1

Re: [R] writing R extension

2006-12-20 Thread ¨Tariq Khan
Some particularly useful links that explain the process in as few pages as possible! Would highly recommend taking a look. http://www.biostat.jhsph.edu/~bcaffo/statcomp/files/rpacks.pdf http://www.bioinf.uni-hannover.de/teaching/fallstudien/schaarschmidt2.pdf On 12/20/06, Sarah Goslee <[EMAIL P

Re: [R] writing R extension

2006-12-20 Thread Sarah Goslee
Hi, > I am dealing with the same issue here and I was wondering whether it would be > possible to just save > the R compliled function objects in a directory and just attach the directory > to the search path. > (I am using R2.4.0+ESS+Xemacs in windows XP). The other method that I'm familiar w

Re: [R] fit sine?

2006-12-20 Thread Greg Snow
Look at the functions: nls, optimize, optim. Also look at the packages circular and CircStats. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAI

Re: [R] Rotating a distribution plot by 90 degrees

2006-12-20 Thread Greg Snow
Does the following code do what you want? x <- c( rnorm(50,10,2), rnorm(30,20,2) ) y <- 2+3*x + rnorm(80) d.x <- density(x) d.y <- density(y) layout( matrix( c(0,2,2,1,3,3,1,3,3),ncol=3) ) plot(d.x$x, d.x$y, xlim=range(x), type='l') plot(d.y$y, d.y$x, ylim=range(y), xlim=rev(range(d.y$y)), type

Re: [R] writing R extension

2006-12-20 Thread ahmad ajakh
Hi all, I am dealing with the same issue here and I was wondering whether it would be possible to just save the R compliled function objects in a directory and just attach the directory to the search path. (I am using R2.4.0+ESS+Xemacs in windows XP). Thanks. AA. - Original Message F

Re: [R] Newbie data organisation/structures question...

2006-12-20 Thread Michael Kubovy
On Dec 20, 2006, at 11:05 AM, Gav Wood wrote: > So my data is in this sort of format: > > P T I > 1 1 (1, 2, 3) > 2 1 (2, 4) > 1 2 (1, 3, 6, 7) > 2 2 (6) Not knowing why you organized the data as you did, let me suggest another approach: iv <- c(1, 2, 3, 2, 4, 1, 3, 6, 7, 6) p <- c(1

Re: [R] Newbie data organisation/structures question...

2006-12-20 Thread Marc Schwartz
On Wed, 2006-12-20 at 16:05 +, Gav Wood wrote: > Howdo folks, > > So my data is in this sort of format: > > P T I > 1 1 (1, 2, 3) > 2 1 (2, 4) > 1 2 (1, 3, 6, 7) > 2 2 (6) > > And I want to be able to quickly get: > > 1: The I when both P and T are given. e.g.: > P = 2, T = 2; I =

Re: [R] writing R extension

2006-12-20 Thread Hans-Peter
2006/12/20, michele de meo <[EMAIL PROTECTED]>: > I'd like to develop a simple library in R in which to save my > particular functions. > I know that it's possible to save functions in the workspace, but I > need this functions in a library. - As an option you could also save the function in a te

Re: [R] how to do muliple comparisons in linear mixed models

2006-12-20 Thread Dieter Menne
liu, jcheng gmail.com> writes: > i want to compare the several main effects in a linear model. i wonder > how to do the multiple comparisons for less confident intervals. Package multcomp by Torsten Hothorn. Dieter __ R-help@stat.math.ethz.ch mai

Re: [R] writing R extension

2006-12-20 Thread Barry Rowlingson
Sarah Goslee wrote: > If that's still too complex, you could also save your function to a file > and load it as needed with source(). That will give the user the > same effect. > > source("/path/to/my/stuff/myfiles.R") > > Since you didn't tell us OS or anything else about your system, > it's ha

Re: [R] help for multinominal logistic regression code

2006-12-20 Thread Prof Brian Ripley
If you meant 'multinomial', help.search("multinomial") and help.search("logistic") both get you to multinom in package nnet. That is support software for a book, where the detailed documentation is. On Wed, 20 Dec 2006, Mehmet Akif ERDOGAN wrote: > Dear all, > > firstly I would like to say I am

Re: [R] RuleFit & quantreg: partial dependence plots; showing an effect

2006-12-20 Thread Ravi Varadhan
Thanks, Roger. These should be very useful tools. Ravi. --- Ravi Varadhan, Ph.D. Assistant Professor, The Center on Aging and Health Division of Geriatric Medicine and Gerontology Johns Hopkins University Ph:

Re: [R] writing R extension

2006-12-20 Thread Sarah Goslee
For a very simple library, you really only need chapter 1 of the manual on writing R extensions, which even describes the "helper functions" that take some of the work out of making a "proper" package. If that's still too complex, you could also save your function to a file and load it as needed w

[R] writing R extension

2006-12-20 Thread michele de meo
I'd like to develop a simple library in R in which to save my particular functions. I have read the manual on "Writing R Extensions" but it's too difficult. Someone could help me? I want only save my personal function (recorded in R-code, not in C) in "myLibrary" and I want to call it with: >lib

[R] help for multinominal logistic regression code

2006-12-20 Thread Mehmet Akif ERDOGAN
Dear all, firstly I would like to say I am a beginner user for R and also naturally a new member of this mail list. This means I have newer read previous mails before. I need a small help for my analysis. I want to perform a multinominal logistic regression for my data set. But as I

Re: [R] Upgrading

2006-12-20 Thread DEEPANKAR BASU
Thanks for the comments. I had got Thomas' message wrong the first time. Now, I have downloaded and re-installed the latest version of R. After reading the "administration and maintenance" manual, I added deb http://cran.R-project.org/bin/linux/ubuntu dapper/ to my sources.list file and used ap

[R] Newbie data organisation/structures question...

2006-12-20 Thread Gav Wood
Howdo folks, So my data is in this sort of format: P T I 1 1 (1, 2, 3) 2 1 (2, 4) 1 2 (1, 3, 6, 7) 2 2 (6) And I want to be able to quickly get: 1: The I when both P and T are given. e.g.: P = 2, T = 2; I = (6) 2: The concatenated vector of Is when P and a subset of T is given, e.g.:

Re: [R] RuleFit & quantreg: partial dependence plots; showing an effect

2006-12-20 Thread roger koenker
On Dec 20, 2006, at 8:43 AM, Ravi Varadhan wrote: > Dear Roger, > > Is it possible to combine the two ideas that you mentioned: (1) > algorithmic > approaches of Breiman, Friedman, and others that achieve > flexibility in the > predictor space, and (2) robust and flexible regression like QR

Re: [R] biocondutor installation problem

2006-12-20 Thread Seth Falcon
I'm pretty sure your subject should be: Help configure R to use web proxy on RHEL4. You don't have a Bioconductor specific problem (yet). "Kesavan Asaithambi" <[EMAIL PROTECTED]> writes: >> source("http://www.bioconductor.org/biocLite.R";) > > Error in file(file, "r", encoding = encoding) : >

Re: [R] RuleFit & quantreg: partial dependence plots; showing an effect

2006-12-20 Thread Ravi Varadhan
Dear Roger, Is it possible to combine the two ideas that you mentioned: (1) algorithmic approaches of Breiman, Friedman, and others that achieve flexibility in the predictor space, and (2) robust and flexible regression like QR that achieve flexibility in the response space, so as to achieve compl

Re: [R] RuleFit & quantreg: partial dependence plots; showing an effect

2006-12-20 Thread roger koenker
They are entirely different: Rulefit is a fiendishly clever combination of decision tree formulation of models and L1-regularization intended to select parsimonious fits to very complicated responses yielding e.g. piecewise constant functions. Rulefit estimates the conditional mean of th

[R] DBI + ROracle problem with parser ??

2006-12-20 Thread Christian Hoffmann
Hi, Database queries using the combination DBI + ROracle are handicapped by quirks in the pipeline between the textual representation of the query and the database engine Oracle. dbGetQuery(conn, query): work: dbGetQuery(conn, "select * from dual") dbGetQuery(conn, "select * from dual /* comment

Re: [R] Replacing values

2006-12-20 Thread Ivailo Stoyanov
On Mon, 2006-12-18 at 10:58 -0800, downunder wrote: > Hi all, > > I have to recode some values in a dataset. for example changing all zeros to > "." or 999 would be also ok. does anybody know how to do this? thanks in > advance. lars Lars, I would perform all such "editing" tasks in a spreadshee

Re: [R] R Version Problem in using write.foreign+SAS

2006-12-20 Thread Shubha Vishwanath Karanth
Hi, If I am using R 2.4.0, and work with write.foreign command, I get the datetime format as 01Dec2006 00:00:00, because the built in function (foreign:::writeForeignSAS) has the fixed format for Datetimes and the format is "%d%b%Y %H:%M:%S". And that's the reason why I get the datetime as 01Dec20

Re: [R] How to write a two-way interaction as a random effect in a lmer model?

2006-12-20 Thread Martin Henry H. Stevens
Hi Eva, A couple questions: Are repeated measurements taken on rnr? Is rnr "subject"? Is stress a continuous variable? See below. On Dec 11, 2006, at 11:04 AM, Fucikova, Eva wrote: > Dear All, > > > > I am working with linear mixed-effects models using the lme4 > package in > R. I created a mod

[R] R Version Problem in using write.foreign+SAS

2006-12-20 Thread Shubha Karanth
Hi experts, I have a problem in Write.foreign command (SAS). I have a data frame called d. >d Datetime 2006-12-01 00:00:00 2006-12-01 00:10:00 2006-12-01 00:20:00 2006-12-01 00:30:00 2006-12-01 00:40:00 >class(d$Datetime) [1] "POSIXt" "POSIXct" Then I tried with, write.foreig

[R] Fwd: Where to find the code in the foreign library?

2006-12-20 Thread Shubha Karanth
Hi R users, If I give foreign:::writeForeignSAS in command prompt of R, I get the code used for it. Now I need to make some changes here and run my usual foreign command. So, where can I find the above code in SAS foreign library? Thanks in advance. [[alternative HTML version deleted

[R] Where to find the code in the foreign library?

2006-12-20 Thread Shubha Vishwanath Karanth
Hi R users, If I give foreign:::writeForeignSAS in command prompt of R, I get the code used for it. Now I need to make some changes here and run my usual foreign command. So, where can I find the above code in SAS foreign library? Thanks in advance. [[alternative HTML version dele

[R] FW: R Version Problem in using write.foreign+SAS

2006-12-20 Thread Shubha Vishwanath Karanth
Hi R experts, I have a problem in Write.foreign command (SAS). I have a data frame called d. >d Datetime 2006-12-01 00:00:00 2006-12-01 00:10:00 2006-12-01 00:20:00 2006-12-01 00:30:00 2006-12-01 00:40:00 >class(d$Datetime) [1] "POSIXt" "POSIXct" Then I tried with, wri

Re: [R] Query regarding linking R with Matlab

2006-12-20 Thread Henrik Bengtsson
Hi. On 12/20/06, Bhanu Kalyan.K <[EMAIL PROTECTED]> wrote: > > Sir, > > I am still new to the R-matlab interfacing. I will explain you the problem > statement more clearly. > > The following is a matlab code. (swissroll.m) > = > > % SWISS R

[R] Rotating a distribution plot by 90 degrees

2006-12-20 Thread Benjamin Otto
Hi, Can I rotate a plot (e.g. a distribution plot) by 90 degrees? The barplot function provides the "horiz" command but that's not availeable for the base package functions. I found an old advice from Paul Murrell on a similar problem suggesting to use viewports (grid package). Yet I couldn't repr

Re: [R] R.matlab question

2006-12-20 Thread Henrik Bengtsson
Hi. On 12/20/06, Aimin Yan <[EMAIL PROTECTED]> wrote: > Does anyone know how to solve this question about R.matlab? > I am in windowsXP, my matlab is matlab 7.0.0 19920(R14) > > thanks, > > Aimin > > > matlab <- Matlab(host="localhost", port=9998) > > if (!open(matlab)) throw("Matlab server is n

[R] Query regarding linking R with Matlab

2006-12-20 Thread Bhanu Kalyan.K
Sir, I am still new to the R-matlab interfacing. I will explain you the problem statement more clearly. The following is a matlab code. (swissroll.m) = % SWISS ROLL DATASET N=2000; K=12; d=2; clf; colordef none; colormap jet; se

Re: [R] RuleFit & quantreg: partial dependence plots; showing an effect

2006-12-20 Thread Mark Difford
Dear List, I would greatly appreciate help on the following matter: The RuleFit program of Professor Friedman uses partial dependence plots to explore the effect of an explanatory variable on the response variable, after accounting for the average effects of the other variables. The plot method

Re: [R] Upgrading

2006-12-20 Thread Petr Pikal
Hi as the addition to Thomas's answer. On 19 Dec 2006 at 15:20, Thomas Lumley wrote: Date sent: Tue, 19 Dec 2006 15:20:18 -0800 (PST) From: Thomas Lumley <[EMAIL PROTECTED]> To: DEEPANKAR BASU <[EMAIL PROTECTED]> Copies to: R-help@s

Re: [R] call by reference

2006-12-20 Thread guillet
> Can you tell us what you want to do with 'pass/call by reference'? > > If you want an R function to alter its argument then it is possible > (KalmanLike is an example), but it should only be possible via C code. > And if you know enough to do that, you probably would not be asking (and > definite

Re: [R] call by reference

2006-12-20 Thread Prof Brian Ripley
Can you tell us what you want to do with 'pass/call by reference'? If you want an R function to alter its argument then it is possible (KalmanLike is an example), but it should only be possible via C code. And if you know enough to do that, you probably would not be asking (and definitely not be