Re: [R] generating symmetric matrices

2007-07-30 Thread Benilton Carvalho
after dat.cor use: Rmat[lower.tri(Rmat)] - dat.cor Rmat - t(Rmat) Rmat[lower.tri(Rmat)] - dat.cor b On Jul 27, 2007, at 11:28 PM, Gregory Gentlemen wrote: Greetings, I have a seemingly simple task which I have not been able to solve today. I want to construct a symmetric matrix of

Re: [R] generating symmetric matrices

2007-07-30 Thread Douglas Bates
On 7/27/07, Gregory Gentlemen [EMAIL PROTECTED] wrote: Greetings, I have a seemingly simple task which I have not been able to solve today. I want to construct a symmetric matrix of arbtriray size w/o using loops. The following I thought would do it: p - 6 Rmat - diag(p) dat.cor -

Re: [R] xtable with vector

2007-07-30 Thread Benilton Carvalho
a - matrix(1:6, nr=1) colnames(a) - paste(col, 1:6) xtable(a) On Jul 28, 2007, at 12:39 PM, Stefan Nachtnebel wrote: Hello, Is there a possibility to use xtable with a vector to generate a latex table? I always get an error, that no applicable method is available. For example: b-1:12

Re: [R] generating symmetric matrices

2007-07-30 Thread John Logsdon
Since a symmetric matrix must be square, the following code does it: X-5 MAT-matrix(rnorm(X*X),X,X) MAT-MAT+t(MAT) MAT [,1] [,2] [,3] [,4] [,5] [1,] 0.1084372 -1.7867366 -0.9620313 -1.0925719 -0.5902326 [2,] -1.7867366 -0.0462097 -1.2707656 0.6112664

Re: [R] generating symmetric matrices

2007-07-30 Thread Gregory Gentlemen
Douglas Bates [EMAIL PROTECTED] wrote: On 7/27/07, Gregory Gentlemen wrote: Greetings, I have a seemingly simple task which I have not been able to solve today. I want to construct a symmetric matrix of arbtriray size w/o using loops. The following I thought would do it: p - 6 Rmat -

Re: [R] Bind together two vectors of different length...

2007-07-30 Thread Gabor Grothendieck
Using the zoo package create zoo objects and merge them. This code assumes that A and B each have unique and ascending values: library(zoo) A - seq(1, 10, 1) B - seq(1, 10, 2) t(merge(A = zoo(A,A), B = zoo(B,B), fill = 0)) 1 2 3 4 5 6 7 8 9 10 A 1 2 3 4 5 6 7 8 9 10 B 1 0 3 0 5 0 7 0 9 0

Re: [R] how to combine data of several csv-files

2007-07-30 Thread Antje
Hello, thank you for your help. But I guess, it's still not what I want... printing df.my gives me df.my v1 v2 v3 v4 v5 v6 v7 v8 1 NA -0.6442149 0.02354036 -1.40362589 -1.1829260 1.17099178 -0.046778203 NA 2 NA -0.2047012 -1.36186952

Re: [R] beta regressions in R

2007-07-30 Thread ronggui
see the package betareg in CRAN. 2007/7/29, Walter R. Paczkowski [EMAIL PROTECTED]: Good morning, Does anyone know of a package or function to do a beta regression? Thanks, Walt Paczkowski _ Walter R. Paczkowski, Ph.D. Data Analytics

Re: [R] Matrix nesting (was Re: Obtaining summary of frequencies of value occurrences for a variable in a multivariate dataset.)

