search for: adaikalavan

...10, main=paste(x, expression(Delta), "values")) plot(1:10, main=paste(x, eval(expression(Delta)), "values")) plot(1:10, main=paste(x, expression(Delta, "values" ))) plot(1:10, main=paste(x, expression(paste(Delta, "values" )))) Many thanks. Regards, -- Adaikalavan Ramasamy ramasamy at cancer.org.uk Centre for Statistics in Medicine http://www.ihs.ox.ac.uk/csm/ Cancer Research UK Tel : 01865 226 677 Old Road Campus, Headington, Oxford Fax : 01865 226 962

...: Object "formu" not found I think this has something to do with formula and environments but I do not know enough to solve it myself. I searched the archive without much help (perhaps I was using the wrong keywords). Any help would be very much appreciated. Thank you. Regards, -- Adaikalavan Ramasamy ramasamy at cancer.org.uk Centre for Statistics in Medicine http://www.ihs.ox.ac.uk/csm/ Wolfson College Annexe Tel : 01865 284 408 Linton Road, Oxford OX2 6UD Fax : 01865 284 424

...PCRE TRUE TRUE TRUE TRUE TRUE TRUE > plot( 1:10 ) Bus Error and R terminates to the shell prompt. I have searched the archives with no success. Can anyone kindly help me fix the problem or at least show me how to get more information about the error. Thank you. -- Adaikalavan Ramasamy ramasamy at cancer.org.uk Medical Statistician Centre for Statistics in Medicine http://www.ihs.ox.ac.uk/csm/ Medical Statistics Group Tel : 01865 226 677 Cancer Research UK Fax : 01865 226 962

Just out of curiosity, what is the difference between the terms for package and library ? Why are we loading a package with the library() command ? If this is a case of RTFM, I would be happy to do so if pointed in the right direction. I have searched the FAQ and mail archives and only came up with http://tolstoy.newcastle.edu.au/R/help/05/02/12162.html but this still does not explain what is

...sp <- ( (n1-1)*v1 + (n2-1)*v2 )/(n1 + n2 - 2) denom2 <- sp * sqrt( 1/n1 + 1/n2 ) num / denom2 # gives 0.5913777 I tested this using R-1.9.1 (21/06/2004) on Redhat Fedora Core 2 and Windows 2000 Professional with the same results. Any suggestions would be kindly appreciated. Regards, -- Adaikalavan Ramasamy ramasamy@cancer.org.uk Centre for Statistics in Medicine http://www.ihs.ox.ac.uk/csm/ Cancer Research UK Tel : 01865 226 677 Old Road Campus, Headington, Oxford Fax : 01865 226 962

thank you to both Adaikalavan and Patrick. on the basis of Adaikalavan example this is an example that point out my problems. >log_raise=function(num, exp, base){return(log(num^exp,base))} I would like to have optional fields, so some settings to be default parameters; if I declare the function as above the operator mus...

Hi there, Suppose mu is constant, and error is normally distributed with mean 0 and fixed variance s. I need to find a statistics that: Y_i = mu + beta1* I1_i beta2*I2_i + beta3*I1_i*I2_i + +error, where I_i is 1 Y_i is from group A, and 0 if Y_i is from group B. It is large when beta1=beta2=0 It is small when beta1 and/or beta2 is not equal to 0 How can I find it by R? Thank you very much

