Re: [R] Problem at subsetting matrix by using dimnames

Thanks a lot. Based on your suggestions, would the following diag(mitest[c(a, b, c), c(b, c, b)]) work in general? (it does for my example). If so, I would prefer this form as I understand it better. Agus Henrique Dallazuanna wrote: Try this: mapply(function(x, y)mitest[x, y], c(a, b, c),

Re: [R] Pull Coefficients from MCMCpack models

Try this: First of all: predict(foo) Error in UseMethod(predict) : no applicable method for predict Therefore, an alternative approach may be taken: As per: bar - apply(foo,2,mean) bar[1] is the intercept bar[2] is the coefficient for x1 bar[2] is the coefficient for x2 Then

[R] Problem while using KLdiv()

I am receiving the following error while executing KLdiv method: kl-KLdiv(y) Error in function (classes, fdef, mtable) : unable to find an inherited method for function KLdiv, for signature data.frame the code is as follows: require(flexmix) KLdiv(y) I tried searching in the R mail archives

[R] Linear Model NA Value Test

Hello, I am deriving near real-time liner relationships based on 5-min precipitation data, sometimes the non-qced data result in a slope of NA. I am trying to read the coefficient (in this example x) to see if it is equal to NA, if it is equal to NA assign it a value of 1. I am having trouble

[R] Plotting waver files

Is there a good tutorial or just sample code on how to plot .wav files in R? I would really appreciate it! -- View this message in context: http://www.nabble.com/Plotting-waver-files-tp25530418p25530418.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Point Density calculation

Sorry for the confusion about my question. Here is the actual data I want to analyze: Test Data: 45 40 42.11 N 111 02 32.79 W 5 45 40 31.30 N 111 02 32.55 W 5 45 40 38.29 N 111 02 52.81 W 5 45 40 45.90 N 111 03 25.57 W 5 45 42 52.26 N 110 58 26.36 W 10 45.43 20.60 N 110 58 05.19 W 15

[R] cox memory

Hi there, I have a rather large data set and perform the following cox model: test1 - list(tstart,tstop,death1,chemo1,radio1,horm1) out1-coxph( Surv(tstart,tstop, death1) ~

Re: [R] Pull Coefficients from MCMCpack models

Try this: apply(foo, 2, mean) or apply(foo, 2, median) Thanks, Deb Noah Silverman n...@smartmediacorp.com 22/09/2009 12:34 pm Hi, I've been testing some models with the MCMCpack library. I can run the process and get a nice model object. I can easily see the summary and even plot

Re: [R] Quadratic Constraints

Thanks Hans for ur help I will try to solve the problem accordingly if there are any problems will contact u surely . -- View this message in context: http://www.nabble.com/Quadratic-Constraints-tp25528480p25530553.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Putting a text box in a plot

Hello, Jim Thank you a lot for suggestions, I will check that package out. It could be the one I used way back then! :-) Best, Sergey On Tue, Sep 22, 2009 at 01:55, Jim Lemon j...@bitwrit.com.au wrote: On 09/21/2009 07:42 PM, Sergey Goriatchev wrote: Hello everyone, I have a plot and I

[R] Markov Switching

DeaR list, I am looking for a descent implementation of a frequentist estimation of a standard Markov Switching Vector Autoregressive Model. Until now I found the following packages: 1) MSBVAR: It seems that this package estaimates Markov-Switching VAR-model only from a Bayesian point of

Re: [R] How to use nls when [selfStart] function returns NA or Inf??

Thanks Gabor, but my problem isn't finding reasonable starting parameter values, it's preventing nls giving up when it tries parameter values resulting in NA or Inf. I know queries are often over-specific and the appropriate response is don't start there, so I'm trying to balance between

Re: [R] Basic function output/scope question

Thanks Erik, Schalk, Patrick, and David for you helpful advice. I hope I'll, at some point, I'll become versed enough in R to return the favor. -- Best regards, David Young Marketing and Statistical Consultant Madrid, Spain +34 913 540 381 http://www.linkedin.com/in/europedavidyoung

[R] snowfall: missing MPI node

