Re: [R] How to get commands history as a character vector instead of displaying them?

Thanks, Wacek and Romain! Your solutions worked very well in RGui and Rterm in interactive mode. My final purpose was to obtain the file name of the postscript device in Rweb (http://pbil.univ-lyon1.fr/Rweb/); now I found savehistory() would not work because Rweb was non-interactive. I didn't

[R] Convert Character to Date

Hello, I have a date in the format Year-Month Name (e.g. 1990-January) and R classes it as a character. I want to convert this character into a date format, but when I try as.Date(1990-January, %Y-%B), I get back NA. The function strptime also gives me NA back. Thanks.

Re: [R] Convert Character to Date

The yearmon class in the zoo package can represent year/months: library(zoo) ym - as.yearmon(1990-January, %Y-%B); ym [1] Jan 1990 as.Date(ym) [1] 1990-01-01 On Tue, Mar 31, 2009 at 3:14 AM, Bob Roberts quagmire54...@yahoo.com wrote: Hello,   I have a date in the format Year-Month Name

[R] Package candisc

Hi listers, I am working on an canonical discriminant analysis, but I am having some trouble to use the CANDISC function... I just installed the last R version and installed the package CANDISC... But, I am getting the following message because about a permission: The downloaded packages are in

[R] [R-pkgs] data.table is on CRAN (enhanced data.frame for time series joins and more)

Dear all, The data.table package was released back in August 2008. This email is to publicise its existence in response to several suggestions to do so. It seems I didn't send a general announcement about it at the time and therefore perhaps, not surprisingly, not many people know about it.

[R] Convert date to integer

Hello, I have a dataframe containing dates, times and other parameters. The times have the format h:m, e.g. 13:00 or 5:30, R classes them as factors. I want to change the times to integers e.g. 13:00 - 1300. I tried to use chron to create a chronological object, but it didn't work for the times

[R] Efficient calculation of partial correlations in R

Hello, I'm looking for an efficient function for calculating partial correlations. I'm currently using the pcor.test () function, which is equivalent to the cor.test() function, and can receive only single vectors as input. I'm looking for something which is equivalent to the cor() function, and

[R] how to increase the limit for max.print in R

Hi All, I am using DNAcopy package in R for copy number analysis of 500K chip. The final output which I get from DNA copy is too big to be printed in a file. So I am getting an error as reached getOption(max.print) -- omitted 475569 rows Can somebody please provide me the pointers with

[R] Row/columns names within 'assign' command

Dear all, I am attempting to add row and column names to a series of tables (120 in total) which have 2 variable parts to their name. These tables are created as follows: # Create table indexes index - expand.grid(year = sprintf(%04d, seq(1986, 1995)), month = sprintf(%02d, 1:12)) #

[R] Does R support double-exponential smoothing?

I want to use double-exponential smoothing to forecast time series datas,but I couldn't find it in the document,does R support this method? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] how to increase the limit for max.print in R

On Tue, 2009-03-31 at 14:47 +0530, pooja arora wrote: Hi All, I am using DNAcopy package in R for copy number analysis of 500K chip. The final output which I get from DNA copy is too big to be printed in a file. So I am getting an error as reached getOption(max.print) -- omitted

Re: [R] summarize logical string

you need rle(), e.g., x - c(0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0) rle(x)$lengths I hope it helps. Best, Dimitris dba...@cnb.csic.es wrote: Hello everyone, I am a newbie, working on a gene clustering problem, and I have problems in summarizing a logical string into number of repeats of

Re: [R] summarize logical string

Thank you very much, rapid and very helpful. Best, Djordje you need rle(), e.g., x - c(0, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0) rle(x)$lengths I hope it helps. Best, Dimitris dba...@cnb.csic.es wrote: Hello everyone, I am a newbie, working on a gene clustering problem, and I

[R] summarize logical string

Hello everyone, I am a newbie, working on a gene clustering problem, and I have problems in summarizing a logical string into number of repeats of each value. In other words, how could I obtain from 0 1 1 1 0 0 0 0 1 1 0 1 0 0 this: 1 3 4 2 1 1 2 so a string that gives me the number of

Re: [R] To save Trellis Plots on A3 size paper (Portrait and Landscape)

Debabrata Midya Debabrata.Midya at commerce.nsw.gov.au writes: I like to save Trellis Plots on A3 size paper (Portrait and Landscape). Since a3 is not among the paper choices, you could give the width and height in inches (b...) explicitly. Dieter

Re: [R] To save Trellis Plots on A3 size paper (Portrait and Landscape)

Debabrata Midya wrote: I like to save Trellis Plots on A3 size paper (Portrait and Landscape). Since a3 is not among the paper choices, you could give the width and height in inches (b...) explicitly. Dieter -- View this message in context:

Re: [R] unicode only works with a second one

Hi Thomas, I get the zodiac symbols okay without the trailing \u3030 on FC9 Linux and R-2.8.1. Perhaps it's only on Windows. Might try it at work tomorrow where I can boot into WinXP. Jim __ R-help@r-project.org mailing list

Re: [R] how to increase the limit for max.print in R

Thanks, it Worked. Do you have any idea how much is the max limit for max.print Thanks and Regards, Pooja Arora | Associate Software Engineer- Washu-caB2B | Persistent Systems Limited pooja_ar...@persistent.co.in | Cell: +91 9923223668 | Tel: +91 (20) 3023 4007 Innovation in software product

[R] Programming objects in R

I am planning a project where an object-oriented approach would be appropriate, and for a number of reasons I would prefer using either Python or R. My problem at the moment is to find out how to do OO programming in R. Are there any introductory texts anywhere ? Tom --

Re: [R] Convert date to integer

thoeb t.hoebinger at gmail.com writes: Hello, I have a dataframe containing dates, times and other parameters. The times have the format h:m, e.g. 13:00 or 5:30, R classes them as factors. Probably you have read in the data from a file with read.table; check stringsAsFactors in the docs to

[R] Nonparametric analysis of repeated measurements data with sm library

Dear all, Does anybody know how to get more evaluation points in performing Nonparametric analysis of repeated measurements data with sm library. The following command gives the estimation for 50 points, but I would like to increase to 100 points But I do not know how to do that.

[R] Nonparametric analysis of repeated measurements data with sm library

Dear all, Does anybody know how to get more evaluation points in performing Nonparametric analysis of repeated measurements data with sm library. The following command gives the estimation for 50 points, but I would like to increase to 100 points But I do not know how to do that.

[R] plot Geotiff

Is there a way to plot Geotiff, combining it with spatial (vector) data? Paulo E. Cardoso [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Programming objects in R

I think that 'Software for Data Analysis' by John Chambers should certainly be on the reading list. Patrick Burns patr...@burns-stat.com +44 (0)20 8525 0696 http://www.burns-stat.com (home of The R Inferno and A Guide for the Unwilling S User) Tom Backer Johnsen wrote: I am planning a project

[R] error during DPpackage compilation

Dear All, I've had trouble compiling DPpackage as a user in one system. It works fine as root in other machines. I can see any clues in error messages My guess is that it is a permissions matter. Any help is appreciated. OS: Linux Kernel: 2.6.27 SMP Arch: Intel 64 bits gfortran not available

Re: [R] Programming objects in R

The third chapter of the book R Programming for Bioinformatics by Robert Gentleman is on object-oriented programming in R. Dan Viar Chesapeake, VA On Tue, Mar 31, 2009 at 6:15 AM, Patrick Burns pbu...@pburns.seanet.com wrote: I think that 'Software for Data Analysis' by John Chambers should

[R] R Language Bloggers and Web Sites Owners

Dear List, Apologies in advance for non bloggers for the spam on the slightly off topic. This is an invitation to all R language bloggers to help spread the world at a forum where leading authors come together. Current only 1 R blogger is there- David Smith, from Revolution Computing. There

Re: [R] which rows are duplicates?

Wacek Kusnierczyk wrote: Michael Dewey wrote: At 05:07 30/03/2009, Aaron M. Swoboda wrote: I would like to know which rows are duplicates of each other, not simply that a row is duplicate of another row. In the following example rows 1 and 3 are duplicates. x - c(1,3,1) y

Re: [R] Programming objects in R

Patrick Burns wrote: I think that 'Software for Data Analysis' by John Chambers should certainly be on the reading list. It seems that that is true. I have a number of R books, but not that one. Thanks for the suggestion. Tom Patrick Burns patr...@burns-stat.com +44 (0)20 8525 0696

[R] Jarque-Bera test and Ljung-Box test for multivariate time series

Hi! I know that there is function in fBasics package for univariate Jarque-Bera test and a funtion for univariate Ljung-Box test in stats package. But I am wondering if there is a function somewhere to do the tests for multivariate time series? Thanks, John [[alternative HTML version

Re: [R] R Language Bloggers and Web Sites Owners

Ajay ohri wrote: Dear List, Apologies in advance for non bloggers for the spam on the slightly off topic. This is an invitation to all R language bloggers to help spread the world at a forum where leading authors come together. Current only 1 R blogger is there- David Smith, from

[R] R commands equivalent for STATA commands

Hello I need help to convert a few commands of STATA to R, can someone help me for these ones: 1) stsum, by(variable) 2) sts list, by(variable) 3) stset Time, failure(Status) Also where can I find R commands equivalent to STATA commands. Thank you in advance Sunita -- View this message in

Re: [R] unicode only works with a second one

Thanks for the feedback. I did now try Vista (2.8.1), XP (2.9.0alpha) and Win2000 (2.8.1) and non did work compeletely, only on Vista/2.8.1 I got some symbols if I add the other sign. I will try the my.symbols later today, when i am at home. Thomas __

Re: [R] Factor Analysis Output from R and SAS

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi, I ran factor analysis using R and SAS. However, I had different outputs from R and SAS. Why they provide different outputs? Especially, the factor loadings are different. I did real dataset(n=264), however, I had an extremely different from

Re: [R] Efficient calculation of partial correlations in R

Dear Schragi, There's a function named partial.cor() in the Rcmdr package, but it's so simple that I'll just reproduce it here: partial.cor - function (X, ...) { R - cor(X, ...) RI - solve(R) D - 1/sqrt(diag(RI)) R - -RI * (D %o% D) diag(R) - 0 rownames(R) - colnames(R)

Re: [R] Mapping in R

Kelsey Scheitlin kns07g at fsu.edu writes: Hi, I am looking for a specific mapping capability in R that I can't seem to find, but think exists. I would like to make a border of a map have alternating black and white squares instead of the common latitude and longitude grid. (example:

[R] R for Mac OS X Search in Help

Hi: A small query about R for Mac on OSX: when looking at the help system screens, it does not seem to be possible to search within them. For example, when looking at a long help page (e.g. for par) I often want to search for a particular parameter. But I don't see how. I just have to scroll

[R] scope of variables in R

I need to allocate (using C nomenclature) a set of global variables, some integer scalars and some double vectors. I have placed the name of such variables in the file containing the main script and the called functions: #

[R] Factor Analysis Output from R and SAS

Dear Users, I ran factor analysis using R and SAS. However, I had different outputs from R and SAS. Why they provide different outputs? Especially, the factor loadings are different. I did real dataset(n=264), however, I had an extremely different from R and SAS. Why this things happened? Which

Re: [R] Factor Analysis Output from R and SAS

Dear TY, Considering that you used different methods -- maximum-likelihood factor analysis in R and principal components analysis in SAS -- the results are quite similar (although the three rotated factors/components come out in different orders). I hope this helps, John -Original

Re: [R] which rows are duplicates?

Wacek Kusnierczyk wrote: Wacek Kusnierczyk wrote: Michael Dewey wrote: At 05:07 30/03/2009, Aaron M. Swoboda wrote: I would like to know which rows are duplicates of each other, not simply that a row is duplicate of another row. In the following example rows 1 and 3 are duplicates.

Re: [R] which rows are duplicates?

Dimitris Rizopoulos wrote: Wacek Kusnierczyk wrote: Wacek Kusnierczyk wrote: Michael Dewey wrote: At 05:07 30/03/2009, Aaron M. Swoboda wrote: I would like to know which rows are duplicates of each other, not simply that a row is duplicate of another row. In the following example

[R] Function rmultiregfp from package bayesm gone?

Hi, I have been using the function rmultiregfp from the bayesm package for a long time. In the current version of the package this function (and the corresponding init.rmultiregfp) is no longer available, does anyone know why? Can someone who has the old bayesm still on his computer please post

Re: [R] scope of variables in R

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi, A) you should give and return the variables through input/output parameter of the functions, e.g. pwtset - function(n, inout) { if (n 100){ inout$ncof - n } else { inout$ncof - n -100 } ... return(inout) } ford -

