Displaying 16 results from an estimated 16 matches for "eruptions".

2005 Oct 21

1

finite mixture model (2-component gaussian): plotting component gaussian components?

...trate. (Note: my real dataset will not have peaks this well separated, but I needed to find a small example.)
##################
# (1) get some data
##################
data(faithful)
##################
# (2) get model
##################
library(mclust)
MyMixtureModel<-summary(EMclust(faithful$eruptions),faithful$eruptions)
##################
# (3) plot mixture model
##################
attach(MyMixtureModel)
mclust1Dplot(data=faithful$eruptions,z=z,mu=mu,sigmasq=sigmasq,pro=pro,ask=FALSE,type=c("density"))
do.call("mclust1Dplot",c(list(data=faithful$eruptions,ask=FALSE,type=c(...

2006 Jul 20

2

function names in a vector used by for (){} character problem ?

...y.default$kernels, but thats not my question
here. my own kernel functions are running and working.
my problem is the following is not working:
dev.off()
par(mfrow=c(3,3))
for(i in 1:length(bw))
{
for(j in 1:length(kernels))
{
par(mfg = c(i, j))
nawaline2(eruptions,waiting,kernels[j],bw[i],1000)
# FYI: the following worked but ofcourse all rows were the same
# nawaline2(eruptions,waiting,gauss,bw[i],1000)
}
}
# here are the standard arguments for nawaline, kern exspects a
function like gaussian, epan or rectangular.
nawaline2(xi,yi,kern,h,N=1...

2006 Apr 24

1

Rumbling Java volcano may erupt

http://www.physorg.com/news65033298.html
Switch to Rails before it''s too late, folks!
Jeff
--
Posted via http://www.ruby-forum.com/.

2005 Oct 22

0

package mclust: cdens, EMclust?

...------------------
SAMPLE CODE
--------------------
Here is a little test script to illustrate. (Note: my
real dataset will not have peaks this well separated,
but I needed to find a small example.)
##################
data(faithful)
library(mclust)
MyMixtureModel<-summary(EMclust(faithful$eruptions),faithful$eruptions)
attach(MyMixtureModel)
mclust1Dplot(data=faithful$eruptions,z=z,mu=mu,sigmasq=sigmasq,pro=pro,ask=FALSE,type=c("density"))
do.call("mclust1Dplot",c(list(data=faithful$eruptions,ask=FALSE,type=c("density")),MyMixtureModel))
# plot components...

2006 Aug 25

3

plot question

Hi everyone,
I have what may appear to be a newbie question, but I have looked
everywhere I can think to look and I cannot find an answer. On page 35
of "An Introduction to R" the following command appears:
plot(ecdf(eruptions), do.points=FALSE, verticals=TRUE). What is the
do.points argument? I know what it does (suppresses printing of the
points) but where can I find help on it? I want to be able to explain it
fully to my students.
Thanks for your help,
Cathy
--
Dr. Cathy Carter
Department of Geography
University...

2017 Oct 09

1

example of geom_contour() with function argument

...ataset in the
> example). Can you please point to the specific example that might help?
>
> Here is what I get:
>
> require(mvtnorm)
> require(ggplot2)
> set.seed(1234)
> xx <- data.frame(rmvt(100, df = c(13, 13)))
>
>
> v <- ggplot(faithfuld, aes(waiting, eruptions, z = drmvt, df = c(13,13)))
> v + geom_contour()
>
> Don't know how to automatically pick scale for object of type function.
> Defaulting to continuous.
> Error: Aesthetics must be either length 1 or the same as the data (5625):
> x, y, z, df
>
That's not what I get:...

2017 Oct 09

1

example of geom_contour() with function argument

...int to the specific example that might help?
> >
> > Here is what I get:
> >
> > require(mvtnorm)
> > require(ggplot2)
> > set.seed(1234)
> > xx <- data.frame(rmvt(100, df = c(13, 13)))
> >
> >
> > v <- ggplot(faithfuld, aes(waiting, eruptions, z = drmvt, df = c(13,13)))
> > v + geom_contour()
> >
> > Don't know how to automatically pick scale for object of type function.
> > Defaulting to continuous.
> > Error: Aesthetics must be either length 1 or the same as the data (5625):
> > x, y, z, df
>...

2017 Oct 09

2

example of geom_contour() with function argument

...elp?
>> >
>> > Here is what I get:
>> >
>> > require(mvtnorm)
>> > require(ggplot2)
>> > set.seed(1234)
>> > xx <- data.frame(rmvt(100, df = c(13, 13)))
>> >
>> >
>> > v <- ggplot(faithfuld, aes(waiting, eruptions, z = drmvt, df = c(13,13)))
>> > v + geom_contour()
>> >
>> > Don't know how to automatically pick scale for object of type function.
>> > Defaulting to continuous.
>> > Error: Aesthetics must be either length 1 or the same as the data (5625):
>&g...

2017 Oct 09

1

example of geom_contour() with function argument

...n which calculates a KDE (that is, something from the dataset in the
example). Can you please point to the specific example that might help?
Here is what I get:
require(mvtnorm)
require(ggplot2)
set.seed(1234)
xx <- data.frame(rmvt(100, df = c(13, 13)))
v <- ggplot(faithfuld, aes(waiting, eruptions, z = drmvt, df = c(13,13)))
v + geom_contour()
Don't know how to automatically pick scale for object of type function.
Defaulting to continuous.
Error: Aesthetics must be either length 1 or the same as the data (5625):
x, y, z, df
Can you please tell me how to use this here? Or is some other...

2009 Jun 14

5

Leave One Out Cross Validation

Hi All,
I have been trying to get this LOO-Cross Validation method to work on R for
the past 3 weeks but have had no luck. I am hoping someone would kindly help
me.
Essentially I want to code the LOO-Cross Validation for the ''Local Constant''
and ''Local Linear'' constant estimators. I want to find optimal h, bandwidth.
Thank you very much!
-M
--
View this

2000 Feb 29

1

No subject

Does anyone know of any comprehensive literature (like a book, collection of
articles, etc.) that walk thorugh example of ARIMA modeling with the ts
package of R? Genarlly if you've just started to learn the R-language and
are all week long with ARIMA, what is a fast-pace source that you would
recommend (to a beginner) ?
Thank you,
Krasi

2010 Apr 07

0

exporting rotating plot into Powepoint

Hi,
I need to export a rotating plot like the one generated from the following function into a PPT slide or something that can be viewed by reviewers of a journal:
library(np)
library(datasets)
data(faithful)
f <- npudens(~eruptions+waiting,data=faithful)
plot(f,view="rotate",neval=100,phi=30,main="",xtrim=-0.2)
Thanks!
Gabriele Zoppoli, MD
Ph.D. Fellow, Experimental and Clinical Oncology and Hematology, University of Genova, Genova, Italy
Guest Researcher, LMP, NCI, NIH, Bethesda MD
Work: 301-451-85...

2011 Apr 19

0

Error message in package:bayesSurv. Why?

Dear folks,
I have been struggling to create what I fondly imagined would be a
straightforward adaptation of the package''s example to my own dataset, which
looks at incidence of depression following the birth of initial and up to 3
subsequent children (4 children in all, with all subjects having a first
child). I''ve regarded these as equivalent to tooth eruption and subsequent

2000 Feb 23

1

Large datasets under R

Hello,
I recall reading a thread months ago on this mailing list about handling
very datasets under R, but I can't seem to find it. This has become
particularly important recently, because I've been playing with a dataset
containing information about every fatal car accident in the U.S. since
1975; in total, the relevant files are about 120 megs. I'd like to load
all of these into R

2017 Oct 09

1

example of geom_contour() with function argument

Hi BFD,
?geom_contour() *does* have helpful examples. Your Google-foo is weak:
Searching for geom_contour brought me:
http://ggplot2.tidyverse.org/reference/geom_contour.html as the first
result.
HTH
Ulrik
On Mon, 9 Oct 2017 at 08:04 Big Floppy Dog <bigfloppydog at gmail.com> wrote:
> Can someone please point me to an example with geom_contour() that uses a
>

2000 Apr 30

5

Help Need with aov()

Hi there,
I'm using R1.0.1 Windows 98.
This file contains some inputs and an aov function code. Can someone
check it for me? Somehow I got completely different answer when typing
them in R and in Splus.
Splus gives me this:
> summary( Turnip.aov )
Error: Blocks
Df Sum of Sq Mean Sq F Value Pr(F)
Residuals 3 163.7367 54.57891
Error: Plots %in% Blocks