On Thursday 24 July 2003 05:02, Dowkiw, Arnaud wrote:
Dear R-helpers,
I have performed a PLS regression with the mvr function from the pls.pcr
package an I have 2 questions : 1- do you know if mvr automatically centers
the data ? It seems to me that it does so...
Yup, it does... common
Hi,
I've made a selfStart function for use with gnls and the
following piece of code works nicely:
check1 - gnls(y ~ spot.shape.fct(xcord, ycord, background, spotintensity,
rho, sigma, delta, mux, muy),
start=getInitial(y ~ spot.shape.fct(xcord,
Dear all,
I wonder if it is possible to obtain marginal tests for effects in generalized linear
models. Indeed, the anova function produces sequential tests and it doesn't have any
type argument to specify that we would like marginal tests instead, as in the
similar anova function for lme
?drop1 for the valid marginal tests.
If you want `type III' tests (which aren't marginal in the usual sense)
see Anova in package car, which also does `type II' tests (as performed by
drop1).
On Thu, 24 Jul 2003, Eve Corda wrote:
I wonder if it is possible to obtain marginal tests for
Dear Kjetil,
As I already mentioned, it appears that there isn't a function
available calculating the quantiles directly (at least, it doesn't appear
in the C source of ctest). So as I already suggested, uniroot (or a similar
C routine which calls the corresponding C code directly) is probably
I only concern about the boundary of the shape and I dont concern gray
level.
U[m] = x[m] + jy[m] m= 0,1,,N-1 is the N points on the boundary.
It is easy to express U as a discrete Fourier series.
But it seems like that wd in the wavethresh package cannot deal with complex
values.
And it is not
OK, I really am struggling with this one! Forgive me if I am being stupid
I am running R 1.7.1 on Suse Linux 8.1. I connect to the internet through a proxy so
I have:
IAHC-LINUX03:~ # echo $http_proxy
wwwcache.bbsrc.ac.uk:8080
IAHC-LINUX03:~ # echo $HTTP_PROXY
wwwcache.bbsrc.ac.uk:8080
One way is to use an object-oriented design and wrap up the reference
functionality in a common superclass. At
http://www.maths.lth.se/help/R/ImplementingReferences/ I have got some
discussions which are in line what you are trying to achieve and that
you might be able to adopt.
Also, note that
Dear R-helpers!
May I kindly ask the pure statistics-experts to help me for a
purpose which first part is not directly concerned with R.
Consider two distribution functions, say f and g. For both, the
median is smaller than a half. Now, the multiplicative or additive
linkage of both distribution
I've made some changes (in bold) following your comments :
# class definition :
setClass(MyClass, representation(mynumber=numeric));
# the initilization method :
setMethod(initialize,MyClass,
function(.Object)
{
[EMAIL PROTECTED] - 10;
return(.Object);
On Thu, 24 Jul 2003, michael watson (IAH-C) wrote:
OK, I really am struggling with this one! Forgive me if I am being stupid
I am running R 1.7.1 on Suse Linux 8.1. I connect to the internet
through a proxy so I have:
IAHC-LINUX03:~ # echo $http_proxy
wwwcache.bbsrc.ac.uk:8080
Hi all,
I tried to produce some kriged surfaces with geoR (latest version). The size
of the grid should be around 900 x 650 cells (what I find is not a very big
grid), and the number of points is around 2500. The command krige.conv
stopped after arround 5 min saying it can not allocate a
Thanks for all the help.
Juergen Gross supplied a program which does just what Belsley
suggested.
Chuck Cleland, John Fox and Andy Liaw all made useful programming
suggestions.
John Fox asked
(1) I've never liked this approach for a model with a constant, where
it
makes more sense to me to
For distribution functions F and G we have the (Frechet) bounds:
max{0, F(x)+G(y) -1} = H(x,y) = min{F(x),F(y)}
where H is the joint df of (X,Y) having marginals F and G. If
X and Y are comonotonic (Schmeidler (Econometrica, 1989)), that is
if there is a random variable Z, such that X =
Dear Salvatore,
Assuming that you mean convolution when you write
additive linkage, the answer is that there is no general
answer. It will depend heavily on the joint distribution
of the two random variables.
Just to give a simple example, let X~f, Y~g, and
P(X=0.4)=P(Y=0.4)=1. Then, X and Y
Hello Professor
If you are suggesting that I am simply missing the http://; part of my cache URL, or
that I am missing a trailing /, then I pre-empted this response and it still doesn't
work.
I have tried setting both http_proxy and HTTP_PROXY to all of:
wwwcache.bbsrc.ac.uk:8080
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Laurent Faisnel
Sent: den 24 juli 2003 14:04
To: John Chambers
Cc: [EMAIL PROTECTED]; Duncan Murdoch
Subject: Re: [R] S3 and S4 classes
[snip]
Thank you for your fast and useful answers. I've
On 23 Jul 2003 09:12:07 -0500
Douglas Bates [EMAIL PROTECTED] wrote:
High-level control of axes in xyplot is implemented by the scales
argument to xyplot. You can include components 'at' and 'labels' in
a list given as the scales argument. See ?xyplot.
Wladimir Eremeev [EMAIL PROTECTED]
Dear Peter,
At 08:24 AM 7/24/2003 -0400, Peter Flom wrote:
(1) I've never liked this approach for a model with a constant, where
it
makes more sense to me to centre the data. I realize that opinions
differ
here, but it seems to me that failing to centre the data conflates
collinearity with
On Wed, 23 Jul 2003, Jason Liao wrote:
Thanks for Prof. Ripley and Andy for your technical explantion. It
seems that that the real CPU speed has not advanced as fast as these
Ghz or other performance indicator suggest.
Yes, my program is totally CPU intensive.
It may not be even if you
You're right, I have tried keeping the double quotes in there too and that doesn't work
What is clear is that R is completely ignoring my http_proxy environment variable, as
it's error says unable to connect on port 80, when every single one of my http_proxy
attempts has stipulates port
Hi,
is there advice how to subset a data.frame, where until R version 1.7.0
it was possible to write
data[[subset]]
which meant the same as
data$subset
(data is a data.frame). From 1.7.1 on this does not seem to work longer.
Could there be a bug in downwards compatibility?
sincerely
Greetings,
Does anyone know how to get an id number in the little header
above each individual plot within a trellis plot? The default
seems to be to print the word id and add a line indicating on
a linear scale where the current id sits.
Thanks in advance for any help you can send,
Suzette
I'm trying to use the source command to run commands from a file. For
instance: source(do.R), where do.R is a file in the same directory in
which I am running R.
The contents of do.R are:
ls()
print(hello)
sum(y1)
mean(y1)
After source(do.R), all I see is:
source(do.R)
[1] hello
I'm
I got the following using the benchmark program at
http://www.sciviews.org/other/benchmark.htm
under Windows 2000 prof., 256 MB of RAM
R Benchmark 2
=
Number of times each test is run__: 3
I. Matrix calculation
-
Creation,
Hi all,
I was wondering whether it is possible to plot math functions, for
example sin, cos or a Gaussian type function, in R, and if so, how to
do it. I have been searching through the archives and the R manual but
had no luck in finding any hints on how to go about this.
Any help is much
Jonck van der Kogel wrote:
Hi all,
I was wondering whether it is possible to plot math functions, for
example sin, cos or a Gaussian type function, in R, and if so, how to do
it. I have been searching through the archives and the R manual but had
no luck in finding any hints on how to go about
Dear John
An interesting discussion!
I would be the last to suggest ignoring such diagnostics as Cook's D;
as you point out, it diagnoses a problem which condition indices do not:
Whether a point is influential.
OTOH, condition indices diagnose a problem which Cook's D does not:
Would shifting
Peter Muhlberger wrote:
I'm trying to use the source command to run commands from a file. For
instance: source(do.R), where do.R is a file in the same directory in
which I am running R.
The contents of do.R are:
ls()
print(hello)
sum(y1)
mean(y1)
After source(do.R), all I see is:
Dirk Repsilber wrote:
Hi,
is there advice how to subset a data.frame, where until R version 1.7.0
it was possible to write
data[[subset]]
which meant the same as
data$subset
(data is a data.frame). From 1.7.1 on this does not seem to work longer.
Could there be a bug in downwards
I was wondering whether it is possible to plot math functions, for example
sin ...
also maybe:
?plot
snipped from the help page
Examples:
...
plot(sin, -pi, 2*pi)
unsnip
bob
-Original Message-
From: Uwe Ligges [mailto:[EMAIL PROTECTED]
Sent: Thursday, July 24, 2003
On Thu, 24 Jul 2003 16:27:28 +0200, Dirk Repsilber
[EMAIL PROTECTED] wrote :
Hi,
is there advice how to subset a data.frame, where until R version 1.7.0
it was possible to write
data[[subset]]
which meant the same as
data$subset
(data is a data.frame). From 1.7.1 on this does not seem
On Thu, 24 Jul 2003 10:42:05 -0400, Peter Muhlberger
[EMAIL PROTECTED] wrote :
I'm trying to use the source command to run commands from a file. For
instance: source(do.R), where do.R is a file in the same directory in
which I am running R.
The contents of do.R are:
ls()
print(hello)
sum(y1)
Suzette Blanchard [EMAIL PROTECTED] writes:
Does anyone know how to get an id number in the little header
above each individual plot within a trellis plot? The default
seems to be to print the word id and add a line indicating on
a linear scale where the current id sits.
See the
Well, we feel that source does work the way that it is documented to
work. Please read the documentation and notice that the entire file
is evaluated as one step in the read-eval-print loop. If you want to
print the results of individual function calls you will need to change
your script to
On Thu, 24 Jul 2003 16:55:55 +0200, Jonck van der Kogel
[EMAIL PROTECTED] wrote :
Hi all,
I was wondering whether it is possible to plot math functions, for
example sin, cos or a Gaussian type function, in R, and if so, how to
do it. I have been searching through the archives and the R manual
Dear all,
I am a relative newcomer to the R language, and am sussing out
the possibilities of using R to do integer programming (which I am
also new to). Is there something along the lines of the NUOPT S-
PLUS package that is available, or on the way, for R? Are there
other options for
Thanks to everyone for their suggestions on getting source to print! It
seems not everyone was aware of a couple options that gets source to print
out everything. I'm now using the following command:
source(do.R, print.eval=TRUE, echo=TRUE)
__
[EMAIL
All:
I am trying to use gam in a scatterplot smoothing problem.
The data being smoothed have greater 1000 observation and have
multiple humps. I can smooth the data fine using a function
something like:
out - ksmooth(x,y,normal,bandwidth=0.25)
plot(x,out$y,type=l)
The problem is when I try
Hi,
I've got a wireframe 3D surface plot, but I don't want a frame around it.
Is there any way to remove the frame, or (worst case)
change the color of the frame to the background color (which looks like
grey).
I'm using ver 1.7.1
I've tried frame.plot = F, but that doesn't seem to work for
Still works for me in both 1.7.1 and 1.8.0 (in development)
data(women)
women[[height]]
[1] 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72
women$height
[1] 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72
Can you provide an example of how it is failing for you?
Dirk Repsilber [EMAIL PROTECTED]
You can specify some options in par.box. Use col=NA to make the frame
transparent. See example below (which was modified from the help file).
See also the scales parameter if you want to remove the arrows as well.
Cheers,
Jerome
library(lattice)
x - seq(-pi, pi, len = 20)
y -
Jerome Asselin wrote:
You can specify some options in par.box. Use col=NA to make the frame
transparent. See example below (which was modified from the help file).
See also the scales parameter if you want to remove the arrows as well.
Cheers,
Jerome
library(lattice)
x -
hi jerome,
thank you for your quick reply.
your advice removes the 3D box in which the 3D plot is generated,
but i like THAT box. (sorry for being unclear earlier)
what i want to do is remove the 2D frame (a box of thin black lines) that
circumscribes all of my plot area. any ideas?
thanks
Hi Y'all,
Alexis, I think Jerome's example works except for one change:
This one has a box but no wires:
wireframe(z ~ x * y, g, drape = TRUE,
perspective = FALSE,
aspect = c(3,1), colorkey = FALSE,
par.box = list(col=1),col=NA)
versus with no
hi ladies and gents:
I am using R 1.7.0. Alas, the pch=. in the postscript device is too
large for me, but apparently not scaleable via cex. Maybe a bug, though
I read a complaing about this in the list archives from 1999, and given
that it is still around, maybe it is a feature.
Now, I
hi jake,
thanks for your e-mail.
what i am actually trying to do is remove the 'picture frame' 2D box (the
super-frame) that circumscribes the entire figure-- for example, the left
vertical side of this box i'm concerned about is to the left of the 'z' axis label.
i hope i'm finally making my
Finally, I think this works. The function trellis.par.get() contains all
default values of the parameters. I had to reset the option axis.line as
below. The bounding (square) box is made transparent. See
trellis.par.get() for the list of all default parameter values.
Cheers,
Jerome
On Thu, 24 Jul 2003, Hotz, T. wrote:
Dear Kjetil,
As I already mentioned, it appears that there isn't a function
available calculating the quantiles directly (at least, it doesn't appear
in the C source of ctest). So as I already suggested, uniroot (or a similar
C routine which calls the
Hello to all, I'm a biologist trying to tackle a fish (Poisson Regression) which is
just too big for my modest understanding of stats!!!
Here goes...
I want to find good literature or proper mathematical procedure to calculate a
confidence interval for an inverse prediction of a Poisson
Hello,
I've done most all the various steps outlined in a
recent posting to the mailing list archives
(and in the help files) to load and run R-WinEdt. I
can get it to run fine but I am not successful in
getting it to interface with RGui (1.6.2, not
minimized). I try to use the R-line,
hi all,
contourplot generates a nice plot, but there are cut lines that don't have
labels in my figure.
i think this is because there's a high ridge stretching
north-west/south-east entirely across the plot, diagonal-to-diagonal,
and contourplot only wants to label each elevation-level only once,
Right click on your WinEdt icon and go to Properties menu-item. In the
Target box look for something like,
[PathToWinEdt]\WinEdt\WinEdt.exe -C=R-WinEdt -e=r.ini
If you don't find something like this, make the change and try it.
* It may be better to have a new icon for R-WinEdt on your desktop
Dear R List,
Is there a function to plot vector fields in R?
I'm looking for something similar to PlotVectorField in Mathematica
(see below), e.g. given functions f(x) and f(y), I would like to plot
the resulting field of vectors (f(x),f(y)) over some range of x and y.
PlotVectorField[{f(x),
1. If you provide a toy data set with, e.g., 5 observations, to
accompany your example, it would be much easier for people to try out
ideas and then give you a more solid response.
2. Have you tried something like log(dose+0.5) or I(log(dose+0.5))
in your model statement in conjunction
If your datasets are small, then it is not difficult to roll your own
using optim(). If you look at
http://www.qimr.edu.au/davidD/sib-pair.R
you will find such a routine at lines 1511-1562. This calls
kinship.rel() that produces NRM etc.
--
| David Duffy (MBBS PhD)
Hi
I am using R 1.7.1 on RH linux 9.0
sum(unlist(lapply(ls(),function(x)object.size(get(x)/1024^2
[1] 2.424263
so I am not using much memory (I have a gig of ram on my machine)
now in nlme
gtest-groupedData(log(X8)~Time|sub,all[,c(names(all)[1:9],X8)],outer=~A*B)
object.size(gtest)/1024
Robin,
This is not a direct answer to your question, but there exists a pretty
good linear, mixted and integer programming package called lp_solve. I
have used it successfully on a couple of projects. As far as I know it is
in the public domain and there is source code available. There is
Hi All,
Is it possible for system.time() to measure the time
it takes for a machine to evaluate more than one R
expression?
For example,
# This I can do:
system.time(x - rnorm(10))
[1] 0.07 0.00 0.13 0.00 0.00
# But this I can't:
system.time(x - rnorm(10); new - sample(x, 10,
On Thu, Jul 24, 2003 at 07:10:27PM -0700, Robert Keefe wrote:
# This I can do:
system.time(x - rnorm(10))
[1] 0.07 0.00 0.13 0.00 0.00
# But this I can't:
system.time(x - rnorm(10); new - sample(x, 10, replace=T))
Error: syntax error
Just use curly braces:
system.time({x -
60 matches
Mail list logo