Re: [R] which rows are duplicates?

Dimitris Rizopoulos wrote: another approach (maybe a bit cleaner) seems to be: data - data.frame(x=sample(1:2, 5, replace=TRUE), y=sample(1:2, 5, replace = TRUE)) vals - do.call('paste', c(data, sep = '\r')) data$class - match(vals, unique(vals)) data I have tried benchmarking it.

[R] How to specify axis interval unit...

Dear R Users: Greetings! Is there any way to specify the axis interval unit for barplots in R.? Any help is deeply appreciated. Thanks in advance! Prasanth VP, Global Manager - Biometrics, Delta Technology Management Services Pvt Ltd, Plot No: 13/2, Sector - I, Third Floor, HUDA

[R] 3d cloud plot with point size reflecting variable value

Hello all, I'm using the cloud function to plot the iris data as per the document: http://www.stat.ucl.ac.be/ISpersonnel/lecoutre/stats/fichiers/_gallery.pdf I'd like to change the point size to reflect a fourth variable, as done here

[R] naming vectors/ matrices in R

I would like to create a series of vectors in a loop and name them with letters. I know the command letters gives me the 26 letters of the alphabet. For example, I have a dataframe, called A, with 6 columns. B-length(A) C-letters(1:B) [1] a b c d e f I would like to extract the first letter

