Hi there,
I want to generate a random point pattern using the Neyman-Scott
cluster process using spatstat package in R. After running the
following procedures, why i can not see any figures?
nclust - function(x0, y0, radius, n) {return(runifdisc(n, radius, x0, y0))}
nclust
function(x0, y0,
I am amazed that the traffic decreased the last days. I am afraid that
there is any problem with the spamfilter of my provider.
There are not messages from march 04rd until march 15th, Two messages
from march 16 and nothing else since march 03rd
Please answer also with PM to the mail address. I
Hi,
there's been plenty of messages sent during that period. The r-help
archive is at https://www.stat.math.ethz.ch/pipermail/r-help/ (you can
find this via http://www.r-project.org/ - Mailing Lists - r-help/web
interface - R-help Archives). There are similar pages for r-devel
and the other
Dear Dr. L. Y Hin,
Thank you very much for your help once again!
Would it possible for the gam() just pick a default df to use? For example,
its default df is 4?
?s
s(x, df=4, spar=1)
-
s(col3, 9)
is constructing the smoothing matrix using df=9,
s() is a function
Problem solved,
it was a home made problem.
My mailbox was flooded with spam. The mailbox was open for all addresses
coming to the URL [EMAIL PROTECTED]
There were no problem the last year but now the spammer found the open
mail account.
We changed to defined addresses and mistyped the r-help
has cox process been impelemented in R?
Linda
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
You can open .ssc files in Windows Notepad, then either cut and paste
the commands or type them out. Typing them out may seem a little
excessive but it is one way to learn the syntax and structure of S/R, or
should that be R/S
HTH
Phineas
xpRt.wannabe [EMAIL PROTECTED] 03/17/06 4:44 AM
Dear
Hi
May the following document for factanal.fit.mle be helpful for
you as you requested. You may try with princomp for doing principal
component factor analysis in R.
Maximum Likelihood Estimate of Factor Analysis Model
DESCRIPTION:
Returns an
Use column names:
myts - ts(cbind(First = 1:10, Second = 11:20))
plot(myts)
On 3/17/06, Tong Wang [EMAIL PROTECTED] wrote:
Hi you guys:
I have been wondering if there is any way to change the labeling in
plot.ts( ), for example , if I plot two sequences,
i always got y labels as
No seems that there is an other problem with the spam filter.
Another guy told me that they have the same problem. Some of the mails
will delivered, most of the mails not
Regards Knut
__
R-help@stat.math.ethz.ch mailing list
Dear all,
could someone suggest some strategies for detecting fitting problems in
neural network estimation?
I'm using the nnet package for fitting standardized simulated data (some
thousands estimations are required).
The estimation is generally ok, but sometimes (about 1-3 every 1000) I found
One other approach is to use zoo:
library(zoo)
x - ts(matrix(1:24, 12))
plot(as.zoo(x), ylab = c(A, B))
On 3/17/06, Gabor Grothendieck [EMAIL PROTECTED] wrote:
Use column names:
myts - ts(cbind(First = 1:10, Second = 11:20))
plot(myts)
On 3/17/06, Tong Wang [EMAIL PROTECTED] wrote:
x - factor(1:3, labels = c(b , f, minus))
x
[1] b f minus
Levels: b f minus
I want to change all minus to b. I know that the simplest way to do this is
levels(x) - c(b, f, b)
and also that
x[x == minus] - b
x - factor(x)
works. But why not
x - ifelse(x == minus, b, x)
x -
The problem is that ifelse strips attributes. For example, consider
this which has nothing to do with factors but illustrates the point
with ifelse:
x - ts(1:12)
y - ifelse(TRUE, x, x)
str(x)
Time-Series [1:12] from 1 to 12: 1 2 3 4 5 6 7 8 9 10 ...
str(y)
int 1
On 3/17/06, Göran
Sorry, here is the code:
x - ts(1:12)
y - ifelse(rep(TRUE, 12), x, x)
str(x)
Time-Series [1:12] from 1 to 12: 1 2 3 4 5 6 7 8 9 10 ...
str(y)
int [1:12] 1 2 3 4 5 6 7 8 9 10 ...
On 3/17/06, Gabor Grothendieck [EMAIL PROTECTED] wrote:
The problem is that ifelse strips attributes. For
I do not know how many user are not able to receive the messages.
It seems that only messages from the list server are not delivered to
the mailserver.
I am able to receive private mails and I talked with the sysadmin of the
mailserver. He told me that he is not able to found any new mail from
Knut == Knut Krueger [EMAIL PROTECTED]
on Fri, 17 Mar 2006 14:53:27 +0100 writes:
Knut I do not know how many user are not able to receive the messages.
Knut It seems that only messages from the list server are not delivered to
Knut the mailserver.
Knut I am able to
Here is a snip of code I used in a program that was looping for a really
long time, maybe this can be helpful. You will need the svMisc package.
print(i)
progress(i)
Sys.sleep(.05)
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Jeffrey Racine
Sent:
You could try this:
eol - if (Platform$OS.type == windows) \n else \r
and then use the eol variable in your cat statement.
On 3/17/06, Jeffrey Racine [EMAIL PROTECTED] wrote:
Hi, and thanks in advance for your time.
Background - I am working on a package and wish to have a routine's
I have noticed a slightly puzzling behaviour exhibited by
smooth.spline(). If I do
sss - smooth.spline(x,y)
for a certain pair of data vectors x and y, and then do
length(sss$x)
I get the result ``18''. However if I do
length(unique(x))
I get ``27''. Trying to
I have a dataframe named temp, and another dataframe
named descriptions.
I wish to rename temp, and to call it the names of
a certain column in the dataframe descriptions.
Is there a good way to do this?
A similar question:
I am using a for loop to create several new
dataframes.
I think you meant
.Platform$OS.type == windows
Wouldn't be useful to add an eol element to the .Platform list storing the
os-specific eol character as was done with the file.sep character.
-Christos Hatzis
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf
For example, I saved several .Rdata files with different names:
a.Rdata, b.Rdata ... But, if I forget the specific file name, how can
I get all files list?
One more question about .R file:
source(~/Desktop/allchapters.R,local=FALSE)
Error in parse(file, n = -1, NULL, ?) : syntax error at
135:
On 3/17/06, jia ding [EMAIL PROTECTED] wrote:
For example, I saved several .Rdata files with different names:
a.Rdata, b.Rdata ... But, if I forget the specific file name, how can
I get all files list?
One more question about .R file:
source(~/Desktop/allchapters.R,local=FALSE)
Error in
Hi,
I have tuples of data in rows of a data.frame, each column is a variable
for the 'items' (one per row).
One of the variables is the 'size' of the item (row).
I would like to cut my data.frame into groups such that each group has
the same *total size*. So, assuming that we order by size,
Dan Bolser wrote:
Hi,
I have tuples of data in rows of a data.frame, each column is a variable
for the 'items' (one per row).
One of the variables is the 'size' of the item (row).
I would like to cut my data.frame into groups such that each group has
the same *total size*. So,
?assign, but _don't_ use it; lists are better.
dfr - list()
for(j in 1:9) {
dfr[[as.character(j)]] - ...
}
Don't try to imitate the limited macro approach of other software
(e.g. SAS). You can do all that in R, but it's much simpler and much
safer to rely on list indexing and functions that
Jean-Pierre GERBAL [EMAIL PROTECTED] wrote in
message
news:[EMAIL PROTECTED]
i'm a mathematic teacher and i have a question for R-developers :
is it possible to have (in the future) a boxplot with whiskers from
the first decile to the ninth decile, as usual in secondary french
schools... by
Earl F. Glynn [EMAIL PROTECTED] wrote in message
news:[EMAIL PROTECTED]
Sorry, I mixed up some code in Example #3. Here is a correction:
boxplot.meansd - list(stats = bp.limits,
n = length(len),
boxplot.meansd - list(stats = boxplot.limits,
n =
Hello.
I have some basic question for you!.
Has calculating R2 sense in mixed model? I think no! But i don't why!
Thank in advance for your help
Angelo
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Is there a way of calculating the derivative of a function returned
by splinefun()? Such a function is a cubic spline, whence it has a
calculable derivative, but is there a (simple) way of getting at it?
One workaround that I have thought of is to take a fine grid of
points, evaluate the
If you are just looking for something simple that may be good enough
then assign the largest one to group 1, the second largest to group 2,
..., the 8th largest to group 8 and then start over again with group 1
and so on.
# test data
set.seed(1)
x - sample(100, 100, rep = TRUE)
xs - sort(x)
g -
You can open them in R. On Windows, File:Open Script,
change Files of type to All Files, then open the .ssc file.
Tim Hesterberg
I have some S-PLUS script files (.ssc). Does there exist an R
function/command that can read such files? I simply want to view the
code and practice in R to help
Dear R list,
I've just installed R and then followed the instructions to install the
package geoR from within R by using the following line of code:
install.packages(geoR, contriburl = http://www.est.ufpr.br/geoR/windows;)
This installed okay under the folder C:\documents and settings\ ...
On Fri, 17 Mar 2006, Terrence Murphy wrote:
Dear R list,
I've just installed R and then followed the instructions to install the
package geoR from within R by using the following line of code:
install.packages(geoR, contriburl = http://www.est.ufpr.br/geoR/windows;)
This installed
Try adding depend=TRUE to install.packages() and see if that helps.
Andy
From: Terrence Murphy
Dear R list,
I've just installed R and then followed the instructions to
install the
package geoR from within R by using the following line of code:
install.packages(geoR, contriburl =
On 3/17/2006 9:44 AM, Jeffrey Racine wrote:
Hi, and thanks in advance for your time.
Background - I am working on a package and wish to have a routine's
progress reported. The routine can take some time, and I would like to
inform the user about the routine's progress. I have scoured the
Hi,
I define a generic function that has many parameters (~20). I then want to
define methods for the generic function with setMethod. The default behavior
for setGeneric is to allow dispatching on any argument in the def function,
expect ...
For example,
# define a generic function
Hi,
It appears that deal does not support missing values (NA), so I need to
remove them (NAs) from my data frame.
how do I do this?
(I am very new to R, so a detailed step-by-step
explanation with code samples would be nice).
Some columns (variables) have quite a few NAs, so I would rather drop
Hello ,
Hey! I have been trying to get in touch with you. I finally got a cam so you
can see me when we talk. http://ca.geocities.com/kellyluvsmhe/cam.html
--
This e-mail was sent through the Seismic OS Files file database. The
webmasters of the Seismic OS Files file database take no
Hi Sam
If you are new to R it will definitively pay off to start from the basics.
Go to the help menu- manuals in pdf and select An Introduction to R.
After you read that document you will be able to answer your questions :-)
Good luck!
Francisco
From: Sam Steingold [EMAIL PROTECTED]
* Francisco J. Zagmutt [EMAIL PROTECTED] [2006-03-17 21:09:48 +]:
Go to the help menu- manuals in pdf and select An Introduction to
R. After you read that document you will be able to answer your
questions :-)
I did. I still need help.
The matter is not so much with getting things
I am trying to make predictions with se's using a nlme (kew11.nlme
below). I get an error indicating levels for a factor are not allowed.
I have searched and read Rnews, MEMSS, MASS, R-Help, and other lists
in Spanish where I found questions similar to mine but not solution.
I do not really
Sam Steingold sds at podval.org writes:
Hi,
It appears that deal does not support missing values (NA), so I need to
remove them (NAs) from my data frame.
how do I do this?
(I am very new to R, so a detailed step-by-step
explanation with code samples would be nice).
If you wanted to
This won't be much help, I'm afraid but ...
The problem is, of course, that the prediction is a **nonlinear** function
of the parameters (including possibly the variance components, depending on
what level of the variance hierarchy you are trying to predict). So you need
to somehow estimate how
If I understand it correctly, something like this should do what you want
x[!apply(x, 1, function(y) any(is.na(y)), ]
where x is the dataframe in question.
Hope that helps.
Kevin
- Original Message -
From: Ben Bolker [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Sent: Friday, March
Haifeng Xie xieh at wmin.ac.uk writes:
If I understand it correctly, something like this should do what you want
x[!apply(x, 1, function(y) any(is.na(y)), ]
where x is the dataframe in question.
Hope that helps.
Kevin
I believe he wants to remove *columns* with NAs, not rows
geoR depends on the package sp
which needs to be installed as well
On Fri, 17 Mar 2006, Terrence Murphy wrote:
Date: Fri, 17 Mar 2006 14:16:06 -0500
From: Terrence Murphy [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Subject: [R] can't load geoR
Dear R list,
I've just installed R and then
OK,
I give up, I search everywhere and tried a few things and I just
cannot seem to figure out how to perform a MANOVA in R on my data.
It is a 2x2 within (repeated measures) x 2 between design. Can
someone help me here? Lets assume AxB within and C between. Someone
send
Hi all,
this might be a very stupid question, but I tried for several hours now
and I'm helpless. I'm running windows and R 2-2-1. Recently, every time
when I have an error in my code a message box pops up telling me the
error message. The code won't continue before I click ok in the box.
Before
?manova, which points you to manova.stats
On Mar 17, 2006, at 7:28 PM, John Christie wrote:
I give up, I search everywhere and tried a few things and I just
cannot seem to figure out how to perform a MANOVA in R on my data.
It is a 2x2 within (repeated measures) x 2 between
Hi,
I have some general questions about statistical analysis for a research
dataset and a request for advice on using R and associated packages for a
valid analysis of this data. I can only pose the problem as how to run
multiple ANOVA tests on time series data, with reasonable controls of the
The residual vector autocorrelation paper you cite looks
interesting, but I haven't used it, and unless you've received replies
to the contrary that I haven't seen, I doubt if many people have much
relevant experience.
Have you tried writing to one or more of the authors of
Hi,
After several days of google searches, checking colormaps, R coding,
sleepless nights, etc., I would really appreciate some advice from the
experts :-)
I have a matrix of p-values. It is possible to use image to display this
matrix, using a custom colormap. I am very cautious about using
Your first model y~x is a special case of y~a*x^n, so it should
the comparison of those models should be straightforward: Does the
confidence inteval for n include 1?
It is not so easy to compare the second model with either the first
or the third. While model 1 is a
The coefficient of deterimination, R^2, is defined as the percent of
the variance explained. For a mixed model, we need to ask, percent of
WHICH variance? For more comments, I tried RSiteSearch with various
kew words. The most relevant comment I found was
56 matches
Mail list logo