[R] Selecting Best Model in an anova.

2010-03-25 Thread Hector Guilarte
Hello, I have a simple theorical question about regresion... Let's suppose I have this: Model 1: Y = B0 + B1*X1 + B2*X2 + B3*X3 and Model 2: Y = B0 + B2*X2 + B3*X3 I.E. Model1 = lm(Y~X1+X2+X3) Model2 = lm(Y~X2+X3) The Ajusted R-Square for Model1 is 0.9 and the Ajusted R-Square for Model2 is

[R] Error in library(gplots) : there is no package called 'gplots'

2010-03-25 Thread Martin Valere
Dear all, I have an issue trying to install new packages (have tried with RODBC_1.3-1, gplots_2.6.1, gtools_2.7.4 packages) and get the same error message : Error in library(gplots) : there is no package called 'gplots' Only clue I have found so far on the Web is related to Perl (Perl

Re: [R] How to use the paste function to create an already used variable

2010-03-25 Thread Biedermann, Jürgen
Thank you very much for the solutions! :-) It's my fault, but I forgot to mention that at the end of the loop, I would like to give Data back to data1, data2, data3, data4 The eval formula only seems to work in one direction (it's the same with the get function) for (i in 1:4) {

[R] help with breaking loops used to fit covariates in nlme model building procedure

2010-03-25 Thread Katya Mauff
Dear All I'm attempting to speed up my model building procedure, but need some help with the loops I've created...please bear with me through the explanation! My basic model call is something like: m0sulf.nlme-nlme(conc~beta0*exp(-beta1*day)+beta2*exp(-beta3*day),

Re: [R] Error in library(gplots) : there is no package called 'gplots'

2010-03-25 Thread Andrew Dolman
Do you mean installing packages as in install.packages(gplots) ? for which you need an internet connection. or attaching packages that are already installed, when you want to use them library(gplots) Or are you trying to install from source code or from downloaded zip packages?

Re: [R] Error in library(gplots) : there is no package called 'gplots'

2010-03-25 Thread Martin Valere
Dear Andrew, thanks for our interesrt. yes I mean for instance install.packages(gplots_2.7.4.tar.gz, repos=NULL, type=source) with the correct path to the package location of course Valère [Martin Valere] -Ursprüngliche Nachricht- Von: Andrew Dolman

Re: [R] Filling a grid based on existing data

2010-03-25 Thread Steve Murray
Sorry, that was a poorly worded question. You're right in that the gaps are in fact NAs and I would be proposing to remove these entirely. So it's really a case of filling up a 720 x 360 grid (by row) based on the data in the 18556 rows x 19 columns data frame. I've tried doing: data_mat -

[R] Stratified random sampling

2010-03-25 Thread wesley mathew
Dear All Could you please tell me, is there any package or function for finding Stratified random samples from the population ? Thanks in advance Kind Regards Wesley C Mathew [[alternative HTML version deleted]] __ R-help@r-project.org

[R] New paper on why not to log-transform count data

2010-03-25 Thread Graziella Iossa
Users of R-Help might be interested in this new article by Bob O'Hara and Johan Kotze published on Methods in Ecology and Evolution on why you should not log-transform count data, including an R code: http://www3.interscience.wiley.com/journal/123328987/abstract Thanks Kind regards, Graziella

[R] Counting a number of elements in an object

2010-03-25 Thread Greg Gilbert
I apologize if this has been answered. I have researched this to the best of my ability, that's not to say the answer isn't in the archives just I am a new user and I don't know the proper terms to search under. I have an object: f - mpr100 ~ time + nhb + hispanic + other +

[R] Error using Rcpp

2010-03-25 Thread Abhisek
Hi, Im not sure if this is the right place to post this. I am using Xubuntu Karmic Koala and am trying to use the Rcpp package. I am testing it using a simple code that takes in a vector and adds 1 to each element: #include Rcpp.h // This file takes in a vector and adds one to each entry

[R] setting sensitivity of r to errors

2010-03-25 Thread Jannis
Dear all, does anyone of you know how to increase Rs sensitivity to errors? I just migrated back from Matlab and really enjoyed there that Matlab just pops up with (really helpful!) error messages as soon as there is anything slightly wrong with my code. This is certainly anoying on the first

Re: [R] translating SQL statements into data.table operations

2010-03-25 Thread Matthew Dowle
Nick, Good question, but just sent to the wrong place. The posting guide asks you to contact the package maintainer first before posting to r-help only if you don't hear back. I guess one reason for that is that if questions about all 2000+ packages were sent to r-help, then r-help's traffic

[R] Returning Data Frame from Function for use Outside Function

2010-03-25 Thread Greg Gilbert
I have a function (see below) that does some bootstrapping (I am happy to expand offline why I could use existing functions.) I put my results into and empty matrix and add a row of results with each iteration. My problem is i am a new user to R and I don't understand data frames, matrices,

[R] Returning Data Frame from Function for use Outside Function

2010-03-25 Thread Greg Gilbert
I have a function (see below) that does some bootstrapping (I am happy to expand offline why I could use existing functions.) I put my results into and empty matrix and add a row of results with each iteration. My problem is i am a new user to R and I don't understand data frames, matrices,

[R] Returning Data Frame from Function for use Outside Function

2010-03-25 Thread Greg Gilbert
I have a function (see below) that does some bootstrapping (I am happy to expand offline why I could use existing functions.) I put my results into and empty matrix and add a row of results with each iteration. My problem is i am a new user to R and I don't understand data frames, matrices,

[R] Returning Data Frame from Function for use Outside Function

2010-03-25 Thread Greg Gilbert
I have a function (see below) that does some bootstrapping (I am happy to expand offline why I could use existing functions.) I put my results into and empty matrix and add a row of results with each iteration. My problem is i am a new user to R and I don't understand data frames, matrices,

Re: [R] Error using Rcpp

2010-03-25 Thread Romain Francois
Le 25/03/10 13:16, Abhisek a écrit : Hi, Im not sure if this is the right place to post this. It is not. The Rcpp-devel mailing list is the right place: https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel I am using Xubuntu Karmic Koala and am trying to use the Rcpp

[R] nls, predict() problem

2010-03-25 Thread Kay Cichini
hello, can anyone help with this: ### ###data: measurments (response = trans) run several times at the same predictor value level (press) por-data.frame(list(structure(list(run = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L,

Re: [R] How to use the paste function to create an already used variable

2010-03-25 Thread Biedermann, Jürgen
Perfekt! So, the code below does, what I want it does. for (m in 1:AnzRat) { Daten_akt - eval(parse(text=paste(Daten,m,_akt,sep=))) ... ... assign(paste(Daten,m,_akt,sep=),Daten_akt) } Patrick Burns schrieb: I think you are looking for 'assign'. On 25/03/2010 09:40,

[R] Do the missing values in random effect structure matter?

2010-03-25 Thread willow1980
Hi, everyone, I met a problem with mixed effect model with lmer. Basically, I want to treat region and mother ID as a nested random effect structure with mother ID nested within a specific region. However, in my dataset, many subjects had no mother ID (about 25% of them). I want to know, if I

Re: [R] setting sensitivity of r to errors

2010-03-25 Thread Barry Rowlingson
On Thu, Mar 25, 2010 at 12:17 PM, Jannis bt_jan...@yahoo.de wrote: Dear all, does anyone of you know how to increase Rs sensitivity to errors? There's a few things you can do with options(): check.bounds warns you extending vectors 'accidentally', 'warn' lets you do stuff like turn warnings

[R] *** caught segfault *** address 0x18, cause 'memory not mapped'

2010-03-25 Thread bernardo lagos alvarez
Hello R Community, I've been run the following codes. However, I've been getting an unusual segfault that I'm unable to trace its origin. Please give me a light to decipher the caught segfault Thanks for you attention. Bernardo. options(STERM='iESS', editor='emacsclient') rm(list = ls())

[R] Creation of a new data frame

2010-03-25 Thread Steve Murray
Dear all, I have a data frame of 18556 rows and 19 columns and wish to create a new grid from these data of dimensions 360 rows and 720 columns. The existing data frame is set up so that every 38 rows makes up one row of the new data frame, with 2 NA values at the end of each 'block' that

Re: [R] Multi-panel Pie Charts.

2010-03-25 Thread Hrishi Mittal
Hi Jim, what's the tab.title function you are using? - Try http://prettygraph.com Pretty Graph , the easiest way to make R-powered graphs on the web. -- View this message in context: http://n4.nabble.com/Multi-panel-Pie-Charts-tp1687026p1690627.html Sent from the R help mailing list

Re: [R] Mosaic Plots

2010-03-25 Thread Petr PIKAL
Hi There is also mosaicplot in base graphics ?mosaicplot Petr r-help-boun...@r-project.org napsal dne 22.03.2010 14:16:58: Hi, maybe vcd package is what you are looking for. Example from vcd library library(vcd) ?mosaic data(Titanic) mosaic(Titanic) ## Formula interface for

[R] I have a question on nomograms.

2010-03-25 Thread 笑啸
Dear volunteer: I am a graduate student of medcine in china.And now,I am devoting myself to constructing a nomograms of bladder cancer.I want to do it with R-project.However, I do not know how to construct a nomograms with R-project.I want to get yours help,thank you! I wish you can tell me

Re: [R] I have a question on nomograms.

2010-03-25 Thread Marc Schwartz
On Mar 25, 2010, at 8:54 AM, 笑啸 wrote: Dear volunteer: I am a graduate student of medcine in china.And now,I am devoting myself to constructing a nomograms of bladder cancer.I want to do it with R-project.However, I do not know how to construct a nomograms with R-project.I want to get

Re: [R] *** caught segfault *** address 0x18, cause 'memory not mapped'

2010-03-25 Thread Henrik Bengtsson
FYI, your more likely to get any response to this if you first check if the issue has been solved in a more recent version of R. Your version is 1.5 years old. /Henrik On Thu, Mar 25, 2010 at 2:03 PM, bernardo lagos alvarez blacert...@gmail.com wrote: Hello R Community, I've been run the

[R] Error in png concerning 'type = cairo'

2010-03-25 Thread Alla Bulashevska
Dear R-users, though i have successfully installed Cairo_1.4-5 package on my Ubuntu machine, i become the following error message: Error in png(plot_name, width = 1600, height = 1200) : 'type = cairo' requires cairo = 1.2 : try 'type = cairo1' By using png(plot_name, width = 1600, height =

Re: [R] I have a question on nomograms.

2010-03-25 Thread Ted Harding
On 25-Mar-10 13:54:38, ЦХ wrote: Dear volunteer: I am a graduate student of medcine in china.And now,I am devoting myself to constructing a nomograms of bladder cancer.I want to do it with R-project.However, I do not know how to construct a nomograms with R-project.I want to get yours

Re: [R] Stratified random sampling

2010-03-25 Thread Girish A.R.
Check out package survey cheers, -Girish -- View this message in context: http://n4.nabble.com/Stratified-random-sampling-tp1690475p1690730.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] Insert .eps files in to an R plot.

2010-03-25 Thread mnstn
Hello Everybody, I have an eps figure an awesome bacteria and a plot (generated using R) also in eps format. Now it looks like there is space for only one figure and I have to insert the picture of the bacteria into the plot. Is there a way to insert figures (eps/png/jpg) in to plots (may be

Re: [R] define F-ratio computations with aov

2010-03-25 Thread Galanidis Alexandros
Many thanks to both of you! I used the summary(obj.aov)$Mean way and I have the correct values now, although slightly different from the analogous SPSS computations. I'll try splplot and Error methods in a calmer period, but I have the feeling that in such cases aov is not the best method

[R] Goodness of Fit for Word Frequency Counts

2010-03-25 Thread Thiemo Fetzer
Dear Mailing list! sorry to bother you - but maybe you can help me out. I have been searching and searching for appropriate tests. I have a huge dataset of loan requests and I have data at portfolio level, with average portfolio size 200 loans. I want to test whether portfolios are randomly

Re: [R] I have a question on nomograms.

2010-03-25 Thread Frank E Harrell Jr
Marc Schwartz wrote: On Mar 25, 2010, at 8:54 AM, 笑啸 wrote: Dear volunteer: I am a graduate student of medcine in china.And now,I am devoting myself to constructing a nomograms of bladder cancer.I want to do it with R-project.However, I do not know how to construct a nomograms with

Re: [R] How to use the paste function to create an already used variable

2010-03-25 Thread Sharpie
quot;Biedermann, Jürgenquot; wrote: Perfekt! So, the code below does, what I want it does. for (m in 1:AnzRat) { Daten_akt - eval(parse(text=paste(Daten,m,_akt,sep=))) ... ... assign(paste(Daten,m,_akt,sep=),Daten_akt) } You may want to consider replacing

Re: [R] Error using Rcpp

2010-03-25 Thread Abhisek
Many thanks Romain! It worked. And thanks for pointing out the errors. best abhisek On Thu, Mar 25, 2010 at 12:45 PM, Romain Francois rom...@r-enthusiasts.comwrote: Le 25/03/10 13:16, Abhisek a écrit : Hi, Im not sure if this is the right place to post this. It is not. The Rcpp-devel

Re: [R] Choosing the larger of row pairs

2010-03-25 Thread James Rome
On 3/24/2010 10:58 PM, Peter Alspach wrote: jim key value 1 1 1 2 2 0 3 2 2 4 3 0 5 4 0 6 5 1 7 6 3 8 6 2 9 7 0 tt - rle(jim$key)$lengths ttJim - jim[cumsum(tt)-tt+tapply(jim$value, jim$key, which.max),] ttJim key value 1

[R] Odp: Creation of a new data frame

2010-03-25 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 25.03.2010 14:32:46: Dear all, I have a data frame of 18556 rows and 19 columns and wish to create a new grid from these data of dimensions 360 rows and 720 columns. The existing data frame is set up so that every 38 rows makes up one row of

Re: [R] Insert .eps files in to an R plot.

2010-03-25 Thread Sharpie
mnstn wrote: Hello Everybody, I have an eps figure an awesome bacteria and a plot (generated using R) also in eps format. Now it looks like there is space for only one figure and I have to insert the picture of the bacteria into the plot. Is there a way to insert figures (eps/png/jpg) in

Re: [R] how to deal with vector[0]?

2010-03-25 Thread Sharpie
Jin wrote: Hi, I have a vector with 4 elements, e.g., tau_i=c(100,200,300,400), but potentially tau_i[0]=0. In a for loop, tau_i=c(100,200,300,400) m=4 tau_i[0]=0 # --- ? P_i=1 for(i in 2:m) { P_i = P_i*(tau_i[i-1]-tau_i[i-2]) } Error

[R] Abstract classes

2010-03-25 Thread Albert-Jan Roskam
Hi R lovers, I'm looking for more information about, and implementations of abstract classes. After reading Head First Design Patterns (O'Reilly). I want to know whether I could use this to build an interface, i.e. a placeholder for all possible methods in the underlying subclasses. The final

[R] Basic graphs: something like groups, but where each plot has independent axis?

2010-03-25 Thread Ben Bimber
I am trying to graph weight data for our colony. We a data frame with Animal, Weight and Date. I'd like to display this data a series of simple xyplots. We only need to display these plots together, we do not need to make comparisons between animals. Each animal has been weighted over a

Re: [R] Abstract classes

2010-03-25 Thread Sharpie
Albert-Jan Roskam wrote: Hi R lovers, I'm looking for more information about, and implementations of abstract classes. After reading Head First Design Patterns (O'Reilly). I want to know whether I could use this to build an interface, i.e. a placeholder for all possible methods in the

[R] Logit/probit model

2010-03-25 Thread Ana De Barros
Deal all, I have a population with the following characteristics: N=1000 X0=rep(1,N) X1=rnorm(N) X2=rnorm(N) I also know that the population distribution is a linear logistic function with parameters alpha0=0 (intercept), alpha1=0.4 and alpha2=1.1. So easily I can get the dependent variable (in

Re: [R] Basic graphs: something like groups, but where each plot has independent axis?

2010-03-25 Thread Sharpie
Ben Bimber wrote: I am trying to graph weight data for our colony. We a data frame with Animal, Weight and Date. I'd like to display this data a series of simple xyplots. We only need to display these plots together, we do not need to make comparisons between animals. Each animal has

Re: [R] I have a question on nomograms.

2010-03-25 Thread Marc Schwartz
On Mar 25, 2010, at 9:51 AM, Frank E Harrell Jr wrote: Marc Schwartz wrote: On Mar 25, 2010, at 8:54 AM, 笑啸 wrote: Dear volunteer: I am a graduate student of medcine in china.And now,I am devoting myself to constructing a nomograms of bladder cancer.I want to do it with

[R] Reading SQL Server Tables using RODBC

2010-03-25 Thread Paul Miller
Hello Everyone,   I'm still quite new to R and am having trouble reading SQL Server Tables using RODBC. I've looked though the RODBC documentation as well as material I found online and in the book Data Manipulation with R, but I just can't quite seem to get things to work.   Right now, I have

[R] how to deal with vector[0]?

2010-03-25 Thread Jin
Hi, I have a vector with 4 elements, e.g., tau_i=c(100,200,300,400), but potentially tau_i[0]=0. In a for loop, tau_i=c(100,200,300,400) m=4 tau_i[0]=0 # --- ? P_i=1 for(i in 2:m) { P_i = P_i*(tau_i[i-1]-tau_i[i-2]) } Error in P_i = P_i * (tau_i[k

[R] Sweave running R chunks for figure generation several times

2010-03-25 Thread Florian Burkart
Hi, following issue: Sweave is running the R code to generate figures three times under standard behaviour, twice if run with eps=FALSE, once if run with eps=FALSE,pdf=FALSE. However, ideally, I'd want it to run once to generate pds with eps=FALSE being set. How come this isn't so? Am I

[R] matrix to write table with new vector after calculation of one vector

2010-03-25 Thread nabiul khan
Dear Friends, I have a data matrix with 500 rows and 5 columns. The 1st row contains integer code (1, 2, 3). Other rows are real data. I want to convert the real data in 5th column based on CODES in 1st column, meaning conversion process are three types (depending on codes 1,2,3. Finally I want to

[R] Read SAS data

2010-03-25 Thread Nicola Sturaro Sommacal
Hi! I need to import in R some SAS dataset (sas7bdat). I found two functions to do it: read.ssd from the package foreign and sas.get from Hmisc. df = read.ssd(libname = path2data, sectionnames = sasSmallDataset, tmpXport = path2data, tmpProgLoc = path2data, sascmd = path2sas) sas.get(libraryName

Re: [R] Reading SQL Server Tables using RODBC

2010-03-25 Thread Bos, Roger
Paul, I have never tried odbcDriverConnect, but I can tell you how I connect to my SQL Server 2005 db. I create a DNS entry using odbcad32 (use startrun). You can specify everything there and test the connection. Once that tests okay, just use the following: library(RODBC) xf -

Re: [R] Reading SQL Server Tables using RODBC

2010-03-25 Thread Marc Schwartz
On Mar 25, 2010, at 11:46 AM, Paul Miller wrote: Hello Everyone, I'm still quite new to R and am having trouble reading SQL Server Tables using RODBC. I've looked though the RODBC documentation as well as material I found online and in the book Data Manipulation with R, but I just can't

Re: [R] Logit/probit model

2010-03-25 Thread Charles Annis, P.E.
Perhaps I am missing something but it appears that because X1 and X2 are random normal, that the influence of X2 is much like a second sampling of X1, and thus you would expect just what you observed, especially with a large (1000) sample size. Try making X2 and X1 different. Charles Annis, P.E.

[R] sample size mixed poisson glm

2010-03-25 Thread caspar hallmann
Dear list members, I am trying to estimate power/sample size based on already collected pilot data. The setting: Five orchards have been sampled with respect to fruit damage. In each orchard two plots were selected (semi-randomly), where one was treated and one left as control. Within each plot

Re: [R] I have a question on nomograms.

2010-03-25 Thread R Heberto Ghezzo, Dr
On the same topic but from a different perspective. A Nomogram or better a Line Aligned Nomogram is a graph with 2 or more scales, maybe linear where by alignig values in each scale you can read values in the other scale. The relationship can be linear, then the scales a straght lines or

[R] Problem with reshaping data from wide to long format

2010-03-25 Thread Lauri Nikkinen
Hi, I have a data.frame in wide format which I would like to reshape into a long format: example (nonsense) data: dput(perus2) structure(list(id = c(30L, 38L, 21L, 12L, 22L, 28L, 31L, 44L, 8L, 47L, 23L, 20L, 41L, 42L, 29L, 50L, 5L, 33L, 4L, 17L, 11L, 1L, 18L, 6L, 9L, 32L, 16L, 14L, 39L, 48L,

Re: [R] Problem with reshaping data from wide to long format

2010-03-25 Thread Gabor Grothendieck
Try new.row.names = 1:150 as an arg to reshape. On Thu, Mar 25, 2010 at 2:04 PM, Lauri Nikkinen lauri.nikki...@iki.fi wrote: Hi, I have a data.frame in wide format which I would like to reshape into a long format: example (nonsense) data: dput(perus2) structure(list(id = c(30L, 38L, 21L,

Re: [R] Problem with reshaping data from wide to long format

2010-03-25 Thread Lauri Nikkinen
Thanks, that solved my problem! -Lauri 2010/3/25 Gabor Grothendieck ggrothendi...@gmail.com: Try new.row.names = 1:150 as an arg to reshape. On Thu, Mar 25, 2010 at 2:04 PM, Lauri Nikkinen lauri.nikki...@iki.fi wrote: Hi, I have a data.frame in wide format which I would like to reshape

Re: [R] Read SAS data

2010-03-25 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Nicola Sturaro Sommacal Sent: Thursday, March 25, 2010 9:16 AM To: r-help@r-project.org Subject: [R] Read SAS data Hi! I need to import in R some SAS dataset (sas7bdat). I

[R] c(), or cat(), or paste(), all cause unwanted reordering

2010-03-25 Thread Jeff Brown
Hi, I would expect the following: paste( as.character( cat( rep( ., 2 ) ) ), a string, as.character( cat( rep( ., 3 ) ) ) ); to yield this string: . . a string . . ., but instead it yields this: . .. . .[1] a string The third argument has been stuck immediately

Re: [R] I have a question on nomograms.

2010-03-25 Thread Ted Harding
On 25-Mar-10 17:44:55, R Heberto Ghezzo, Dr wrote: On the same topic but from a different perspective. A Nomogram or better a Line Aligned Nomogram is a graph with 2 or more scales, maybe linear where by alignig values in each scale you can read values in the other scale. The relationship can

Re: [R] Sweave running R chunks for figure generation several times

2010-03-25 Thread Marc Schwartz
On Mar 25, 2010, at 10:14 AM, Florian Burkart wrote: Hi, following issue: Sweave is running the R code to generate figures three times under standard behaviour, twice if run with eps=FALSE, once if run with eps=FALSE,pdf=FALSE. However, ideally, I'd want it to run once to

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

2010-03-25 Thread Erik Iverson
Hello, Jeff Brown wrote: I would expect the following: paste( as.character( cat( rep( ., 2 ) ) ), a string, as.character( cat( rep( ., 3 ) ) ) ); to yield this string: . . a string . . ., but instead it yields this: . .. . .[1] a string cat is writing its

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

2010-03-25 Thread Henrique Dallazuanna
Maybe: gsub(^(.*)$, ..\\1..., a string) On Thu, Mar 25, 2010 at 3:53 PM, Jeff Brown dopethatwantsc...@yahoo.com wrote: Hi, I would expect the following: paste(        as.character( cat( rep( ., 2 ) ) ),        a string,        as.character( cat( rep( ., 3 ) ) ) ); to yield this

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

2010-03-25 Thread Erik Iverson
Erik Iverson wrote: Hello, Jeff Brown wrote: I would expect the following: paste( as.character( cat( rep( ., 2 ) ) ), a string, as.character( cat( rep( ., 3 ) ) ) ); to yield this string: . . a string . . ., but instead it yields this: . .. . .[1] a string cat is

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

2010-03-25 Thread Ted Harding
On 25-Mar-10 19:07:23, Erik Iverson wrote: Hello, Jeff Brown wrote: I would expect the following: paste( as.character( cat( rep( ., 2 ) ) ), a string, as.character( cat( rep( ., 3 ) ) ) ); to yield this string: . . a string . . ., but instead it yields this: .

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

2010-03-25 Thread Sarah Goslee
cat() isn't doing what you think, and as.character() is unnecessary here because paste() takes care of that. A couple of more elegant solutions have already appeared, but here's one that parallels your version: paste(c( rep( ., 2 ), a string, rep( ., 3 )), collapse= ) On

Re: [R] Filling a grid based on existing data

2010-03-25 Thread David Winsemius
On Mar 25, 2010, at 7:01 AM, Steve Murray wrote: Sorry, that was a poorly worded question. You're right in that the gaps are in fact NAs and I would be proposing to remove these entirely. So it's really a case of filling up a 720 x 360 grid (by row) based on the data in the 18556 rows

[R] RODBC : reading binary data from a TXT field belonging to a PostgeSQL table

2010-03-25 Thread Laurent Rhelp
Dear R-List, I am working with binary data that I want to store in a PostgreSQL DataBase. I decided to use a TXT field. I read my binary file with readBin function, I succeed in my data storage in the database but I have some trouble to extract the data : the correct amount of bytes is

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

2010-03-25 Thread Sarah Goslee
On Thu, Mar 25, 2010 at 3:19 PM, Ted Harding ted.hard...@manchester.ac.uk wrote: One needs to be very circumspect with this sort of thing! For instance, experimenting with simplifications of Jeff's expression:  paste(          rep( ., 2 ),          a string,          rep( ., 3 )  )  #

Re: [R] Multi-panel Pie Charts.

2010-03-25 Thread Jim Lemon
On 03/26/2010 12:41 AM, Hrishi Mittal wrote: Hi Jim, what's the tab.title function you are using? It's in the plotrix package. The panes function will be in the next version. Jim __ R-help@r-project.org mailing list

Re: [R] lattice grob

2010-03-25 Thread baptiste auguie
Hi, On 24 March 2010 23:22, Paul Murrell p.murr...@auckland.ac.nz wrote: Hi baptiste auguie wrote: Thanks Felix and Paul. I had overlooked grid.grabExpr, assuming that one had to draw on a device before grabbing the output. Now I'm not sure if there's any difference between either

Re: [R] c(), or cat(), or paste(), all cause unwanted reordering

2010-03-25 Thread Jeff Brown
Wow, you guys are awesome. Thanks! -- View this message in context: http://n4.nabble.com/c-or-cat-or-paste-all-cause-unwanted-reordering-tp1691133p1691198.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

[R] RODBC : reading binary data from a TXT field belonging to a PostgeSQL table

2010-03-25 Thread Laurent Pantera
Dear R-List, I am working with binary data that I want to store in a PostgreSQL DataBase. I decided to use a TXT field. I read my binary file with readBin function, I succeed in my data storage in the database but I have some trouble to extract the data : the correct amount of bytes is

Re: [R] how to deal with vector[0]?

2010-03-25 Thread Jin
Hello Charlie, Thanks for your valuable help. You are right. So I changed tau_i=c(100,200,300,400) to tau_i=c(0,100,200,300) because 400 can be defined by another term. I appreciate for your time and knowledge. Best, JinNational Renewable Energy Lab.Golden CO, 80401(303)

[R] 3 levelplots and 1 colorbar

2010-03-25 Thread Joaquin Rapela
I want to create a simple plot containing three levelplots with one colorbar. I used the Three levelplots code below, but the third levelplot is drawn smaller than the first two. However, if I try the Two levelplots code below it works well. Can anybody tell me how could I draw three levelplots

[R] Creating dataframe of all possible variable combinations

2010-03-25 Thread Hosack, Michael
Hello, I need to create a dataframe containing all possible combinations of three variables: SITE (101,102,103,104), WDAY (MON,TUE,WED,THR,FRI), and TOD (MORN, AFTN). There should be a total of 40 unique combinations in my dataframe. I used expand.grid() successfully(?) to create my dataframe,

[R] Creating dataframe of all possible variable combinations

2010-03-25 Thread Hosack, Michael
Hello, I need to create a dataframe containing all possible combinations of three variables: SITE (101,102,103,104), WDAY (MON,TUE,WED,THR,FRI), and TOD (MORN, AFTN). There should be a total of 40 unique combinations in my dataframe. I used expand.grid() successfully(?) to create my dataframe,

Re: [R] Creating dataframe of all possible variable combinations

2010-03-25 Thread Jorge Ivan Velez
Hi Mike, Try this: SITE - c(101,102,103,104) WDAY - c('MON','TUE','WED','THR','FRI') TOD - c('MORN', 'AFTN') out - expand.grid(SITE, WDAY, TOD) out HTH, Jorge On Thu, Mar 25, 2010 at 4:21 PM, Hosack, Michael wrote: Hello, I need to create a dataframe containing all possible combinations

[R] Convert number to Date

2010-03-25 Thread anna
Hello, I have a date value in excel: 1/4/2010 which in number format gives me 40182. When I read this with read.xls from R I get same 40182 so what I do is that I use the as.Date() function but here is what the as.Date() function returns me: as.Date(40182) [1] 2080-01-06 Why don't I get the same

Re: [R] Creating dataframe of all possible variable combinations

2010-03-25 Thread Rolf Turner
On 26/03/2010, at 9:22 AM, Hosack, Michael wrote: Hello, I need to create a dataframe containing all possible combinations of three variables: SITE (101,102,103,104), WDAY (MON,TUE,WED,THR,FRI), and TOD (MORN, AFTN). There should be a total of 40 unique combinations in my dataframe. I

Re: [R] Convert number to Date

2010-03-25 Thread Rolf Turner
On 26/03/2010, at 9:33 AM, anna wrote: Hello, I have a date value in excel: 1/4/2010 which in number format gives me 40182. When I read this with read.xls from R I get same 40182 so what I do is that I use the as.Date() function but here is what the as.Date() function returns me:

Re: [R] Creating dataframe of all possible variable combinations

2010-03-25 Thread Stephan Kolassa
Hi Mike, the following works for me: SITE - ordered(c(101,102,103,104)) WDAY - ordered(c(MON,TUE,WED,THR,FRI),levels=c(MON,TUE,WED,THR,FRI)) TOD - ordered(c(MORN,AFTN),levels=c(MORN,AFTN)) foo - expand.grid(SITE=SITE,WDAY=WDAY,TOD=TOD) foo[order(foo$SITE),] If this doesn't solve your

Re: [R] 3 levelplots and 1 colorbar

2010-03-25 Thread RICHARD M. HEIBERGER
Joaquin, p -levelplot(frame, colorkey=TRUE) print(update(p, legend=NULL), position=c(0/4,0,1/4,1), more=TRUE) print(update(p, legend=NULL), position=c(1/4,0,2/4,1), more=TRUE) print(update(p, legend=NULL), position=c(2/4,0,3/4,1), more=TRUE) draw.colorkey(p$legend$right$args$key, draw=TRUE,

Re: [R] Convert number to Date

2010-03-25 Thread Joshua Wiley
Dear Anna, Rolf's explanation not withstanding, it has to do with differences in how R and Excel treat dates. If you use as.Date(40182, origin=1899-12-30) you will get the same date as Excel. You can look at:

Re: [R] Read SAS data

2010-03-25 Thread Mike Reese
Nicola, I found its easier if you convert your SAS data to the XPORT format (a SAS transport format). Here's instructions for creating the transport file in SAS - use xport not cport. http://support.sas.com/documentation/cdl/en/movefile/59598/HTML/default/a002575816.htm I also found the

[R] A file with extension .sdb in a codebook section of a large database from a survey?

2010-03-25 Thread Douglas Bates
The TIMSS2007 database http://timss.bc.edu/TIMSS2007/idb_ug.html seems to provide both kinds of universal data formats - either SPSS saved data sets or SAS saved data sets. (Yes, I am being sarcastic.) These, of course, are accompanied by massive codebooks explaining the nature of each of the

Re: [R] Insert .eps files in to an R plot.

2010-03-25 Thread Greg Snow
The subplot function in the TeachingDemos package will create a new plot inside of the existing one, so if you can plot your one image, the subplot function will help (using base graphics). The EBImage package may help with the reading and plotting of the figure. (there is also grImport, but

[R] Manually calculate SVM

2010-03-25 Thread Noah Silverman
Hi, I'm learning more about SVMs and kernels in general. I've gotten used to using the svm function in the e1071 package. It works great. Now, I want to do/learn some more interesting stuff. (Perhaps my own kernel and/or scoring system). So I want to better understand 1) how calculation of

[R] Frequencies from x/y data into a 2d table (for 3d histogram or heatmap)

2010-03-25 Thread Ralf B
Hi all, I have simple x/y data from screen recording in a sequence: number,x,y 1,10,30 1,20, 1,43,110 1,74,18 1,88,112 and would like to create a 3d histogram data structure that i can use to create a 3d histogram or, more likely a heatmap. The unterlying data structure

Re: [R] Convert number to Date

2010-03-25 Thread Joshua Wiley
You are quite welcome. After poking around a bit more, I can offer a more detailed explanation on Excel. It does treat 1 January 1900 as the origin. However, while R treats origin as 0, Excel treats it as 1. This explains 1 of the two day change needed for R to get the same results as Excel.

Re: [R] Interleaving elements of two vectors?

2010-03-25 Thread Jeff Brown
I just had to solve this problem for myself, after not having luck with the code posted above. I'm posting in case others need a completely general function. riffle - function (a,b) { # Interleave a b, starting with a, without repeating. x - NULL; count = 1;

[R] degree index application in network

2010-03-25 Thread reccos
Hi Could somebody help me how to get degree index in network, where I have weighted edges? Thanks -- View this message in context: http://n4.nabble.com/degree-index-application-in-network-tp1691304p1691304.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Convert number to Date

2010-03-25 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Joshua Wiley Sent: Thursday, March 25, 2010 1:48 PM To: anna Cc: r-help@r-project.org Subject: Re: [R] Convert number to Date Dear Anna, Rolf's explanation not

Re: [R] Creating dataframe of all possible variable combinations

2010-03-25 Thread Erik Iverson
Hello, Hosack, Michael wrote: Hello, I need to create a dataframe containing all possible combinations of three variables: SITE (101,102,103,104), WDAY (MON,TUE,WED,THR,FRI), and TOD (MORN, AFTN). There should be a total of 40 unique combinations in my dataframe. I used expand.grid()

[R] Precision level

2010-03-25 Thread Dimitri Liakhovitski
Hello! I am wondering at what point does R consider a numeric value to be equal to zero - for statements of the type x==0 and x %in% 0. Thank you very much! -- Dimitri Liakhovitski Ninah.com dimitri.liakhovit...@ninah.com __ R-help@r-project.org

[R] Absolutely No idea how to plot my Spatial Data

2010-03-25 Thread Babaorumi
Hi, I have a data set of points which are represented by 3 variables x,y,z where x is the position of the point on the x-absciss and y on the y-absciss. Each of my points has a value z, which I want to be displayed as follows: the more z is high, the more the color on the map is dark. How can I

  1   2   >