[R] Unexpected result of read.dbf

Hi there, Unexpected result is given with read.dbf from foreign 0.8-9, which is reproducible in R 2.1.1 with the following sample data: test.dbf KEYCODE,N,19,0 422010010 42201002101 42201002102 42201002103 42201002104 422010060 422010071 422010072 42201008001 42201008002 422010100

[R] Problem with building R packages under Windows

Dear all, I am coming to the guru for advise here. I am a native user of Windows, and S-plus, and the process of migrating my libraries from S to R has been more than a painstaking task... I am currently using R version 2.1.1 in Windows XP SP2. I have read the Writing R extensions, the FAQ in

Re: [R] trouble with wilcox.test

On Wed, 17 Aug 2005, Greg Hather wrote: I'm having trouble with the wilcox.test command in R. Are you sure it is not the concepts that are giving 'trouble'? What real problem are you trying to solve here? To demonstrate the anomalous behavior of wilcox.test, consider wilcox.test(c(1.5,5.5),

Re: [R] Problem with building R packages under Windows

Dr L. Y Hin wrote: Dear all, I am coming to the guru for advise here. I am a native user of Windows, and S-plus, and the process of migrating my libraries from S to R has been more than a painstaking task... I am currently using R version 2.1.1 in Windows XP SP2. I have read the Writing

Re: [R] trouble with wilcox.test

Prof Brian Ripley wrote: On Wed, 17 Aug 2005, Greg Hather wrote: I'm having trouble with the wilcox.test command in R. Are you sure it is not the concepts that are giving 'trouble'? What real problem are you trying to solve here? To demonstrate the anomalous behavior of

[R] Binary kernel discrimination

Hello, Could you tell me if a package exists to perform a binary kernel discrimination using a training set compose of molecules represented by binary fingerprint. This method was first described by Harper in J. Chem. Inf. Comput. Sci 2001 41 1295 and is also described in recent papers

Re: [R] line/bar through median in lattice's bwplot?

M. K. mk_lists at yahoo.ca writes: Is there a way to render a line through the median point in the boxplot generated by the Lattice command bwplot? The line basically bisects the bar at the median point... bwplot(height~voice.part , pch='|', data=singer, xlab=Height (inches)) How to find

Re: [R] accesing slots of S4 class in C code

