[R] Rmetrics New Built

2004-05-31 Thread Diethelm Wuertz
*www.Rmetrics.org Rmetrics - new Built 190.10053 *The new built has now implemented my 'timeDate' and 'timeSeries' classes which became part of the fBasics package. Furthermore, MS Windows specifics were removed from the packages, so we can try to build Rmetrics on Linux and on Mac OSX.

Re: [R] Getting the same values of adjusted mean and standard errors as SAS

2004-05-31 Thread David J. Netherway
Thanks for the help. Both the Design package and the effects package look as though they are what I need although it will probably take me a while to get on top of both. I have had a brief go at the Design package and the contrast function is particularly useful. A question on the Design

Re: [R] glmm?

2004-05-31 Thread Douglas Bates
Spencer Graves [EMAIL PROTECTED] writes: Is there an easy way to get confidence intervals from glmm in Jim Lindsey's library(repeated)? Consider the following slight modification of an example from the help page: df - data.frame(r=rbinom(10,10,0.5), n=rep(10,10),

Re: [R] bar plot patterns

2004-05-31 Thread Osman
Thank you for your answer. I have about 7 stacks it is not very appealing to have just stripes with changing angles. I was wondering if there is a way to have varying patters in black and white. Osman - Original Message - From: Petr Pikal To: Osman Cc: [EMAIL PROTECTED]

Re: [R] Rmetrics New Built

2004-05-31 Thread elijah wright
What are the next steps? Maybe somebody is around and can try to build the packages for (Debian) Linux and Mac OSX. This would be a great help for me! The *.tar.gz files are availalble on If the Linux and Mac OSX builds are successfully done, I will submit the packages to the CRAN server.

[R] Calculating distances between points in a data frame?

2004-05-31 Thread Sander Oom
Dear list, I would like to calculate the distance between consecutive points in a data frame. Of course the first point in the data frame does not have a point of origin, and should get a value NA. I have tried two different loops, which both result in error: num - seq(0,10,1) X -

[R] [OT] plot y against x

2004-05-31 Thread Ted Harding
Hi Folks, I'd be grateful for some views on the following. When I say plot y against x I mean that y is on the vertical axis and x is on the horizontal axis. I acquired this usage so long ago that I can no longer remember how I acquired it, and therefore can not cite my authority for my usage.

Re: [R] Calculating distances between points in a data frame?

2004-05-31 Thread Wolski
Hi! It may bee that the function dist can be of some use to you? ?diff I have something like this in mind. (you do not need a loop.) XY$DistXY - sqrt(diff(X)^2+diff(Y)^2) Have fun trying. Sincerely Eryk *** REPLY SEPARATOR *** On 5/31/2004 at 5:25 PM Sander Oom wrote: Dear

RE: [R] [OT] plot y against x

2004-05-31 Thread Charles Annis, P.E.
Engineers have been plotting stress (or strain, or log of either) as y against log(fatigue lifetime) as x for 100 years or more. Thus it appears that allowable stress is a function of required cycles, and that is how the plot is interpreted. The underlying regression, however, correctly treats

[R] Question about building library and BLAS

2004-05-31 Thread Zhu Wang
Dear helpers, I am trying to create a library which uses some Fortran source files and Lapack and Blas subroutines. The Fortran source files from the original author contain subroutines isamax.f, sgefa.f and sgesl.f, which are part of BLAS subroutines on my Linux computer, but maybe different

Re: [R] Rmetrics New Built

2004-05-31 Thread Dirk Eddelbuettel
On Mon, May 31, 2004 at 09:54:21AM -0500, elijah wright wrote: If the Linux and Mac OSX builds are successfully done, I will submit the packages to the CRAN server. all four of the packages successfully build on Debian unstable - no errors, nor warnings. this is a Good Thing. I would

Re: [R] Question about building library and BLAS