[R] Doubt(load large amount of data) -reg

Dear Developers, I want to import the large size (566.639 kb) of data in R. I will try through txt, excel,xml,in this way I successes only small amount of data .but I want to load or import large amount of data in R .so please give me Suggestion Shanthi Rangasamy Ph. D. Student Electronics

Re: [R] which rows are duplicates?

Dimitris Rizopoulos wrote: Wacek Kusnierczyk wrote: Wacek Kusnierczyk wrote: Michael Dewey wrote: At 05:07 30/03/2009, Aaron M. Swoboda wrote: I would like to know which rows are duplicates of each other, not simply that a row is duplicate of another row. In the following example

Re: [R] Mapping in R

It is relatively easy to do a small function for this in R: a) select the desired width for your border b) through a fix number or through the use of pretty(), establish the number of rectangles you want along your border. Then e.g., you can use rect() you build the sequence of black-white

Re: [R] RMySQL compile

Try this: PKG_CPPFLAGS=-I/usr/mysql/5.0/include/mysql -I/opt/csw/include -I/usr/sfw/include -I/usr/X11/include-I/usr/lib/gtk-2.0/include -I/usr/local/include PKG_LIBS=-L/usr/sfw/lib -L/usr/lib -L/usr/local/lib -L/opt/csw/lib -lmysqlclient R CMD INSTALL RMySQL I had to

