HI,
I wonder does anyone have experience with doing sequential gaussian
simulation with krige() function in gstat?
I find it VERY slow compared to use krige() to achieve kriging function
itself.. I wonder why, is that because it has to model the variogram, and
do the kriging separately for each
Thomas Lumley <[EMAIL PROTECTED]> writes:
> While the rules for constructing the tree are part of the language
> definition, the order of evaluation might not be. In C, for example, the
> order of evaluation is not specified except to the extent that precedence
> constrains it (and for a few spec
> On Thu, 26 Jun 2003 14:43:41 +0100,
> Xu Yun (XY) wrote:
> Dear All:
> I'm a newbie to R and chemometrics.
> Now I'm trying apply bclust on fuzzy c-means like this:
>> bc1 <- bclust(iris[,1:4], 3, base.centers=20,iter.base=100,
> base.method="cmeans")
> Committee Member:
>
Dear R-helpers,
I'd like to be able to post-process contours coming from contour(). Does
anyone have a version of contour() (or something similar) which will return
the contour coordinates?
In searching the archives I've come across a message in Nov 01 which had
this on a wish-list, but I ca
Dear all,
i have a question about the "sample" function used in R, does it work as a
pseudo-dandom number generator programmed with C, like it is described in Modern
Applied Statics with S-Plus 3d edition chapter 5 section 2?
Thanks a lot
-
[[alte
Hi,
I would like to plot a graph on the jpeg device for example and to write
a table associated to this plot where I have the pixel coordonates of
each plotted point so that I could include the jpeg image onto a html
page and get all the information about each point when the mouse button
is pr
R is not S-PLUS, and you need Modern Applied Statistics in S (4th ed) for a
description including R.
sample in R used a PRNG: see ?RNG in R for the details of PRNGs in R.
On Fri, 27 Jun 2003, [iso-8859-1] Ramzi Feghali wrote:
> i have a question about the "sample" function used in R, does it wor
Philippe Hupé wrote:
I would like to plot a graph on the jpeg device for example and to write
a table associated to this plot where I have the pixel coordonates of
each plotted point so that I could include the jpeg image onto a html
page and get all the information about each point when the mo
Hello to everybody:
Im triying to develop a GAM model using Grasper. However, I cant obtain
any result.
Im obtaining always the same result:
Error in "[.data.frame"(thindata, , cols) :
invalid subscript type
when I use grasp.cormat
or
Error in "[.data.frame"(XXX, , c(1, sX)) :
I have been having similar problems trying to compile R 1.7.1 in OS X
10.2.6 with no success. If you really want to use R just use the Fink
package installer and install the "r-base" package ( R 1.7.0). You
might like to get in touch with the package maintainer (check package
info at the Fink
Hello:
I am running MCD robust method on R and Splus 2000, I have two
questions:
1. Why is it slower the cov.mcd R version than the S-Plus version? (it
takes quite a long time more the R version)
2. The R version shows only the "cov", but the S-Plus version shows both
the "cov" and "raw.cov" cova
On 6/27/03 8:03, Hector L. Ayala-del-Rio wrote
>I have been having similar problems trying to compile R 1.7.1 in OS X
>10.2.6 with no success. If you really want to use R just use the Fink
>package installer and install the "r-base" package ( R 1.7.0). You
>might like to get in touch with the
I'm not even close to expert at solving these kinds of problems, so
I'll just offer some information from my successful build of R 1.7.1
in OS X 10.2.6 (on a G4 desktop machine), and hope it helps.
I used this configure command:
../source/configure --prefix=/Users/macq/R/R-1.7.1 --enable-R-shli
Hello All,
Is there a package that contains the PDF and CDF of the bivariate
lognormal distribution?
Thanks,
John.
--
--
Dr. John Janmaat
Department of Economics, Acadia University, Wolfville, NS, B4P 2R6
E-mail: [EMAIL PRO
Are you familiar with :
@article{box,:hunt:1962,
Author = {Box, G. E. P. and Hunter, William G.},
Title = {A Useful Method for Model-building},
Year = 1962,
Journal = {Technometrics},
Volume = 4,
Pages = {301--318},
Keywords = {Adaptive modeling; Empirical model; Kinetic
Dear All,
is there any problem using the PPC version of Mandrake to run R?
Cheers,
Federico Calboli
=
Federico C.F. Calboli
Department of Biology
University College London
Room 327
Darwin Building
Gower Street
London
WClE 6BT
Tel: (+44) 020 7679 4395
Fax (+44) 020 76
Webcast - Integration with R: Extending Statistical
Analysis in DecisionSite
July 24, 2003 11:00 am EDT
Duration: 1 hour
For more information and to register, go to:
http://www.spotfire.com/webcasts
Statisticians and statistical programmers often
face the challenge of how to deploy the metho
Hi,
Is there an easy way to read matlab .mat file from R? There are R
packages (R.matlab, R.io from google) around but I might just missed
some easy functions.
Many thanks in advance.
--
Bo Peng
__
[EMAIL PROTECTED] mailing list
https://www.stat.m
Gurus,
I used
predict(loess(Y~X));
where Y and X are of the same length. But there are same NA's in both Y and
X. Those NA's are in the same locations in Y and X. The following is the
error messageI got:
Error in "[<-"(*tmp*, , i, value = predict(loess(Y~X))) :
number of items t
Hi
Is there is there an equivalent of S-plus "feffects" argument in
model.tables() in R for calculating effects of a factorial design? Or is
there another way they can be calculated from an aov fit.
Thanks
Peter
__
[EMAIL PROTECTED] mailing list
ht
On Fri, 27 Jun 2003, Yao, Minghua wrote:
> Gurus,
>
> I used
> predict(loess(Y~X));
>
> where Y and X are of the same length. But there are same NA's in both Y and
> X. Those NA's are in the same locations in Y and X. The following is the
> error messageI got:
>
> Error in "[<-"(*tmp*, ,
On Fri, 27 Jun 2003, Peter Ho wrote:
> Is there is there an equivalent of S-plus "feffects" argument in
> model.tables() in R for calculating effects of a factorial design?
No (as reading the help page would have told you).
> Or is
> there another way they can be calculated from an aov fit.
Hi evryone,
I would like to access to DB2 on R on a linux machine (R and DB2 are on 2 pc
linux )
I have read on the R-archive that RODBC work with DB2, but it is not clear
for me, if it's only for R on a windows machine or if it's work also for R
under linux ?
My data base administrator say that
On Fri, 27 Jun 2003, laurent buffat wrote:
> I would like to access to DB2 on R on a linux machine (R and DB2 are on 2 pc
> linux )
> I have read on the R-archive that RODBC work with DB2, but it is not clear
> for me, if it's only for R on a windows machine or if it's work also for R
> under lin
Hello-
I'm trying to run R from the command line using a function as an input; it
seems to run ok when I use "Rterm --vanilla < name_of_R_function.R" (by ok,
I mean that it will say that it's reading in the correct items) but the
graphics won't display on screen. Do I need to explicitly define
I have data on a number of patients. Essentially, for each patient I know
his/her age and whether he/she exhibits certain symptoms:
age symptom1 symptom2
50 0 1
53 0 0
70 1 1
...
I have started off by fitting simple models with forms like
>-Original Message-
>From: [EMAIL PROTECTED]
>[mailto:[EMAIL PROTECTED] On Behalf Of Prof
>Brian Ripley
>Sent: Friday, June 27, 2003 1:35 AM
>To: Marc Schwartz
>Cc: [EMAIL PROTECTED]; 'Nick Bond'
>Subject: RE: [R] dropping factor levels in subset
>
>
>A more transparent solution is
>
>old
On 27 Jun 2003 at 0:38, Peter Ho wrote:
> Hi all,
>
> I have been trying to reproduce an analysis from Douglas Montgomery´s
> book on design and analysis of experiments.
If you are learning design of experiments, there are better books
around. Montgomery's sems very populat, though. Anybody k
Re: [, drop=TRUE} for factors
It's been in S-PLUS (but not S I believe) for a long time, probably since
before 1994: it is in S+3.4, 1996 vintage.
It appears to have been added to R around August 1998.
Yes, Frank Harrell argues for the default to be true and I believe his
Hmisc package overrides
On Fri, 27 Jun 2003 11:57:02 -0500
Marc Schwartz <[EMAIL PROTECTED]> wrote:
> >-Original Message-
> >From: [EMAIL PROTECTED]
> >[mailto:[EMAIL PROTECTED] On Behalf Of Prof
> >Brian Ripley
> >Sent: Friday, June 27, 2003 1:35 AM
> >To: Marc Schwartz
> >Cc: [EMAIL PROTECTED]; 'Nick Bond'
>
On Fri, 27 Jun 2003, Nick Helseth wrote:
So, you are using Windows: please say so up front.
> I'm trying to run R from the command line using a function as an input; it
> seems to run ok when I use "Rterm --vanilla < name_of_R_function.R" (by ok,
> I mean that it will say that it's reading in the
hello,
I only want to get the slope of a linear regression of ca. 100 variables
against time.
I can do for each response (100 times)
var1.lm <- lm(response~predictor)
but I thought that there might be an easier way of doing this. If I am
including more variables it is doing a multiple regres
>-Original Message-
>From: [EMAIL PROTECTED]
>[mailto:[EMAIL PROTECTED] On Behalf Of Prof
>Brian Ripley
>Sent: Friday, June 27, 2003 12:11 PM
>To: Marc Schwartz
>Cc: [EMAIL PROTECTED]; 'Nick Bond'
>Subject: RE: [R] dropping factor levels in subset
>
>
>Re: [, drop=TRUE} for factors
>
>It'
Consider the following:
> tstdat <- data.frame(x=1:8, y1=rep(1:4, each=2), y2=0.01*rnorm(8))
>
> ys <- names(tstdat)[2:3]
>
> slopes <- rep(NA,2)
> names(slopes) <- ys
>
> for(i in 1:2){
+ mdl <- formula(paste(ys[i], "~ x"))
+ slopes[i] <- coef(lm(mdl, tstdat))[2]
+ }
> slopes
y1
Martin -
My recollection is that if the left hand side in a model formula
is a matrix (in your case an [n x 100] matrix) then either lm()
or glm() will return a matrix of coefficients. These are the point
estimates for a multivariate regression (meaning, multivariate
response). I hunted just a
Great suggestions, Thomas. Consider the following:
> y <- as.matrix(data.frame(y1=rep(1:4, each=2), y2=0.01*rnorm(8)))
> x <- 1:8
> coef(lm(y~x))
y1 y2
(Intercept) 0.3571429 0.011181184
x 0.4761905 -0.001528493
hth. spencer graves
Thomas W Blackwell wrote:
Thank you, Prof. Ripley. By adding na.action=na.exclude to loess function,
the program runs great.
I did omit some trivial things. But I still think it was NA's that caused
troubles in predict or loess.
Minghua
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent:
On Fri, 27 Jun 2003, Yao, Minghua wrote:
> Thank you, Prof. Ripley. By adding na.action=na.exclude to loess function,
> the program runs great.
>
> I did omit some trivial things. But I still think it was NA's that caused
> troubles in predict or loess.
*NO* the errors message is from your usag
try the following:
first cbind all the design vectors,
data<-cbind()
response<-x
apply(data, 2, function(x) lm(response~x)$coef[2])
It should give you all the slope coef (not multivariate but with an
intercept try lm(response~x-1)$coef if you do not want the intercept).
- Original Messag
During the installation I do my './configure' but the 'make'
command does not work. I'm using the Sun native make, full
path /usr/xpg4/bin/make and we have a gcc compiler. I could
really use some help trying to understand what the make com-
mand is doing. (or not doing in my case)
Hi:
This is not a typical R posting, but I was quite surprised to read
Prof. Ripley's comment about the inappropriate use of AIC to
compare "non-nested" models. As he says, While it is indeed true that
Akaike's (1973) develops AIC for nested models, i.e. models which can
be obtained by variou
I believe the recommended `make' command is GNU make. You may want to
install/try this and see if it works. It may already be installed on
your system but called `gmake'.
-roger
George Caunt wrote:
During the installation I do my './configure' but the 'make'
command does not work. I'm using
As the R-admin manual says, you can use either the GNU make or the Sun
native make which is in /usr/ccs/bin.
The phrase `does not work' is dreaded by all technical support personnel:
it is maximally uniformative. So please try again to explain the actual
problem!
On Fri, 27 Jun 2003, Roger D.
Please try to read the whole paragraph. When you have done that, read my
exposition in my PRNN book.
Hint 1: AIC is about maximized likelihood, and one of the models being
compared is not being fitted by ML.
Hint 2: Differences in AIC for non-nested models are subject to large
sampling fluctua
Please, I am having problems in the process of
installing the R-packages in the R_HOME/library.
The following message is the error that I am getting:
# R CMD INSTALL -l /usr/local/bin/R/library base
WARNING: package 'base' does not exist
ERROR: no packages specified
# R CMD INSTALL -l /usr/local/b
Hi,
Newbie here, just joined the list.
I discovered something interesting when performing a mixed between and
within (repeated measure) ANOVA in R. The error terms (SS actually I
guess) are inflated by the number of factors of the within measure.
This just seems wrong. I think I'm doing it
Please forgive my ignorance on grapics. I'm trying to make a relatively
simple plot with two line plots, same axes, by mean over a series of
dates. I can make the plot well like this:
plot(sort(tapply(first.anti.auth.sum,date,mean), partial=1), type="l",
col="yellow",ylim=c(0,2.0))
par(new=TRUE)
p
How does one get the character that corresponds to a number? Also
the inverse?
The %c format code is not supported by the R sprintf function and I
could not find another way to do it.
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mai
How does one get a list of all the color names supported by
R and the mapping between color names and hex codes? Color
names are things like "red", "salmon", etc.
I am using Windows 2000 Professional and R 1.7.1.
__
[EMAIL PROTECTED] mailing list
http
On Fri, 27 Jun 2003, Gabor Grothendieck wrote:
> How does one get a list of all the color names supported by
> R and the mapping between color names and hex codes? Color
> names are things like "red", "salmon", etc.
>
colors() gives the list, and as its help page tells you, col2rgb converts
colo
Gabor Grothendieck wrote:
How does one get a list of all the color names supported by
R and the mapping between color names and hex codes? Color
names are things like "red", "salmon", etc.
I am using Windows 2000 Professional and R 1.7.1.
See
?colors
?col2rgb
?rgb
Regards,
Sundar
__
plot(sort(tapply(first.anti.auth.sum,date,mean), partial=1), type="l",
col="yellow",ylim=c(0,2.0))
par(new=TRUE)
\\ par(xaxt='n')
plot(sort(tapply(first.pro.auth.sum,date,mean), partial=1), type="l",
col="red",ylim=c(0,2.0))
check up ?mtext()
but I would like the X axis to be dates, not in
> check up ?mtext()
Thanks, but that's not what I'm looking for. I don't want to change the
text string but the axis labels themselves.
ap
--
Andrew J Perrin - http://www.unc.edu/~aperrin
Assistant Professor of Sociology, U of N
Why does the following code generate a Histogram of 7 cells
instead of 5?
Thanks for your help.
--Tommy
> data(islands)
> str(hist(sqrt(islands), br = 5, col="lightblue", border="pink"))
List of 7
$ breaks : num [1:8] 0 20 40 60 80 100 120 140
$ counts : int [1:7] 40 1 1 2 2 1 1
$ inte
Is there anyone who would like to give me some examples of plots or data
frames on clustering anaylis?
if so, great thanks in advance!
Files can be sent to my big mail box as [EMAIL PROTECTED]
I want t operform cluster analysis on a set of data, the data is composed of
time-evolution rms deviatio
Is there anyone who would like to give me some examples of plots or data
frames on clustering anaylis?
if so, great thanks in advance!
Files can be sent to my big mail box as [EMAIL PROTECTED]
I want t operform cluster analysis on a set of data, the data is composed of
time-evolution rms deviatio
Hi everyone,
I'm looking for an R package that is capable of _write_ as well as read
in netcdf format. Could someone offer advice on this?
I have some data in netcdf format and I only need to make simple
modifications to it.
I thought doing that via R could be handy and would save me reading u
?hist says
breaks: one of:
...
* a single number giving the number of cells for the
histogram,
...
In the last three cases the number is a suggestion only.
On Fri, 27 Jun 2003, Tommy E. Cathey wrote:
> Why does the following code generate a Histogram
58 matches
Mail list logo
