*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. Please
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 pac
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
David J. Netherway wrote:
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.
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(1
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]
S
> 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 s
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 <- seq(0,3
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
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 l
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
ob
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 (o
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 wou
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
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' :-)
>
>
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 =
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
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 the
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
[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
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 <- dat
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.
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
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%
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.
-Orig
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 `http://cran.r-project.org/bin/windows/contrib/1.9
> > > 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
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:
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
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 class(n
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
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),
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 foll
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,
19
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)
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 so
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
>>
>> 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
mak
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
Rmetri
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 INSTA
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 d
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 Functionalit
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 b
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
> 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...
45 matches
Mail list logo