Hello, I don't know if the question pertains to Rmpi, snow or snowfall. I run my job by: mpirun -np N -hostfile $PBS_NODEFILE RMPISNOW -f my-script.r --slave In the snowfall sfInit call I have to specify one less CPU respect to the mpirun call sfInit(parallel=TRUE, cpus=N-1, type=MPI)

Re: [R] matching pairs regardless of order,multiple matches

Dear Jim and Henrique, thank you both for your help. I have done this but run into another problem:   In the example below loc1,loc2 occurs in the (now correct, thanks to your advice) list twice.     trips=(loc1,loc2,loc2,loc3)   DF$listoftrips=(loc1,loc2, loc1,loc3, loc2,loc3,loc1,loc2).   I 

Re: [R] problem using KLdiv - flexmix

Hai Can someone please tell me which package is good for CART in R? I have found something called PARTYKIT is that a good package for CART?? Thanks and Regards Naveena -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Robert Hall

Re: [R] problem using KLdiv - flexmix

Ohh i forgot to mention one thing the value 'y' is obtained as follows: s1p5-read.table(s1p5,header=TRUE) s2p5-read.table(s2p5,header=TRUE) y-cbind(s1p5,s2p5) I know the problem is that the value Y is not compatible with what KLdiv requires .. can anyone please tell how to build the

[R] trying to install plotrix

