[R] Howto Compute Pairwise Similarity/Correlation Matrix from a Data Frame

Hi, I have the following 5 vectors. I wish to compute the pairwise Pearson Correlation matrix with this data.frame. Is there a compact way to do it? At the end I hope to create a heatmap out of this correlation matrix. __BEGIN__ > data <- read.table("mydata.txt") > print(data) V1

Re: [R] read.spss {foreign} doesn't work over network?

This looks like a restriction on your network (assuming the file is named correctly). I can read SPSS files with read.spss over ours, and the code in read.spss is just standard C I/O (fopen) that knows nothing about Windows network mappings. On Tue, 17 Jun 2008, Farley, Robert wrote: I'm un

[R] matrix

Hello everyone,I am having a difficult time understanding what commands are used to add and remove elements froma matrix or data frame.For ex ample if I want to remove element=[1,50] or elements [1,50:63] or elements  [2:4,20:35]from a matrix and replace with NA or any number what command must I

Re: [R] x labels out of Quartz canvas

On Tue, 17 Jun 2008, MeMooMeM wrote: Thanks a lot! I looked at the document. It shows how to set the size of the canvas, I don't believe it does. but not how to change it *after* plotting. Now I start with a bigger canvas, but the plot is scaled into it, so I am having the same problem ag

Re: [R] Using the shape () function

hadley wickham wrote: You might try using the reshape package instead: last <- function(x) x[length(x)] names(d) <- c("value", "person", "time") cast(d, person ~ time, last) The first and the last line I think is clear, although I will have to experiment more to understand the call on cast () b

[R] Problem in Binning of a data set

Hello, I am having problem with binning the data. I have a 50X3 matrix and I binned the data for all the 3 columns. Using table command I got the total no. of elements in a particular bin. Could you please tell me how to see that what all elements are there in a particular bin and then create a di

Re: [R] Differential Equations

On Tue, Jun 17, 2008 at 6:25 PM, David Arnold <[EMAIL PROTECTED]> wrote: > All, > > I've found odesolve and lsoda. > > Any other packages for differential equations? > There's Rsundials, which gives you an algebraic ode solver. > > Any good tutorials on using R and solving differential and part

Re: [R] Sweave: problem with usepackage{C:/PROGRA~1/R/R-27~1.0/share/texmf/Sweave}

From the NEWS file for 2.7.0: o The default for 'stylepath' in Sweave's (default) RweaveLatex driver can be set by the environment variable SWEAVE_STYLEPATH_DEFAULT: see ?RweaveLatex. You want to set 'stylepath' to be false, and the easiest way to do this is to set that e

Re: [R] longest common substring

See: http://www.omegahat.org/Rlibstree/ Binds R to libstree for suffix tree operations. Libstree is included with the package, so don't worry about building it separately. Michael On Tue, Jun 17, 2008 at 10:28 PM, Daren Tan <[EMAIL PROTECTED]> wrote: > > i need to compute the longest common su

[R] longest common substring

i need to compute the longest common substring of two strings, can R do that ? _ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

Re: [R] Sweave: problem with usepackage{C:/PROGRA~1/R/R-27~1.0/share/texmf/Sweave}

Thanks. Actually I didn't get the point from Duncan's post but I think I found a workaround: \usepackage{C:/progra\string~1/R/R-27\string~1.0/share/texmf/Sweave} Cheers, Lauri 2008/6/17, Gabor Grothendieck <[EMAIL PROTECTED]>: > See: > > https://stat.ethz.ch/pipermail/r-devel/2008-May/049609.h

Re: [R] keeping original order in factor()