2004-05-31 Thread Douglas Bates
Zhu Wang [EMAIL PROTECTED] writes: I am trying to create a library which uses some Fortran source files Someone named Martin Maechler will shortly be sending you email regarding the distinction between 'library' and 'package' :-) (You are creating a package, not a library, despite the

Re: [R] Question about building library and BLAS

2004-05-31 Thread Douglas Bates
Douglas Bates [EMAIL PROTECTED] writes: Zhu Wang [EMAIL PROTECTED] writes: I am trying to create a library which uses some Fortran source files Someone named Martin Maechler will shortly be sending you email regarding the distinction between 'library' and 'package' :-) (You are

Re: [R] Calculating distances between points in a data frame?

2004-05-31 Thread Gabor Grothendieck
Try using running from the gregmisc package with pad = TRUE: require(gregmisc) XY - data.frame(num = seq(0,10), X = seq(0,30,3), Y = seq(0, 40, 4) ) DistXY - function(idx) { i - idx[2] with(XY, sqrt( (X[i]-X[i-1])^2 + (Y[i]-Y[i-1])^2 ) ) } XY$Dist - running( 1:nrow(XY), width=2, fun =

[R] Putting referenced titles on plots

2004-05-31 Thread grr
I have a data frame ctx and an array names, where names[i] is the column name for ctx[i], and am making histograms for each column of ctx: for (i in 2:ncol(ctx)){hist(ctx[,i], br=100, main=names[i])} The titles don't come out like I expect. Each names[i] is something like 1098_s_at and R

[R] Contrasts

2004-05-31 Thread Kimberly Ann Fernandes
Hello, I am trying to figure out how to conduct a t-test on a specific contrast for my data. I have four factors in my data and would like to conduct a t-test on the average of the data from the first two factors against the average of the data on the second two factor (i.e. is the average of

Re: [R] Regression model type II

2004-05-31 Thread Angel Lopez
Hi Benjamin, Maybe this link is useful to you: http://eeb37.biosci.arizona.edu/~brian/splus.html It has a function 'slope' that calculates different type II regressions and a link to a paper comparing them. Although it was written for S-plus it works in R too. If you get any better solutions, let

Re: [R] Putting referenced titles on plots

2004-05-31 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: I have a data frame ctx and an array names, where names[i] is the column name for ctx[i], and am making histograms for each column of ctx: for (i in 2:ncol(ctx)){hist(ctx[,i], br=100, main=names[i])} Works for me. Are you sure names is a vector containing the stuff you

Re: [R] Putting referenced titles on plots

2004-05-31 Thread Spencer Graves
1. That seems strange. The following simplification of your function produced sensible results for me: hist(1:4, main= 1098_s_at) I got similar results from your exact statement after first defining names and ctx as follows: names - letters[1:3] ctx -

Re: [R] Question about building library and BLAS

2004-05-31 Thread Zhu Wang
On Mon, 2004-05-31 at 11:39, Douglas Bates wrote: I should have read your message more carefully. Those three Fortran routines that you mention are single precision and it would be silly to use them on modern computers. R only provides single precision floating point for compatibility. All

Re: [R] Calculating distances between points in a data frame?

2004-05-31 Thread Sander Oom
Hi Gabor, Thanks for your suggestion. However when installing the package gregmisc, I get the following error: local({a - CRAN.packages() + install.packages(select.list(a[,1],,TRUE), .libPaths()[1], available=a)}) trying URL `http://cran.r-project.org/bin/windows/contrib/1.9/PACKAGES' Content

Re: [R] Contrasts

2004-05-31 Thread Spencer Graves
I hope you get a reply from someone who knows more about this than I. However, in the spirit that a quick hack is sometimes better than an elegant answer later, consider the following: DF - data.frame(a=rep(letters[1:4], 2), y=1:8) DF$b - ((DF$a %in% letters[1:2])-(DF$a%in%

RE: [R] [OT] plot y against x

2004-05-31 Thread Bill Vinyard
It would not matter if the relationship you are trying to plot is monotonic and therefore invertible. If the relationship is non-monotonic and therefore not invertible, then it does matter which variable you call dependent, since in this latter case you have a multi-valued relationship.

Re: [R] Calculating distances between points in a data frame?

2004-05-31 Thread Uwe Ligges
Sander Oom wrote: Hi Gabor, Thanks for your suggestion. However when installing the package gregmisc, I get the following error: local({a - CRAN.packages() + install.packages(select.list(a[,1],,TRUE), .libPaths()[1], available=a)}) trying URL

Re: [R] Rmetrics New Built

2004-05-31 Thread elijah wright
If the Linux and Mac OSX builds are successfully done, I will submit the packages to the CRAN server. all four of the packages successfully build on Debian unstable - no errors, nor warnings. this is a Good Thing. I would guesstimate that Ah, thanks, good to know. Did you try 'R

Re: [R] Putting referenced titles on plots

2004-05-31 Thread grr
Thank you both for the help. The below 1. works for me, too, both cases. I am not actually using names as an object name. I am writing: for (i in 2:nrow(ctxheadlogtrans)){hist(ctxheadlogtrans[,i], br=100, main=gene.names.head[i])} which gives me histograms titled 1, 2, 3, etc. But:

Re: [R] Putting referenced titles on plots

2004-05-31 Thread Roger D. Peng
Is `names' by any chance a factor? What is class(names)? If `names' is a factor, then you're getting the underlying numeric representation rather than the factor levels. -roger [EMAIL PROTECTED] wrote: I have a data frame ctx and an array names, where names[i] is the column name for ctx[i],

Re: [R] Putting referenced titles on plots

2004-05-31 Thread grr
Yes! Thank you. It is a factor (I don't know how it became a factor). I created an array based on the 'names' factor, and can get the titles by referencing the new array. Thanks so much, Graham On May 31, 2004, at 11:44 AM, Roger D. Peng wrote: Is `names' by any chance a factor? What is

RE: [R] Converting data frame to array?

2004-05-31 Thread TAPO (Thomas Agersten Poulsen)
Dear John, Thank you for your helpful answer. I was obviously being stupid, as I have, as you point out, more predictors than observations. What I was hoping to get was some sort of an explaining linear combination of my predictors: which predictors are important for the results

[R] ffnet problem

2004-05-31 Thread v . demartino2
Context:Linux debian testing, compiled R 1.9.0 from source. I've just installed the contributed ffnet package wit no problem at all. But when loading the library the following error message is popping up and no ffnet command seems to work: library(ffnet) Error in dyn.load(x, as.logical(local),

Re: [R] ffnet problem

2004-05-31 Thread Dirk Eddelbuettel
On Mon, May 31, 2004 at 09:38:54PM +0200, [EMAIL PROTECTED] wrote: Context:Linux debian testing, compiled R 1.9.0 from source. I've just installed the contributed ffnet package wit no problem at all. No problem at all seems unlikely in light of ... But when loading the library the following

[R] Doubts on anova and use of contrasts in multcomp package

2004-05-31 Thread José Cláudio Faria
Dear list, I have been studying R and I would like the aid of more experienced to solve the problems of the analysis below: r = gl(3, 8, label = c('r1', 'r2', 'r3')) e = rep(gl(2, 4, label = c('e1', 'e2')), 3) y = c(26.2, 26.0, 25.0, 25.4, 24.8, 24.6, 26.7, 25.2, 25.7, 26.3, 25.1, 26.4,

[R] Help on parameters

2004-05-31 Thread Ronaldo Reis Jr.
Hi, I have a follow analysis Trat1 = quantitative variable Trat2 = qualitative variable with 3 levels (A, B, C) Trat3 = qualitative variable with 3 levels (D, E, F) Resp = Response I try to get the parameters to compare with zero, so I make this model: glm(Resp~Trat1*Trat2+Trat1*Trat3-Trat1-1)

Re: [R] Converting data frame to array?

2004-05-31 Thread Spencer Graves
stepAIC in library(MASS) or step? hope this helps. spencer graves TAPO (Thomas Agersten Poulsen) wrote: Dear John, Thank you for your helpful answer. I was obviously being stupid, as I have, as you point out, more predictors than observations. What I was hoping to get was

RE: [R] Converting data frame to array?

2004-05-31 Thread John Fox
Dear Thomas, I doubt whether there's anything useful that you can do with 300 predictors and only 10 observations. A naïve application of variable selection will likely allow you to account perfectly for the variation in the response variable just by capitalizing on chance. John -Original

Re: [R] ffnet problem

2004-05-31 Thread v . demartino2
I've just installed the contributed ffnet package wit no problem at all. No problem at all seems unlikely in light of ... Dirk, You're right! Here you are what happened: desktop:/tmp# R INSTALL -l /usr/local/lib/R/library ffnet/ * Installing *source* package 'ffnet' ... ** libs make:

Re: [R] Rmetrics New Built

2004-05-31 Thread Diethelm Wuertz
Hello Elijah, Hello Dirk First of all many thanks for your help and support. 1) Sorry, I used non-standard kyewords and thus I added to my doc/KEYWORDS.db database the following lines: Rmetrics: Rmetrics Rmetrics|winRmetrics: winRmetrics Rmetrics|fBasics: Basics Rmetrics|fSeries: Series

Re: [R] ffnet problem

2004-05-31 Thread Dirk Eddelbuettel
On Mon, May 31, 2004 at 10:38:37PM +0200, [EMAIL PROTECTED] wrote: I've just installed the contributed ffnet package wit no problem at all. No problem at all seems unlikely in light of ... Dirk, You're right! Here you are what happened: desktop:/tmp# R INSTALL -l

RE: [R] ffnet problem

2004-05-31 Thread Liaw, Andy
Given that those files are dated 27/06/2000, they are unlikely to work with the current version of R. (R has changed quite a bit since 2000!!) Don't know why/how the code depends on NR codes, as one can find better quality code in many instances. It might not be too difficult to get rid of the

[R] LAPACK and ScaLAPACK new functionality survey

2004-05-31 Thread Liaw, Andy
Apologies for the off-topic post, and apologies for those already seen this elsewhere... From NA-digest last week: --- From: Jack Dongarra [EMAIL PROTECTED] Date: Fri, 28 May 2004 16:17:17 -0400 Subject: LAPACK and ScaLAPACK New Functionality

[R] Tree() and confidence intervals

2004-05-31 Thread Briggs, Meredith M
Hello I'm currently using Monte Carlo techniques to estimate prices (variable not static) from the following type of data: 1,22,40,22,33,5,2000 3,45,33,6,7,0,3000 22,22,33,44,55,66,7 Each row is a record from group A and the cells in all

[R] rcmd with a libgsl.a

2004-05-31 Thread lamack lamack
Dear all, I'm new to R and am trying to create a dll in order to be able to use the dyn.load. I work with some examples and its works fine. Now, I would like to create a dll that use a lib (libgsl.a). I have linked its in a main c program and its works fine too. How can I instruct rcm ... to

Re: [R] ffnet problem

2004-05-31 Thread Jason Turner
Anyway, what steps should I take now? should isn't quite what I'm telling you; just free advice. :) I use nnet from the nnet package (VR bundle), and find it very good. And it doesn't require any additional libraries. If you've got a binary installation, you've probably got it already...