Here are the error messages: install.packages('plotrix') Warning in install.packages(plotrix) : argument 'lib' is missing: using '/Users/dbae/Library/R/library' --- Please select a CRAN mirror for use in this session --- Error in .install.macbinary(pkgs = pkgs, lib = lib, contriburl =

Re: [R] Printing column names before each block of output (fwd)

Hi Jim, I checked vcdf for names() , (this is my test data set), and it does have names: names(vcdf) [1] male.parent family offspring.id P1L55P1L73 [6] P1L74P1L77P1L91P1L96P1L98 [11] P1L100 P1L114 P1L118 peg.no ec.length

Re: [R] problem using KLdiv - flexmix

Got it! read.matrix(file) does the job! Rob Hall On Tue, Sep 22, 2009 at 3:27 PM, Robert Hall hall.robert@gmail.comwrote: Ohh i forgot to mention one thing the value 'y' is obtained as follows: s1p5-read.table(s1p5,header=TRUE) s2p5-read.table(s2p5,header=TRUE) y-cbind(s1p5,s2p5) I

Re: [R] cox memory

On Sep 21, 2009, at 7:27 PM, Λεωνίδας Μπαντής wrote: Hi there, I have a rather large data set and perform the following cox model: test1 - list(tstart,tstop,death1,chemo1,radio1,horm1) out1-coxph( Surv(tstart,tstop, death1) ~ chemo1+chemo1:log(tstop

Re: [R] trying to install plotrix

On 09/22/2009 08:13 PM, David Epstein wrote: Here are the error messages: install.packages('plotrix') Warning in install.packages(plotrix) : argument 'lib' is missing: using '/Users/dbae/Library/R/library' --- Please select a CRAN mirror for use in this session --- Error in

Re: [R] How to use nls when [selfStart] function returns NA or Inf??

With good starting values it often won't need to try anything outside feasible region. If after trying better starting values you find that that it is still attempting to move outsie the feasible region then another approach is to use the boundary value plus some monotonic function of the

Re: [R] matching pairs regardless of order,multiple matches

On Sep 22, 2009, at 5:52 AM, Juliane Struve wrote: Dear Jim and Henrique, thank you both for your help. I have done this but run into another problem: In the example below loc1,loc2 occurs in the (now correct, thanks to your advice) list twice. trips=(loc1,loc2,loc2,loc3)

[R] snowfall: sfExport apparently harmless error

I'm running my script using mpirun -mp 4 and using snowfall+Rmpi on Linux 64bits. I receive the following error, but apparently without consequences on the results. Any idea? I'm able to reproduce it with a minimal script (below). Seems the critical issue is the for loop. Without it no error.

[R] problem in CART (RPART)

Hai Can someone tel me whether it is possible to force in Variables at various levels of the tree in CART. I have been using RPART for CART which comes with PARTYKIT package. I am facing a problem with the tree it generates, it picks up variables which don't make much sense. I want it to pick up

Re: [R] matching pairs regardless of order,multiple matches

Dear list, apologies for not including a working example. match() does not seem to be required for my second question. Using the not-working example from below I managed to do something like aggregate(DF$Data,by=list(Trip=listoftrips),FUN=mean) This yields means by trip, i.e. the values I

[R] duncan test in R

Hello. Sorry if the question I'm asking is trivial but I'm unable to find the answer. How can I make a duncan test of multiple rank in R? That is, the designed test to see which mean levels are different in an one factor ANOVA. Thanks in advance, Arnau.

[R] Manipulate directory-help

Dear R-users, I have a problem: I have many files in the directory:C:/Documents and Settings/lma/Desktop/FamilyAEntrepreneurs/Entrepreneurs/Juha I want to copy those files into the new directory: C:/Documents and Settings/lma/My Documents/Juha/book (-use paste), I first create the directory

[R] R execution by VBA Macro

Hello everyone I am using RExcel, an add-in to use R interface within MS Excel workbook. The add-in is freely distributed by Statconn. Pls refer http://www.statconn.com/ I need help to start writing VBA Macro to invoke and execute R to read .csv file. This is what I made so far… But Compilation

[R] use of class variable in r as in Proc means of sas

Hi,everyone i need to calculate quartile values of a variable grouped by the other variable . same as in aggregate function(only median,mean or functions is possible-i think so) Could you please help me to achieve the same for other quartile values(5,10,25,75,90) as for median using aggregate.

[R] Semi continous variable- define bounds using lpsolve

How to define bounds for a semi continous variable in lp_solve. Min 5x1 +9x2 +7.15x3 +0.1x4 subject to x1+x2+x3+x4=6.7 x1+x4 = 6.5 And x3 can be 0 or greater than 3.6 hence x3 is a semi continous variable how to define bounds as well as semicontinous function because using set.semicont and

Re: [R] Manipulate directory-help

On Sep 22, 2009, at 4:59 AM, Tammy Ma wrote: Dear R-users, I have a problem: I have many files in the directory:C:/Documents and Settings/lma/ Desktop/FamilyAEntrepreneurs/Entrepreneurs/Juha I want to copy those files into the new directory: C:/Documents and Settings/lma/My

Re: [R] Quadratic Constraints

Hi Hans, I am very much thankful to you for helping me to solve the problem. I went through the link provided by you and tried to solve the problem. But I think there's some problem with the code. So it is not arriving at the optimal solution.In Some cases it shows optimum solution while in some

[R] converting a character vector to a function's input

Hi all, I have been trying to solve this problem and have had no luck so far. I have numeric vectors VAR1, VAR2, and VAR3 which I am trying to cbind. I also have a character vector VAR1,VAR2,VAR3. How do I manipulate this character vector such that I can input a transformed version of the

Re: [R] Manipulate directory-help

Try this: sourceDir - C:\\temp\\ destDir - C:\\new\\new if(is.na(file.info(destDir)$isdir)) dir.create(destDir, recursive = TRUE) file.copy(dir(sourceDir, full.names = TRUE), destDir) On Tue, Sep 22, 2009 at 5:59 AM, Tammy Ma metal_lical...@live.com wrote: Dear R-users, I have a

Re: [R] converting a character vector to a function's input

first part: VAR1 - 1:10 VAR2 - 11:20 VAR3 - 21:30 input - VAR1,VAR2,VAR3 # split the input i.s - unlist(strsplit(input, ',')) # create matrix do.call(cbind, lapply(i.s, get)) [,1] [,2] [,3] [1,]1 11 21 [2,]2 12 22 [3,]3 13 23 [4,]4 14 24 [5,]

[R] Odp: use of class variable in r as in Proc means of sas

Hi r-help-boun...@r-project.org napsal dne 22.09.2009 11:51:18: Hi,everyone i need to calculate quartile values of a variable grouped by the other variable . same as in aggregate function(only median,mean or functions is possible-i think so) Could you please help me to achieve the

[R] Odp: converting a character vector to a function's input

Hi r-help-boun...@r-project.org napsal dne 22.09.2009 14:57:59: Hi all, I have been trying to solve this problem and have had no luck so far. I have numeric vectors VAR1, VAR2, and VAR3 which I am trying to cbind. I also have a character vector VAR1,VAR2,VAR3. How do I manipulate this

[R] Programming insight

Hi all, Most of the time in my programming I use for loop however there will always be some way to use other loop like while, if etc to implement same query. My question is among those different loops which is theoretically fastest and why? Thanks -- View this message in context:

Re: [R] use of class variable in r as in Proc means of sas

Is this what you want: x zip price 1 6567000 2 60001478654 3 60004485647 4 60001 2783958 5 60005 97845848 6 60006 378383478 7 60002 397895735 8 60001 487587575 9 60002478848 10 60003 49847874 11 60004467648 12 60005567489 13 60006 4776746 14

[R] Function similar to cumsum/cumprod

Hello, everyone I wonder if there is in R somewhere a function similar to cumsum(). The function calculates a statistic (say mean or standard deviation) buy adding consequtively one more data point. So, say I have a timeseries of 100 observations. I start by calculating mean of first 30

[R] Singular model.matrix of nested designs

Hi, I want to do ANOVA for nested designs like following. I don't understand why the matrix (t(X) %*% X) is singular. Can somebody help me understand it? Regards, Peng a=2 b=3 n=4 A = as.vector(sapply(1:a,function(x){rep(x,b*n)})) B = as.vector(sapply(1:(a*b), function(x){rep(x,n)}))

Re: [R] Function similar to cumsum/cumprod

Try this; set.seed(123) x - rnorm(100) sapply(30:length(x), function(i)mean(x[1:i])) On Tue, Sep 22, 2009 at 11:12 AM, Sergey Goriatchev serg...@gmail.com wrote: Hello, everyone I wonder if there is in R somewhere a function similar to cumsum(). The function calculates a statistic (say mean

Re: [R] Function similar to cumsum/cumprod

Henrique, thank you! On Tue, Sep 22, 2009 at 16:32, Henrique Dallazuanna www...@gmail.com wrote: Try this; set.seed(123) x - rnorm(100) sapply(30:length(x), function(i)mean(x[1:i])) On Tue, Sep 22, 2009 at 11:12 AM, Sergey Goriatchev serg...@gmail.com wrote: Hello, everyone I wonder if

Re: [R] trying to install plotrix

David Epstein wrote: Here are the error messages: install.packages('plotrix') Warning in install.packages(plotrix) : argument 'lib' is missing: using '/Users/dbae/Library/R/library' --- Please select a CRAN mirror for use in this session --- Error in .install.macbinary(pkgs = pkgs,

Re: [R] Function similar to cumsum/cumprod

For the case of the mean, you can use the following: cummean - function (x) cumsum(x) / seq_len(length(x)) for more a general function, one option is cumFUN - function (x, FUN = mean) { sapply(seq_len(length(x)), function (i) FUN(x[1:i])) } # Examples: x - rnorm(100) cummean(x) cumFUN(x)

Re: [R] Function similar to cumsum/cumprod

Sergey Goriatchev sergeyg at gmail.com writes: Hello, everyone I wonder if there is in R somewhere a function similar to cumsum(). The function calculates a statistic (say mean or standard deviation) buy adding consequtively one more data point. So, say I have a timeseries of 100

Re: [R] converting a character vector to a function's input

Here are a few alternatives: VAR1 - VAR2 - VAR1 - 1:3 v - VAR1,VAR2,VAR3 # 1 eval(parse(text = sprintf(cbind(%s), v))) # 2 do.call(cbind, lapply(strsplit(v, ,)[[1]], get)) # 3 library(gsubfn) strapply(v, \\w+, get, combine = list, simplify = ... ~ do.call(cbind, ...)) On Tue, Sep 22, 2009 at

[R] (correction) snowfall+Rmpi: sfExport error

Sorry for wasting your time with an incorrect diagnosis in my previous mail! I'm running my script using mpirun -mp 4 and using snowfall+Rmpi on Linux 64bits. I receive the following error, that terminates the run. Any idea? I'm able to reproduce the problem with a minimal script (below). If n

[R] microarray analysis in R without replicates

Dear all, I have received a microarray data set in standard Affymetrix CEL-format consisting of only six samples without any replicates (same organism and cell type, but different individuals and different biological conditions for each individual; the same Affymetrix GeneChip platform was

[R] run R script automatically by double-clicking WinXP desktop icon

I've written a simple script that does some surveillance analysis on daily counts of walk-in clinic visits, for our county health department. (Actually, Michael Hohle's surveillance package does all the work; I just customized it a little to work with the way our data are recorded.) The output

[R] No parametric methods

Hello I am interested in finding out a method of power analysis (effect size and sample size calculation ) using R in non parametric methods? I am running R 2.8.1 running on linux open SUSE Any libraries or documentation , I was not bale to google up any. Thanks in Advance, Ben-Ari Alon, MD

Re: [R] run R script automatically by double-clicking WinXP desktop icon

Here's something I use (in a batch file): Rterm --no-restore --file=EveningStartup.r Change EveningStartup.r to your particular file. When you create the shortcut, make sure to set the working directory to where your R script is located. Then in your file, you could have the graphs write out

[R] How to get the current script file name in the script?

Hi, Right now, I hardcode the file name in the script. But when I change the script name, I have to change the file name in the script to make sure the file name inside the script is synchronized with the script name. I am wondering if there is an automatic way to get the script file name in a

Re: [R] How to get the current script file name in the script?

Try this one line R script: print(sub(.*=, , commandArgs()[4])) Also see the geopt package for fancier handling of arguments. On Tue, Sep 22, 2009 at 11:47 AM, Peng Yu pengyu...@gmail.com wrote: Hi, Right now, I hardcode the file name in the script. But when I change the script name, I have

[R] Problem in using BATCH command

I tried to use BATCH command to run an input file (a simple text file) which contained a few R commands. I first saw the BATCH help page using ?BATCH The help page mentioned the Usage as: Batch Execution of RR CMD BATCH [options] infile [outfile] The Detailed section mentioned the

Re: [R] How to handle missing values for the GeneMeta package?

I am not sure if this is the solution you want, but complete.cases will remove rows that contain any NA values at all from a matrix. These lines get rid of NA values in my little function: predmat1 - as.matrix(predmat[complete.cases(predmat,respmat),]) respmat1 -

[R] Latex using expression

Hi all, I would like to write the following latex code in a plot x=\sqrt[4]{x[0]} using x0 - rgamma(1000,2,1) y0 - rgamma(100,2,1) x - x0^0.25 y - y0^0.25  plot(x, y,xlab=expression(x=sqrt[4]{x[0]}),ylab=expression(y=sqrt[4]{y[0]})) Many thanks, Amor [[alternative HTML version

[R] tkmessageBox font size

Hi, I´m using the tkmessageBox widget of package tcltk. I would like to change the font size of the message. How can I do it? Thanks in advance, Srpd Partilhe as suas memórias com quem quer que deseja. Com quem quer que deseja.

[R] Coordinates of the current plot

Hi, Is is possible to obtain/ get coordinates of the current graphics? I am using plot function for hclust object and I want to add some additional graphics but I cannot figure out what ylim and xlim is used for actual ploting of tree. I know that I can probably used ploting of dendrogram instead

Re: [R] Coordinates of the current plot

On 9/22/2009 11:04 AM, p...@orbit.umbr.cas.cz wrote: Hi, Is is possible to obtain/ get coordinates of the current graphics? I am using plot function for hclust object and I want to add some additional graphics but I cannot figure out what ylim and xlim is used for actual ploting of tree. I know

[R] cut and re-factor data

Hello R-users, I have a data frame with a factor of ages in 5 year increments, and various count data for each age group. I only have this summary information in R at the moment. I want to create a new factor that aggregates the age factors if the existing factors have insufficient counts. Then

Re: [R] Problem in using BATCH command

You probably want to use the source() command in R. R CMD BATCH is used from the command line - it's rare that people use it (unless they really know what they are doing). /Henrik On Tue, Sep 22, 2009 at 9:01 AM, Manuj Sharma smanuj1...@yahoo.in wrote: I tried to use BATCH command to run an

[R] cox memory

If I use the following newdata=data.frame(chemo1=0, horm1=0, age1=mean(age1), grade1=0, positive1=1, size1=mean(size1)   ) then I get Error in model.frame.default(Terms, newdata, ...) :   variable lengths differ (found for 'log(tstop + 1)') In addition: Warning message: 'newdata' had 1 rows

[R] glm analysis repeated for 900 variables

Dear R users, Could you help my with the following problem? I want to repeat a glm analysis with 2 independent variables for all 900 variables (snps) in my data set. So, I want to check whether snp1 has a different effect on my outcome variable in patients and controls(phenotype). And repeat that

[R] R PNG graph width limitation

Hi folks, I am trying to do a clustering and generate a long dendrogram in R on Linux server: = data-read.table(file=mean_ratio.txt.noheader,row.names=1,sep=\t,quote=\) hc-hclust(dist(data),ward) png(file=mean_ratio.txt.noheader.ward.png,bg=white,pointsize=8,width=32767,height=1536)

Re: [R] Latex using expression

Try this: plot(x, y, xlab=expression(x=sqrt(x[0], 4)),ylab=expression(y=sqrt(x[0], 4))) On Tue, Sep 22, 2009 at 1:07 PM, amor Gandhi amorigan...@yahoo.de wrote: Hi all, I would like to write the following latex code in a plot x=\sqrt[4]{x[0]} using x0 - rgamma(1000,2,1) y0 -

[R] How to combine matrices?

Dear All, Let a, b and c are three matrices with same no. of column but different no. of row. a - matrix(1, 1, 2) b - matrix(2, 2, 2) c - matrix(3, 3, 2) Could someone help me to combine these matrices together as a single matrix? Thank you Fir

Re: [R] Create directory and copy files in R

You assigned values to 'fn' and 'dpath' outside the function body and use these names as the arguments of your function, so actually 'cyfun' does not know what 'fn' and 'path' are. The solution is to put them in the arguments as default values: cyfun - function(fn = C:/Documents and

Re: [R] How to combine matrices?

is rbind(a,b,c) what you mean? Daniel - cuncta stricte discussurus - -Ursprüngliche Nachricht- Von: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Im Auftrag von FMH Gesendet: Tuesday, September 22, 2009 12:40 PM An:

Re: [R] How to combine matrices?

abc - rbind(a, b, c) On Tue, Sep 22, 2009 at 1:40 PM, FMH kagba2...@yahoo.com wrote: Dear All, Let a, b and c are three matrices with same no. of column but different no. of row. a - matrix(1, 1, 2) b - matrix(2, 2, 2) c - matrix(3, 3, 2) Could someone help me to combine these matrices

Re: [R] Create directory and copy files in R

See copyDirectory() in the R.utils package. /hb On Tue, Sep 22, 2009 at 9:43 AM, Yihui Xie xieyi...@gmail.com wrote: You assigned values to 'fn' and 'dpath' outside the function body and use these names as the arguments of your function, so actually 'cyfun' does not know what 'fn' and 'path'

Re: [R] How to combine matrices?

rbind(a, b, c) another function you might be interested in is cbind. cheers, On Tue, Sep 22, 2009 at 9:40 AM, FMH kagba2...@yahoo.com wrote: Dear All, Let a, b and c are three matrices with same no. of column but different no. of row. a - matrix(1, 1, 2) b - matrix(2, 2, 2) c -

Re: [R] cox memory

I did not test to see if the design of the data.frame you offered was valid input to your model, since you did not provide a means to test it within R. I only meant to say that creating a very long vector with identical entries should not be necessary if the newdata argument were

Re: [R] R PNG graph width limitation

Given the DPI=72, do you really need a graph that's wider than 450 inches? Or can you really read a picture that is so wide? Regards, Yihui -- Yihui Xie xieyi...@gmail.com Phone: 515-294-6609 Web: http://yihui.name Department of Statistics, Iowa State University 3211 Snedecor Hall, Ames, IA On

Re: [R] Problem in using BATCH command

BATCH is not an R object, and it is not for use in R console. As Henrik mentioned, you should use it in the command line. Open a terminal window/command window (depending on your OS), and type R CMD BATCH ... there. If you are using Windows, cd to the 'bin' directory of R first (or put the 'bin'

[R] How to sort the elements in a matrix?

Dear All, Let ab is the combination of matrices a and b. a - matrix(c(1,2,3,4), nrow = 2, ncol = 2) b - matrix(c(4,3,1,2), nrow = 2, ncol = 2) ab - rbind(a, b)   From matrix ab, could someone give some advice on the way to sort all elements in column 2, by the sequence of the elements in 

Re: [R] How to sort the elements in a matrix?

Hi, On Sep 22, 2009, at 1:20 PM, FMH wrote: Dear All, Let ab is the combination of matrices a and b. a - matrix(c(1,2,3,4), nrow = 2, ncol = 2) b - matrix(c(4,3,1,2), nrow = 2, ncol = 2) ab - rbind(a, b) From matrix ab, could someone give some advice on the way to sort all elements in

Re: [R] How to sort the elements in a matrix?

ab[order(ab[,1]),] - cuncta stricte discussurus - -Ursprüngliche Nachricht- Von: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Im Auftrag von FMH Gesendet: Tuesday, September 22, 2009 1:20 PM An: r-help@r-project.org

Re: [R] run R script automatically by double-clicking WinXP desktop icon

that helps, thanks. I can put a final line in my batch file that opens a viewer for the png graph. I was hoping to find a way to do make it run in the Rgui (picky, I know.) My graph is clearer in the default graph device that pops up in the Rgui; it's blurry in the Windows viewer for png

Re: [R] Create directory and copy files in R

Tammy Try avoid directories with spaces... may be this solve your problem. Case not, write us again. bests milton On Mon, Sep 21, 2009 at 10:05 AM, Tammy Ma metal_lical...@live.com wrote: HI, All R users, My problem is: fn [1] C:/Documents and

Re: [R] Plotting waver files

blahbleh wrote: Is there a good tutorial or just sample code on how to plot .wav files in R? I would really appreciate it! See for example packages tuneR or seewave Best wishes, Uwe Ligges __ R-help@r-project.org mailing list

Re: [R] How to sort the elements in a matrix?

On Sep 22, 2009, at 1:20 PM, FMH wrote: Dear All, Let ab is the combination of matrices a and b. a - matrix(c(1,2,3,4), nrow = 2, ncol = 2) b - matrix(c(4,3,1,2), nrow = 2, ncol = 2) ab - rbind(a, b) From matrix ab, could someone give some advice on the way to sort all elements in column

Re: [R] R PNG graph width limitation

Well, if you don't care about its width, I'd suggest you use the pdf() device instead, e.g. # 200 inches! pdf(hugeplot.pdf, width = 200, height = 200) par(mar = rep(0, 4)) plot(rnorm(1), pch = 19) dev.off() Regards, Yihui -- Yihui Xie xieyi...@gmail.com Phone: 515-294-6609 Web:

[R] any way to make it work faster (deleting rows that contain certain values)

Hello, dear R'ers, index-expand.grid(1:7,1:4,1:4,1:4,1:4,1:4,1:4,1:4,1:4,1:4,1:4) In this case, dim(index) is 7,340,032 (!) and 11. I realize it's huge. Then, I am trying to get rid of the undesired combinations of columns. They should not contain identical values in any 2 columns. Also if

[R] how to force R to evaluate variable?

Hi, I want to construct a list as Lst - list(name_1=object_1, ..., name_m=object_m) If name_1 is a variable with value NAME1, how can I ask R to use NAME1 instead of 'name_1' as the name of the list element? Thanks! Yupu [[alternative HTML version deleted]]

Re: [R] how to force R to evaluate variable?

Try this: as.list(c(object_1, object_2)) On Tue, Sep 22, 2009 at 3:23 PM, YUPU LIANG li...@cbio.mskcc.org wrote: Hi, I want to construct a list as Lst - list(name_1=object_1, ..., name_m=object_m) If name_1 is a variable with value NAME1, how can I ask R to use NAME1 instead of 'name_1'

Re: [R] how to force R to evaluate variable?

On 9/22/2009 2:23 PM, YUPU LIANG wrote: Hi, I want to construct a list as Lst - list(name_1=object_1, ..., name_m=object_m) If name_1 is a variable with value NAME1, how can I ask R to use NAME1 instead of 'name_1' as the name of the list element? Lst - list(object_1, ..., object_m)

Re: [R] microarray analysis in R without replicates

Try asking bioconductor list http://www.bioconductor.org/docs/postingGuide.html -Johannes 2009/9/22 Rainer Tischler rainer_...@yahoo.de: Dear all, I have received a microarray data set in standard Affymetrix CEL-format consisting of only six samples  without any replicates (same organism

Re: [R] cut and re-factor data

Try e.g. function recode from package car. 2009/9/22 Chris Hane christopher.a.h...@gmail.com: Hello R-users, I have a data frame with a factor of ages in 5 year increments, and various count data for each age group. I only have this summary information in R at the moment. I want to create a

[R] Evaluating expresssions as parameter values

I need to play games with an expression similar to the following one: print(xyplot(DepVar ~ Group|Covar, groups=Othergroup, data=mydf, pch = 18 ,main=Testcase,auto.key = TRUE)) The problem is that the formula argument (the first argument) an the groups argument are passed over from

[R] how to convert character string with only month and year into date

Dear R helpers. I am new to plotting time data using R. wonder how to convert character time info into date in R. I searched over the web but did not find answer. the input character string is something like 03_1993 or 03-1993, so the precision is at month level. I tried the following but

Re: [R] R PNG graph width limitation

You *have to* report sessionInfo() when you ask for this kind of help. It might be that the PNG is actually generated, but it is the viewer that cannot display it. You can check with other viewers or editors. You can also try to downscale using ImageMagick, e.g. convert -geometry 10% large.png

Re: [R] how to convert character string with only month and year into date

Try this: strptime(paste(1, c(03_1993)),%d %m_%Y) strptime(paste(1, c(03-1993)),%d %m_%Y) On Tue, Sep 22, 2009 at 4:03 PM, Sean Zhang seane...@gmail.com wrote: Dear R helpers. I am new to plotting time data using R. wonder how to convert character time info into date in R. I searched over

Re: [R] Evaluating expresssions as parameter values

Try this: f - function(fo, data, groups) { g - do.call(xyplot, list(as.formula(fo), groups = as.name(groups), data)) print(g) } f(yield ~ variety | site, data = barley, groups = year) On Tue, Sep 22, 2009 at 3:03 PM, Erich Neuwirth erich.neuwi...@univie.ac.at wrote: I need

Re: [R] how to convert character string with only month and year into date

Try this: library(zoo) # as.yearmon as.yearmon(03_1993, %m_%Y) [1] Mar 1993 as.Date(as.yearmon(03_1993, %m_%Y)) [1] 1993-03-01 On Tue, Sep 22, 2009 at 3:03 PM, Sean Zhang seane...@gmail.com wrote: Dear R helpers. I am new to plotting time data using R. wonder how to convert character

Re: [R] how to convert character string with only month and year into date

On Sep 22, 2009, at 3:03 PM, Sean Zhang wrote: Dear R helpers. I am new to plotting time data using R. wonder how to convert character time info into date in R. I searched over the web but did not find answer. the input character string is something like 03_1993 or 03-1993, so the

Re: [R] how to convert character string with only month and year into date

David, Gabor, and Henriuqe, Thanks a lot for help! Another (inelegant) way is to use ts() and then supply the start and end time. this inelegant way works (I guess at least for equally spaced data.) . -Sean On Tue, Sep 22, 2009 at 3:19 PM, David Winsemius dwinsem...@comcast.netwrote: On Sep

[R] Help with output height

Hello everyone I have a problem using the R console. I have written 2000 lines in a notepad, this lines are commands for reading dbf files. When I paste these lines in the R console I can only visualize the end of the output, this is because the R console has a default number of rows which is

  1   2   >