Dear all, I am trying to create a web interface using Perl-CGI to call R plots and to display them. The following codes works perfectly fine when I copy and paste into the console directly or if I save it into script.file and then R --no-save < script.file producing the graphs. jpeg("graph.jpeg", width=400, height=400) plot(rnorm(100)) dev.off() Now, I put the line system("R

Hi all! I 'm beginner and i develop a bio-application with VB and i need some statistic functions! could i calculate StdError, CoeffOfVariance, SumSquared with R langage? if yes, what are functions to use? I need also to use ANOVA and t-test... Thanks for your help! Laurent Houdusse Analyste Programmeur

Hello guys, I am wondering the default way of transferring arguments in R. Is it by value or by ref in default case, or could that be changed explicitly? Cheers, Xiaofan --- Xiaofan Li Department of Applied Mathematics and Theoretical Physics University of Cambridge

Hi all, If the question below as been answered before I apologize for the posting. I would like to get the frequencies of occurrence of all values in a given variable in a multivariate dataset. In short for each variable (or field) a summary of values contained with in a value:frequency pair, there can be many such pairs for a given variable. I would like to do the same for several such variables.

Hello, I am trying to use self created functions in other scripts than the one where they are stored. For the moment I am using the following structure of commands to do that: 1. Load the text file with the functions in the current script: x=parse("path") 2. transform the tex in a function: f1=eval(x[1]), f2=eval(x[2]) if more than one function is stored in the text file 3. use the

...ncorrelated). Even then you still have 80 variances to estimate. Besides, I don't think lda() does that. If you really have to discrinimate 40 data points with 80 variables, use methods that do not rely on the estimation of covariance matrix. Andy > -----Original Message----- > From: Adaikalavan Ramasamy [mailto:ramasamy at stats.ox.ac.uk] > Sent: Friday, July 12, 2002 4:08 PM > To: r-help at r-project.org; allstat at jiscmail.ac.uk > Subject: [R] meaning of error message about collinearity > > > Just a quick question. I am trying to fit an LDA model with a > restri...

Just out of curiosity, can some please explain the following return NA. x <- 6 var(x) y <- c( NA, NA, 10000 ) var(y, na.rm=T) Unless I am seriously misguided, I believe that the variance of a single number (i.e. a constant) should be zero. Thanks. Regards, Adai.

Thanks Adaikalavan, however the problem remains. Considering AIC() as applied to the linear model in AIC() help documentation: > data(swiss) > lm1 <- lm(Fertility ~ . , data = swiss) > AIC(lm1) [1] 326.0716 Clearly this includes the estimation of the residual standard error as an estimated parameter...

Dear all, I am trying to compile the package e1071 (version 1.3-11) with R CMD INSTALL. I tried with R 1.7.0 on Redhat Linux 2.4.7-10 and R 1.6.2 on Linux 2.4.9-34smp but keep getting the same error message during configure : WARNING: g++ 2.96 cannot reliably be used with this package. Please use a different compiler. Can anyone help me with this or at least point me in the right direction ?

Sorry if this is a question more on regular expressions. I am dealing with several files which have been badly named. For example the files are given either the extensions txt, TXT or Txt. I wish to select all those files ending with 'txt' ignoring case. Here is how I would do it in bash (Redhat FC2) : touch a1.txt a2.TXT a3.Txt txt.control TXT.control ls -1 | grep -i

I am comparing two different algorithms in terms of speed and memory usage. I can calculate the processing time with proc.time() as follows but am not sure how to calculate the memory usage. ptm <- proc.time() x <- rnorm(1000000) proc.time() - ptm I would like to be within R itself since I will test the algorithm several hundred times and in batch mode. So manually looking up

...ifference(mat) A.vs.B A.vs.C A.vs.D B.vs.C B.vs.D C.vs.D g1 -6 7 -5 13 1 -12 g2 13 6 -1 -7 -14 -7 g3 3 -1 -6 -4 -9 -5 g4 12 8 3 -4 -9 -5 g5 16 -3 8 -19 -8 11 Regards, -- Adaikalavan Ramasamy ramasamy at cancer.org.uk Centre for Statistics in Medicine http://www.ihs.ox.ac.uk/csm/ Cancer Research UK Tel : 01865 226 677 Old Road Campus, Headington, Oxford Fax : 01865 226 962

I really appreciate all the helpful answers from Richard Müller, Marc Schwartz, Adaikalavan Ramasamy, Vito Ricci. I tried Marc's simple suggestion which was to copy only the data area and paste to a new excel sheet, and that solved my problem. I'll try other suggestions while I'm learning more about R. Unedited answers are provided below followed by my question. Thanks again....