[R] digits in round()

Hi Folks, Compare print(1234567890,digits=4) # [1] 1.235e+09 print(1234567890,digits=5) # [1] 1234567890 Granted that digits: a non-null value for 'digits' specifies the minimum number of significant digits to be printed in values. how does R decide to switch from the

Re: [R] two monitors

Hello, I work with Windows XP and do what you are hoping to do. Instead of using the built in script dialogue editor I use TINN-R. I can drag either this editor or the R Console to different monitors. I often have to make sure that the TINN-R editor panel is not maximized as that will hide any

Re: [R] R Language Bloggers and Web Sites Owners

Ajay ohri ohri2007 at gmail.com writes: Apologies in advance for non bloggers for the spam on the slightly off topic. This is an invitation to all R language bloggers to help spread the world at a forum where leading authors come together. Since the author has a notorious history of

Re: [R] Factor Analysis Output from R and SAS

As it was already pointed out by others, you used different methods (principal components in SAS vs. factor analysis in R). When you use the same method (+ varimax rotation) in both programs, there may still be a *small* difference: this comes from (possibly) different stopping criteria. In R, the

Re: [R] Jarque-Bera test and Ljung-Box test for multivariate time series

John, Tsay has code for multivariate Ljung-Box test statistic on http://faculty.chicagobooth.edu/ruey.tsay/teaching/bs41202/sp2009/ -Hannu On Tue, Mar 31, 2009 at 1:47 PM, John Seppänen john.seppane...@gmail.comwrote: Hi! I know that there is function in fBasics package for univariate

Re: [R] digits in round()

On 3/31/2009 8:59 AM, (Ted Harding) wrote: Hi Folks, Compare print(1234567890,digits=4) # [1] 1.235e+09 print(1234567890,digits=5) # [1] 1234567890 If printing in scientific notation takes up as much space or more than printing in regular notation, then regular notation will be used.

