() in the spatstat package
(after appropriately reconfiguring your polygons as objects of class
owin; something like w - owin(poly=plygn) where plygn is your
polygon).
You could also use the undocumented function inside.xypolygon().
cheers,
Rolf Turner
In respect of fitting piecewise linear regressions, have you looked at
the segmented package?
cheers,
Rolf Turner
On 18/01/12 04:30, crimsonengineer87 wrote:
Dear Forum,
I have been wracking my head over this problem for the past few days. I have
a dataset of (x,y). I have been
)
points(t,var)
indicates a fairly reasonable fit. One wonders however about the
validity of the model
considering the lack of fit at t=0.
R has great and easy-to-use graphics capabilities. Use them!
cheers,
Rolf Turner
__
R-help@r
columns of dataset have extra
attributes tagging
along. E.g. the column could actually be a numeric matrix of zeroes ---
in which case
it wouldn't get dropped.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch
- list()
j - 0
for(i in 2:length(xxx)) {
tmp - read.csv(textConnection(xxx[i]),header=FALSE)
if(ncol(tmp)==ncol(yyy)) yyy - rbind(yyy,tmp) else {
j - j+1
bad[[j]] - tmp
}
}
closeAllConnections()
HTH
cheers,
Rolf Turner
1.000951 with absolute error 0.00011, rather than giving
exactly 1, so there's a bit of slop in the system.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
have no idea why nls() is throwing an error.
* Your parameterisation is bad.
* A better parameterisation can be readily fitted to your data
using optimize().
* The model is probably too complicated and inappropriate for these
data.
HTH
cheers,
Rolf Turner
On 01/02/12
,
Rolf Turner
On 07/02/12 08:54, Berend Hasselman wrote:
On 06-02-2012, at 20:22, Abhishek Pratap wrote:
Hey Guys
I want to divide(numerically) all the columns of a data frame by
different numbers. Here is what I am doing but getting a weird error.
I don't think you are getting
No. This does not do anything like what the OP wanted.
Did you ***try*** it, for crying out loud?
cheers,
Rolf Turner
On 07/02/12 08:29, ilai wrote:
duplicated(Data..)
On Mon, Feb 6, 2012 at 11:34 AM, LCOG1jr...@lcog.org wrote:
Hi all
For the data below, I would like
* brackets (``[[ ]]'') when assigning
a value to the i-th entry of your pre-defined list.
(4) For efficiency it is probably best to use the estimates returned
at step i as the starting values for step i+1.
cheers,
Rolf Turner
__
R-help@r
)
V - tiles(tess(image=U))
lapply(V,centroid.owin)
lapply(V,area.owin)
You should read the help on the functions invoked above. It might be
illuminating to do some plotting along the way (of X, Y, Z, U, and V).
cheers,
Rolf Turner
It has been pointed out to me (in a private email) that there was
a typo in my previous post;
X - rpoisspp(100)
should read
X - rpoispp(100)
(i.e. just one s). Sorry 'bout that! :-)
cheers,
Rolf
__
R-help@r-project.org
and ones.)
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
(3) brute force:
x - readLines(filename)
n - length(x)
Having determined n, do:
y - read.csv(filename,nrows=n-2)
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r
David:
Why do you say the OP wants a 503rd or 504th degree polynomial?
He/she wants an interpolating spline.
But that is precisely what spline() or splinefun() would give him/her.
So I don't understand the problem either!!!
To the OP:
What did you actually *do* to create your spline?
and effort into providing helpful
advice on this list, but is understandably put off by silly questions or
failures to follow the posting guide.
His efforts should be appreciated and applauded, not whined about.
cheers,
Rolf Turner
__
R-help@r
--- and it's pretty clear that at this stage of
your development, you don't.
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
$post.f.crwn.length
However, the height of the two curves are about 1/3 to 1/4 the height that
they should be compared to the histogram. Any ideas?
Yes. Read the help on hist! (Hint: Pay particular attention to the
freq and/or probability arguments.)
cheers,
Rolf Turner
On 14/05/11 02:28, Duncan Murdoch wrote:
SNIP
One way to end up with a corrupted workspace is to save it with a
dependency on a particular package, then try to load it when that
package is unavailable. I'm going to see if R can handle that case
more gracefully.
SNIP
Pardon me for
On 14/05/11 14:09, David Winsemius wrote:
On May 13, 2011, at 9:09 PM, Rolf Turner wrote:
On 14/05/11 02:28, Duncan Murdoch wrote:
SNIP
One way to end up with a corrupted workspace is to save it with a
dependency on a particular package, then try to load it when that
package is unavailable
On 15/05/11 13:41, Dan Abner wrote:
Hello everyone,
Is there an R function that returns an object's search path position?
?find
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
.
If that doesn't work, please provide a minimal *reproducible* (no one
but you
has the ``rwb'' data object) example of the problem that you are having
(as the posting guide requests).
cheers,
Rolf Turner
On 16/05/11 17:01, Benjamin Caldwell wrote:
Hmm; still missing something - hist
is the sky blue? Why are fire engines red? Why don't they teach
linear algebra these days?]
Because if A%*%x = lambda*x then
A%*%(-x) = lambda*(-x).
cheers,
Rolf Turner
P. S. Why is a duck? Answer: Because the higher you go, the
vest has no sleeves
value that you really want.
You might try (the appropriate analogue of) x +
get(annotIndex,envir=.GlobalEnv).
That should satisfy the package checker. (I haven't tested it, but.)
Dunno if this has an adverse impact on the efficiency of your code.
HTH
cheers,
Rolf Turner
checker! :-)
Might be more efficient than my suggestion of using ``get()''.
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
Please ignore this message.
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
.
Could someone suggest how to do this? I have come close to what I need with
loops and such, but there must be a less clumsy way...
L.new - L[tapply(L,nrow) 1]
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch
On 19/05/11 10:26, bill.venab...@csiro.au wrote:
SNIP
Most of [the Google style guide's] advice is very good (meaning I agree with
it!) but some is a bit too much (for example, the blanket advice never to use
S4 classes and methods - that's just resisting progress, in my view).
SNIP
I must
Sorry for the noise.
cheers,
Rolf
__
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,
(1) What has this to do with recursion?
(2) You probably need to use ifelse(). I believe that this is
(in effect) an FAQ.
cheers,
Rolf Turner
On 20/05/11 07:42, Tremblay, Pierre-Olivier wrote:
Hi,
I created a function for obtaining the normal cumulative distribution (I
know
a *vector* (of length 3 in your setting), the last
entry of which is the error (or residual) mean square, which is
probably what you want since you refer the ``mean square value''
(singular).
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
suggests.
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
the command
line will get the lm() in whichever package (teaser or stats) is closer
to the start of the search path. And there's really nothing much that
can be done about this. Have I got that right?
cheers,
Rolf Turner
On 20/05/11 18:26, Jari Oksanen wrote:
Duncan
On reflection, it seems to me that what we really need here
is Prof. Ripley's mind_read() function that was foreshadowed
some years ago (see fortune(mind_read)) so that R could
determine just *which* lm() (for example) function the user
has in mind when he or she types ``lm(...)'' at the
and not Type III
SNIP
You ***really should*** have a look at
www.stats.ox.ac.uk/pub/MASS3/Exegeses.pdf
cheers,
Rolf Turner
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r
to apologize. Your English is better than my Deutsch.
(The ``.at'' address *does* mean ``Austria'', nicht wahr?)
HTH
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
.
cheers,
Rolf Turner
[[alternative HTML version deleted]]
__
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
I think cognizance should be taken of fortune(very uneasy).
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
of the appropriate terms
to search on.
Can anyone remember/point out where this discussion can be found?
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
.
However if you want to use a sledge-hammer to crack a peanut,
install the spatstat package and then do:
require(spatstat)
W - owin(poly=your ellipse)
inside.owin(x1,x2,W) # Where x1 and x2 are the x and y coordinates of
the points you wish to test.
cheers,
Rolf Turner
in the future dependant on the conditions
during the observation.
Any advice would be much appretiated.
Have a look at the CRAN packages:
* msm
* HiddenMarkov
* hmm.discnp
cheers,
Rolf Turner
__
R-help@r-project.org
(1) Use an informative subject line.
(2) Do not use html mail.
Then you might get a response.
cheers,
Rolf Turner
On 01/06/11 02:12, mustafabinar wrote:
Hi R people.
nbsp;
I have a problem. What can I create by using functions the
combinations ofnbsp;distances in multiple
any sense either. You
present the values of only one variable. For a regression you
need to have a y-variable and at least one x-variable. It would
appear that you're not thinking very clearly.
cheers,
Rolf Turner
On 01/06/11 11:32, J S wrote:
Dear forum members,
How can I force
to specify height and
width.
cheers,
Rolf Turner
On 02/06/11 07:17, Michael Friendly wrote:
[Env: R 2.12.2, Win XP]
I'm creating figures using MASS::eqcsplot to provide equal scaling of
the axes. My figures work OK
when I plot to the screen, but when I try to do the same
().
cheers,
Rolf Turner
On 03/06/11 03:04, Graves, Gregory wrote:
I have 2 datafiles 'target' and 'observed' as shown below (I will gladly
email these 2 small files to whomever). X25. And X75. Indicate the
value of 25th and 75th-percentile of the target ('what should
think that
plot(1,main=paste(,
))
might do what you want.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
--- but ***not***
in the output! Your advice is the antithesis of what is required.
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 is *not*
the source of the OP's problems.
Sarah Goslee has already pointed out what at least one of the sources
of his problems is.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman
On 03/06/11 12:16, Joshua Wiley wrote:
Okay, that might have been a little strong. screams bloody murder
is a warning, not technically an error, and does not occur when simply
running R CMD build. That said, the OP did mention using R CMD check
and pdflatex is not an issue when only building
around with
calls to .Internal() !!!) and is pretty straightforward, so the hack won't
be too hard to implement.
Good luck.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
in later graphics.
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
the object. I didn't find anything like this in the R docs.
?comment
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
On 04/06/11 00:02, Barry Rowlingson wrote:
SNIP
. someone on
StackOverflow was wondering why 0= x= 1 isn't a valid
expression, neither in the context of what he expected (TRUE if x is
between 0 and 1) nor how I explained it would function (as (0=x)=
1, and then comparing a TRUE/FALSE
(13,7,9,3,6,6,3,1,11,5,15,5)
p372- data.frame(trt,blk,res)
fit- lm(res ~ trt + blk,data=p372) # DON'T use attach()! Use the data
argument for lm().
summary(fit)
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch
).
cheers,
Rolf Turner
On 09/06/11 16:14, Dennis Murphy wrote:
Hi:
You can try something like this: assuming the factor variables of
interest and the response variable are in a data frame named df,
ivset- c(comma separated vector of factor names)
myaovs- lapply(ivset, function(x
: Histogram
I think the command you want is barplot
x = rbinom(10,15,0.65)
y = rbinom(10,15,0.25)
barplot(rbind(x,y),beside=TRUE)
RTFM. I.e. execute
?rbinom
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https
gone in circles on wildcard, lapply, and names threads and
haven't managed to get anything to work.
I think that
outmap - do.call(merge,myobjects)
should work for you.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https
-code-style-guide/
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
learn about factors; they are
important and useful.
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
(). As has been discussed recently on
this list, the syntax for curve() is a bit delicate.
A workaround for your problem is:
plot(function(t){B1(T1,t)},0,3650)
HTH
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https://stat.ethz.ch
(1) Repeatedly sending the same message to this list is not a good idea.
(2) It would probably be most efficacious for you to contact the maintainer
of the mhsmm package. Type
maintainer(mhsmm)
to find out who this is!
cheers,
Rolf Turner
Extensions -- Writing R Documentation Files -- Mathematics
You need eqn{} and/or deqn{}.
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
intervals, is to be found in
Theoretical comparison of bootstrap confidence intervals,
Peter Hall, Annals of Statistics vol. 16 no. 3 1988, pp. 927 -- 953.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
https
(1) You really ought to do your own homework.
(2) What has this to do with R?
cheers,
Rolf Turner
On 29/06/11 14:03, Lisa wrote:
Dear all,
I just want to get the derivative of a function that looks like:
y = exp(x1*b) / (exp(x1*b) + exp(x2*b))
where y is a scalar, x1, x2
does not do automatically.
cheers,
Rolf Turner
D(expression(exp(x1*b) / (exp(x1*b) + exp(x2*b))), b)
exp(x1 * b) * x1/(exp(x1 * b) + exp(x2 * b)) - exp(x1 * b) *
(exp(x1 * b) * x1 + exp(x2 * b) * x2)/(exp(x1 * b) + exp(x2 *
b))^2
See ?D and also note deriv on the same help
The simplest way is:
xxx - with(clyde,tapply(CONTTIME,SCRNO,sum))
You could also do:
xxx - by(clyde,clyde[[SCRNO]],function(x){sum(x[[CONTTIME]])})
but this gives somewhat messy output; the aforesaid output may be
convenient for some purposes, not for others.
cheers,
Rolf
Homework? We don't do people's homework for them.
cheers,
Rolf Turner
On 01/08/2013 05:13 AM, Agnes Ayang wrote:
Hello R-helpers,
I need to generate standard variates normal to 'create' chi-squared variates.
To make you more understand,
(1) a-rnorm(3,0,1)
*after do (1), I
,
Rolf Turner
On 01/10/2013 02:33 PM, John Sorkin wrote:
windows 7, R 2.12
I am trying to run a piecewise linear regression with a single knot, i.e. a regression composed of two straight lines where the two lines intersect at an x value given by the variable knot. I wish to estimate the slope
available.
I also tried
sudo yum install texlive-fonts-extra
with a similar result.
Can anyone give me a simple recipe as to how to get the inconsolata font
and the associated inconsolata.sty?
Thanks.
cheers,
Rolf Turner
__
R-help@r
, and when I did
sudo yum install libX11-devel.x86_64
I was told that this package was also already installed.
What do I do to make the headers/libs available please?
cheers,
Rolf Turner
P. S. Is there a way to get a binary version of R-patched for Fedora 17,
e.g. by using yum
. Also associated functions like tapply() and sapply(). Very
useful gadgets. And if you really get hooked on such functionality, you
might want to have a look at the plyr package.
cheers,
Rolf Turner
__
R-help@r-project.org mailing list
On 01/12/2013 12:41 AM, Prof Brian Ripley wrote:
On 11/01/2013 11:28, Milan Bouchet-Valat wrote:
Le vendredi 11 janvier 2013 à 22:45 +1300, Rolf Turner a écrit :
I am trying to build R-patched from source on a (newly installed)
Fedora 17 system on a new laptop.
When I do the usual ./configure
#TeX_Live_2012_2
I haven't actually *tried* that, given that I had got things working and
didn't dare do anything that might put my system back into bozo mode! :-)
Again, thanks to all who proffered advice.
cheers,
Rolf Turner
__
R-help@r
methods base
other attached packages:
[1] spatstat_1.30-0 deldir_0.0-21 mgcv_1.7-22 misc_0.0-15
loaded via a namespace (and not attached):
[1] grid_2.15.2 lattice_0.20-10 Matrix_1.0-10 nlme_3.1-105
[5] tools_2.15.2
cheers,
Rolf Turner
On 01/13/2013 04:36 AM, Marc Schwartz wrote:
On Jan 11, 2013, at 10:35 PM, Rolf Turner rolf.tur...@xtra.co.nz wrote:
Tried to do a plot just now and got an error:
Error in text.default(2, 6, main, cex = cex) :
X11 font -adobe-helvetica-%s-%s-*-*-%d-*-*-*-*-*-*-*, face 1 at size 16 could
mean when A is at level 1 and B is at level 1, and so
on.
Suppose mu_11 = 1, mu_12 = -1, mu_21 = -1, and mu_22 = 1.
Then there are no main effects; A averages to 0, as does B.
But there is an elephant-ful of interaction.
cheers,
Rolf Turner
cheers,
Rolf Turner
On 01/13
Hi Erin,
What you're missing is that the TukeyHSD function is generic. There seems
to be only one method for it, namely TukeyHSD.aov. Try typing:
methods(TukeyHSD)
On my system if I type
TukeyHSD.aov
I see a splendiferous amount of lovely --- and (mirabile dictu!) reasonably
=confidence,level=0.90)
plot(x,y)
lines(x,pfit[,fit])
lines(x,pfit[,upr],col=red)
You are then, for any given x value, 95% confident that the true mean of
Y lies *below* the
corresponding y-value on red curve that was plotted.
cheers,
Rolf Turner
On 01/15/2013 09:17 AM, li li wrote:
Hi all,
For the simple linear regression, I want to find the input x value so
that the
lower confidnece limit is a specific number, say 0.2.
In other words, I want to find the value of x so that the lower
confidence bound crosses the horizontal line
On 01/15/2013 08:04 AM, Duncan Murdoch wrote:
SNIP
I've never understood why Microsoft thinks it is good to display
incomplete filenames. This practice causes so many problems.
SNIP
Personally, I never understand Why Microsoft?.
cheers,
Rolf Turner
, but
all of these produce warnings or errors.
The \usage{} section appears to be treated verbatim or as if it were
preformated.
So just put the line break literally in your *.Rd and it will
automagically appear in the
resulting help.
cheers,
Rolf Turner
The help facility is applicable to functions and data sets. It is not
designed
or intended to give help with respect to R syntax (with the exception of
the basic syntax of the operators --- unary and binary --- and the
associated
rules of precedence).
cheers,
Rolf Turner
There does not appear to be any such function as barplot2 in
the current version (3.4-5) of the plotrix package. Moreover
I can find no reference to such a function in the NEWS for
plotrix.
cheers,
Rolf Turner
On 01/23/2013 07:28 AM, Martin Batholdy wrote:
Hi,
is there any way
Given that your labels are no and yes, what do you expect R to
do? To quote a well-known fortune, R is lacking a mind_read() function!
cheers,
Rolf Turner
On 01/23/2013 10:58 PM, Francesco Sarracino wrote:
Thanks,
this works! but I am surprised that R has such a strange behavior
as to whether the result returned is sensible.
It is also possible that the fix I have suggested will induce other
problems
in other contexts; I don't really understand what's going on well enough to
be certain.
I hope that this is of some help, but.
cheers,
Rolf Turner
On 01/12/2010
,col=red)
cheers,
Rolf Turner
On 01/27/2013 10:12 AM, Giuseppe Amatulli wrote:
Hi,
I'm using the R library(car) to draw confidence/prediction ellipses in a
scatterplot.
From what i understood the ellipse() function return an ellipse based
parameters: shape, center, radius .
If i
I just saw a message from David Winsemius, responding to an inquiry
from Carol White:
On Jan 28, 2013, at 9:06 PM, carol white wrote:
Should I understand that this message was received?
It's always possible to check the Archives for this question.
This prompted me to ask about a problem
I was going to attempt to answer your question, but I refuse to
respond to anyone who writes wanna when they mean want to.
The use of wanna is unacceptable in written English unless you
are *trying* to be funny, and I don't think you are.
cheers,
Rolf Turner
On 01/30/2013 09:48 AM
,
Rolf Turner
On 01/29/2013 11:26 PM, Larissa Modica wrote:
Hello everybody,
I am sorry if my questions are too simple or not easily understandable. I’m
not a native English speaker and this is my first analysis using this
function.
I have a problem with a cross correlation function and I
On 01/30/2013 10:28 AM, Sarah Goslee wrote:
The searchable archives may lag, and apparently do. The main list
archive is here:
https://stat.ethz.ch/pipermail/r-help/
and is complete. That's the one to check if you wish to know whether
something made it to the list.
If you go to the r-help
that it's difficult to say what your problem is without having
access to your actual data.
cheers,
Rolf Turner
On 01/30/2013 09:29 AM, Abby Rudolph wrote:
Hello,
I am using the following code to create ppp files from csv data and map shape
files, but I am getting some errors which I
Just look at the code of predict.lm(). It is reasonably perspicuous.
In particular look at res.var.
cheers,
Rolf Turner
On 01/31/2013 05:50 AM, Kurt Rinehart wrote:
For a given linear regression, I wish to find the 2-tailed t-dist
probability that Y-hat = newly observed values
more detail.
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
This smacks of homework to me.
cheers,
Rolf Turner
On 02/02/2013 12:45 PM, Jin Choi wrote:
I am trying to figure out how to calculate p-values for the difference in
prevalence of a risk factor between men and women. For example, I find that
277 out of 710 male patients and 125
,
., Y-hat_{N+10} = Y-hat_{N+9} + Z-hat_10.
In R, let your forecast values be the vector Zhat (a vector of length 10).
Then do:
Yhat - cumsum(c(Y[N],Zhat))[-1]
Xhat - exp(Yhat)
Get error bounds on the forecasts is more problematic.
cheers,
Rolf Turner
On 02/05/2013 11:49 PM
, minus infinity.
cheers,
Rolf Turner
On 02/05/2013 04:59 AM, francesca casalino wrote:
I am sorry I have confused you, the logs are all base e:
ln(a) = 1347
ln(b) = 1351
And I am trying to solve this expression:
exp( ln(a) ) - exp( ln(0.1) + ln(b) )
Thank you.
2013/2/4 francesca
Try:
n - length(mlist)
result - lapply(1:n,function(i,m,v){m[[i]]%*%v[[i]]},m=mlist,v=vlist)
These days the use of lapply is unlikely to be much, if at all, faster than
the use of a for loop.
cheers,
Rolf Turner
On 02/06/2013 06:50 PM, David Romano wrote:
Hi everyone,
I'd like
.
HTH
cheers,
Rolf Turner
P. S. Does anyone know why or whither Vito Muggeo has disappeared?
R. T.
On 02/06/2013 11:30 AM, oli tills wrote:
Hi,
I would like to calculate the area under segmented regression lines (single
breakpoints). I had thought that I could do
Perhaps you want to look at
?plotmath
cheers,
Rolf Turner
On 02/07/2013 11:41 AM, David Arnold wrote:
Hi,
I'd like to add alpha in latex code to my image. Any suggestions?
library(UsingR)
BagA = c(rep(10,6),rep(20,5), rep(30,4),rep(40,3),rep(50,2),60,70)
BagA
BagB = c(10
pattern
object; too easily confused with the name of the function points().
cheers,
Rolf Turner
On 02/08/2013 02:57 AM, Pavel_K wrote:
Hello,
I am writing cause I have a problem when try to create confidence envelopes
in spatstat.
I have a point data representing firms in my study
901 - 1000 of 2038 matches
Mail list logo