search for: eruptions

...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(...

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

Both the document and R itself are at 2.6.2 (2008-02-08). When I run in the Windows R console the command lines(density(eruptions, bw=0.1)) found in ? 8.2, I get: Error in plot.xy(xy.coords(x, y), type = type, ...) : plot.new has not been called yet Regards, Olivier Lefevre

Hi, I am trying to create a Cumulative Frequency graph and I am using the following example: http://www.r-tutor.com/elementary-statistics/quantitative-data/cumulative-frequency-graph When I plot this data, how do I put in "real values" on the x-axis, rather than the values that are used for the breaks Thanks -- Shane [[alternative HTML version deleted]]

...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...

...------------------ 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...

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...

...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:...

...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...

...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 >...

...", na.rm = FALSE, contour = TRUE, n = 100, ...) de esta función requiero ver el cálculo númerico de la estimación de la densidad, cuento con el siguiente codigo y sólo puedo ver el grafico, pero no, los valores de la estimacion de la densidad (valores numericos). p <-ggplot(faithful,aes(x=eruptions,y=waiting)) p +stat_density2d(aes(fill=..density..),geom="raster",contour=FALSE) las lineas anteriores muestran el grafico, pero no sé, de que forma puedo ver los valores numericos de la estimacion de la densidad De que forma podria ver el valor numerico de la estimacion de densidad?...

...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...

I'm trying to shade part of a density plot. The code I'm trying (using the Old Faithful data as an example) is something like this: # The Old Faithful geyser data data(faithful) d <- density(faithful$eruptions, bw = "sj") plot(d) polygon(d[d$x>4], col = "wheat") I expected that the part of the curve to the right of 4 on the x axis should be shaded, but nothing gets shaded at all. In fact, when I just try to print d[d$x>4] I get many lines of $"NA" NULL I must be mist...

...ll, In many cases, I need a plotting region much bigger than the screen (e.g. for maps or for graphs with many labels). A. MS-Windows if I try windows(width=25, height=25, rescale="fixed") it seems to be OK (a screen with scrollbars, exactly what I need) but if I try then plot(faithful$eruptions, faithful$waiting) I receive Error in plot.new() : Outer margins too large (fig.region too small) B. Linux I try X11(width=25, height=25) (and the same) No Error message, but no scrollbars !! Is there a solution ? with many thanks in advance Alain Guerreau directeur de recherche au CNRS...

I'm trying to plot a series of densities using/comparing differing weights. I see the reference to weights and subscripts, but I don't understand how to implement that. My data are of the form: I, J, Actual, Distance, Subset, Weight1, Weight2, ... I'm trying to see the effect of the distance distribution (Actual by Distance) compared to the various weighted distributions

...ally Simple Syndication feeds. including renewable energy research and development. We all want a higher standard of living for ourselves and our children. Aceasta este vestea rea. Cross off those plans to tool to the shore in a convertible. They think it is too hard to change. In fact, spectacular eruptions like that of Mount Pinatubo are demonstrated to contribute to global cooling through the injection of solar energy reflecting ash and other small particles. It is important for us to have a concrete, shared understanding of where we want to go. Educator Mary Jean LeTendre said it best when she obse...

I'm using Asterisk 1.4.22.1 When I'm on active call it happens many times the call gets interrupted by music-on-hold without my pressing any button. MOH just kicks in and int erupt the call and I have no way of getting the call back. Did anybody experienced anything like this? -- Joseph

What is a good way to enter a very long model formula. For example: y ~ Input.2 + Input.3 + ... + Input.1000 (assuming the corresponding dataframe has many other columns). Is there a way to convert a character string to a formula? Are there command line expansions in R besides the simple '.'? Thanks. [[alternative HTML version deleted]]

...t) [1] 35.65342 > mean(~age, data=HELPrct) [1] 35.65342 > mean(age ~ sex, data=HELPrct) female male 36.25234 35.46821 > mean(age ~ sex & treat, data=HELPrct) female.no male.no female.yes male.yes 37.56364 35.90173 34.86538 35.03468 > interval(binom.test( ~ eruptions > 3, faithful)) probability of success lower upper 0.6433824 0.5832982 0.7003038 > pval(binom.test( ~ eruptions > 3, faithful)) p.value 2.608528e-06 > xchisq.test(phs) # physicians health study example (da...

...t) [1] 35.65342 > mean(~age, data=HELPrct) [1] 35.65342 > mean(age ~ sex, data=HELPrct) female male 36.25234 35.46821 > mean(age ~ sex & treat, data=HELPrct) female.no male.no female.yes male.yes 37.56364 35.90173 34.86538 35.03468 > interval(binom.test( ~ eruptions > 3, faithful)) probability of success lower upper 0.6433824 0.5832982 0.7003038 > pval(binom.test( ~ eruptions > 3, faithful)) p.value 2.608528e-06 > xchisq.test(phs) # physicians health study example (da...