further developments since we last communicated
about this problem?
cheers,
Rolf Turner
P. S. BTW Chris: Does plot.im() screw up with other images? E.g. what
happens if you do
with(bei.extra,plot(elev))
???
R. T.
__
R-help@r-project.org mailin
On 16/06/14 12:35, Jim Lemon wrote:
On Mon, 16 Jun 2014 08:54:41 AM Rolf Turner wrote:
On 15/06/14 15:34, Paul Rydelek wrote:
Relatively new user with a question regarding abline.
I want to draw a st. line fit to some data but I want the line to span
the
range of the x-data and NOT the plot
xmax (of x data).
Can this but done with abline or do I have to use lines and construct the
line from a+bx?
A function with this capability was posted by Remko Duursma on 22/1/09. See:
http://article.gmane.org/gmane.comp.lang.r.general/137344/match=ablinepiece
cheers,
Rolf Turner
P. S. Jim
very much like homework. The r-help list has a "no homework"
policy. Do your own homework. That way you might actually learn something!
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
alues, e.g. values can't exceed the interval /[-3.7; 3.7]/. How
can I implement such an interval?
What you are asking for amounts to drawing a sample from a *truncated*
normal distribution.
Try
RSiteSearch("truncated normal")
to find appropriate functions to effect this
egards
Zakir
On Tuesday, June 3, 2014 11:32 PM, Rolf Turner
wrote:
You might find it helpful to look at the sim.hmm() function in the
"hmm.discnp" package, or the simHMM() function in the "HMM" package.
cheers,
Rolf Turner
On 03/06/14 21:17, Bukar Alhaji wrote:
> De
You might find it helpful to look at the sim.hmm() function in the
"hmm.discnp" package, or the simHMM() function in the "HMM" package.
cheers,
Rolf Turner
On 03/06/14 21:17, Bukar Alhaji wrote:
Dear R buddies,
Sorry for this silly question but am new to R. I am tryin
tax is easy to learn and quite powerful. Get off your
duff and learn it.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
a
On 28/05/14 00:48, Jim Lemon wrote:
Hi all,
I thought it was all over, but when I started my R help (help.start()) to
check something, none of the help files in plotrix were available, each
displaying an error message that some 2-3 Gb of memory could not be
allocated. The same message appeared wh
people just use the "system" library; e.g.
/usr/local/lib64/R/library; you will need to assume superuser privileges
(use sudo, perhaps) to copy the package to that location.
It's as easy as that.
cheers,
Rolf Turner
__
R-help@r-
ex=0.2)
But this sets the font size to 0.5 for both the x and y axis? How can I
make the y-axis to 1.5?
How about:
> boxplot(x,cex=0.2,axes=FALSE)
> axis(1,cex.axis=0.5)
> axis(2,cex.axis=1.5)
Does that get you what you want?
cheers,
Rolf Turner
__
)
Try doing something like:
newbm <- bm[!duplicated(bm),]
and then hit "newbm" with both deldir() and delaunayn().
cheers,
Rolf Turner
On 23/05/14 03:36, Raphael Päbst wrote:
Hello again,
I have found further depths of confusion concerning delaunay
triangulations to explore.
Her
nt to do, and/or what you want to do
makes no sense at all.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and pr
On 23/05/14 06:59, Claire wrote:
Thanks Bert,
Will post on r-sig-mixed-models list. Can't help it being in html
though as i sent the query via -email.
Huh? What do you mean by "-email"? Do you really mean "email" (without
the minus sign)? We *all* send our queries via email. This *is* emai
On 21/05/14 23:34, Raphael Päbst wrote:
I believe you are right. A night of sleep has done wonders for my
understanding of the problem.
Oh sleep it is a gentle thing
Beloved from pole to pole!
Thank you for your patience and help!
Patience? Moi? This must be some new use of the word "pati
output of each function (RTFM!) and you will find that there is no
problem.
cheers,
Rolf Turner
(author of "deldir")
On 21/05/14 09:33, Raphael Päbst wrote:
Thanks for the answer!
I'll post a sample tomorrow, I have however found the following:
triang.list() gives me the coo
colour given by ccc[1] again. The
tile corresponding to (x[4],y[4]) gets coloured with the colour given by
ccc[2] again. And so on.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
uld have inferred from context.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Dear Patrick,
OK. I think I have figured out what you are trying to do.
You want:
* each "type" of point (where type is the corresponding
level of "Concept") to be represented by a different symbol,
* the size of the symbol to be determined by the value of "New".
The "maxsize" ar
essary and nonsensical here.)
But path[5] has not been set yet --- so it is NA. Crash!!!
Think through what you are doing more carefully.
(4) Talk to your instructor! I am *fairly* confident that he or she
will not bite.
cheers,
Rolf Turner
On 16/05/14 05:03, chantel777 wrote:
The q
On 15/05/14 19:27, Patrick Connolly wrote:
On Thu, 15-May-2014 at 03:00PM +1200, Rolf Turner wrote:
|>
|> Dear Patrick,
|>
|> (a) Your attachment did not come through; only attachments with
|> extensions belonging to a limited class are passed through by the
|> list softwa
p users to accomplish the tasks that they need to do, not to hinder,
impede and stone-wall them.
The tell-tale word here is "administrator" --- those of that breed seem
to get their kicks from frustrating those who do the actual work.
cheers,
Rolf Turner
Dear Patrick,
(a) Your attachment did not come through; only attachments with
extensions belonging to a limited class are passed through by the list
software.
Perhaps you could email your data set to me privately,
(b) I have only had a quick glance at your code. A few comments:
* In yo
.
RSiteSearch("Gumbel") turns up lots of hits.
cheers,
Rolf Turner
On 14/05/14 20:42, Frede Aakmann Tøgersen wrote:
Oh, I see now. I'm blushing.
Yours sincerely / Med venlig hilsen
Frede Aakmann Tøgersen
Specialist, M.Sc., Ph.D.
Plant Performance & Modeling
Technology &
bout the Gumbel
distribution with a pointer to the help on the *Weibull* distribution? I
was very puzzled by this and thought that I must be
missing/misunderstanding something.
cheers,
Rolf Turner
Yours sincerely / Med venlig hilsen
Frede Aakmann Tøgersen
Specialist, M.Sc., Ph.D.
Plant Perform
CRAN repository.
Formerly available versions can be obtained from the archive.
Archived on 2014-05-01 as undefined-behaviour errors were not
corrected (and previously archived on 2014-03-07).
You could try retrieving the package from the archive and fixing it
yourself. Good luck!
e(c(0.01,0.10),n,TRUE).
You really need to learn something about R if you are going to use R.
Start with "An Introduction to R" available (under "Manuals") from the R
web site.
You also need to learn to express yourself clearly and unambiguously.
Do not expect your readers to b
On 05/05/14 17:05, Ragia Ibrahim wrote:
Dear group,
How to generate uniform probability choosing p to be 2% and 5%, in separate
trials for 100 times.
No idea WTF you are talking about. Can you formulate a question that is
comprehensible to the human mind?
cheers,
Rolf Turner
s own work.
Well put.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
On 04/05/14 10:16, Dr Eberhard W Lisse wrote:
Google Pressdram :-)-O
el
On 2014-05-03, 23:42 , Rolf Turner wrote:>
On 04/05/14 00:05, Dr Eberhard W Lisse wrote:
Thank you very much, Mr Arkell.
I don't get it. Can anyone explain the (joke? allusion?) ?
Thank you.
cheers,
Rol
On 04/05/14 00:05, Dr Eberhard W Lisse wrote:
Thank you very much, Mr Arkell.
I don't get it. Can anyone explain the (joke? allusion?) ?
cheers,
Rolf Turner
On 2014-05-03, 07:11 , Bert Gunter wrote:
By making the effort to learn R?
See e.g. the "Introduction to R" tuto
o
All of the foregoing may be a red herring, off the point and have
nothing to do with the problem that you are having. I don't really know
what I'm doing. I just follow recipes. :-) Good luck.
cheers,
Rolf Turner
On 29/04/14 07:26, Fong Chun Chan wrote:
Hi,
I am trying to get
, Apr 25, 2014 at 12:50 AM, Rolf Turner mailto:r.tur...@auckland.ac.nz>> wrote:
As usual I am too lazy to fight my way through the rather convoluted
code presented, but it seems to me that you just want to calculate a
log likelihood. And that is bog-simple:
The log likelihoo
plnorm().
cheers,
Rolf Turner
On 25/04/14 07:27, Jacob Warren (RIT Student) wrote:
I'm trying to figure out if there is a way in R to get the loglikelihood of
a distribution fit to a set of data where the parameter values are fixed.
For example, I want to simulate data from a given alte
about R.
(4) Don't expect others to do your work for you.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and pr
Sorry, your communication *was* on-list. I guess my eyes are going!
Please ignore the first sentence of my previous post.
cheers,
Rolf
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guid
window if you
don't have or know the window a priori.
cheers,
Rolf Turner
On 16/04/14 01:23, Monaly Mistry wrote:
Hi Rolf,
Thank you very much, I just have one other question, do I need to
change the data frame into ppp format? I've tried using as.ppp but get
the following error me
The behaviour that you get is exactly the behaviour that I, at least,
would expect, and it seems to me to be exactly the correct behaviour.
I do not understand what you are complaining about.
cheers,
Rolf Turner
On 11/04/14 06:31, ivo welch wrote:
I just spent about an hour bug-tracking
On 10/04/14 09:28, thanoon younis wrote:
hi
i want to simulate multivariate dichotomous data matrix with categories
(0,1) and n=1000 and p=10.
Nobody's stopping you! :-)
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
data frames are saved (using the
save() function).
Attachment did not come through.
What does all.equal() say about your two data frames?
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE d
) but in R?
Have you tried Google? Have you tried RSiteSearch()?
These would have led you fairly quickly to the NHPoisson package which
*might* be able to accomplish what you want.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https
have *installed* spatstat!
W <- ripras(ddd) # Where "ddd" is the name of your data frame.
X <- ppp(x=ddd[,1],y=ddd[,2],window=W)
n5 <- nnwhich(X,k=1:5)
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mail
e1 <- df1[1:5,1:2]
e2 <- df2[1:5,1:2]
ee <- rbind(e1,e2)
If that throws the error, include dput(e1) and dput(e2) in your posting.
If it *doesn't* then this might give you some insight into just what is
triggering the error.
Look at rownames(df1) and rownames(df2) as
On 29/03/14 01:34, peter dalgaard wrote:
On 28 Mar 2014, at 02:37 , Rolf Turner wrote:
So there you are. Feel enlightened?
Somewhat, actually, but not to such an extent as to have reached nirvana.
"Promises" blow me away.
Here's the most useful part of the post: to ge
On 28/03/14 12:49, Duncan Murdoch wrote:
On 27/03/2014, 7:17 PM, Rolf Turner wrote:
I was under the impression that
do.call(foo,list(x=x,y=y))
should yield the same result as
foo(x,y).
However if I do
x <- 1:10
y <- (x-5.5)^2
do.call(plot)
I get the expected pl
ed and
subject to change.
Am I being bitten by an instance of this phenomenon? Seems strange.
I would be grateful for enlightenment.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the p
On 26/03/14 12:51, David Winsemius wrote:
On Mar 25, 2014, at 9:52 AM, Luo Weijun wrote:
Dear Robert and R project team,
I notice that the Google search function on the R mail list archives page has
stopped working for quite a while, http://tolstoy.newcastle.edu.au/R/.
Is there any solution o
On 24/03/14 08:37, David Winsemius wrote:
On Mar 22, 2014, at 3:22 PM, Tham Tran wrote:
Dear R users,
Given two vectors x and y
a=1 2 3
b=4 5 6
i want to combine them into a single vector z as 1 4 2 5 3 6
One way:
c( matrix(c(a,b), nrow=2, byrow=TRUE) )
It is more perspicuous to use
On 19/03/14 10:59, Prof Brian Ripley wrote:
On 18/03/2014 21:41, Rolf Turner wrote:
I am currently having, uh, difficulties, with the latest version of
lme4, which I did not have with an elderly version that I was using
previously. To check things out I installed the elderly version in a
rsion!!!
Can anyone explain WTF is going on? What is hanging around in .RData
that causes library() to ignore the "lib.loc" argument? How can I keep
library() from ignoring the "lib.loc" argument?
Thanks.
cheers,
Rolf Turner
__
can figure that out, then create an appropriate argument (a
list) for bxp() and call bxp() with that argument.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http
cle on quadrilaterals and
outlines a procedure for finding the area centroid (I presume that's
what you actually want) of a convex quadrilateral.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
P
Sys.getenv("LD_LIBRARY_PATH"),sep=":"))
You could put that line into your .Rprofile so that you wouldn't need to
enter it on each startup.
HTH
cheers,
Rolf Turner
On 12/03/14 16:00, Ross Boylan wrote:
I am trying to ensure that a directory is searche
xis is numbered only
relatively"; this makes no sense at all. What *do* you want? The
y-axis labelling that you get will be/is the appropriate labelling.
The arrows will go where you tell them to go, so if they are "going
diagonal" you are telling them to go diagonal.
chee
On 08/03/14 10:22, Yuanzhi Li wrote:
Hello, everyone,
Do you have any idea to get a set of random points within a regular
pentagon? Thanks in advance!
You can easily do this using the spatstat package.
cheers,
Rolf Turner
__
R-help@r-project.org
one that does things right, but doing it right is still the
right thing to do.
cheers,
Rolf Turner
Thanks very much,
Christian
From: Duncan Murdoch [murdoch.dun...@gmail.com]
Sent: Saturday, March 01, 2014 10:50 AM
To: Christian De Santis; r-help@r-project
ta>(letterR)
plot <http://inside-r.org/r-doc/graphics/plot>(quadrats(letterR, 5, 7))
but my environment does not know what square() is.
You have to install and load the spatstat package!!! Sigh.
cheers,
Rolf Turner
__
R-help@r-p
:
dfX <- data.frame(i=as.vector(row(X)),j=as.vector(col(X)),
value=as.vector(X))
where X is your matrix. The fact that X is non-square is of course
completely irrelevant.
cheers,
Rolf Turner
__
R-help@r-project.org mailing l
heers,
Rolf Turner
On 01/03/14 13:40, Bill wrote:
Hi. I have a dataframe called radSamp
which looks like this:
head(radSamp[,c(3,4,5)])
Latitude Longitude Value
1 39.16094 140.488345
2 32.84428 -117.224047
3 35.80605 139.378928
4 35.07816 -106.612350
5 35.83174 136.2027
On 27/02/14 05:43, varin sacha wrote:
Many thanks Rolf,
These codes below are ok :
LinearModel.1 <- lm(GDP.per.head ~ Competitivness.score +
Quality.score, data=Dataset)
summary(LinearModel.1)
predict(LinearModel.1, se.fit = FALSE, scale = NULL, df = Inf,interval =
c("prediction"),level = 0.95
argument from your call and don't clutter things up with
irrelevancies.
But here you *don't* want the default (which is "none"). You want
"prediction". So set interval="prediction" in your call. Actually you
can just set interval="p" because of
ere are no further progress reports.
I have attached a tar-ball and a zip archive both containing the
necessary code to demonstrate the forgoing. Unpack, source all the *.R
files and then source "demoScript" to see it happen before your very
eyes. :-)
cheers,
Rolf Turner
Please see at end.
On 26/02/14 01:13, Duncan Murdoch wrote:
On 14-02-25 3:08 AM, Rolf Turner wrote:
I have a function that makes use of the ode() function from the
"deSolve" package. I am trying to find a way of getting it to put out a
"progress report" every &q
ign matrix. I.e. replace
each factor column by k-1 columns of dummy variables (where k is the
number of levels of the given factor). Note that "x" should really be a
*matrix*, not a data frame although it seems that data frames (all of
whose columns are numeric) get coerced to matri
t? I'm going mad!
***MAD*** I tell you! :-)
Suggestions as to a better way of accomplishing my desired goal of
producing progress reports would also be welcome.
I am not at all sure that assigning "tdone" in parent.env(environment())
is the right thing to do.
On 22/02/14 11:53, Greg Snow wrote:
Why are you testing your data for normality? For large sample sizes
the normality tests often give a meaningful answer to a meaningless
question (for small samples they give a meaningless answer to a
meaningful question).
Fortune!!!
cheers,
Rolf
___
is generally considered that
testing for normality is a waste of time and a pseudo-intellectual
exercise of academic interest at best.
cheers,
Rolf Turner
Hope this helps,
Rui Barradas
Em 21-02-2014 15:59, Gonzalo Villarino Pizarro escreveu:
Dear R users,
Please help with with this maybe b
hing comes before the superscript, replace the nothing by {}.
The following seems to work:
plot(1:10,xlab="",ylab=expression(Temperature~{}^o*C))
The degree symbol however may be rendered better by "degree" than by "{}^o".
cheers,
Rolf Turner
On 15/02/14 06:
testsum <- sum(testcsv[2,2:4])
should give what you want. The use of with() is uncalled for in this
context. The with() function allows you to refer to (e.g.) columns of
a data frame by name, as if these columns were objects in your workspace
("global environment"). That is *not
e same approach. I.e. Mod(z)^2 is
equal to Re(z)^2 + Im(z)^2.
Check your code carefully.
Also it would be a good idea to try a "toy example", something much less
complicated than your real problem, and see what happens there.
cheers,
Rolf Turner
__
s() cannot handle what you want handled. You
need to code up the function to be minimized yourself, using Mod() and
then use optim() or nlm() to minimize it that function with respect to
"theta".
HTH
cheers,
Rolf Turner
On 10/02/14 11:45, Andrea Graziani wrote:
Hi everyone,
rote:
On Feb 4, 2014, at 4:57 PM, Rolf Turner wrote:
If I have a character such as "£" stored in a object called "xxx", how can I obtain the
hex code representation of this character? In this case I know that the hex code is "\u00A3", but
if I didn't, ho
rable (id est, have the same weights) this is not a
problem.
As I say, this is off the top of my head. Others older (???) and wiser
than I may correct me.
cheers,
Rolf Turner
On 05/02/14 11:56, IamRandom wrote:
I am trying to do weighted Poisson regression. I have count data.
Simple ex
;? I thought I could do
xxx <- paste("\u",yyy,sep="")
but R won't let me use "\u" "without hex digits". How can I get around
this?
Thanks.
cheers,
Rolf Turner
__
R-help@r-project.org mai
t some insight into what's going on if you were to read
Annette Dobson's book "An Introduction to Generalized Linear Models"
(Chapman and Hall, 1990).
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch
For what it's worth, I would like to say that I concur completely with
Don and Bert. (Also I would like second Bert's vote of thanks to Don
for expressing the position so clearly.)
cheers,
Rolf Turner
On 04/02/14 09:56, Bert Gunter wrote:
Don:
First, I apologize if this is
able=TRUE)
HTH
cheers,
Rolf Turner
On 03/02/14 10:11, David Parkhurst wrote:
Don't worry. I'm also plotting Y against X, which is what really
matters. But I want the plots I'm asking about to see how much
persistence each of those variable has with time.
Thanks for your warnin
l tell you how. But I
repeat: ***DON'T***!!!
cheers,
Rolf Turner
On 03/02/14 08:09, David Parkhurst wrote:
I've tried to figure out how to do this from what I read, but haven't
been successful. Suppose I have a dataframe with variables Date, X, and
Y (and maybe U, V, and
On 02/02/14 15:37, mohammad javad Azadchehr wrote:
Hi
I want R code for bayesian CIF competing risks.
best regards
I want wealth, wisdom, good looks, and to play full forward for the
Sydney Swans.
cheers,
Rolf Turner
__
R-help@r-project.org
On 01/02/14 08:13, Duncan Murdoch wrote:
On 31/01/2014 1:46 PM, Paul A. Steckler wrote:
I'm fairly new to R, and have checked the R FAQ and done an
RSiteSearch for help
on this topic, to no avail.
I want to write some R code that has functions at the top-level that
are not visible when
the code
source: forecast_5.0.tar.gz
MacOS X binary: forecast_5.0.tgz
Windows binary: forecast_5.0.zip
Old sources:forecast archive
so the MacOS X binary is definitely there.
Dunno what you're doing, but you are doing something weird to mess
yourself up like this.
cheers,
dput().
That being said, you should probably (in the first instance) be
contacting the maintainer(s) of the ETAS package.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
) function loads ***packages*** from a
library (= collection of packages). Correct terminology is essential to
effective communication.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read t
On 24/01/14 21:30, PIKAL Petr wrote:
Why difference in the result? rbind and cbind function are something
wrong?
Or
I, am something wrong?
you are something wrong.
Fortune?
cheers,
Rolf
__
R-help@r-project.org mailing list
https://sta
s why.
Probably because it is ***extremely*** rare that any sane person would
be interested in the mode of the data.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
On 19/01/14 06:20, peter dalgaard wrote:
Offhand, I'd say that the issue isn't so much the starting values as the risk of
hitting the upper bound on gamma. If you get gamma set to min(x) at any point during
the iteration, the density at the corresponding point is zero, the log-likelihood is
-I
Can you please tell us (me!) how you chose starting values?
Out of curiosity I tried the following:
set.seed(42)
x <- rlnorm(100,1,2) + 3
require(MASS)
strt <- list(mu=1,sigma=2,gamma=3)
fit <- fitdistr(x,densfun=function(x,mu,sigma,gamma)
{dlnorm(x-gamma,mu,s
x27;t. Nor should you. It would be unethical. Think before you
post, because once you've posted, that's it. "The moving finger writes,
and having writ/Moves on, nor all thy piety nor wit/Shall lure it back
to cancel half a line/Nor all thy tears wash out a word of it.&quo
view on CRAN.
There *may* be ways to avoid the problematic specification of search
intervals to be found in this task view Good luck to the OP!
cheers,
Rolf Turner
On Wed, Jan 8, 2014 at 8:55 PM, Aurélien Philippot
wrote:
Dear R experts,
I want to use numerical methods to solve a
s
exactly the same.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
See in-line below.
cheers,
Rolf Turner
On 04/01/14 12:48, Stefan Mühlbauer wrote:
Hello!
Â
I am having two variables x and y (whereas y is a set of raster images) and
want to
quantify the correlation between x and y by calculating the Pearson
Correlation Coefficient. In
creating S; then do set.seed(123)
again before creating T. If you do so, S and T will be
identical.
cheers,
Rolf Turner
P. S. "T" is not a good name for an object; too easy to confuse
with "TRUE". Not an egregious sin, but to be avoided.
R. T.
On 0
On 25/12/13 06:23, Jeff Newmiller wrote:
In short, if you ask us why one divided by zero doesn't give 3, we have to
wonder if you don't belong in some other educational forum.
Fortune!
cheers,
Rolf Turner
_
As far as I can discern, your question makes no sense at all.
Suppose you *know* that y = 2 + 3*x1 + 4*x2.
Now what should you predict when x1 = 6 (with x2 "missing"/unknown)?
See fortune("magic").
On 19/12/13 07:18, Chris Wilkinson wrote:
I would like to predict a new response from a fitte
On 16/12/13 20:59, Jeff Newmiller wrote:
This is not a statistics theory forum, but posting a solution to this
nonsensical problem would be irresponsible.
Fortune?
cheers,
Rolf
__
R-help@r-project.org mailing list
https://stat.et
On 16/12/13 12:23, Duncan Murdoch wrote:
[After a number of other "Don'ts"]
Don't trust the preview to tell you the quality of the graph, try
printing the document.
Word isn't quite as bad as it appears.
Don't use Word.
Fortune?
See inline below.
On 12/11/13 11:28, Bert Gunter wrote:
This is not really an R question -- it is statistics.
In any case, you should do better posting this on the
R-Sig-Mixed-Models list, which concerns itself with matters like this.
However, I'll hazard a guess at an answer: maybe. (Vague q
Uh, no. You are forgetting to take the square root of 10, and to divide
by the square root of 12.
The variance of Y is (exactly) (56^2 - 1)/12, so the variance of Y-bar
is this quantity over 10,
so the standard deviation of Y-bar is sqrt((56^2 - 1)/12)/sqrt(10).
Which is approximately
(ig
pe(amacrine,fun=foo,simulate=expression(rmh(fit,verbose=FALSE)),i="on",j="off",spar=0.7)
plot(E)
I have cc-ed this message to Prof. Baddeley who will doubtless correct
any errors
or misconceptions in the forgoing.
cheers,
Rolf Turner
On 12/05/13 20:35, DingYi wro
?RSiteSearch
cheers,
Rolf Turner
On 12/05/13 03:22, Iman Chambari wrote:
hi all,
i have a question about R
this code for me not work: predict.arfima(x)
error: could not find function predict.arfima
i test other same code that error again : pred.arfima predict.arima
pred.arima
i
701 - 800 of 2101 matches
Mail list logo
