Dear all,
I noticed the following in the call of lme using msVerbose.
fm1 - lme(distance ~ age, data = Orthodont, control =
lmeControl(msVerbose=T))
9 318.073: -0.567886 0.152479 1.98021
10 318.073: -0.567191 0.152472 1.98009
11 318.073: -0.567208 0.152473
And qplot(x=Categorie,y=Total,data=mydata,geom=bar,fill=Part) + coord_flip()
makes it a bit easier to read the labels.
Hadley
On Dec 5, 2007 8:33 AM, Domenico Vistocco [EMAIL PROTECTED] wrote:
This command works:
qplot(x=Categorie,y=Total,data=mydata,geom=bar,fill=Part)
for your data.
Hans W Borchers wrote:
Dear R-help,
I am looking for a function that will plot error bars in x- or y-direction
(or
both), the same as the Gnuplot function 'plot' can achieve with:
plot file.dat with xyerrorbars,...
Rsite-searching led me to the functions 'errbar' and 'plotCI'
8rino-Luca Pantani wrote:
Dear R-users.
I eventually bought myself a new computer with the following
characteristics:
Processor AMD ATHLON 64 DUAL CORE 4000+ (socket AM2)
Mother board ASR SK-AM2 2
Ram Corsair Value 1 GB DDR2 800 Mhz
Hard Disk WESTERN DIGITAL 160 GB SATA2 8MB
I'm a newcomer to
Dear R-users.
I eventually bought myself a new computer with the following
characteristics:
Processor AMD ATHLON 64 DUAL CORE 4000+ (socket AM2)
Mother board ASR SK-AM2 2
Ram Corsair Value 1 GB DDR2 800 Mhz
Hard Disk WESTERN DIGITAL 160 GB SATA2 8MB
I'm a newcomer to the Linux world.
I started
Note that Ottorino has only 1GB of RAM installed, which makes a 64-bit
version of R somewhat moot. See chapter 8 of
http://cran.r-project.org/doc/manuals/R-admin.html
I would install a i386 version of R on x86_64 Linux unless I had 2Gb or
more of RAM. I don't know how easily that works on
I'm coming late to this, but this *does* need a correction
just for the archives !
MS == Marc Schwartz [EMAIL PROTECTED]
on Sat, 01 Dec 2007 13:33:21 -0600 writes:
MS On Sat, 2007-12-01 at 18:40 +, David Winsemius wrote:
David Winsemius [EMAIL PROTECTED] wrote in
Hi everyone,
Has anyone written a parser in Java for either the ASCII or binary format
produced by save()? I need to parse a single large 2D array that is
structured like this:
list(
32609_1 = c(-9549.39231289146, -9574.07159324482, ... ),
32610_2 = c(-6369.12526971635, -6403.99620977124,
Prof Brian Ripley wrote:
Note that Ottorino has only 1GB of RAM installed, which makes a 64-bit
version of R somewhat moot. See chapter 8 of
http://cran.r-project.org/doc/manuals/R-admin.html
Only somewhat. The Opteron actually has 1GB too (Hey, it was bought in
2004! And the main point was
On Wed, Dec 05, 2007 at 06:11:40PM +0100, Peter Dalgaard wrote:
One oddity about Ubuntu is that there are no CRAN builds for 64bit.
Volunteers would be welcomed with open arms.
Dirk,
--
Three out of two people have difficulties with fractions.
__
Prof Brian Ripley wrote:
You mean the package nls2 at
http://w3.jouy.inra.fr/unites/miaj/public/AB/nls2/install.html
and not the unfortunately named nls2 that has just appeared on CRAN?
The first is not really a 'package' in the R sense.
Actually, it is one (sort of), but it is broken. The
I am relatively new to R and object oriented programming. I have relied on
SAS for most of my data analysis. I teach an introductory undergraduate
forecasting course using the Diebold text and I am considering using R in
addition to SAS and Eviews in the course. I work primarily with univariate
On Wed, 05-Dec-2007 at 06:11PM +0100, Peter Dalgaard wrote:
[]
| One oddity about Ubuntu is that there are no CRAN builds for 64bit.
| Presumably, the Debian packages work, or you can get the build
| script and make your own build. This is not really within my
| domain, though.
I've always
On Wed, 5 Dec 2007, David Coppit wrote:
Hi everyone,
Has anyone written a parser in Java for either the ASCII or binary format
produced by save()? I need to parse a single large 2D array that is
structured like this:
list(
32609_1 = c(-9549.39231289146, -9574.07159324482, ... ),
8rino-Luca Pantani wrote:
Dear R-users.
I eventually bought myself a new computer with the following
characteristics:
Processor AMD ATHLON 64 DUAL CORE 4000+ (socket AM2)
Mother board ASR SK-AM2 2
Ram Corsair Value 1 GB DDR2 800 Mhz
Hard Disk WESTERN DIGITAL 160 GB SATA2 8MB
I'm a
On Wed, 5 Dec 2007, Prof Brian Ripley wrote:
[...]
Thanks I'll read it more carefully.
Perhaps if you told us what you are trying to achieve we might be able to
help you achieve it.
I have a function which takes a date as an argument. I've tested it, and
I'd like to run it over a range
Hi I have a matrix with NA value that I would like to convert these to a value
of 0.
any suggestions
Kind Regards
Amit Patel
___
ttp://uk.promotions.yahoo.com/forgood/
[[alternative HTML version deleted]]
Jeff Delmerico wrote:
I'm working on a shared library of C functions for use with R, and I want
to create a matrix in R and pass it to the C routines. I know R computes
and supposedly stores numerics in double precision, but when I create a
matrix of random numerics using rnorm(), the
Peter Dalgaard wrote:
The toolchain availability tends to get in the way. Linux-based gadgets
could prove easier. I do wonder from time to time whether there really
is a market for R on cellphones...
As soon as someone writes library(ringtone) there might be :)
And I think you'd have to
On 05-Dec-07 18:57:58, Barry Rowlingson wrote:
Peter Dalgaard wrote:
The toolchain availability tends to get in the way. Linux-based
gadgets could prove easier. I do wonder from time to time whether
there really is a market for R on cellphones...
As soon as someone writes
I just read the description in ?Krig in the package fields which says:
Fits a surface to irregularly spaced data.
Yes, that is the most general case. Regular data location is a subset
of irregular. Anyway, kriging, just one g, after the name of Danie
Krige, the south african statistician who
Hi all,
I hope this question is not too trivial. I can't find an explanation
anywhere (Stats and R books, R-archives) so now I have to turn to the R-list.
Question:
If you have a factorial design with two factors (say A and B with two
levels each). What does the intercept coefficient with
I am not sure to understand your problem, but it seems to me that you
can use directly the function on the range of the dates:
x=as.Date(c('2007-01-01','2007-01-02'))
fff=function(x){y=x+1;return(y)}
fff(x)
[1] 2007-01-02 2007-01-03
class(fff(x))
[1] Date
Perhaps your function use a
hello,
I have this simple question. This is my dataset
size
1 57
2 97
3 105
4 123
5 136
6 153
7 173
8 180
9 193
10 202
11 213
12 219
13 224
14 224
15 248
16 367
17 496
18 568
19 618
20
On Wednesday 05 December 2007, Simone Gabbriellini wrote:
hello,
I have this simple question. This is my dataset
size
1 57
2 97
3 105
4 123
5 136
6 153
7 173
8 180
9 193
10202
11213
12219
13224
14224
15248
16
On Wed, 2007-12-05 at 18:42 +0100, Martin Maechler wrote:
I'm coming late to this, but this *does* need a correction
just for the archives !
MS == Marc Schwartz [EMAIL PROTECTED]
on Sat, 01 Dec 2007 13:33:21 -0600 writes:
MS On Sat, 2007-12-01 at 18:40 +, David Winsemius
You may want to look at package dlm.
Giovanni
Date: Wed, 05 Dec 2007 12:05:00 -0600
From: Alexander Moreno [EMAIL PROTECTED]
Sender: [EMAIL PROTECTED]
Precedence: list
Hi,
I'm trying to use the kalman filter to estimate the variable drift of a
random walk, given that I have a vector
Hi Ottorino,
I have been using R on 64-bit Ubuntu for about a year without
problems, both Intel and AMD CPUs. Installing and using several
packages (e1071, svmpath, survival) also works. However, I had to
install R from source:
$ gunzip -c R-2.6.1.tar.gz | tar xvf -
$ cd R-2.6.1
$ ./configure
Jeff Delmerico wrote:
Thanks Ben, that fixed the display within R. However, even after changing
the display settings, the matrix elements still appear to be exported in
single precision. The matrix object is being passed into my C routines as
an SEXP Numeric type, and somewhere along
x[is.na(x)] - 0
On 05/12/2007, Amit Patel [EMAIL PROTECTED] wrote:
Hi I have a matrix with NA value that I would like to convert these to a
value of 0.
any suggestions
Kind Regards
Amit Patel
___
Gustaf Granath wrote:
Hi all,
I hope this question is not too trivial. I can't find an explanation
anywhere (Stats and R books, R-archives) so now I have to turn to the R-list.
Question:
If you have a factorial design with two factors (say A and B with two
levels each). What does the
Hi --
I just noticed the following (R 2.6.1 on OSX)
lapply(c(as.Date('2007-01-01')), I)
[[1]]
[1] 13514
This is a bit surprising.. Why does lapply unclass the object? Sorry for
such a basic question, I don't seem able to produce the right google keywords.
Ranjan
You estimate a model with the Factors A or B either present (1) or not
present (0) and with an intercept. Thus you would predict:
For both A and B not present: Intercept
For A only present: Intercept+coef(A)
For B only preseent: Intercept+coef(B)
For both present: Intercept+coef(A)+coef(B).
I'm wanting to run R scripts non-interactively as part of a
technology independent framework.
I want control over the behaviour of these processes by specifying
various global variables in a configuration file that would be
passed as a command line argument.
I'm wondering if you know of any R
David Coppit wrote:
Hi everyone,
Has anyone written a parser in Java for either the ASCII or binary format
produced by save()?
You might want to consider using the hdf5 package to save the array in
HDF5 format. There are HDF5 libraries for Java as well
Hello,
I am running R 2.6.1 on windows xp
I am trying to fit a cox proportional hazard model with a shared
Gaussian frailty term using coxme
My model is specified as:
nofit1-coxme(Surv(Age,cen1new)~ Sex+bo2+bo3,random=~1|isl,data=mydat)
With x1-x3 being dummy variables, and isl being the
Hello,
I am running R 2.6.1 on windows xp
I am trying to fit a cox proportional hazard model with a shared
Gaussian frailty term using coxme My model is specified as:
nofit1-coxme(Surv(Age,cen1new)~ Sex+bo2+bo3,random=~1|isl,data=mydat)
With x1-x3 being dummy variables, and isl being the
On 2007-December-05 , at 16:47 , Scionforbai wrote:
I just read the description in ?Krig in the package fields which
says:
Fits a surface to irregularly spaced data.
Yes, that is the most general case. Regular data location is a subset
of irregular. Anyway, kriging, just one g, after
Dear R-Helpers
(1)
After a night's sleep, I realized why the other helpers think differently
from me.
I agree with others that it may be better to use multi-stratum model but I
was a bit surprised since they seem to think 'block' variable should *not*
be a fixed effect.
A. Others seemed to
On 12/5/07, Jeff Delmerico [EMAIL PROTECTED] wrote:
Does anyone know if there is a way to change the display or storage settings
so that the values will be displayed to their full precision? Or does
rnorm only produce values to single precision?
Any assistance would be greatly appreciated.
This command works:
qplot(x=Categorie,y=Total,data=mydata,geom=bar,fill=Part)
for your data.
domenico vistocco
Stéphane CRUVEILLER wrote:
Hi,
the same error message is displayed with geom=bar as parameter.
here is the output of dput:
dput(mydata)
structure(list(Categorie =
Dear R-Users,
I was wondering if there is a way to test the significance of the
clarkevans statistic in spatstat package?
I did not find any related function or the related values to calculate
it by hand.
does someone has any ideas?
thank you,
Jens
--
I'm getting unexpected results from the coxph function when using
weights from counter-matching. For example, the following code
produces a parameter estimate of -1.59 where I expect 0.63:
I agree with Thomas' answer wrt using offset instead of weight. One way to
understand this is to
On Wed, 5 Dec 2007, bernardo lagos alvarez wrote:
Hi,
Know anyone where to find the package rgdal for R.2.4.0?
On CRAN: the current version has
Package: rgdal
Title: Bindings for the Geospatial Data Abstraction Library
Version: 0.5-20
Date: 2007-11-07
Depends: R (= 2.3.0), methods, sp
Or
Thank you
As per advice from several R users I have set
r-project.org, stat.math.ethz.ch,fhcrc.org, stat.ethz.ch, math.ethz.ch,
hypatia.math.ethz.ch
all to be safe domains
But still some R emails go to Junk and require to be found manually
I have explored the issue with Univ Wash computing
Hi, Salut,
You should use the package nsl2 (only for Linux distribution)
Vous pouvez utiliser le package nls2 (Linux seulement)
Regards,
Souleymane Date: Tue, 4 Dec 2007 16:07:57 +0100 From: [EMAIL PROTECTED] To:
[EMAIL PROTECTED] CC: [EMAIL PROTECTED] Subject: Re: [R] confidence intervals
G'day Serge,
On Wed, 5 Dec 2007 11:25:41 +0100
de Gosson de Varennes Serge (4100)
[EMAIL PROTECTED] wrote:
I am using the quadprog package and its solve.QP routine to solve and
quadratic programming problem with inconsistent constraints, which
obviously doesn't work since the constraint
Hi, I am trying to get 95% CI s around a quantile growth curve, but
my result looks strange to me.
Here is the function I defined myself
boot.qregress-function(mat1, group, quantile, int, seed.1){
boot.fit-NULL
set.seed(seed.1)
for (i in 1:int){
but unfortunately it did not work.
What did not work? Provide some information more...
By the way, isn't 'gfortran' the new GNU fortran compiler which
replaced 'g77'? Or not on Ubuntu?
I immediately realize that an amd64 is not the right processor to make
life easy.
??? Example of bad
[if you get this twice: it seems to have not made it through, yesterday]
Earl == Earl F Glynn [EMAIL PROTECTED]
on Mon, 3 Dec 2007 13:26:11 -0600 writes:
Earl affy snp [EMAIL PROTECTED] wrote in message
Earl news:[EMAIL PROTECTED]
For example, it should go from very
On Wed, 5 Dec 2007, Ranjan Bagchi wrote:
Hi --
I just noticed the following (R 2.6.1 on OSX)
lapply(c(as.Date('2007-01-01')), I)
[[1]]
[1] 13514
This is a bit surprising.. Why does lapply unclass the object? Sorry for
such a basic question, I don't seem able to produce the right google
Hello everyone,
The thread title speaks for itself. Here's the code that worked for me:
numSummary(finance[,Employees], statistics=c(mean, sd, quantiles))
mean sd 0% 25% 50% 75% 100% n NA
11492.92 29373.14 1777 3040 4267 6553 179774 53 5
Hello,
I'm searching for a method which maps variables of this kind of table, see
below, to 2-dimensional space, like in multidimensional scaling. However, this
table is asymmetric: for example, variable T1 affects T2 more than T2 affects
T1(0.41 vs. 0.21).
DEPTABLE
T1T2 T3
On 2007-December-04 , at 21:38 , Scionforbai wrote:
- krigging in package fields, which also requires irregular spaced
data
That kriging requires irregularly spaced data sounds new to me ;) It
cannot be, you misread something (I feel free to say that even if I
never used that package).
Charles C. Berry wrote:
On Wed, 5 Dec 2007, Gad Abraham wrote:
Hi,
The following error looks like a bug to me but perhaps someone can shed
light on it:
library(splines)
library(survival)
s - survreg(Surv(futime, fustat) ~ ns(age, knots=c(50, 60)),
data=ovarian)
n -
Hi,
I need to compute a prediction R-squared for a linear regression. I
have figured out how to compute Allen's PRESS statistic (using the PRESS
function in the MPV library), but also want to compute the R-squared
that goes along with this statistic. I have read that this is computed
like an
I'm working on a shared library of C functions for use with R, and I want to
create a matrix in R and pass it to the C routines. I know R computes and
supposedly stores numerics in double precision, but when I create a matrix
of random numerics using rnorm(), the values are displayed in single
Michael Bibo wrote:
I am running R on a corporate Windows XP SP2 machine on which I do not
have
administrator privileges or access to most settings in Control Panel.
R is
installed from my limited user account. The version of the JVM I have
installed is perhaps best described as
Hi,
I am trying to figure out how to get the correlation coefficient for a
QQ plot (residual plot). So to be more precise, I am creating the plot
like this:
qq.plot(rstudent(regrname), main = rformula, col=1)
But want to also access (or compute) the correlation coefficient for
that plot.
On 12/5/07, Romain Francois [EMAIL PROTECTED] wrote:
Hello,
My guess is that you are actually talking about the numSummary function in
Rcmdr, not in abind. In that case, you can look how the structure of the
output is like:
str( numSummary( iris[,1:4] ) )
List of 4
$ type :
Henrique Dallazuanna wwwhsd at gmail.com writes:
x[is.na(x)] - 0
On 05/12/2007, Amit Patel amitpatel_ak at yahoo.co.uk wrote:
Hi I have a matrix with NA value that I would like to convert these to a
value of 0.
any suggestions
also
library(gdata)
x - matrix(rnorm(16), nrow=4, ncol=4)
Hi Thanks for your suggestion, I'm trying to install this package in Ubuntu
(7.10) but unsuccessfully. Also tried in MacOSX, and no success too.
_
De: Ndoye Souleymane [mailto:[EMAIL PROTECTED]
Enviado el: miércoles, 05 de diciembre de 2007 13:38
Para: [EMAIL PROTECTED];
Peter Dalgaard p.dalgaard at biostat.ku.dk writes:
Has something changed from 2.6.0 to 2.6.1 that may require JVM 1.4.1?
If
so, I can use that information to request an upgrade of my JVM.
Hmm, could be. They got rebuilt on my system and committted at some
point in the
See the following (the email seems to have wrapped 2 lines onto one
at one point but it should be obvious):
http://tolstoy.newcastle.edu.au/R/e2/help/07/06/18853.html
On Dec 5, 2007 5:05 PM, Thomas Allen [EMAIL PROTECTED] wrote:
I'm wanting to run R scripts non-interactively as part of a
Hello,
how is, for example, the Schwarz criterion is defined for kmeans? It should
be something like:
k - 2
vars - 4
nobs - 100
dat - rbind(matrix(rnorm(nobs, sd = 0.3), ncol = vars),
matrix(rnorm(nobs, mean = 1, sd = 0.3), ncol = vars))
colnames(dat) - paste(var,1:4)
(cl -
Dear Jens,
I've submitted a new version (0.7-4) of the polycor package to CRAN. The
hetcor() function now uses your nearcor() in sfsmisc to make the returned
correlation matrix positive-definite if it is not already.
I know that quite some time has elapsed since you raised this issue, and I
Dear R-help,
I am looking for a function that will plot error bars in x- or y-direction (or
both), the same as the Gnuplot function 'plot' can achieve with:
plot file.dat with xyerrorbars,...
Rsite-searching led me to the functions 'errbar' and 'plotCI' in the Hmisc,
gregmisc, and plotrix
I don't think this is referring to the nls2 package on CRAN but
rather something else.
On Dec 5, 2007 7:37 AM, Ndoye Souleymane [EMAIL PROTECTED] wrote:
Hi, Salut,
You should use the package nsl2 (only for Linux distribution)
Vous pouvez utiliser le package nls2 (Linux seulement)
Regards,
Dear friends :
using the glm function and setting family=binomial, I got a list of
coefficients.
The coefficients reflect the effects of predicted variables on the
probability of the response to be 1.
My response variable consists of A and D . I don't know which level of
the response was
On Wed, 2007-12-05 at 18:06 -0800, Bin Yue wrote:
Dear friends :
using the glm function and setting family=binomial, I got a list of
coefficients.
The coefficients reflect the effects of predicted variables on the
probability of the response to be 1.
My response variable consists of A
Dear Marc Schwartz:
When I ask R2.6.0 for windows, the information it gives does not contain
much about family=binomial .
You said that there is a detail section of ?glm. I want to read it
thoroughly. Could you tell me where and how I can find the detail section
of ?glm.
Thank you very much
anscombe is built into R already so you don't need to read it in.
An intercept is the default in lm so you don't have to specify it.
opar - par(mfrow = c(2,2))
plot(y1 ~ x1, anscombe)
reg - lm(y1 ~ x1, anscombe)
reg
abline(reg)
...etc...
par(opar)
Note that plot(anscombe[1:2]) and
Is the rmpi package (or rpvm) needed to exploit multiple cores on a
single unix box using the snow package. The documentation of the package
does not provide info about setting up a single machine with multiple
cores. Also, if how effective is it to run a bayesian simulation on
parallel (or
Hello all,
I have 3 time series (tt) that I've fitted segmented regression models
to, with 3 breakpoints that are common to all, using code below
(requires segmented package). However I wish to specifiy a zero
coefficient, a priori, for the last segment of the KW series (green)
only. Is this
Check the Rowv, Colv options to heatmap.2
data(mtcars)
x - as.matrix(mtcars)
heatmap.2(x, Rowv=FALSE, dendrogram=column)
-Ashoka
Scientist - Pacific Northwest National Lab
On Dec 5, 2007 4:20 PM, affy snp [EMAIL PROTECTED] wrote:
Dear list,
I am using heatmap.2(x) to draw a heatmap.
Dear all:
By comparing glmresult$y and model.response(model.frame(glmresult)), I
have found out which one is
set to be TRUE and which FALSE.But it seems that to fit a logistic
regression , logit (or logistic) transformation has to be done before
regression.
Does anybody know how to
On Thu, 6 Dec 2007, Emmanuel Charpentier wrote:
Prof Brian Ripley a écrit :
Note that Ottorino has only 1GB of RAM installed, which makes a 64-bit
version of R somewhat moot. See chapter 8 of
http://cran.r-project.org/doc/manuals/R-admin.html
Thank you for this reminder|tip ! I didn't
Michael Kubovy kubovy at virginia.edu writes:
Dear r-helpers,
How do I tell xYplot() and Key() that I want the labels in italic?
Key(x = 0.667, y = 0.833, other = list(title = expression(italic(v)),
cex.title = 1,
labels = c(expression(italic(b)), expression(italic(c)),
78 matches
Mail list logo