2007-07-30 Thread Allan Kamau
Success, thanks Patrick. Below is the final matrix construction code. x=list() x[length(myVariableNames)]-NA names(x)-names(x.val) for (i in myVariableNames){ residues=names(x.val[[i]]) residuesFrequencies=as.vector(x.val[[i]]) someList=list() names(residuesFrequencies)=residues

Re: [R] the large dataset problem

2007-07-30 Thread Bernzweig, Bruce \(Consultant\)
Hi Eric, I'm facing a similar problem. Looking over the list of packages I came across: R.huge: Methods for accessing huge amounts of data http://cran.r-project.org/src/contrib/Descriptions/R.huge.html I haven't installed it yet so I don't know how well it works. I probably

Re: [R] how to combine data of several csv-files

2007-07-30 Thread 8rino-Luca Pantani
Ok. I missed the grouping factor Try this. You can modify my factor to fit your needs. As to avoid list, I cannot help, sorry I use them only when I have to collect different classes of objects. v1 - NA v2 - rnorm(6) v3 - rnorm(6) v4 - rnorm(6) v5 - rnorm(6) v6 - rnorm(6) v7 - rnorm(6) v8 -

Re: [R] text() and vector arguments like adj

2007-07-30 Thread Prof Brian Ripley
On Sat, 28 Jul 2007, [EMAIL PROTECTED] wrote: Hello, I remarked that the function ## Default S3 method: text (x, y = NULL, labels = seq(along = x), adj = NULL,pos = NULL, offset = 0.5, vfont = NULL,cex = 1, col = NULL, font = NULL, ...) accepts vectors of arguments (of the same length)

Re: [R] error in using R2WinBUGS on Ubuntu 6.10 Linux

2007-07-30 Thread J. Patrick Meyer
I'm using Wine 0.9.41, and the patched version of WinBUGS. Everything seems to run correctly. I get the right output from WinBUGS. I'm not sure to what the error message refers. Thanks, Patrick Thomas Harte wrote: what version of Wine are you running? and is this the patched (1.4.1) version

Re: [R] beta regressions in R

2007-07-30 Thread apjaworski
Walter, There is a betareg package in CRAN. Cheers, Andy __ Andy Jaworski 518-1-01 Process Laboratory 3M Corporate Research Laboratory - E-mail: [EMAIL PROTECTED] Tel: (651) 733-6092 Fax: (651) 736-3122

Re: [R] generating symmetric matrices

2007-07-30 Thread Charles C. Berry
On Fri, 27 Jul 2007, Gregory Gentlemen wrote: Greetings, I have a seemingly simple task which I have not been able to solve today. I want to construct a symmetric matrix of arbtriray size w/o using loops. The following I thought would do it: p - 6 Rmat - diag(p) dat.cor -

[R] how to install rattle() in R-GUI2.5.1

2007-07-30 Thread j.joshua thomas
Dear Group, I have installed R-GUI 2.5.1, thereby i tried install install.packages(RGtk2) install.packages (rattle) where i have problems such as below: install.packages(RGtk2) Warning: unable to access index for repository http://cran.au.r-project.org/bin/windows/contrib/2.5 Warning message:

Re: [R] Looping through all possible combinations of cases

2007-07-30 Thread Charles C. Berry
Lots of ways. Here is a simpler one. start by reading ?combn ?subset ?Subscript ?is.element ?apply ?paste - note the 'collapse' arg ?names - note the 'names(x) - value' usage ?list ?unlist then write a function that calc's

Re: [R] Matrix Multiplication, Floating-Point, etc.

2007-07-30 Thread Charles C. Berry
7.31 Why doesn't R think these numbers are equal? On Fri, 27 Jul 2007, Talbot Katz wrote: Hi. I recently tried the following in R 2.5.1 on Windows XP: ev2-c(0.8,-0.6) ev1-c(0.6,0.8) ev1%*%ev2 [,1] [1,] -2.664427e-17 sum(ev1*ev2) [1] 0 (I got the same result with R

Re: [R] beta regressions in R

2007-07-30 Thread Ben Bolker
Walter R. Paczkowski dataanalytics at earthlink.net writes: Does anyone know of a package or function to do a beta regression? RSiteSearch(\\\beta regression\\\) (all those extra s seem necessary to get the equivalent of beta AND regression ...) betareg package on CRAN

[R] duplicate DATE at in lattice scale handled differently from base graphics OR lattice numeric scales

2007-07-30 Thread Alex Brown
When at, label chains contain duplicate at values, axis ticks are dropped. I believe this is handled incorrectly for Date ats in lattice 0.15-4 when compared to how it is handled for numeric, or for dates in base plot. This can result in mis-labelled axes in some circumstances. I have

Re: [R] the large dataset problem

2007-07-30 Thread Ben Bolker
Eric Doviak edoviak at earthlink.net writes: Dear useRs, I recently began a job at a very large and heavily bureaucratic organization. We're setting up a research office and statistical analysis will form the backbone of our work. We'll be working with large datasets such the SIPP as well

[R] how to install rattle() with R-GUI2.5.1 ?

2007-07-30 Thread j.joshua thomas
i am using R-Gui2.5.1 where can i find rattle() ? previous version of R-Gui2.4.1 with rattle() works fine i tried with the new version and i could not locate the rattle() i follow the steps 1.install R-Gui-2.5.1 2. install gtk-2.10.7-win32-1 3. install gtk-dev-2.10.7-win32-1 4. install

Re: [R] generating symmetric matrices

2007-07-30 Thread Douglas Bates
On 7/30/07, Gregory Gentlemen [EMAIL PROTECTED] wrote: Douglas Bates [EMAIL PROTECTED] wrote: On 7/27/07, Gregory Gentlemen wrote: Greetings, I have a seemingly simple task which I have not been able to solve today. I want to construct a symmetric matrix of arbtriray size w/o using

Re: [R] getting the name of variables passed to a function

2007-07-30 Thread Liaw, Andy
Here's one possibility: R f - function(...) { call - match.call(); sapply(as.list(call[-1]), deparse) } R f(x, y) [1] x y R f(x=x, y=y) x y x y You basically need to know how to manipulate call objects. The relevant section in the R Language Definition should help. Andy From: Horace

[R] stop criteria when L-BFGS-B needs finite values of 'fn' in optim

2007-07-30 Thread Dae-Jin Lee
Hi all! I'm running some simulations and I need to estimate some paramaters with optim( ), in some cases optim stops with the next message: L-BFGS-B needs finite values of 'fn' I would like to know how to include and if condition when this happen, could it be

Re: [R] Slightly OT - use of R

2007-07-30 Thread Bert Gunter
Why? You might receive more useful replies from a relevant subset of users if you specify the purpose you have in mind. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of John Logsdon Sent: Monday, July 30,

[R] Release of an Ecological Modelling Software using R (with Rserve) and ImageJ

2007-07-30 Thread Bio7
Dear R users, With the hope that this application can also be useful for some R users, I like to announce the first release of Bio7. This application is an Integrated Development Environment for ecological modelling with a main focus on individual based modelling and spatially explicit models

Re: [R] lattice grayscale theme

2007-07-30 Thread Patrick Drechsler
Patrick Drechsler [EMAIL PROTECTED] writes: is there a grayscale setting for lattice plots? Solved using ## Set background color of strips to grayscales: strip.background - trellis.par.get(strip.background) trellis.par.set(strip.background = list(col = grey(7:1/8))) ## Set color of

Re: [R] the large dataset problem

2007-07-30 Thread Ted Harding
On 30-Jul-07 11:40:47, Eric Doviak wrote: [...] Sympathies for the constraints you are operating in! The Introduction to R manual suggests modifying input files with Perl. Any tips on how to get started? Would Perl Data Language (PDL) be a good choice? http://pdl.perl.org/index_en.html

[R] Piecewise Regression with a known slope

2007-07-30 Thread Jarrett Byrnes
Hey, all. I'm working on a data set with a broken stick linear regression where I know one of the two slopes. It is a negative linear function until the line intersects with the x-axis, at which point it becomes 0. It is not a nonlinear asymptotic function, and, indeed, using negative

Re: [R] generating symmetric matrices

2007-07-30 Thread Bert Gunter
See ?dist for an object oriented approach that may be better. Directly, you can do something like (see ?row ?col): x - matrix(NA, 10,10) ## Lower triangular : x[row(x) = col(x) ] - rnorm(55) x[row(x) col(x)] - x[row(x) col(x)] ## or you could have saved the random vector and re-used it.

Re: [R] Call R program from C++ code

2007-07-30 Thread Kisas
Hi Vladimir: Thank you very much! I'm very interested in this sentence: R can be compiled as a shared library object, that you can dynamically load from your application and use its functions. That's would be great if R programs can be compiled as shared library object(.lib? or .dll ? )

Re: [R] Matrix Multiplication, Floating-Point, etc.

2007-07-30 Thread Talbot Katz
Thank you for responding! I realize that floating point operations are often inexact, and indeed, the difference between the two answers is within the all.equal tolerance, as mentioned in FAQ 7.31 (cited by Charles): (as.numeric(ev1%*%ev2))==(sum(ev1*ev2)) [1] FALSE

Re: [R] 2nd R Console

2007-07-30 Thread Ted Harding
On 28-Jul-07 08:19:10, Michael Janis wrote: Hi, I was reading a thread: [R] 2nd R console and had a similar question regarding having more than one R console open at a time. However, my question differs from that of the thread: Is it possible, or is there a wrapper that will allow one,

Re: [R] Matrix Multiplication, Floating-Point, etc.

2007-07-30 Thread Doran, Harold
Talbot The general advice on this list is to read the following http://docs.sun.com/source/806-3568/ncg_goldberg.html -Original Message- From: Talbot Katz [mailto:[EMAIL PROTECTED] Sent: Monday, July 30, 2007 1:55 PM To: [EMAIL PROTECTED] Cc: r-help@stat.math.ethz.ch; Doran,

Re: [R] how to combine data of several csv-files

2007-07-30 Thread jim holtman
This should do it: v1 - NA v2 - rnorm(6) v3 - rnorm(6) v4 - rnorm(6) v5 - rnorm(6) v6 - rnorm(6) v7 - rnorm(6) v8 - rnorm(6) v8 - NA list - list(v1,v2,v3,v4,v5,v6,v7,v8) categ - c(NA,cat1,cat1,cat1,cat2,cat2,cat2,NA) # create partitioned list list.cat - split(list, categ) #

Re: [R] generating symmetric matrices

2007-07-30 Thread Cleber Borges
hello try ?upper.tri # example p - 6 Rmat - diag(p) dat.cor - rnorm(p*(p-1)/2) Rmat[upper.tri(Rmat)]- dat.cor Rmat[lower.tri(Rmat)]- dat.cor Cleber Borges Greetings, I have a seemingly simple task which I have not been able to solve today. I want to construct a

[R] Overlaying a single contour from a new data array in levelplot

2007-07-30 Thread Jenny Barnes
Dear Deepayan Thank you for your response - it has proved very very helpful, I can't thank you enough! I have another question for you if you have time to reply. I know you have been asked about the colour of the polygon outline before (27 April 2007) and you replied that is a bug and the

Re: [R] generating symmetric matrices

2007-07-30 Thread Ted Harding
On 28-Jul-07 03:28:25, Gregory Gentlemen wrote: Greetings, I have a seemingly simple task which I have not been able to solve today. I want to construct a symmetric matrix of arbtriray size w/o using loops. The following I thought would do it: p - 6 Rmat - diag(p) dat.cor -

[R] zoo, ts and data.frames

2007-07-30 Thread Edna Bell
Hi R gurus: I have some zoo objects that I have put into a data.frame. However, when I try to plot the objects from the data frame, the x axis is now Index rather than time Is there a sort of zoo data frame or a multiple zoo object, please? Thanks for any help! Sincerely, Edna [EMAIL

[R] filenames

2007-07-30 Thread Dong GUO 郭东
Dear all, I want to create filename from a loop, say, i=(1:10), and the final names will be file1.csv, file2.csv, file3.csv in Python, it seems easer as - file+str(i) +'.csv , but how can i do in R list = as.character(1:10) list [1] 1 2 3 4 5 6 7 8 9 10 [[alternative HTML

Re: [R] generating symmetric matrices

2007-07-30 Thread Ben Bolker
Gregory Gentlemen gregory_gentlemen at yahoo.ca writes: Greetings, I have a seemingly simple task which I have not been able to solve today. I want to construct a symmetric matrix of arbtriray size w/o using loops. The following I thought would do it: [snip] p - 6 Rmat - diag(p) vals

Re: [R] beta regressions in R

2007-07-30 Thread Tobias Verbeke
Walter R. Paczkowski wrote: Good morning, Does anyone know of a package or function to do a beta regression? http://cran.r-project.org/src/contrib/Descriptions/betareg.html HTH, Tobias -- Tobias Verbeke - Consultant Business Decision Benelux Rue de la révolution 8 1000 Brussels -

Re: [R] Constructing correlation matrices

2007-07-30 Thread Bert Gunter
See ?dist for an object oriented approach that may be better. Directly, you can do something like (see ?row ?col): x - matrix(NA, 10,10) ## Lower triangular : x[row(x) = col(x) ] - rnorm(55) x[row(x) col(x)] - x[row(x) col(x)] ## or you could have saved the random vector and re-used it.

[R] zoo, ts and data frames - SOLVED

2007-07-30 Thread Edna Bell
Sorry...there is a merge command that solves this. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal,

Re: [R] Constructing correlation matrices

2007-07-30 Thread Horace Tso
Greg, I take it that you're trying to generate a random correlation matrix, so first create a covariance matrix, p = 6 v = matrix(rnorm(p*p), ncol=p) cov = t(v) %*% v Then convert it to a correlation matrix, cov2cor(cov) HTH. Horace Gregory Gentlemen [EMAIL PROTECTED] 7/29/2007 7:31:36

Re: [R] array writing and their filenames

2007-07-30 Thread Dong GUO 郭东
Many thanks to Edna and Roland. I followed the suggestions, and it worked well. Regards, Dong On 7/30/07, Edna Bell [EMAIL PROTECTED] wrote: #Original 3x4x2 array xb , , 1 [,1] [,2] [,3] [,4] [1,] 0.4 0.9 5.6 0.1 [2,] 0.3 2.3 3.3 0.7 [3,] 0.6 0.8 0.2 0.7 , , 2

Re: [R] the large dataset problem

2007-07-30 Thread Roland Rau
Eric Doviak wrote: I need to find some way to overcome these constraints and work with large datasets. Does anyone have any suggestions? I might be not the most authoritative person on this subject but I put all my large datasets[1] into an SQLite database and extract/summarize data from it

Re: [R] Looping through all possible combinations of cases

2007-07-30 Thread jim holtman
Here is how to do it for 2; you can extend it: # test data n - 100 x - data.frame(id=sample(letters[1:4], n, TRUE), values=runif(n)) # get combinations of 2 at a time comb.2 - combn(unique(as.character(x$id)), 2) for (i in 1:ncol(comb.2)){ + cat(sprintf(%s:%s %f\n,comb.2[1,i],

[R] array loop

2007-07-30 Thread Dong GUO 郭东
Dear all, here are two arrays: region(26,31,8), nation(8) I tried to get a new array, say, giGi(26,31,8) giGi - array(0,dim = c(region_dim)) for (i in (1:region_dim[3])) { giGi[,,i] = region[,,i]-nation[,i] } As the above is part of function, but results shows only giGi[,,1] has the right

Re: [R] write.csv

2007-07-30 Thread Dong GUO 郭东
Thanks Ted for the help. I will try and see if my case will get the expected results. Dong On 7/30/07, Ted Harding [EMAIL PROTECTED] wrote: On 29-Jul-07 17:41:58, Dong GUO ¹ù¶« wrote: Hi, I want to save an array(say, array[6,7,8]) write a cvs file. How can I do that??? can I write in

[R] Odp: Looping through all possible combinations of cases

2007-07-30 Thread Petr PIKAL
Hi I am not sure about using loops but something like colSums(combn(DATA[,1],2)) gives you sums of all pairs and colSums(combn(DATA[,1],3)) sums of all triplets etc. if you want to know which ones they are just use combn on names. e.g. y X1.5 a1 b2 c3 d4 e5

Re: [R] array writing and their filenames

2007-07-30 Thread Roland Rau
Dong GUO 郭东 wrote: Hi, I want to save a array (say, array[6,7,8]) write a cvs file. How can I do that??? can I write in one file? For array[6,7,8], you don't need a csv(!) file since it is only a scalar. If this is what you want, check ?write.table But what you probably meant is how to

Re: [R] Order by the columns

2007-07-30 Thread jim holtman
?order You could do something like: mat[order(mat[,1], mat[,2], mat[,3]),] On 7/29/07, Am Stat [EMAIL PROTECTED] wrote: Dear useR, I have a data matrix, it has n columns, each column is a two-level variable with entires -1 and +1. They are randomly generated, now I want to order them like

Re: [R] write.csv

2007-07-30 Thread Dong GUO 郭东
the dim of my results is (26,31,8) -(years, regions and variables). so, if i save each (years, regions) in 8 csv files, later, I could connect the (26,31) to dbf file in ArcGIS to show in a map. This is what I intend to do. I dont know a better way to do it directly in R... On 7/31/07, jim

Re: [R] write.csv

2007-07-30 Thread jim holtman
Then you can just write a 'for' loop to write out each submatrix: for (i in 1:dim(x)[3]){ write.csv(x[,,i], paste(x, i, .csv, sep=)) } On 7/30/07, Dong GUO 郭东 [EMAIL PROTECTED] wrote: the dim of my results is (26,31,8) -(years, regions and variables). so, if i save each (years, regions)

Re: [R] Creating an instance of R from MS Access?

2007-07-30 Thread D L McArthur
Felipe Carrillo mazatlanmexico at yahoo.com writes: Hi all: Does anyone know if it's at all possible to create a connection to R from MS access? ... See function odbcConnectAccess in package RODBC. --D L McArthur dmca at ucla.edu __

Re: [R] write.csv

2007-07-30 Thread Dong GUO 郭东
Many thanks, Jim... On 7/31/07, jim holtman [EMAIL PROTECTED] wrote: Then you can just write a 'for' loop to write out each submatrix: for (i in 1:dim(x)[3]){ write.csv(x[,,i], paste(x, i, .csv, sep=)) } On 7/30/07, Dong GUO ¹ù¶« [EMAIL PROTECTED] wrote: the dim of my results is

[R] 'non-standard' folder names in R package

2007-07-30 Thread Tao Shi
Hi list, I use a .xml file for a function's demo in the R package I'm creating. Since it doesn't belong to any of the 'standard' folders, i.e. those mentioned in the 'Writing R Extension', I put it in a folder call myXML, much like the 'iris.xl' file in 'xls' folder from 'gdata' package, for

Re: [R] getting the name of variables passed to a function

2007-07-30 Thread Horace Tso
Thanks to Prof. Brian Ripley, Marc, and Andy. match.call() is what i need. After spending a weekend away from email, I stumbled on a thread on exactly the same topic this morning, http://finzi.psych.upenn.edu/R/Rhelp02a/archive/101445.html Horace Liaw, Andy [EMAIL PROTECTED] 7/30/2007

Re: [R] Overlaying a single contour from a new data array in levelplot

2007-07-30 Thread Deepayan Sarkar
On 7/30/07, Jenny Barnes [EMAIL PROTECTED] wrote: Dear Deepayan Thank you for your response - it has proved very very helpful, I can't thank you enough! I have another question for you if you have time to reply. I know you have been asked about the colour of the polygon outline before

[R] Error message in lmer

2007-07-30 Thread Gang Chen
I'm trying to run a simple one-way ANCOVA with the lmer function in R package lme4, but have encountered some conceptual problem. The data file MyData.txt is like this: Group Subj Cov Resp A 1 3.90 4.05 A2 4.05 4.25 A3 4.25 3.60 A4 3.60

Re: [R] write.csv

2007-07-30 Thread Ted Harding
On 29-Jul-07 17:41:58, Dong GUO ¹ù¶« wrote: Hi, I want to save an array(say, array[6,7,8]) write a cvs file. How can I do that??? can I write in one file? if I could not write in one file, i want to use a loop to save in different files (in the array[6,7,8], should be 8 csv files), such

[R] add custom strip to lattice plot

2007-07-30 Thread Patrick Drechsler
Hi, what is the recommended way of adding a strip to a lattice plot? In the example below I would like to add the value of mean(y) to a new strip.: --8---cut here---start-8--- library(lattice) ## Small sample data set: p0 - xyplot(uptake ~ Type | Treatment,

Re: [R] lattice grayscale theme

2007-07-30 Thread Deepayan Sarkar
On 7/28/07, Patrick Drechsler [EMAIL PROTECTED] wrote: Hi, is there a grayscale setting for lattice plots? I like the default color settings. I also like the settings that are available for setting black and white with something like this: --8---cut

Re: [R] problems saving and loading (PLMset) objects

2007-07-30 Thread Quin Wills
Ah, didn't realize that I couldn't re-assign in one step... I was trying to load in various data (on separate occasions), using a common object name to run through some template code. Many thanks! -Original Message- From: jim holtman [mailto:[EMAIL PROTECTED] Sent: 30 July 2007 23:51

[R] simple coding question

2007-07-30 Thread Kirsten Beyer
I have a list of ICD9 (disease) codes with various formats - 3 digit, 4 digit, 5 digit. The first three digits of these codes are what I am most interested in. I would like to either add zeros to the 3 and 4 digit codes to make them 5 digit codes or add decimal points to put them all in the

Re: [R] Tabs in PDF documents

2007-07-30 Thread Sundar Dorai-Raj
Dennis Fisher said the following on 7/30/2007 6:25 AM: Colleagues, I am using R 2.5.1 on an Intel Mac (OS 10) to create PDF outputs using pdf(); same problem exists in Linux (RedHat 9) While adding text to the document with text() and mtext(), I encounter the following problem:

Re: [R] how to combine data of several csv-files

2007-07-30 Thread jim holtman
Is this what you want: x - read.table(textConnection( v1 v2 v3 v4 v5 v6 v7 v8 + 1 NA -0.6442149 0.02354036 -1.40362589 -1.1829260 1.17099178 -0.046778203 NA + 2 NA -0.2047012 -1.36186952 0.13045724 2.1411553 0.49248118 -0.233788840 NA

Re: [R] Extract random part of summary nlme

2007-07-30 Thread Andrew Robinson
Hi Rense, Try: fm1 - lme(distance ~ age, data = Orthodont) VarCorr(fm1) Subject = pdSymm(age) Variance StdDevCorr (Intercept) 5.41508675 2.3270339 (Intr) age 0.05126947 0.2264276 -0.609 Residual1.71620460 1.3100399 I hope that this helps, Andrew

Re: [R] lattice grayscale theme

2007-07-30 Thread Deepayan Sarkar
On 7/30/07, Patrick Drechsler [EMAIL PROTECTED] wrote: The Gmane interface seems to have some lag at the moment... Deepayan Sarkar [EMAIL PROTECTED] writes: On 7/28/07, Patrick Drechsler [EMAIL PROTECTED] wrote: is there a grayscale setting for lattice plots? I like the default color

Re: [R] 2nd R Console

2007-07-30 Thread Michael Janis
Oh now that *is* something interesting! Thank you Greg, I'll have to give this a try. -Original Message- From: Greg Snow [mailto:[EMAIL PROTECTED] Sent: Monday, July 30, 2007 3:44 PM To: [EMAIL PROTECTED]; r-help@stat.math.ethz.ch Subject: RE: [R] 2nd R Console Have you looked at the

Re: [R] Tabs in PDF documents

2007-07-30 Thread ONKELINX, Thierry
I think you need to escape the '\' character in your string by adding an extra '\' in front of it. So try pdf(junk.pdf) par(family=mono) plot(1,1) text(1,1, \\txx) mtext(\\txx) dev.off() HTH, Thierry -Oorspronkelijk bericht- Van: [EMAIL

Re: [R] the large dataset problem

2007-07-30 Thread Greg Snow
Check out the biglm package for some tools that may be useful. -Original Message- From: Eric Doviak [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch r-help@stat.math.ethz.ch Sent: 7/30/07 9:54 AM Subject: [R] the large dataset problem Dear useRs, I recently began a job at a very large and

Re: [R] lattice grayscale theme

2007-07-30 Thread Patrick Drechsler
Deepayan Sarkar [EMAIL PROTECTED] writes: On 7/30/07, Patrick Drechsler [EMAIL PROTECTED] wrote: The Gmane interface seems to have some lag at the moment... Deepayan Sarkar [EMAIL PROTECTED] writes: On 7/28/07, Patrick Drechsler [EMAIL PROTECTED] wrote: is there a grayscale setting for

Re: [R] problems saving and loading (PLMset) objects

2007-07-30 Thread jim holtman
you just need to say: load(expr.RData) You should not be assigning it to 'expr' since it is already 'load'ed On 7/30/07, Quin Wills [EMAIL PROTECTED] wrote: Hi I'm running the latest R on a presumably up to date Linux server. 'Doing something silly I'm sure, but can't see why my saved

Re: [R] duplicate DATE at in lattice scale handled differently from base graphics OR lattice numeric scales

2007-07-30 Thread Deepayan Sarkar
On 7/30/07, Alex Brown [EMAIL PROTECTED] wrote: When at, label chains contain duplicate at values, axis ticks are dropped. I believe this is handled incorrectly for Date ats in lattice 0.15-4 when compared to how it is handled for numeric, or for dates in base plot. This can result in

Re: [R] deriv; loop

2007-07-30 Thread jim holtman
for question 1, is this what you want (BTW allocate 'result' to the size you want - the example keeps extending it which is OK for small numbers, but for larger size preallocate): result - numeric(0) for (i in 1:6) result[i] - i result [1] 1 2 3 4 5 6 prod(result) [1] 720 On 7/29/07,

Re: [R] the large dataset problem

2007-07-30 Thread jim holtman
FYI. I used your script on a Windows machine with 1.5GHZ and using the CYGWIN software that has the UNIX utilities. The field as 1000 lines with 10,000 fields on each line. Here is what it reported: gawk 'BEGIN{FS=,}{print $(1) , $(1000) , $(1275) , $(5678)}' tempxx.txt newdata.csv real

Re: [R] Constructing correlation matrices (follow up)

2007-07-30 Thread Horace Tso
Greg, in light of Doug Bates' question, what i have suggested a little early in response to your question is known as a Wishart matrix with n degree of freedom, which is guarenteed to be positive definite. If this is not what you want, you have to be more specific about the property of this

Re: [R] filenames

2007-07-30 Thread Dirk Eddelbuettel
On Mon, Jul 30, 2007 at 09:58:27PM +0200, Dong GUO 郭东 wrote: I want to create filename from a loop, say, i=(1:10), and the final names will be file1.csv, file2.csv, file3.csv in Python, it seems easer as - file+str(i) +'.csv , but how can i do in R list = as.character(1:10) list [1]

Re: [R] lattice grayscale theme

2007-07-30 Thread Patrick Drechsler
The Gmane interface seems to have some lag at the moment... Deepayan Sarkar [EMAIL PROTECTED] writes: On 7/28/07, Patrick Drechsler [EMAIL PROTECTED] wrote: is there a grayscale setting for lattice plots? I like the default color settings. I also like the settings that are available for

Re: [R] how to install rattle() in R-GUI2.5.1

2007-07-30 Thread Prof Brian Ripley
1) What is R-GUI? And what is R-GUI2.4.4? 2) In R 2.5.1 under Windows this works if you choose a properly populated CRAN mirror. I think http://cran.au.r-project.org has access problems. On Tue, 31 Jul 2007, j.joshua thomas wrote: Dear Group, I have installed R-GUI 2.5.1, thereby i tried

Re: [R] Matrix Multiplication, Floating-Point, etc.

2007-07-30 Thread jim holtman
One thing to realize is that although it appears that the operations are the same, the code that is being executed is different in the two cases. Due to the different sequence of instructions, there may be round-off errors that are then introduced On 7/30/07, Talbot Katz [EMAIL PROTECTED] wrote:

Re: [R] Creating an instance of R from MS Access?

2007-07-30 Thread Prof Brian Ripley
On Mon, 30 Jul 2007, D L McArthur wrote: Felipe Carrillo mazatlanmexico at yahoo.com writes: Hi all: Does anyone know if it's at all possible to create a connection to R from MS access? ... See function odbcConnectAccess in package RODBC. That is the other way, from R to an Access

Re: [R] Matrix Multiplication, Floating-Point, etc.

2007-07-30 Thread Moshe Olshansky
After multiplication by 10 you get 6*8 = 48 - the result is an exact machine number so there is no roundoff, while 0.6*0.8 = 0.48, where neither of the 3 numbers (0.6, 0.8, 0.48) is an exact machine mumber. However, (-0.6)*0.8 should be equal EXACTLY to -(0.6*0.8), and in fact you get that

Re: [R] simple coding question

2007-07-30 Thread Peter Alspach
Kirsten One way to do this: kirsten - c(123, 1234, 12345) 100*as.numeric(paste(substring(kirsten, 1, 3), substring(kirsten, 4, 5), sep='.')) HTH Peter Alspach -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Kirsten Beyer Sent:

<    1   2