Thanks for everyone's suggestions. I think factor(foo, levels = unique(foo)) works best for my needs. (By the way, I'm still trying to figure out how to use the ordered option in factor().) Andrew On Wed, Jun 18, 2008 at 12:32 AM, Peter Alspach <[EMAIL PROTECTED]> wrote: > Andrew > > levels(fac

Re: [R] keeping original order in factor()

Andrew levels(factor(foo, levels=c('b','a'))) should work. You can make foo an ordered factor too, but that is not necessary. HTH ... Peter Alspach > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Andrew Yee > Sent: Wednesday, 18 June 2008

Re: [R] keeping original order in factor()

> foo <- sample(letters, 10, replace = TRUE) > foo [1] "i" "j" "u" "f" "a" "k" "m" "k" "i" "w" > unique(foo) [1] "i" "j" "u" "f" "a" "k" "m" "w" > bar <- factor(foo, levels = unique(foo)) > bar [1] i j u f a k m k i w Levels: i j u f a k m w > Is this what you are trying to do? HTH Steve McKi

[R] keeping original order in factor()

Apologies for the naieveness of this question, but I'm having trouble figuring out to have factor() maintain original ordering. For example, foo <- c("b","b","a","a") levels(factor(foo, ordered=T)) #I'd like this to return as "b" "a" #not "a" "b" I thoug

Re: [R] try catch block

This is a useful reference on the tryCatch() function: http://www1.maths.lth.se/help/R/ExceptionHandlingInR/ I think something like this should work: xml <- tryCatch(xmlTreeParse(xmlTxt, useInternal=TRUE), error=function(err) xmlMalFormed()) - Original Message From: ppatel3026 <[E

Re: [R] combining two data frames

Hi, haibing . you should try rbind or cbind. Type >help(rbind) or >? rbind in the command line . regards. On 2008-6-18, at 上午6:33, calundergrad wrote: i know this is a very simple question but i have two data frames (one of which is a continuation of the first data frame) and i was just

[R] combining two data frames (different question)

i have two data frames. One data frame contains one column with the identificication number of a geographic place. The other columns are just some data about the corresponding geographic place. The other data frame also contains one column with the id number of a geographic place. the other co

[R] combining two data frames

i know this is a very simple question but i have two data frames (one of which is a continuation of the first data frame) and i was just wondering how do you combine those? any simple solution would be helpful. thank you very much -- View this message in context: http://www.nabble.com/combinin

[R] Deming Regression

Hi all. Has anyone ever done a Deming Regression in R? I'm wondering if there's a simple way to do it. Thanks for all your help! -Ed -- View this message in context: http://www.nabble.com/Deming-Regression-tp17949318p17949318.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] Differential Equations

hi , David. I think odesolve is what you need , do you ever try it ? Solving differential equations is not a strong function of R, that means, maybe you should try to write your own procedure using your own numerical method for your own specific equation . regards. On 2008-6-18, at 上午9:25,

Re: [R] subarray

If you are trying to subset the matrix, the correct format is: Y <- dat[, 46:63] You probably need to reread the introduction to R to understand how indexing works. On Tue, Jun 17, 2008 at 10:01 PM, Paul Adams <[EMAIL PROTECTED]> wrote: > Hello everyone, > I am trying to subset a matrix with the

Re: [R] Capturing coxph warnings and errors

On Tue, Jun 17, 2008 at 9:47 AM, Daniel Brewer <[EMAIL PROTECTED]> wrote: > Hi, > > I have a script that takes a subset of genes on a microarray and tries > to fit a coxph model to the expression values for each gene. This seems > to work fine but in some cases it produces warnings and/or errors.

[R] subarray

Hello everyone, I am trying to subset a matrix with the following code: dat<-read.table(file="C:\\Documents and Settings\\Owner\\My Documents\\Yeast\\Yeast.txt",header=T,row.names=1) file.show(file="C:\\Documents and Settings\\Owner\\My Documents\\Yeast\\Yeast.txt") Z.matrix<-as.matrix(dat) Y<-ar

Re: [R] printing in binary

I think that this one finds a binary representation of integer but does not print data in binary format (creating binary file, like write(unit,*) in FORTRAN). --- On Wed, 18/6/08, jim holtman <[EMAIL PROTECTED]> wrote: > From: jim holtman <[EMAIL PROTECTED]> > Subject: Re: [R] printing in bina

[R] Differential Equations

All, I've found odesolve and lsoda. Any other packages for differential equations? Any good tutorials on using R and solving differential and partial differential equations? D. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/lis

Re: [R] printing in binary

check out the as.character.binmode in the R.util package: > ?as.character.binmode > intToBin(12345) [1] "1100111001" > On Tue, Jun 17, 2008 at 8:05 PM, Edna Bell <[EMAIL PROTECTED]> wrote: > Hi r Gurus: > > Is there a way to print numbers in binary format, please? > > I know that you can use

[R] printing in binary

Hi r Gurus: Is there a way to print numbers in binary format, please? I know that you can use sprintf for hex, but I'm not sure how to get binary. Thanks Edna __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] PCA analysis

I am not entirely sure after reading your email, but I thought you wanted to do something like this: ###Start of example ###create random data for the example x=rnorm(100,100,10) ##create Xs e=rnorm(100,0,5) ##create Errors y=x+e##create Ys ###plot plot(y~x,pch=NA) ##plo

Re: [R] error in casting matrix

Paul Adams wrote: Hello everyone, I am trying to cast a data frame to a matrix I have used the following code: dat<-read.table(file="C:\\Documents and Settings\.txt",header=T,row.names=1) file.show(file="C:\\Documents and Settingt.txt") Z.matrix<-as.matrix(dat.data.frame) but I keep g

[R] error in casting matrix

Hello everyone, I am trying to cast a data frame to a matrix I have used the following code: dat<-read.table(file="C:\\Documents and Settings\.txt",header=T,row.names=1) file.show(file="C:\\Documents and Settingt.txt") Z.matrix<-as.matrix(dat.data.frame) but I keep getting the error error

Re: [R] Measuring dispersion

Thanks for the suggestion, however I'm looking for a score since my goal is to rank thousands of distributions. For instance, given a large text, I would like to rank all terms according to their distribution (dispersion) within the text. Terms evenly distributed in the text should have a low scor

Re: [R] Measuring dispersion

You could also look at the difference between your empirical distribution and the uniform distribution (something like Kolmogorov-Smirnov test). --- On Tue, 17/6/08, S. Nunes <[EMAIL PROTECTED]> wrote: > From: S. Nunes <[EMAIL PROTECTED]> > Subject: [R] Measuring dispersion > To: [EMAIL PROTECT

Re: [R] Reshape or Stack? (To produce output as columns)

>> It sounds like something is going wrong with the melting. Could you >> please include the output of str(original data frame), and >> str(melted)? (Or even better a small version of your data created >> with dput) > > And this is a str output of the original data frame (first few rows of that >

Re: [R] array elements incorrect

If ind is the vector of rows of dat you are interested in then dat[ind,] is what you need. --- On Wed, 18/6/08, Paul Adams <[EMAIL PROTECTED]> wrote: > From: Paul Adams <[EMAIL PROTECTED]> > Subject: [R] array elements incorrect > To: r-help@r-project.org > Received: Wednesday, 18 June, 2008, 8

Re: [R] how to fit a curve of form Y = X^Z

I see two possibilities: 1) Taking logarithm yields log(Y) = log(X)*Z and this is the regular linear regression with intercept = 0, and in this case Z = Sum(log(Xi)*log(Yi))/Sum(log(Xi)^2). This is very simple but not necessarily what you want (but this solution can be used as a starting point

[R] array elements incorrect

Hello everyone, I have a question as to what code should be used if one wanted to subset from a dataframe of 7000rows by 38 columns.For example if you wanted to generate a array of 100 by 38 would you not use the following: Z<-(dat,dim=c(100,38)) where dat is the dataframe of 7000by 38 Whenever I

Re: [R] try catch block

ppatel3026 <[EMAIL PROTECTED]> writes: > How can I use the try catch block such that if this statement fails > xml <- xmlTreeParse(xmlTxt, useInternal=TRUE) > > then this statement is executed > > xml <- xmlMalFormed() library(XML) xml <- tryCatch({ # other code, then... xml

Re: [R] How to control height of abline

Thanks Dimitris Rizopoulos wrote: > > check this: > > x <- rnorm(200) > dd <- density(x) > plot(dd) > ind <- seq(100, 400, len = 6) > arrows(dd$x[ind], 0, dd$x[ind], dd$y[ind] - 0.015, length = 0.2) > > > I hope it helps. > > Best, > Dimitris > > > Dimitris Rizopoulos > Biostatistical

Re: [R] Simultaneous Confidence/Prediction Bands

On 18 Jun 2008 at 9:51, Rolf Turner wrote: > > I'm sure that others will correct me if I am wrong ... :-) ... but > are these pointwise rather than simultaneous bands? > > I'm not sure how one would go about getting simultaneous bands. > Isn't the Working-Hotelling approach one such possibilit

[R] color2D.matplot axis names

Hi, I finally came up with a nice colored matrixplot, using the color2D.matplot function of the plotrix package. But I can't assign xtics and ytics to this plot. I made sure that the matrix has correct colnames() and rownames(). Here's what I do: a=matrix(1:16, 4, 4) colnames(a) <- c("X1", "X2"

Re: [R] Simultaneous Confidence/Prediction Bands

I'm sure that others will correct me if I am wrong ... :-) ... but are these pointwise rather than simultaneous bands? I'm not sure how one would go about getting simultaneous bands. cheers, Rolf Turner On 18/06/2008, at 9:32 AM, Jorge Ivan Velez wrote: Dear Tom, Se

Re: [R] Simultaneous Confidence/Prediction Bands

Dear Tom, See the examples in ?predict.lm HTH, Jorge On Tue, Jun 17, 2008 at 5:29 PM, Tom La Bone <[EMAIL PROTECTED]> wrote: > > Is there a built-in function in R that will generate simultaneous > confidence > and prediction bands for linear regression? > > Tom > -- > View this message in cont

[R] Simultaneous Confidence/Prediction Bands

Is there a built-in function in R that will generate simultaneous confidence and prediction bands for linear regression? Tom -- View this message in context: http://www.nabble.com/Simultaneous-Confidence-Prediction-Bands-tp17941537p17941537.html Sent from the R help mailing list archive at Nabb

[R] try catch block

How can I use the try catch block such that if this statement fails xml <- xmlTreeParse(xmlTxt, useInternal=TRUE) then this statement is executed xml <- xmlMalFormed() ? This code does not work but assuming its somewhere along these lines: tryCatch(xml <- xmlTreeParse(xmlTxt, useInternal=TRU

Re: [R] Row and Col length of a matrix?

Oh, and don't forget: NROW and NCOL url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Champaign, IL 61820 On Jun 17, 2008, a

Re: [R] Reshape or Stack? (To produce output as columns)

> > Dear all, > > > > Many thanks for the suggestions put forward. I've decided to go with the > > 'melt' command from the 'reshape' library, as this seems to run the > > quickest. > > > > I do have a couple of questions however, regarding the use of the 'melt' > > command. Below are the last

Re: [R] Row and Col length of a matrix?

See ?dim, ?nrow, ?ncol. MeMooMeM wrote: Hi, Is there a direct command to get the row and col length of a matrix? My ugly solution is: length(mtx[,1]) and length(bcf[1,]) Thanks! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

Re: [R] Row and Col length of a matrix?

See ?dim HTH, Jorge On Tue, Jun 17, 2008 at 4:50 PM, MeMooMeM <[EMAIL PROTECTED]> wrote: > > Hi, > > Is there a direct command to get the row and col length of a matrix? > > My ugly solution is: length(mtx[,1]) and length(bcf[1,]) > > Thanks! > > -- > View this message in context: > http://www

[R] Row and Col length of a matrix?

Hi, Is there a direct command to get the row and col length of a matrix? My ugly solution is: length(mtx[,1]) and length(bcf[1,]) Thanks! -- View this message in context: http://www.nabble.com/Row-and-Col-length-of-a-matrix--tp17936978p17936978.html Sent from the R help mailing list archiv

[R] echo Rprofile.site on startup

Hello, I have not been able to find how to have R echo the Rprofile.site contents to the console on start-up. Changing the windoze shortcut to include the '--verbose' option does not show what I am looking for. I would like to see the options and commands listed, in addition to the library() loadi

Re: [R] Using the shape () function

>> You might try using the reshape package instead: >> >> last <- function(x) x[length(x)] >> names(d) <- c("value", "person", "time") >> cast(d, person ~ time, last) > > The first and the last line I think is clear, although I will have to > experiment more to understand the call on cast () better

Re: [R] Sweave: problem with usepackage{C:/PROGRA~1/R/R-27~1.0/share/texmf/Sweave}

See: https://stat.ethz.ch/pipermail/r-devel/2008-May/049609.html On Tue, Jun 17, 2008 at 3:08 PM, Lauri Nikkinen <[EMAIL PROTECTED]> wrote: > R users, > > I'm at a loss with a problem considering running .tex files produced > by Sweave. When I run (R 2.7.0): >

[R] Build a R into a single static library on windows?

Hi R users, I've been developing a C++ library that depends on R through R.dll. The current way is to have R directories somewhere and have R_HOME pointing to it. This works. But it's inconvenient when deploying this library since the whole R distribution has to be deployed with it. It would be v

[R] Error

My code seems to break out with error below for every 1000 files it processes. Then I re-run from the last file where it errored out and it runs without any bugs. Any ideas what might cause error below? Error in match(x, table, nomatch = 0) : formal argument "nomatch" matched by multiple act

Re: [R] Reshape or Stack? (To produce output as columns)

On Tue, Jun 17, 2008 at 1:47 PM, Steve Murray <[EMAIL PROTECTED]> wrote: > > Dear all, > > Many thanks for the suggestions put forward. I've decided to go with the > 'melt' command from the 'reshape' library, as this seems to run the quickest. > > I do have a couple of questions however, regarding

[R] some R code of linear mixed model

Hello everyone: I have some quesions about the R code for linear mixed model,hope some one can give me some hints,thanks a lot~ Say: we have A B C three factors (i)A is fixed ,B and C are random,and B is nested in C,so the R code for this case would be: case1<-lme(y~A+B+C+..,rando

Re: [R] constrOptim with method = "L-BFGS-B"

I believe that 'optim' will not accept equality constraints. However, you do not need the generality of 'optim' to "minimize a quadratic function with boundary conditions and one equality condition". This type of problem is called "quadratic programming", and RSiteSearch("quadratic

Re: [R] x labels out of Quartz canvas

Thanks a lot! I looked at the document. It shows how to set the size of the canvas, but not how to change it *after* plotting. Now I start with a bigger canvas, but the plot is scaled into it, so I am having the same problem again :( I keep on reading tough. Thanks a lot for your help! -Memo

Re: [R] Sweave: problem withusepackage{C:/PROGRA~1/R/R-27~1.0/share/texmf/Sweave}

Thank you for your answer. Unfortunately I don't have administrator rights to my computer (company's PC), so I'm looking for some other solution. Cheers, Lauri 2008/6/17 Scillieri, John <[EMAIL PROTECTED]>: > Try moving your R directory out of the C:\Program Files directory. > > I ended up settin

Re: [R] How to control height of abline

I'm sorry Sata, I typed your name wrongly in the my previous email. My apologizes. Thanks, Jorge On Tue, Jun 17, 2008 at 3:19 PM, Jorge Ivan Velez <[EMAIL PROTECTED]> wrote: > > Dear Satam > > Just take a look at Dimitris Rizopoulos' suggestion at > http://www.nabble.com/How-to-control-height-

Re: [R] How to control height of abline

Dear Satam Just take a look at Dimitris Rizopoulos' suggestion at http://www.nabble.com/How-to-control-height-of-abline-tp17932528p17932528.html HTH, Jorge On Tue, Jun 17, 2008 at 2:15 PM, sata pinal <[EMAIL PROTECTED]> wrote: > I use matplot to get the density curve and then I use > > ablin

Re: [R] Sweave: problem withusepackage{C:/PROGRA~1/R/R-27~1.0/share/texmf/Sweave}

Try moving your R directory out of the C:\Program Files directory. I ended up setting up a C:\programs\ and putting all my unix-happy programs (R, latex, etc) there to ensure path spaces don't screw me up when running under windows. -Original Message- From: [EMAIL PROTECTED] [mailto:

[R] How to control height of abline

I use matplot to get the density curve and then I use abline(v=g$V2, col = 3 ) to get the vertical line. Goal : I want very small lines at the bottom on the x axis , if possible in the arrow forms instead of vertical lines on the whole graph. Thanks a lot. [[alternative

[R] Sweave: problem with usepackage{C:/PROGRA~1/R/R-27~1.0/share/texmf/Sweave}

R users, I'm at a loss with a problem considering running .tex files produced by Sweave. When I run (R 2.7.0): --- #Taken from ?Sweave testfile <- system.file("Sweave", "Sweave-test-1.Rnw", package = "utils") ## enforce par(ask=FALSE) options(device

Re: [R] How to control height of abline

check this: x <- rnorm(200) dd <- density(x) plot(dd) ind <- seq(100, 400, len = 6) arrows(dd$x[ind], 0, dd$x[ind], dd$y[ind] - 0.015, length = 0.2) I hope it helps. Best, Dimitris Dimitris Rizopoulos Biostatistical Centre School of Public Health Catholic University of Leuven Address: K

[R] read.spss {foreign} doesn't work over network?

I'm unable to open an SPSS file over my network. If I copy it to my local C:/ drive I can read it. I saved the command (in a "crib sheet" text file) in order to avoid all the typing, so I'm pretty sure I've done it before. I verified that the file I'm trying to read is OK. This is what happens:

Re: [R] How to control height of abline

mogra wrote: > I use plot to get the density curve and then I use > > abline(v=g$V2, col = 3 ) to get the vertical line for specific point on x > axis. > > Goal : I want very small lines at the bottom on the x axis , if possible in > the arrow forms instead of vertical lines on the whole graph. >

Re: [R] Reshape or Stack? (To produce output as columns)

Dear all, Many thanks for the suggestions put forward. I've decided to go with the 'melt' command from the 'reshape' library, as this seems to run the quickest. I do have a couple of questions however, regarding the use of the 'melt' command. Below are the last few lines of the 'melted' data.

Re: [R] Problems with Cochrane-Orcutt procedures

Aha ! I think this is it. You are using a data frame, and I was using zoo for the regrCMS object. When I cast it into a data frame, it then works ! I guess it would be cool to understand why it doesn't work as a zoo object, but no matter. Thanks John, Tolga > regrCMS<-as.data.frame(regrCMS) >

[R] How to control height of abline

I use plot to get the density curve and then I use abline(v=g$V2, col = 3 ) to get the vertical line for specific point on x axis. Goal : I want very small lines at the bottom on the x axis , if possible in the arrow forms instead of vertical lines on the whole graph. Thanks a lot. -- View th

Re: [R] Problems with Cochrane-Orcutt procedures

Dear Tolga, I'm afraid that your data work fine for me: > regrCMSlm <- lm(regrCMS[,1] ~ regrCMS[,2]) > cochrane.orcutt.lm(regrCMSlm) $coefficients (Intercept) regrCMS[, 2] 23.5679065 -0.1784187 $cov (Intercept) regrCMS[, 2] (Intercept)60.449366 -2.14491371 regrCMS[, 2]

Re: [R] Using the shape () function

hadley wickham wrote: On Tue, Jun 17, 2008 at 9:28 AM, Tom Backer Johnsen <[EMAIL PROTECTED]> wrote: In a research project we are using a web-based tools for collecting data from questionnaire. The system generates files that are simple to read as a data frame in the "long" format, which are si

Re: [R] Problems with Cochrane-Orcutt procedures

Sure, of course. And thanks for looking John. Here is the entire data: > regrCMS a b 09/20/07 26.084 28.40 09/21/07 22.458 28.90 09/24/07 21.297 29.25 09/25/07 21.733 29.40 09/26/07 21.319 28.75 09/27/07 22.507 28.85 09/28/07 19.571 28.90 10/01/07 21.961 29.00 10/02/07 21.729 28

Re: [R] OpenGL and Linux

Floyd poole wrote: Anyone know how to get OpenGL for redhat fedora 8/9? Just install these: [EMAIL PROTECTED] R]$ rpm -qa | grep mesa mesa-libGLU-7.1-0.31.fc9.i386 mesa-libGL-7.1-0.31.fc9.i386 mesa-libGL-devel-7.1-0.31.fc9.i386 mesa-libGLU-devel-7.1-0.31.fc9.i386 -- O__ Peter Dalg

Re: [R] Problems with Cochrane-Orcutt procedures

Dear Tolga, That's a little more information, but because the code seems to work for me on other data (though no longer the message dispatch), I can't say what produces the error. I guess that if you can't debug this yourself, you'll have to share the data (generally a good idea in any event). >

Re: [R] A new task view on survival analysis

And here's the link: http://cran.r-project.org/web/views/Survival.html Hadley On Tue, Jun 17, 2008 at 12:33 PM, Arthur Allignol <[EMAIL PROTECTED]> wrote: > Dear all, > > A new task view on survival analysis > is now online. > It attempts to deal with all the R-packages > that permit to analyze t

[R] OpenGL and Linux

Anyone know how to get OpenGL for redhat fedora 8/9? thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] A new task view on survival analysis

Dear all, A new task view on survival analysis is now online. It attempts to deal with all the R-packages that permit to analyze time-to-event data. Any comments or suggestions to improve the task view are very welcome. Best regards, Arthur Allignol Freiburg Center for Data Analysis and Modeli

Re: [R] Problems with Cochrane-Orcutt procedures

Sure, I can imagine GLS is a much better way to deal with this. I guess I was looking at this because I did try GLS but got exactly the same results as LM and I just wanted to be sure. I did "debug" the code in https://stat.ethz.ch/pipermail/r-help/2002-January/017774.html and the offending li

Re: [R] Scan document including "\n"

Hi there, Try this: your.file=read.table(textConnection(" "),header=FALSE) paste(your.file$V1,"\n",collapse="",sep="") [1] "\n\n\n\n" HTH, Jorge On Tue, Jun 17, 2008 at 1:01 PM, ppatel3026 <[EMAIL PROTECTED]> wrote: > > How do you read in a whole file while preserving end of line "\n" > c

[R] Scan document including "\n"

How do you read in a whole file while preserving end of line "\n" characters? Basically, read in a whole file as one string. Ex: After this file is read into a variable, it should really look like "\n\n\n\n" -- View this message in context: http://www.nabble.com/Scan-document-including-

Re: [R] Z test and proportions

Michael Pearmain wrote: > Yes my mistake, > > I looked at the pwr.2p2n.test but i cannot place both n's and both p > values to determine the sig value > e,g *pwr.2p2n.test(h = , n1 = , n2 = , sig.level = , power = ) > > or am i missing someting obvious? > > * Not quite obvious, but h is the effect

Re: [R] Problems with Cochrane-Orcutt procedures

Dear Tolga, I'm afraid that I don't see an error. (I expect in any event that the Cochrane-Orcott and Prais estimators are now only of historical interest.) Regards, John -- John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada web

[R] useR! 2008: program online

Dear useRs, we have finally put together the program for the useR! 2008 conference, it is available online at http://www.R-project.org/useR-2008/program.html We think the contributions provide an exciting program with a rich variety in the kaleidoscope sessions and many classical as well as

Re: [R] how to fit a curve of form Y = X^Z

> I have two variables X and Y, and think that Y is related > to X by a function of the form : Y = X^Z, where Z is < 1. > However, I'm not sure how to find the best-fit equation to > fit my data to a curve of this form using R. Have you any ideas? You can use nlm() to fit a non-linear model. An

Re: [R] alternative to matching/merge?

Jim, I understand that that do.call will combine all the lists into a matrix. But I have a list,x, which is a list of a list as in x[[1]] and x[[2]] which are themselves lists. Then using do.call(rbind, x) would combine both x[[1]] and x[[2]] into a matrix. I want to keep x[[1]] and x[[2]] sep

[R] The assign function in R

Hello, I want to convert assign("a", b, where =1 ) from SPLUS to R. Is it safe to assume that the equivalent of where=1 is pos=1 in R? Thanks for help! -- View this message in context: http://www.nabble.com/The-assign-function-in-R-tp17918416p17918416.html Sent from the R help mailing list a

[R] RWeka handlers

Dear RWeka users, I was wondering how I can set Weka classifier options in make_Weka_classifier(). What I tried to do is > classifiers[[1]] <- make_Weka_classifier("weka/classifiers/lazy/IBk", handlers = Weka_control(K = 2)) but this seems to have no effect - IBk still uses K=1 (default set

Re: [R] Decision Trees RWeka

On Tue, 17 Jun 2008, [EMAIL PROTECTED] wrote: Hello, I have a question concerning decision trees coming from RWeka : library(RWeka) m =J48(Species~.,data=iris) How could such a decision tree be transferred into a matrix, pretty much in the same fashion, as it is done by getTree()

Re: [R] Z test and proportions

Yes my mistake, I looked at the pwr.2p2n.test but i cannot place both n's and both p values to determine the sig value e,g *pwr.2p2n.test(h = , n1 = , n2 = , sig.level = , power = ) or am i missing someting obvious? i did the sam ein SPSS using a macro and the following code: COMPUTE n1 = Contr

[R] R error using Survr function with gcmrec

Would someone be able to help with this question? I'm using the Gcmrec, Survrec, and Design packages to do a power analysis on simulated data. I'm receiving an error after using the Survr function that all data must have a censoring time even after using the gcmrec function: newdata<-add

[R] stop unwanted boxes/dialogs in tcl/tk

Hello, I have a question about tcl/tk: is there a way to stop more messagers/listchoice/etc. that are set up earlier, but are unwanted later? for example, require(tcltk) ttMain <- tktoplevel() tkwm.title(ttMain,"Question") f.fcn <- function(){ t1 <- modalDialogOK("Elicitation","What's the

Re: [R] Z test and proportions

Michael Pearmain wrote: > Hi All, > > I have a table based on ordial data and i want to compare proportions and > i've seen in the pwr package i can use > power.prop.test > > however i want to find out what the sig. value is based on n1,n2,p1,p2 and > this package doesn't contain this.. > Does anyo

[R] constrOptim with method = "L-BFGS-B"

Hi, i need to minimize a quadratic function with boundary condidtions and one equality condition. In order to do that i converted the equality constraint into 2 inequality constaints and passed everything cia constrOptim, as the manual said: everything included in the ... will be passed to Optim

[R] Z test and proportions

Hi All, I have a table based on ordial data and i want to compare proportions and i've seen in the pwr package i can use power.prop.test however i want to find out what the sig. value is based on n1,n2,p1,p2 and this package doesn't contain this.. Does anyone know of a package that does or is it

[R] Decision Trees RWeka

Hello, I have a question concerning decision trees coming from RWeka : library(RWeka) m =J48(Species~.,data=iris) How could such a decision tree be transferred into a matrix, pretty much in the same fashion, as it is done by getTree() in library(ofw) library(ofw) data(sr

[R] Help fCalendar holidayNYSE for Regressors

I am working with weekly time series data as in: tsData=ts(data,start=c(2004,1),freq=52) I have a table of regression variables that matchs called cReg (loaded from an xls sheet). I would like to append to the cReg table dummy variables for all the holidays as calculated from the fCalendar packa

[R] Fitting Multiple Univariate Distributions to Data

I am looking for procedure that allow one to fit multiple distributions to a variable. For example, based on analysis of the data we suppose that the data can be represented by 3-5 normal distributions added together. I would like to be able to determine the mean, sd, and weight associated with

[R] R error using Survr function with gcmrec

Would someone be able to help with this question? I'm using the Gcmrec, Survrec, and Design packages to do a power analysis on simulated data. I'm receiving an error after using the Survr function that all data must have a censoring time even after using the gcmrec function: newdata<-add

  1   2   >