RE: [R] how to test whether two slopes are sign. different?

Or better yet skip the whole significantly different all together, and figure out if a model with 2 slopes explains the data better than a model with 1-- AIC's! Mark -Original Message- From: Brett Magill [mailto:[EMAIL PROTECTED] Sent: Sun 7/20/2003 7:12

RE: [R] how to test whether two slopes are sign. different?

On Sun, 20 Jul 2003, Herzog, Mark wrote: Or better yet skip the whole significantly different all together, and figure out if a model with 2 slopes explains the data better than a model with 1-- AIC's! That's not what AIC is designed to do: it is about `prediction' not `explanation', as

[R] R commands from a text file ?

Hello I was wondering if it was possible to enter R commands from an external text file. If it is possible, it will be easy for repetetive tasks. Does anyone have an idea ? thanks in advance Ahmet Temiz TURKEY __ __

RE: [R] RODBC and Oracle: error table does not exist

I am re-trying a question I asked 12 days ago, to which unfortunately I got no answer so far. Maybe someone who has succesfully established ODBC connections between R and Oracle can give a hint what I am doing wrong? -Ursprüngliche Nachricht- Von: RINNER Heinrich [mailto:[EMAIL PROTECTED]

AW: [R] R commands from a text file ?

You can use source(): For example source(C:\temp\foo.R) wher foo.R is a text fiel containing R commands. Regards, Heinrich. -Ursprüngliche Nachricht- Von: orkun [mailto:[EMAIL PROTECTED] Gesendet: Montag, 21. Juli 2003 08:47 An: R-help Betreff: [R] R commands from a text file ?

RE: [R] R commands from a text file ?

Or you may be interested in ?BATCH Simon. Simon Blomberg, PhD Depression Anxiety Consumer Research Unit Centre for Mental Health Research Australian National University http://www.anu.edu.au/cmhr/ [EMAIL PROTECTED] +61 (2) 6125 3379 -Original Message- From: RINNER Heinrich

Re: [R] R commands from a text file ?

orkun wrote: Hello I was wondering if it was possible to enter R commands from an external text file. If it is possible, it will be easy for repetetive tasks. Does anyone have an idea ? thanks in advance Ahmet Temiz TURKEY __

[R] grid and gtkDevice package clipping bug.

Hello! When I draw points on grid's viewport with clip set to TRUE clipping does not have effect on some symbols. It happens only when drawing over gtkDevice (GTK ver. 1.2) (with postscript device clipping works OK - take a look at this screenshot

RE: [R] RODBC and Oracle: error table does not exist

Dear Marc, thanks very much for your answer!! Adding quotes to the table names didn't change anything (I had tried that before), but creating synonyms in Oracle did the trick! Everything is working fine now, so thanks again. -Heinrich. -Ursprüngliche Nachricht- Von: Marc Mamin

[R] correlated residuals in gls: Coefficient matrix not invertible

Dear Rers, I have threes series, x, y, z and I want to fit a model z ~ x + y. First of all, I fit a lm. I found the residuals are correlated, by looking at the acf() and pacf(). Then I tried to fit a gls model allowing residuals to be correlated (correlation = corARMA(p=5, q=1)): y.na -

RE: [R] Problem with packages...

Have you tried to do: library( the library you want) before using the function? Federico = Federico C.F. Calboli Department of Biology University College London Room 327 Darwin Building Gower Street London WClE 6BT Tel: (+44) 020 7679 4395 Fax (+44) 020 7679 7096

Re: [R] R commands from a text file ? THANKS

M.Kondrin wrote: orkun wrote: Hello I was wondering if it was possible to enter R commands from an external text file. If it is possible, it will be easy for repetetive tasks. Does anyone have an idea ? thanks in advance Ahmet Temiz TURKEY __

[R] C compiler to R

I'd like to know what the more appropriate C compiler to use with R. Thanks, C. Cezar Freitas (ICQ 109128967) IMECC - UNICAMP Campinas, SP - Brasil __ [EMAIL PROTECTED] mailing list

RE: [R] how to test whether two slopes are sign. different?

Dear Stoet This can be handled well in using a mixed-effects model, library (nmle). You can use the lmList option to check whether the slopes differ across populations. -- Harold C. Doran Director of Research and Evaluation New American Schools 675 N. Washington Street, Suite 220

[R] calling R from C

Hi All, We'd like to use functions provided in R in our application. Our application is written in C/C++ and currently runs on win32, Linux and Mac. We'd be happy to attach the whole R ( i.e. not just transfer some function by hand). The important detail is that we deal with big amount

[R] calling R from C

Hi All, We'd like to use functions provided in R in our application. Our application is written in C/C++ and currently runs on win32, Linux and Mac. We'd be happy to attach the whole R ( i.e. not just transfer some function by hand). It is important that we deal with big amount of

Re: [R] help on barplot

Please relieve me my ignorance on one point: How does rbind(d[, 1], d[, 2]) differe from t(d)? Thanks, Spencer Graves Marc Schwartz wrote: On Sun, 2003-07-20 at 21:16, Murad Nayal wrote: Hello, I am trying to compare two histograms using barplot. the idea is to plot the histograms as pairs

[R] Confidence Band for empirical distribution function

Hi, I was trying to draw an empirical distribution function with uniform confidence bands. So I tried to find a way to calculate values of the Kolmogorov-Smirnov Distribution but failed. I guess it must be hidden somewhere (since the ks-test is implemented), but I was unable to find it. Is

Re: [R] C compiler to R

On Mon, 21 Jul 2003, Cezar Augusto de Freitas Anselmo wrote: I'd like to know what the more appropriate C compiler to use with R. If you have an already-compiled R then you should use the compiler that was used to compile R. If you want to compile R then any C compiler should work, but gcc

Re: [R] help on barplot

No practical difference in this case. Either way you end up with the matrix rotated 90 degrees. Presumably t(d) for a large d would be faster, however at 11:30 last night, using rbind() for some reason was the first approach that came to mind... ;-) Marc On Mon, 2003-07-21 at 08:30, Spencer

[R] Lattice: how to draw some text outside the panel?

Dear r-help I draw plots with xyplot() function. Each plot contains also a line of regression. I want to write the trend value and its significance (obtained with lm()) below each panel. I use ltext() for this. But the text is cut, when it comes outside a panel. Moreover

[R] Changing the labels on a regression tree (repeat post - withadded clarity)

Hello I posted a very similar question last week, but the responses I received indicated that my post was unclear I have a regression tree created in rpart with tr.logypsx - rpart(log(YPSX + 1) ~AGE+drugfact+sexfact+as.numeric(OBSX) +WINDLE + EABUSED + PABAU + positive.par + control.par

[R] problem in pot with subset selection

Hi, I try to make a selective plot, but it make an error. The subset only work with length 4, it is correct? plot((ocorrencia/isca)~frag,subset=especieama==c(grupo1,grupo2,grupo3,AnoplotermesSp1,NeocapritermesOpacus),pch=c(1,2,3,4,5),xlab=Área (ha),ylab=Proporção de iscas ocupadas por

[R] doubt about graphics

Hi, I have an data like this x - c(1:11) y - seq(245,445,20) I make a plot plot(x,y) The first y value in y axis is 250 and the last y valeu is 450, the spacing between values are 50, so the y values showed are: 250, 300, 350, 400, 450 I need that values in y axis are: 240, 280, 320, 360,

Re: [R] doubt about graphics

x - c(1:11) y - seq(245,445,20) plot(x,y, axes=F, ylim=c(240, 480)) axis(1) axis(2, at=seq(240, 480, 40)) Worked for me on R 1.7.1 and S-Plus 6.1 under Win2000. hope this helps. spencer graves Ronaldo Reis Jr. wrote: Hi, I have an data like this x - c(1:11) y - seq(245,445,20) I make a plot

[R] how to calculate Rsquare

Hi, I have something like this: x - 1:10 y2 - 30+5*x+rnorm(x,sd=3) y - c(y1,y2) x - c(x,x) plot(x,y) x - 1:10 y1 - 1+5*x+rnorm(x,sd=2) y2 - 30+5*x+rnorm(x,sd=5) y - c(y1,y2) x - c(x,x) f - factor(rep(c(a,b),c(10,10))) m - lm(y~x+f) anova(m) Analysis of Variance Table Response: y

Re: [R] problem in pot with subset selection

On Mon, 21 Jul 2003, Ronaldo Reis Jr. wrote: Hi, I try to make a selective plot, but it make an error. The subset only work with length 4, it is correct? plot((ocorrencia/isca)~frag, subset=especieama==c(grupo1,grupo2,grupo3,AnoplotermesSp1,NeocapritermesOpacus),

Re: [R] grid and gtkDevice package clipping bug.

Duncan Temple Lang wrote: Hi there. What commands did you use to get this effect? I can't seem to reproduce it on my machine with some simple commands to get the same display. Thanks, D M.Kondrin wrote: Hello! When I draw points on grid's viewport with clip set to TRUE clipping does not

[R] RODBC: problem saving a new table in an Excel database

Hi I am using package RODBC version 1.0-1 under R version 1.7.1 on Windows XP Pro. I am having problems writing a new table to an (Excel) database using sqlSave. I connect to an empty Excel spreadsheet using odbcConnectExcel (which, I believe, uses the Microsoft Excel Driver DSN). Then I try and

Re: [R] how to calculate Rsquare

It should be easy to get a separate R^2 from lm for each level of a factor: Just split the data and run lm once for each level of f. I've done this with a for loop something like the following: dat - data.frame(x=x, y=y, f=f) for(i.f in 1:2){ sel - (f == c(a,

Re: [R] RODBC: problem saving a new table in an Excel database

The driver does not support writing to Excel: did you check? On Mon, 21 Jul 2003, David Khabie-Zeitoune wrote: Hi I am using package RODBC version 1.0-1 under R version 1.7.1 on Windows XP Pro. I am having problems writing a new table to an (Excel) database using sqlSave. I connect to

[R] Rosetta Inpharmatics position announcement

Apologies for cross-posting to both R-Help and S-News. I am not a subscriber to this list so please direct all inquires to [EMAIL PROTECTED] Senior Research Scientist - Data Analysis The Senior Research Scientist Data Analyst position will be responsible for mathematical modeling and analysis

[R] Inconsistent handling of character NA?

[R 1.7.1 on Windows XP Pro] Since R allows missing values for character variables, why are NA's not propagated by character manipulation functions? For example: temp - c(a, NA) temp [1] a NA is.na(temp) [1] FALSE TRUE paste(temp[1], temp[2]) [1] a NA substr(temp, 1, 1) [1] a N

Re: [R] [OT] Modeling strategies

Douglas Bates [EMAIL PROTECTED] writes: I happened to look up Frank Harrell's book Regression Modeling Strategies on Amazon.com today. I was surprised to see that in addition to the typical links to related books they had sponsored links to sites about How to become a model, Try out for

Re: [R] Contents of R-help digest.-contouring

Comments inline: john lewis wrote: R- Users: Can someone indicate what I am during wrong? This is a script essentially from Venerable Ripley's text on interpolating a surface with loess function but I can not get it to run. Thanks. John Lewis Professor McGill University Montreal library(MASS)

Re: [R] Contents of R-help digest.-contouring

The 'z' argument to contour must be a matrix with the appropriate dimensions -- lo.2 as defined originally has dimensions 4356 by 1. You can create a matrix of appropriate dimensions as follows: lo.2 - array(lo.1$fit, dim=sapply(topo.mar, length)) It's also probably worth changing topo.mar so

Re: [R] Inconsistent handling of character NA?

On Mon, 21 Jul 2003, Raubertas, Richard wrote: [R 1.7.1 on Windows XP Pro] Since R allows missing values for character variables, why are NA's not propagated by character manipulation functions? They are in the development version. For example: temp - c(a, NA) temp [1] a NA

Re: [R] [OT] Modeling strategies

On 21 Jul 2003 08:32:33 -0500 Douglas Bates [EMAIL PROTECTED] wrote: I happened to look up Frank Harrell's book Regression Modeling Strategies on Amazon.com today. I was surprised to see that in addition to the typical links to related books they had sponsored links to sites about How to

Re: [R] Contents of R-help digest.-contouring

On Mon, 21 Jul 2003 13:45:39 -0400, john lewis [EMAIL PROTECTED] wrote : R- Users: Can someone indicate what I am during wrong? This is a script essentially from Venerable Ripley's text on interpolating a surface with loess function but I can not get it to run. That's Venables... topo.lo -

Re: [R] Contents of R-help digest.-contouring

On Mon, 21 Jul 2003, Duncan Murdoch wrote: On Mon, 21 Jul 2003 13:45:39 -0400, john lewis [EMAIL PROTECTED] wrote : R- Users: Can someone indicate what I am during wrong? This is a script essentially from Venerable Ripley's text on interpolating a surface with loess function but I

[R] bold AND italic as font in text()

Dear all, Is it possible to somshow plot text as italic AND bold. I tried font=c(2,3) in text(), but it doesn't work. It seems like the latter value is used. Thanks in advance! Sincerely, Tord --- Tord Snäll Avd. f

Re: [R] Contents of R-help digest.-contouring

The reference to Venerable and Ripley reminds me that in the late 1940s, George Box worked for Imperial Chemicals Industries, which appeared in at least one brochure I saw as Empirical Chemicals Industries. George invented response surface methods there, so the typo was not as inappropriate

[R] generate a series of fucntion

Hi there, I want to generate a large amount of functions, say f=function(x,t) exp(-t[1]-t[2]*g_1(x)-t[3]*g_2(1+x)) where g_1(x) and g_2(x) are from a long list of moments, such as x, x^2, log(x), log(1+x) .. and so on. Any suggestions on how to do this efficiently? thanks a lot. x.w

[R] Random Numbers

Thank you to all who replied to my previous question regarding the generation of random numbers. I have read up on the literature and learnt a whole lot more. In my research I found code for many generators which I have translated into java for my project. I'm now at the stage of testing

RE: [R] generate a series of fucntion

Does this do sort of what you want? g1 - list(function(x) log(x), function(x) log(x+1), function(x) x, function(x) x^2) g1 [[1]] function(x) log(x) [[2]] function(x) log(x+1) [[3]] function(x) x [[4]] function(x) x^2 lapply(g1, function(g) {function(x,t) exp(-t[1]-t2*g(x)-t[3]*g(1+x))})

Re: [R] generate a series of fucntion

On Mon, 21 Jul 2003 15:54:24 -0400, ximing wu [EMAIL PROTECTED] wrote : Hi there, I want to generate a large amount of functions, say f=function(x,t) exp(-t[1]-t[2]*g_1(x)-t[3]*g_2(1+x)) where g_1(x) and g_2(x) are from a long list of moments, such as x, x^2, log(x), log(1+x) .. and so on.

Re: [R] Setting name attributes to a vector - join?

I assume that cnd.names is a vector of character strings and landcov[ ,2] is a vector of integers in c(1:length(cnd.names). Then foo - cnd.names[landcov[ ,2]] should be what you want. - tom blackwell - u michigan medical school - ann arbor - On Mon, 21 Jul 2003, Andy Bunn

Re: [R] bold AND italic as font in text()

use text(, , expression(bold(what you want to write))) Mahbub. --- Tord Snall [EMAIL PROTECTED] wrote: Dear all, Is it possible to somshow plot text as italic AND bold. I tried font=c(2,3) in text(), but it doesn't work. It seems like the latter value is used. Thanks in advance!

Re: [R] Contents of R-help digest.-contouring

Spencer Graves [EMAIL PROTECTED] writes: The reference to Venerable and Ripley reminds me that in the late 1940s, George Box worked for Imperial Chemicals Industries, which appeared in at least one brochure I saw as Empirical Chemicals Industries. George invented response surface methods

[R] Re: [Rd] calling R from C

On Mon, 21 Jul 2003, Khamenia, Valery wrote: Hi All, We'd like to use functions provided in R in our application. Our application is written in C/C++ and currently runs on win32, Linux and Mac. We'd be happy to attach the whole R ( i.e. not just transfer some function by hand).

[R] Analysis of Complex Survey Data

Hi all I would like to perform a logistic regression analysis on some complex survey data with R, but am not sure if there are functions within R that will enable me to do so. Also, are there any extensions of the cor function that would enable me to incorporate survey weights when calculating

Re: [R] Analysis of Complex Survey Data

On Tue, 22 Jul 2003 [EMAIL PROTECTED] wrote: Hi all I would like to perform a logistic regression analysis on some complex survey data with R, but am not sure if there are functions within R that will enable me to do so. svyglm() in the survey package will fit generalised linear models to

RE: [R] generate a series of fucntion

ximing wu [EMAIL PROTECTED] wrote: I want to generate a large amount of functions, say f=function(x,t) exp(-t[1]-t[2]*g_1(x)-t[3]*g_2(1+x)) where g_1(x) and g_2(x) are from a long list of moments, such as x, x^2, log(x), log(1+x)

[R] R and C++ compared with only C++

My computer is a pentium 4 running at 2.4 GHz. My R is 1.7.1 I have written a program in R that calls C++. The program spends most of the time in C++ ( 90% ). R basically deals with output and input. How slower can this be compared with the program I would get from rewriting everything in C++?

[R] function default

Greetings, Is there a way to create default arguments for functions? If there is, then please explain or forward me to the instructions. I could not find anything in the archives or in the manual under function or args. Thanks -- __

RE: [R] function default

Simply declare the defaults where you specify the function arguments. e.g. test - function (x=4, y=5) x*y test() [1] 20 Simon. Simon Blomberg, PhD Depression Anxiety Consumer Research Unit Centre for Mental Health Research Australian National University http://www.anu.edu.au/cmhr/ [EMAIL

Re: [R] function default

Most functions in R have default arguments, and these are specified in the function definition. For instance, the following function has a default argument of 10 which is used if none is specified: do.it - function(n=10) { rnorm(n) } do.it() Regards, Andrew C. Ward CAPE Centre Department of

Re: [R] R and C++ compared with only C++

As always, it depends on your priorities and motivation. If that extra 10% is important and if you can be bothered writing the C++, then putting the whole thing into C++ may be worth it. If not then you could let R deal with the input and output. Regards, Andrew C. Ward CAPE Centre Department

Re: [R] function default

fun - function(x=2, y=3){ c(x, y) } fun() should return (2, 3); fun(1) should return (1, 3); fun(y=9) should return (2, 9), etc. hope this helps. spencer graves Nicholas Croglio wrote: Greetings, Is there a way to create default arguments for functions? If there is, then please

Re: [R] Lattice: how to draw some text outside the panel?

Hi Wladimir Eremeev wrote: Dear r-help I draw plots with xyplot() function. Each plot contains also a line of regression. I want to write the trend value and its significance (obtained with lm()) below each panel. I use ltext() for this. But the text is cut, when it comes

Rcmdr problems (was:Re: [R] Minor nuisance with rw1071)

On 18 Jul 2003 at 9:54, John Fox wrote: Sorry for being so late coming back on this, I was travelling the weekend. I tried to remove as much as possible from the workspace (.RData file), but nothing created by Rcmdr. This did'nt change anything. Then I downloaded the latest Rcmdr (0.9-0,

Re: [R] R and C++ compared with only C++

Francisco J Molina [EMAIL PROTECTED] asked: I have written a program in R that calls C++. The program spends most of the time in C++ ( 90% ). ... How slower can this be compared with the program I would get from rewriting everything in C++? Suppose the

Re: [R] bold AND italic as font in text()

Hi Tord Snall wrote: Dear all, Is it possible to somshow plot text as italic AND bold. I tried font=c(2,3) in text(), but it doesn't work. It seems like the latter value is used. Looks like you want font=4. Try the following: plot(0:5, 0:5, type=n) text(1:4, 1:4, paste(font =