Re: [R] digits in round()

TH == Ted Harding ted.hard...@manchester.ac.uk on Tue, 31 Mar 2009 13:59:41 +0100 (BST) writes: TH Hi Folks, TH Compare TH print(1234567890,digits=4) TH # [1] 1.235e+09 TH print(1234567890,digits=5) TH # [1] 1234567890 TH Granted that TH digits: a non-null

Re: [R] RMySQL compile

Jeff, Thank you so much for responding. I am ready to do a lot of tweaking. From your message it looks like you are asking me to provide assignments on the same line as R CMD but precede the command. I am using bash and the usual thing would be to export those variables and then

Re: [R] R Language Bloggers and Web Sites Owners

Hm, this feels reminiscent of the Wellsphere free blog content scam: http://neurocritic.blogspot.com/2009/02/wealthcentral.html -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University Looking to arrange a meeting? Check my public calendar:

[R] labeling panels in lattice plots

I am using windows XP with R 2.8.1 I am generating a lattice plot of annual rain patterns using the following function: xyplot(rain.stats$min+ rain.stats$max + rain.stats$ave ~ rain.stats$month |rain.stats$year, lty = 1, data = rain.stats, type = c(l,l, l), col = c(red, blue,

[R] help about: anova and population no normal

Dear R-Helpers, Parametric statistics are statistics where the population is assumed to fit any parametrized distributions (most typically the normal distribution). My problem is: 1) if my polulation is no normal 2) if the sample data of all replications and treatments were well fitted from the

[R] Upgrade through the R interface?

Why can we not update R through a update.version() command from within R? This may be my ignorance about how software distribution works, but I was interested. -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do

Re: [R] labeling panels in lattice plots

would you mind creating a dummy data.frame or maybe dput() a smaller subset of the data frame? thanks Stephen Sefick On Tue, Mar 31, 2009 at 9:28 AM, steve_fried...@nps.gov wrote: I am using windows XP with R 2.8.1 I am generating a lattice plot of annual rain patterns using the following

Re: [R] labeling panels in lattice plots

Try converting year to a factor xyplot(min + max + ave ~ month | factor(year), data = rain.stats, ...) Also, notice the inclusion of the data argument. HTH, --sundar On Tue, Mar 31, 2009 at 6:28 AM, steve_fried...@nps.gov wrote: I am using windows XP with R 2.8.1 I am generating a

Re: [R] Upgrade through the R interface?

On 3/31/2009 9:37 AM, stephen sefick wrote: Why can we not update R through a update.version() command from within R? This may be my ignorance about how software distribution works, but I was interested. Because nobody has contributed such a function. There are a few tricky issues in writing

Re: [R] labeling panels in lattice plots

Sundar, Your suggestion to convert year to a factor worked like a charm. Thank you very much. Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282

[R] target of assignment expands to non-language object

I'm running the following code numbers - 1:50 for (i in 1:50) assign(paste(model,numbers[i]),i)-(lm(temp$Overall.Scaled.Score~temp$raw.score)) where I want R to create 50 different models-1:50, but get the following error message target of assignment expands to non-language object. I've tried

Re: [R] Upgrade through the R interface?

I was just curious. For a couple of times a year building from the sources doesn't make this a super high priority for me to write such a function (assuming I have the skills to do this anyway) as the core team makes it quite easy to just download whatever I need. Thanks for all of the help,

Re: [R] How to generate natural cubic spline in R?

Hello: If B-splines will suffice, there are many capabilities in R for that. My favorite is the 'fda' package, but 'splines' and other packages are also good. The splinefun function in the base package returns a function to compute spline interpolations optionally using a

Re: [R] Upgrade through the R interface?