I am trying to use a custom S4 object in my C code and I cannot get the access to its slots working. The following is a toy example, but even this crashes. In R I have: setClass(pd, representation(data=numeric)) x - new(pd, data=1:5) test - function(pd.obj) { res -

[R] Error messages using LMER

Dear All, After playing with lmer for couple of days, I have to say that I am amazed! I've been using quite some multilevel/mixed modeling packages, lme4 is a strong candidate for the overall winner, especially for multilevel generzlized linear models. Now go back to my two-level poisson model

Re: [R] trouble with wilcox.test

P Ehlers wrote: Prof Brian Ripley wrote: On Wed, 17 Aug 2005, Greg Hather wrote: I'm having trouble with the wilcox.test command in R. Are you sure it is not the concepts that are giving 'trouble'? What real problem are you trying to solve here? To demonstrate the anomalous

Re: [R] newbie- how to add a row to a matrix

thank you all rbind solved my problems!!! simone Il giorno 18/ago/05, alle ore 00:34, Martin Lam ha scritto: # row bind a - matrix(1:5) a a - rbind(a, 6) a # column bind b - matrix(1:5) b b - cbind(b, 6:12) b b - cbind(b, 13) b Hope this helps, Martin --- Simone

[R] matrix indexing

Dear R-users, I was wondering for the following: Let 'x' be a matrix and 'ind' and indicator matrix, i.e., x - array(1:20, dim = c(4, 5)) ind - array(c(1:3, 3:1), dim = c(3, 2)) I'd like to get (as a vector) the elements of 'x' which are not indexed by 'ind'. Since negative indices are not

[R] GLMM - Am I trying the impossible?

Dear all, I have tried to calculate a GLMM fit with lmer (lme4) and glmmPQL (MASS), I also used glm for comparison. I am getting very different results from different functions, and I suspect that the problem is with our dataset rather than the functions, but I would appreciate help in

[R] R: stepwise procedures

hi all i found step(), stepAIC() and some other functions in leaps and subselect. is there a package/function that does the traditional stepwise regression procedures using F in and F out values? i know that step does the selctions based on AIC. /

Re: [R] trouble with wilcox.test

If this is stack overflow (and I don't know that yet: when I tried this on Windows the traceback was clearly corrupt, referring to bratio), the issue is that it is impossible to catch such an error, and it is not even AFAIK portably possible to find the stack size limit (or even the current

Re: [R] RMySQL not loading on Mac OS X

Bill, thanks a lot for your answer. I did not know about the sig-Mac list. I will post there next time if I do not find a solution. Concerning your suggestion: The system default compiler is gcc 4.0, but RMySQL seems to be built using gcc-3.3 regardless if I switch to 3.3. or not. Would

Re: [R] matrix indexing

Hi, you might not consider this more elegant, but how about this x[-apply(ind, 1, function(i) (i[1]-1)*nrow(x) + i[2])] Ben On 8/18/05, toka tokas [EMAIL PROTECTED] wrote: Dear R-users, I was wondering for the following: Let 'x' be a matrix and 'ind' and indicator matrix, i.e., x -

Re: [R] How to assess significance of random effect in lme4

Yes, it is a different issue. ranef() extracts the empirical Bayes estimates, which are the empirical posterior modes. The bVar slot holds the corresponding posterior variances of these modes. Technically, (according to D. Bates) the values in the bVar slot are the the diagonal elements of

Re: [R] GLMM - Am I trying the impossible?

It is not supported to call anova() on a glmmPQL fit. For the glmmPQL fit you show, you have very large parameter estimates for a logistic and have partial separation (as you comment on for the control group): in that case PQL is not a reasonable method. Try fit - glm(dead ~ Parasite *

Re: [R] RMySQL not loading on Mac OS X

On 18/08/2005, at 7:46 PM, Georg Otto wrote: Concerning your suggestion: The system default compiler is gcc 4.0, but RMySQL seems to be built using gcc-3.3 regardless if I switch to 3.3. or not. Would it be a solution to force RMySQL to use gcc 4.0 during built? (Might be a naive idea,

Re: [R] power of a matrix

Hi, Gabor: Thanks. spencer graves Gabor Grothendieck wrote: Its expm. On 8/17/05, Spencer Graves [EMAIL PROTECTED] wrote: Hi, Peter: I couldn't find mexp in the Matrix package, but I did find it in fMultivar and in Lindsey's rmutil. These are different functions, but

Re: [R] How to assess significance of random effect in lme4

Hi, Harold: Thanks for the clarification. I thought I had read the original post. Obviously, I had misread it. Thanks again. spencer graves Doran, Harold wrote: Yes, it is a different issue. ranef() extracts the empirical Bayes estimates, which are the empirical posterior modes. The

Re: [R] matrix indexing

x - array(1:20, dim = c(4, 5)) ind - array(c(1:3, 3:1), dim = c(3, 2)) # instead of using ind (pairs of coordinates) for getting the items in the matrix, you can convert it to a list of single coordinates to point to the item in the matrix: # t = transpose # nrow = get the number of rows indices

Re: [R] [SPAM] - Re: How to assess significance of random effect in lme4 - Bayesian Filter detected spam

Actually, I re-read the post and think it needs clarification. We may both be right. If the question is I am building a model and want to know if I should retain this random effect? (or something like that) then the LRT should be used to compare the fitted model against another model. This would

[R] lme model: Error in MEEM

Hi, We have data of two groups of subjects: 32 elderly, 14 young adults. for each subject we have 15 observations, each observation consisting of a reaction-time measure (RT) and an activation maesure (betadlpcv). since we want to analyze the influence of (age-)group and RT on the activation, we

[R] select previous date

Hi everybody, could anyone help me in finding a way for selecting from a dataframe each row that comes before another that matches a condition? EXAMPLE: df raw.number name date temp 1aaa 2001-04-15 15 2aaa 2001-01-15 12 3aaa 2001-03-15 13

Re: [R] lme model: Error in MEEM

On 8/18/05, Christoph Lehmann [EMAIL PROTECTED] wrote: Hi, We have data of two groups of subjects: 32 elderly, 14 young adults. for each subject we have 15 observations, each observation consisting of a reaction-time measure (RT) and an activation maesure (betadlpcv). since we want to analyze

[R] problems when installing R in Fedora core 4

Hi, I got a problem when installing R in Fedora core 4. When I ran .configure, it gave the following error message: configure: error: --with-x=yes (default) and X11 headers/libs are not available Could anyone tell me what's wrong? Am I missing some package in Fedora? Thanks a lot for your help.

[R] help with unknown function

Hello I am working on radio tracking data, with a short programme someone gave me and ... which should, supposedly, work ... In this programme, there is the function : getareahr(kern, levels = 95). But i cannot find any 'getareahr' in R ... could anyone help me? thanks! Agnès

Re: [R] accesing slots of S4 class in C code

On 8/17/05, Aniko Szabo [EMAIL PROTECTED] wrote: I am trying to use a custom S4 object in my C code and I cannot get the access to its slots working. The following is a toy example, but even this crashes. In R I have: setClass(pd, representation(data=numeric)) x - new(pd,

[R] display of a loess fitted surface

Good morning, I am Marta Colombo,student at Politecnico,Milan. I am studying local regression models and I am using loess function. My problem is that when I have a loess object I don't know how to display the fitted surface; in fact, while in S when you have a loess object you can see it

[R] display of a loess fitted surface

Good morning, I am Marta Colombo,student at Politecnico,Milan. I am studying local regression models and I am using loess function. My problem is that when I have a loess object I don't know how to display the fitted surface; in fact, while in S when you have a loess object you can see it

Re: [R] problems when installing R in Fedora core 4

You are probably missing the 'devel' package for x11 which includes header files. -roger Peter Yang wrote: Hi, I got a problem when installing R in Fedora core 4. When I ran .configure, it gave the following error message: configure: error: --with-x=yes (default) and X11 headers/libs are

[R] how to draw an ellipse

That's probably a stupid question, but I'm looking for a low-level command which plots ellipse, specifying only center and deformation axes. The purpose is to illustrate bivariates gaussians with 2D .95 confident levels without using any specific library. Thanxs for your help, Guillaume

Re: [R] accesing slots of S4 class in C code

Hi Aniko, On 17 Aug 2005, [EMAIL PROTECTED] wrote: I am trying to use a custom S4 object in my C code and I cannot get the access to its slots working. The following is a toy example, but even this crashes. res - .C(TestFun, pd.obj) I'm pretty sure you want to use the .Call interface for

Re: [R] problems when installing R in Fedora core 4

On 08/18/05 09:31, Peter Yang wrote: Hi, I got a problem when installing R in Fedora core 4. When I ran .configure, it gave the following error message: configure: error: --with-x=yes (default) and X11 headers/libs are not available Could anyone tell me what's wrong? Am I missing some

Re: [R] Error messages using LMER

On 8/18/05, Shige Song [EMAIL PROTECTED] wrote: Dear All, After playing with lmer for couple of days, I have to say that I am amazed! I've been using quite some multilevel/mixed modeling packages, lme4 is a strong candidate for the overall winner, especially for multilevel generzlized

Re: [R] problems when installing R in Fedora core 4

On Thu, 2005-08-18 at 09:31 -0400, Peter Yang wrote: Hi, I got a problem when installing R in Fedora core 4. When I ran .configure, it gave the following error message: configure: error: --with-x=yes (default) and X11 headers/libs are not available Could anyone tell me what's wrong? Am I

Re: [R] how to draw an ellipse

Guillaume Allain guillaume.allain at cbconseil.com writes: That's probably a stupid question, but I'm looking for a low-level command which plots ellipse, specifying only center and deformation axes. The purpose is to illustrate bivariates gaussians with 2D .95 confident levels without

Re: [R] lme model: Error in MEEM

sorry, RT had an error in raw data and was treated as a factor. after correction of the raw data (RT is numeric) now it works fine. thanks a lot christoph --- Ursprüngliche Nachricht --- Von: Douglas Bates [EMAIL PROTECTED] An: Christoph Lehmann [EMAIL PROTECTED] Kopie:

Re: [R] trouble with wilcox.test

I think your guess about stack overflow is probably correct and I definitely don't think it's worth wasting effort recoding. Peter Ehlers Prof Brian Ripley wrote: If this is stack overflow (and I don't know that yet: when I tried this on Windows the traceback was clearly corrupt, referring

Re: [R] help with unknown function

Agnes Gault a écrit : Hello I am working on radio tracking data, with a short programme someone gave me and ... which should, supposedly, work ... In this programme, there is the function : getareahr(kern, levels = 95). But i cannot find any 'getareahr' in R ... could anyone help me?

Re: [R] help with unknown function

Agnes Gault wrote: Hello I am working on radio tracking data, with a short programme someone gave me and ... which should, supposedly, work ... In this programme, there is the function : getareahr(kern, levels = 95). But i cannot find any 'getareahr' in R ... could anyone help me?

Re: [R] how to draw an ellipse

Le 18.08.2005 15:45, Guillaume Allain a écrit : That's probably a stupid question, but I'm looking for a low-level command which plots ellipse, specifying only center and deformation axes. The purpose is to illustrate bivariates gaussians with 2D .95 confident levels without using any specific

[R] axTicks and window resizing

Dear listers, I have written a function to facilitate the drawing of altitude profiles with x (distance), y (altitude) and a z parameter (altitude magnification). profplot-function(x,y,z=10,...){ op - par()$mai par(mai=c(0.95625,0.76875,0.76875,0.95625)) plot(x,y*z,

[R] How do I make a Sweave + latex table out of this ?

Dear list, I have a table that I would like to convert to latex for inclusion into a Sweave file. round(ftable(prop.table(xtabs(~agemF + votcat + Type , data=work),margin=2))*100,1) Type Voiced Voiceless unaspirated Voiceless aspirated agemF votcat

Re: [R] do glm with two data sets

Thanks for your help. # read the two data sets e - as.matrix(read.table(file1.txt, header=TRUE,row.names=1)) g - as.matrix(read.table(file2.txt, header=TRUE,row.names=1)) # solution d1-data.frame(g[1,], e[1,]) fit-glm(e[1,] ~ g[1,], data=d1) summary(fit) I am not sure that is the best

Re: [R] lme model: Error in MEEM

Christoph Lehmann christoph.lehmann at gmx.ch writes: Hi, We have data of two groups of subjects: 32 elderly, 14 young adults. for each subject we have 15 observations, each observation consisting of a reaction-time measure (RT) and an activation maesure (betadlpcv). since we want to

Re: [R] do glm with two data sets

Hu, Ying (NIH/NCI) wrote: Thanks for your help. # read the two data sets e - as.matrix(read.table(file1.txt, header=TRUE,row.names=1)) g - as.matrix(read.table(file2.txt, header=TRUE,row.names=1)) # solution d1-data.frame(g[1,], e[1,]) fit-glm(e[1,] ~ g[1,], data=d1) summary(fit)

Re: [R] do glm with two data sets

On Thu, 2005-08-18 at 10:38 -0400, Hu, Ying (NIH/NCI) wrote: Thanks for your help. # read the two data sets e - as.matrix(read.table(file1.txt, header=TRUE,row.names=1)) g - as.matrix(read.table(file2.txt, header=TRUE,row.names=1)) # solution d1-data.frame(g[1,], e[1,]) This is

Re: [R] how to draw an ellipse

Romain Francois wrote: Why don't you want to use any specific library ? You can't reinvent the wheel !! True, but Madonna is always reinventing herself... There is a package ellipse on CRAN which will do what you are looking for. Have you tried RSiteSearch(ellipse) its simple

Re: [R] how to draw an ellipse

Congratulations guys, you made my day better and funnier! Le 18 août 05, à 17:31, Barry Rowlingson a écrit : Romain Francois wrote: Why don't you want to use any specific library ? You can't reinvent the wheel !! True, but Madonna is always reinventing herself... There is a package

[R] How to do log normal regression?

I want to fit a Log-Normal CDF function between two variables, and estimate the parameters. Is there any package/functions designed for this purpose? Basically, I have data for Y and X, and I suspect the relationship between Y and X is Y = CDF Log-Normal (X), and I want to run this regression

Re: [R] display of a loess fitted surface

On Thu, 18 Aug 2005, Marta Colombo wrote: I am Marta Colombo,student at Politecnico,Milan. I am studying local regression models and I am using loess function. My problem is that when I have a loess object I don't know how to display the fitted surface; in fact, while in S when you have a

[R] standard errors for expression intensities

Hello! When I use the functions rma() or justrma() in the Bioconductor affy package, the standard errors for the expression estimates computed by the function se.exprs() is a matrix of NA's. I am wondering if any of you have encountered this problem and if there is a solution. Any help would be

[R] Regular expressions sub

Dear all, I am struggling with the use of regular expression. I got as.character(test$sample.id) [1] 1.11 10.11 11.11 113.31 114.2 114.3 114.8 and need [1] 11 11 11 31 2 3 8 I.e. remove everything before the . . TIA, Bernd __

Re: [R] select previous date

You can try this: dates - df$date[df$temp==15] one.month.before - sapply(strsplit(dates, -), function(x) paste(x[1], sprintf(%02d, as.numeric(x[2])-1), x[3], sep=-)) df[df$date %in% one.month.before,] Ben On 8/18/05, alessandro carletti [EMAIL PROTECTED] wrote: Hi everybody, could anyone help

Re: [R] axTicks and window resizing

On Thu, 18 Aug 2005, Patrick Giraudoux wrote: Dear listers, I have written a function to facilitate the drawing of altitude profiles with x (distance), y (altitude) and a z parameter (altitude magnification). profplot-function(x,y,z=10,...){ op - par()$mai

Re: [R] do glm with two data sets

You are right. # read the two data sets e - as.matrix(read.table(file1.txt, header=TRUE,row.names=1)) g - as.matrix(read.table(file2.txt, header=TRUE,row.names=1)) # solution 2 summary(glm(e[1,] ~ g[1,])) summary(glm(e[1,] ~ g[2,])) ... They work very well. If I put it in the loop, such as for

Re: [R] Regular expressions sub

x - scan(clipboard, what=) Read 7 items x [1] 1.11 10.11 11.11 113.31 114.2 114.3 114.8 gsub([0-9]*\\., , x) [1] 11 11 11 31 2 3 8 Bernd Weiss wrote: Dear all, I am struggling with the use of regular expression. I got as.character(test$sample.id) [1] 1.11 10.11

Re: [R] Regular expressions sub

One solution is test - c(1.11,10.11,11.11,113.31,114.2,114.3) id - unlist(lapply(strsplit(test,[.]),function(x) {x[2]})) -Original Message- From: Bernd Weiss [mailto:[EMAIL PROTECTED] Sent: Thursday, August 18, 2005 12:10 PM To: r-help@stat.math.ethz.ch Subject: [R] Regular

Re: [R] do glm with two data sets

Hu, Ying (NIH/NCI) wrote: You are right. # read the two data sets e - as.matrix(read.table(file1.txt, header=TRUE,row.names=1)) g - as.matrix(read.table(file2.txt, header=TRUE,row.names=1)) # solution 2 summary(glm(e[1,] ~ g[1,])) summary(glm(e[1,] ~ g[2,])) ... They work very well.

Re: [R] line/bar through median in lattice's bwplot?

On 8/18/05, Dieter Menne [EMAIL PROTECTED] wrote: M. K. mk_lists at yahoo.ca writes: Is there a way to render a line through the median point in the boxplot generated by the Lattice command bwplot? The line basically bisects the bar at the median point... bwplot(height~voice.part

Re: [R] problems when installing R in Fedora core 4

On Thu, 18 Aug 2005, Jonathan Baron wrote: ... Whether you have the needed files already depends on what kind of installation you did. Some of the packages are devel and others are compat. Here is my list of compat rpms that I have installed, and I think I installed all of these just to

Re: [R] Regular expressions sub

Bernd Weiss bernd.weiss at uni-koeln.de writes: I am struggling with the use of regular expression. I got as.character(test$sample.id) [1] 1.11 10.11 11.11 113.31 114.2 114.3 114.8 and need [1] 11 11 11 31 2 3 8 I.e. remove everything before the . . Define the dot as

Re: [R] line/bar through median in lattice's

Dieter Menne dieter.menne at menne-biomed.de writes: Check the code for something like points (which is the default). Find a mysterious '|' and pch Dieter Yes, that was it! Thanks Dieter! The mysterious code is right near the bottom of the function. Too bad this is not documented, at

[R] Set R 2.1.1. in English

(whish R 2.1.1 could be parametrise 'English' even with a French Windows XP) If I understand you correctly, it can. Just add LANGUAGE=en to the shortcut. Wonderful hope but not sure to catch what you term shortcut. I tried to add this command in C:\R\rw2011\etc\Rprofile, the

Re: [R] Welcome to the R-help mailing list (Digest mode)

Hi, thank you for cooperation. --- Dr. (UA) Tetyana Stepanchuk Department of Electronic Commerce Johann Wolfgang Goethe-University Mertonstrasse 17 D-60054, Frankfurt/Main Germany phone: +49 069 798 22379 http://www.ecommerce.wiwi.uni-frankfurt.de

Re: [R] axTicks and window resizing

OK. Things work now and the window can be resized easy after adding at = axTicks() in axis() explicitely. This makes: profplot-function(x,y,z=10,...){ op - par()$mai par(mai=c(0.95625,0.76875,0.76875,0.95625)) plot(x,y*z, type=l,asp=1,las=1,xlab=,ylab=,yaxt=n,...)

Re: [R] Error messages using LMER

Dear Professor Bates, Here is output R 2.1.1 produced with control = list(EMverbose = TRUE, msVerbose = TRUE). I am getting the new devel version and see what will hapen there: EM iterations 0 85289.766 ( 5407.13:

Re: [R] LyX and Sweave

I just wanted to point out that I was there first :) on the Lyx List (Nov 2004): http://www.mail-archive.com/lyx-users@lists.lyx.org/msg36262.html Perhaps somebody who is trying to put all of this together can benefit from both sets of explanations. pj [EMAIL PROTECTED] wrote: On Mon, 25

Re: [R] Error messages using LMER

Here is what happened using R-devel: - EM iterations 0 85289.766 ( 5407.13: 0.0815) ( 26134.4: 0.00387) 1 84544.322 ( 333.732: 0.137) ( 1462.32: 0.00934) 2 84515.108 ( 129.506: 0.0270) ( 446.306: 0.00481) 3 84514.519 (

[R] ?

Hello, I have a small problem with developing design matrix X, which I use in estimation the log-likelihood of a multinomial logit model. I have the data: number of observation - 289 number of choice alternative- 3 number of choice specific variables in matrix X -4 matrix X =289x4

Re: [R] Error messages using LMER

One thing to be noted: after switching to R-devel, even the simplest model can not converge. I always get this: - Warning messages: 1: optim or nlminb returned message See PORT documentation. Code (27) in: LMEopt(x = mer,

Re: [R] Error messages using LMER

Thanks for including all of that output. I believe that in this version the parameters are the relative variances. This would indicate that somehow you are getting fits with very low residual sums of squares in the weighted least squares problem. It could be that you have too many fixed

Re: [R] Set R 2.1.1. in English

On Thu, 18 Aug 2005, Patrick Giraudoux wrote: (whish R 2.1.1 could be parametrise 'English' even with a French Windows XP) If I understand you correctly, it can. Just add LANGUAGE=en to the shortcut. Wonderful hope but not sure to catch what you term shortcut. I tried to add this

Re: [R] standard errors for expression intensities

Please ask Bioconductor Qs on their mailing list! PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html This is indeed covered there, as is the (mis)use of HTML mail. On Thu, 18 Aug 2005, Devarajan, Karthik wrote: Hello! When I use the functions rma() or justrma()

[R] problem with repeated formal arguments and ...

I want to add an argument if it is not present. Following the Green Book, p. 337: test - function(x, ...){ dots - list(...) if (!hasArg(from)) from - 0 else from - dots$from curve(x, from=from, ...) } test(sin) test(sin, from=4) Error in curve(x, from = from, ...) : formal

Re: [R] How do I make a Sweave + latex table out of this ?

Hi Fredrik, You can do this nicely with latex() if you massage this a little so that you end up with a data frame that includes a column for agemF (to indicate which agemF each row belongs to) and then use the rgroup and n.rgroup options of latex(). Dave Fredrik Karlsson wrote: Dear list,

Re: [R] Set R 2.1.1. in English

OK. ¨Perfectly clear. This written for other listers: Shortcut is the shortcuts (alias) everybody use to launch applications (in French: raccourci), and for suckers (as I am): one just have to right-click the shortcut used to launch Rgui.exe, select properties (propriétés in French) and add

Re: [R] problem with repeated formal arguments and ...

Ross Boylan wrote: I want to add an argument if it is not present. Following the Green Book, p. 337: test - function(x, ...){ dots - list(...) if (!hasArg(from)) from - 0 else from - dots$from curve(x, from=from, ...) } test(sin) test(sin, from=4) Error in

Re: [R] trouble with wilcox.test

Ok, I will think more about the appropriateness of the Wilcoxon test here. I was using R version 2.1.1, 2005-06-20 Windows XP SP2 512MB RAM --Greg - Original Message - From: Prof Brian Ripley [EMAIL PROTECTED] To: Greg Hather [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch Sent:

Re: [R] Regular expressions sub

Dirk Eddelbuettel [EMAIL PROTECTED] writes: Bernd Weiss bernd.weiss at uni-koeln.de writes: I am struggling with the use of regular expression. I got as.character(test$sample.id) [1] 1.11 10.11 11.11 113.31 114.2 114.3 114.8 and need [1] 11 11 11 31 2 3 8

[R] Console

I am at my first steps with R... and I already notice that the console has a quite limited number of lines. Can anyone tell me how to visualise all the information, which is actually present? I only see the last part of the output, which obviosly exceeds the maximum number of rows in the

[R] How to put factor variables in an nls formula ?

Hello, I want to fit a Gompertz model for tree diameter growth that depends on a 4 levels edaphic factor (‘Drain’) and I don’t manage to introduce the factor variable in the formula. Dinc is the annual diameter increment and D is the Diameter. treestab Dinc D Drain [1,] 0.03

Re: [R] problem with repeated formal arguments and ...

Sundar Dorai-Raj wrote: Ross Boylan wrote: I want to add an argument if it is not present. Following the Green Book, p. 337: test - function(x, ...){ dots - list(...) if (!hasArg(from)) from - 0 else from - dots$from curve(x, from=from, ...) } test(sin) test(sin, from=4)

[R] Set R 2.1.1. in English

Another way is to create a file called Renviron.site with the following single line: LANGUAGE=en and put it in your ...\R\rw2011\etc folder. In your case, echo LANGUAGE=en c:\R\rw2011\etc\Renviron.site would do it. This would result in R being in English not only if you (1) used the

[R] A. Mani : Avoiding loops

Hello, I want to avoid loops in the following situation. There is a 5-col dataframe with col headers alone. two of the columns are non-numeric. The problem is to calculate statistics(scores) for each element of one column. The functions depend on matching in the other non-numeric

[R] Use of contains in S4 classes

setClass(B, representation=representation(B, extra=numeric)) setClass(B, representation=representation(extra=numeric), contains=B) Are these the same? If not, how do they differ? What about setClass(B, representation=representation(B, extra=numeric), contains=B) ? As far as I

[R] R equivalent to `estimate' in SAS proc mixed

Example: I have the following model model - lmer(response ~ time * trt * bio + (time|id), data = dat) where time = time of observation trt = treatment group (0-no treatment / 1-treated) bio = biological factor (0-absent / 1-present) and I would like to obtain an

Re: [R] Console

Quoting Daniela Salvini [EMAIL PROTECTED]: I am at my first steps with R... and I already notice that the console has a quite limited number of lines. Can anyone tell me how to visualise all the information, which is actually present? I only see the last part of the output, which obviosly

Re: [R] Use of contains in S4 classes

On Thu, 18 Aug 2005, Ross Boylan wrote: setClass(B, representation=representation(B, extra=numeric)) setClass(B, representation=representation(extra=numeric), contains=B) Are these the same? If not, how do they differ? What about setClass(B, representation=representation(B,

Re: [R] Use of contains in S4 classes

Oops, the second class should have been A in the examples. Corrected version: setClass(B, representation=representation(A, extra=numeric)) setClass(B, representation=representation(extra=numeric), contains=A) Are these the same? If not, how do they differ? What about setClass(B,

[R] 0/0, R segfaults

Hi, I noticed that when I was conducting some calculation involving finding correlation coeficients, R stopped abnormally. So I did some research, and find out that 0/0 was the culprit. For sure 0/0 is not a valid expression, but R should give a warning, an error msg or NaN instead of

[R] axis label justified

Hi, I am trying to make my axis labels left justified, and have used adj=0 in the axis() without success. Can anyone have a suggestion? axis(2,at=1:50,labels=paste('a',1:50,sep=''),las=2,cex.axis=0.5,adj=0,tck=0,mgp=c(3,0.5,0)) Thanks __

Re: [R] axis label justified

I note that the axis help seems to refer to padj. After playing around it is obvious that I don't know what is meant by this argument, so maybe I'm doing something wrong. My practical soultion is plot(1:50,axes = FALSE,ylab = ) axis(2,at = 1:50,labels = rep(,50),las = 2,padj = 0)

Re: [R] 0/0, R segfaults

On 18 August 2005 at 16:01, Xing Qiu wrote: | Hi, | | I noticed that when I was conducting some calculation involving | finding correlation coeficients, R stopped abnormally. So I did some | research, and find out that 0/0 was the culprit. For sure 0/0 is not | a valid expression, but R

[R] Problem with get.hist.quote() in tseries

When using get.hist.quote(), I find the dates are broken. This is with R 2.1.1 on Mac OS X `panther'. library(tseries) Loading required package: quadprog 'tseries' version: 0.9-27 'tseries' is a package for time series analysis and computational finance. See

Re: [R] Regular expressions sub

On 18 Aug 2005 at 21:17, Peter Dalgaard wrote: Dirk Eddelbuettel [EMAIL PROTECTED] writes: Bernd Weiss bernd.weiss at uni-koeln.de writes: I am struggling with the use of regular expression. I got as.character(test$sample.id) [1] 1.11 10.11 11.11 113.31 114.2 114.3 114.8