Hi. Em Terça 31 Março 2009, stephen sefick escreveu: I was just curious. I use a bash script [http://asaguiar.med.br/modules.php?name=Downloadsd_op=getitlid=41] that checks for availability of upgrades and, if it is the case, downloads, compiles and install the newest R version. It could be

[R] CV and GCV for finding smoothness parameter

I received an assignment that I have to do in R, but I'm absolutely not very good at it. The task is the following: http://www.nabble.com/file/p22804957/question8.jpg To do this, we also get the following pieces of code (not in correct order): http://www.nabble.com/file/p22804957/hints.jpg

[R] Problem with mathematical expression and loop

Hi to all I use a loop to plot 9 different histograms of 9 different transformations of one dataset (x). I want to label the histograms with the mathematical expression of each transformation (e.g. x^3). For that I've prepared a vector with the labeling names for expression. trans.expr -

[R] HELP WITH SEM LIBRARY AND WITH THE MODEL'S SPECIFICATION

Dear John I want to thank you so much for the help. I've solved my problem following your instructions. The matrix S of the cov was singular because i had few complete observation. I have included 2 wrongs item, because (in the majority of cases) when one is inputed, the other one is NA, so i

[R] portable R editor

Several messages last week debated about a portable R with Tinn-R. I made such an app just a few weeks ago. It's called R-TinnR-portables and is available here: http://otan.ecoledelapaix.org/spip.php?article102 Simply unzip the file downloadable under the name R-2.8.1TinnR-portables. You'll have

Re: [R] RMySQL compile

sten...@go.com wrote on 03/31/2009 08:11 AM: Jeff, Thank you so much for responding. I am ready to do a lot of tweaking. From your message it looks like you are asking me to provide assignments on the same line as R CMD but precede the command. I am using bash and the usual thing would be

[R] Connecting C# to R via rscproxy.dll does not return all results

Hi there, I am using http://www.nabble.com/Using-R.dll-in-.NET-IPC-td13265477.html this solution to hook a C# Application to R via rscproxy.dll. So far it all works fine, except for the fact that whenever a line in the result vector contains more than 1 return parameter (e.g. ttest´s

Re: [R] help about: anova and population no normal

Those with more formal statistical backgrounds may provide better advice, but in my own informal training I've come to wonder why parametric stats persist in the face of modern computing power. As I understand it, Fisher developed ANOVA as low-computation method of approximating the Randomization

Re: [R] Problem with mathematical expression and loop

On 3/31/2009 9:15 AM, Denis Aydin wrote: Hi to all I use a loop to plot 9 different histograms of 9 different transformations of one dataset (x). I want to label the histograms with the mathematical expression of each transformation (e.g. x^3). For that I've prepared a vector with the

Re: [R] naming vectors/ matrices in R

emj83 stp08emj at shef.ac.uk writes: I would like to create a series of vectors in a loop and name them with letters. I know the command letters gives me the 26 letters of the alphabet. For example, I have a dataframe, called A, with 6 columns. B-length(A) C-letters(1:B) [1] a b

Re: [R] scope of variables in R

On Tue, Mar 31, 2009 at 7:10 AM, mau...@alice.it wrote: I need to allocate (using C nomenclature) a set of global variables, some integer scalars and some double vectors. ... My question is: how can I have R interpreter allocate globalvariables visible and accessible by all the functions in

Re: [R] target of assignment expands to non-language object

On Tue, Mar 31, 2009 at 10:05 AM, Alina Sheyman alina...@gmail.com wrote: I'm running the following code numbers - 1:50 for (i in 1:50) assign(paste(model,numbers[i]),i)-(lm(temp$Overall.Scaled.Score~temp$raw.score)) a) What is the purpose of numbers? Why not write paste(model,i) b)

Re: [R] Nonparametric analysis of repeated measurements dat a with sm library

Alphonse Monkamg amonkamg at yahoo.fr writes: Does anybody know how to get more evaluation points in performing Nonparametric analysis of repeated measurements data with sm library. Try ... paramter ngrid: ngrid the number of divisions of the above interval to be considered. Default:

Re: [R] help about: anova and population no normal

Oh, and to answer your question more directly, the randomization test permits testing hypotheses using any metric, so scale shape are definitely testable. Typically one is interested in means, so on each iteration of the test loop one computes the group/condition means. However, it's simple to

[R] 'for Loop'

Hello, I'm trying to create a for loop for a data set, I have a list of results in this data set and I want to take the 1st two add them together and divide by the mean of the 1st to, then do the same for the 3rd and 4th values in the list and so on and each time return a value for the

[R] Deriving Samples from specific, not implemented PDF for a QQ-Plot

Dear All, this is my first post to the R-help, so please don't be too strict. My problem concerns a QQ-Plot: I want to show how well empirical samples match with a theoretical distribution. The theoretical distribution has got several parameters, but I made it to fit via ML.

Re: [R] 'for Loop'

Al, Is there any ID index for the pairs? For example, if the first pair can be labeled a, and second pair labeled b etc., then you can add an index column or you may already have such a column in your list. Then run aggregate(your.data.column, by=index.column, FUN=mean). Or you can just add an

Re: [R] 3d cloud plot with point size reflecting variable value

R User R User ruser2008 at googlemail.com writes: I'm using the cloud function to plot the iris data as per the document: http://www.stat.ucl.ac.be/ISpersonnel/lecoutre/stats/fichiers/_gallery.pdf I'd like to change the point size to reflect a fourth variable, as done here

Re: [R] Factor Analysis Output from R and SAS

At 8:17 AM -0400 3/31/09, John Fox wrote: Dear TY, Considering that you used different methods -- maximum-likelihood factor analysis in R and principal components analysis in SAS -- the results are quite similar (although the three rotated factors/components come out in different orders). I

[R] Using apply to get group means

Hi all, I'm trying to improve my R skills and make my programming more efficient and succinct. I can solve the following question, but wonder if there's a better way to do it: I'm trying to calculate mean by several variables and then put this back into the original data set as a new

[R] Reshape: 'melt' numerous objects

Dear R Users, I'm trying to use the reshape package to 'melt' my gridded data into column format. I've done this before on individual files, but this time I'm trying to do it on a directory of files (with variable file names) - therefore I have to also use the 'assign' command. I have come up

[R] USING MULTILEVEL PACKAGE AND WABA FUNCTION

Dear friends, this time I have a problem with using waba function. Firstly, I'll explain you my situation. In the survey a gruop of supervisors judge the dipendents of a company. One supervisor reported on more than one subordinate. Thus, I need to show that lack of independence is not a

[R] Bug in col2rgb?

col2rgb(#0079, TRUE) [,1] red 0 green0 blue 0 alpha 121 col2rgb(#0080, TRUE) [,1] red255 green 255 blue 255 alpha0 col2rgb(#0081, TRUE) [,1] red 0 green0 blue 0 alpha 129 Any ideas? Thanks, Hadley -- http://had.co.nz/

Re: [R] Using apply to get group means

Not exactly the output you asked for, but perhaps you can consider, library(doBy) summaryBy(x3~x2+x1,data=x,FUN=mean) x2 x1 x3.mean 1 1 A 1.5 2 1 B 2.0 3 1 C 3.5 4 2 A 4.0 5 2 B 5.5 6 2 C 6.0 the plyr package also provides similar functionality, as do

Re: [R] Using apply to get group means

On Tue, Mar 31, 2009 at 11:31 AM, baptiste auguie ba...@exeter.ac.uk wrote: Not exactly the output you asked for, but perhaps you can consider, library(doBy) summaryBy(x3~x2+x1,data=x,FUN=mean)  x2 x1 x3.mean 1  1  A     1.5 2  1  B     2.0 3  1  C     3.5 4  2  A     4.0 5  2  B    

Re: [R] Reshape: 'melt' numerous objects

On Tue, Mar 31, 2009 at 11:12 AM, Steve Murray smurray...@hotmail.com wrote: Dear R Users, I'm trying to use the reshape package to 'melt' my gridded data into column format. I've done this before on individual files, but this time I'm trying to do it on a directory of files (with variable

Re: [R] Using apply to get group means

A different solution (using aggregate for the table of means and merge for adding it to the dataframe): x1-rep(c(A,B,C),3) x2-c(rep(1,3),rep(2,3),1,2,1) x3-c(1,2,3,4,5,6,2,6,4) x-data.frame(x1,x2,x3) #here using data.frame the x1 variable is directly converted to factor x3means -

Re: [R] 'for Loop'

Al, Say, your data file is 'test', execute the following in sequence, aggregate(test[1],test[2],mean)-inter names(inter)[2]='mean' merge(test,inter,all=T)-inter2 inter2$RSV=inter2$Result/inter2$mean The column 'RSV' in inter2 should be what you want. Jun On Tue, Mar 31, 2009 at 11:11 AM,

Re: [R] Bug in col2rgb?

On 3/31/2009 12:29 PM, hadley wickham wrote: col2rgb(#0079, TRUE) [,1] red 0 green0 blue 0 alpha 121 col2rgb(#0080, TRUE) [,1] red255 green 255 blue 255 alpha0 col2rgb(#0081, TRUE) [,1] red 0 green0 blue 0 alpha 129 Any

  1   2   >