search for: xmin

...just build up a piecewise linear model to fit some data, including some interaction and i am not sure of how to interpret the summary:. here it is: -------------------------------------------------------------------------------- Call: lm(formula = weightedDiff ~ angleNoise + (reflection < Break[xMin]) * reflection + (reflection >= Break[xMin]) * reflection + angleNoise:(reflection < Break[xMin]) * reflection + angleNoise:(reflection >= Break[xMin]) * reflection) Residuals: Min 1Q Median 3Q Max -1.073e-03 -1.749e-04 -5.913e-06 1.650e-04 1...

Hi, I want to plot grouped bars to compare 95% confidence interval estimates from two models. Each bar represents a 95% confidence interval estimate of a coefficient from one of the two models. Each group represents confidence interval estimates of the same coefficient from the two models. I think such a bar plot will nicely present whether 95% confidence interval estimates of the same

...ues represent the duration of the cycle (in months). In other words, all rectangles have coordinates (start, duration) and (end, duration). rr I've spent hours trying to figure out the documentation and pouring over Google and RSeek searches and am at an impasse. The documentation refers to xmin, xmax, ymin, and ymax but doesn't say anything about them. The only example gives them both as vectors, so I assume they refer to a sequence of coordinates in which each rectangle's vertices is given by (xmin[i],ymin[i]), (xmin[i],ymax[i]), (xmax[i],ymax[i]), and (xmax[i],ymin[i]). But...

Dear R users, Consider the following example function: f = function(a,b) { g = function(x) a*x + b h = function(x) g(x)^2 + x^2 opt = optimize(h,lower = -1, upper = 1) x.min = opt$minimum h.xmin = opt$objective g.xmin = g(x.min) return(c(x.min, h.xmin, g.xmin)) } In my real problem the function that plays the role of "g" is costly to compute. Now, to minimize "h", "optimize" calls "h" with different values of x. In particular, at the end o...

signif(.) is a <primitive> function. Unfortunately, I couldn't even find WHERE in the source, signif(.) is defined. Here are the symptoms: xmin <- .Machine $ double.xmin signif(xmin,3) #--> NA umach <- unlist(.Machine)[paste("double.x", c("min","max"), sep='')] for(dig in 1:10) {cat("dig=",dig,": "); print(signif(umach, digits = dig))} ------- Why do I "complain&quo...

I use a 32 bit machine. For those using 64 bit machines, what is the .Machine$double.xmin for for machines? regards. Stephen. ----- Original Message ----- From: Achim Zeileis <Achim.Zeileis at wu-wien.ac.at> Date: Tuesday, July 12, 2005 10:51 am Subject: Re: [R] exact values for p-values - more information. > On Tue, 12 Jul 2005, S.O. Nyangoma wrote: > > > >...

...32L, 254622L, 263114L, 274370L, 282955L, 282352L, 278092L)), .Names = c("year", "federal.ps"), class = "data.frame", row.names = c(NA, -13L)) p <- ggplot(fcs, aes(year, federal.ps )) + geom_line() p + geom_rect(data=rect, aes(xmin=xmin, xmax = xmax, ymin=-Inf, ymax = Inf), fill='red', alpha=0.2) ##=============End Code========================== sessionInfo() R version 2.15.3 (2013-03-01) Platform: i686-pc-linux-gnu (32-bit) locale: [1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C [3] L...

...<?xml version="1.0" encoding="UTF-8" ?> - <layout> <name>Layout 2</name> <pagesByX>1</pagesByX> <pagesByY>1</pagesByY> - <elements> <legend refParent="r3" xmax="0.54156171284634758" xmin="0.020151133501259435" ymax="0.31480598161051165" ymin="0.042461931948864044" /> <component background="auto" border="none" borderDegMinSec="false" borderEachPage="true" borderFont="Tahoma, 8" borderFore=&quot...

...panel = panel.superpose, log10(states) ~ size, groups=category, data=data, type='b', lwd = 2, par.settings = list(superpose.symbol=mysymbols), ylim=c(y_min, y_max), scales = list(tck=c(1, 0), axs="r", x=list(tick.number=(xmax - xmin + 1), at=xmin:xmax, labels=xmin:xmax, cex=1.75), y=list(axs="r", rot=c(90, 0), labels=y_labels, at=y_at, cex=1.75 ) ), key = list ( text = list(levels(data$category)), lines = li...

I want to make a plot, but instead of showing _numerical_ values, I would like to show _symbolic_ values. For example, I want to plot a function y = a x + b, where x varies between Xmin and Xmax. I would like the plot to show, in the x-axis, the strings Xmin and Xmax, instead of their numeric values. Is it possible? Alberto Monteiro

...se<-std/sqrt(n) lowb<-mn-se upb<-mn+se data.frame(tx=x$tx[1],day=x$day[1],mn=mn,std=std,lowb=lowb,upb=upb,se=se) }) summ<-do.call("rbind",summ) #Definining x axis range xmax<-unique(max(summ$day,na.rm=TRUE)) xmin<-unique(min(summ$day,na.rm=TRUE)) yaxmin<-unique(min(summ$lowb)) yaxmax<-unique(max(summ$upb)) plot(1,1,type='n',xlab='Day',xlim=c(xmin,xmax),ylim=c(yaxmin,yaxmax), ylab=k, las=1,cex.lab=1,xaxp=c(xmin,xmax,diff(range(c(xmin,xmax))))) points(su...

All: Are there any functions out there for selecting all the points in a region of a plot. I envision something like the identify() function except one could circle a cloud of points (and perhaps a vector would be returned of the same length as the points plotted indicating logical membership in the circled cloud). Perhaps someone has done something with the locator() function that would

...ta" module, I used: > rsaga.get.usage("io_grid", 4) SAGA CMD 2.0.3 library path: C:/Progra~1/saga_vc/modules library name: io_grid module name : Import Binary Raw Data Usage: 4 [-GRID <str>] [-FILE_DATA <str>] [-NX <num>] [-NY <num>] [-DXY <str>] [-XMIN <str>] [-YMIN <str>] [-UNIT <str>] [-ZFACTOR <str>] [-NODATA <str>] [-DATA_OFFSET <num>] [-LINE_OFFSET <num>] [-LINE_ENDSET <num>] [-DATA_TYPE <num>] [-BYTEORDER_BIG <num>] [-TOPDOWN <num>] -GRID:<str> Grid Data Object (optio...

...o this mapping inside somewhere, but I'd like to choose the colors for plotting a set of polygons. Is there a pre-existing function that does something like this well? i.e., is capable of using 'breaks'? quantize<-function(x,n=10, breaks=NULL){ # bin the variable x into n levels xmin<-min(x) xmax<-max(x) 1+floor(n*(x-xmin)/(xmax-xmin)*.999) } x<- -10:10 cbind(x,quantize(x,2),quantize(x),quantize(x,21)) quantize(x,breaks=c(5,7)) # Thanks for your time, Dave -- Dr. David Forrest drf at vims.edu (804)684-7900w drf5n at mapl...

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 limits. In short the line spans xmin to xmax(of plot axis) but I want it to span xmin to xmax (of x data). Can this but done with abline or do I have to use lines and construct the line from a+bx? Paul Rydelek, Ph.D. 5232 Dorset Dr. Memphis, TN 38117

In https://en.wikipedia.org/wiki/Pareto_distribution, it is clear what the parameters are for the pareto distribution: *xmin *the scale parameter and *a* the shape parameter. I am using rmutil to generate random deviates from a pareto distribution. It says in the documentation that the probabilty density of the pareto distribution The Pareto distribution has density f(y) = s (1 + y/(m (s-1)))^(-s-1)/(m (s-1)) where m...

Hello, I wish to count how often zero (0) appears in the vector test. test [1] 1 1 1 1 1 1 2 1 1 1 0 2 0 1 1 0 0 0 1 1 1 0 1 2 1 1 1 1 1 1 I think of something like ... > sapply (test, function (x) if (x==0 ... ... but actually I dont know how to carry on. Could anybody give me a hint? Thanks Hermann [[alternative HTML version deleted]]

Dear list, Does anyone know how to fit the power law distribution? I have the empirical distribution and would like to check whether it fits power law (with the power estimated from the data). Any hints are appreciated Best regards Galina [[alternative HTML version deleted]]

...out of bounds This always occurs after the second lines command, and doesn't happen with all of my data points (some do not have errors). Any ideas? Thanks, Sean function(x,annot,rat1,rat2,rf,...) { par(las=2) wh <- which(annot[,5]==x) xmax <- max(annot[wh,4]) xmin <- min(annot[wh,3]) chr <- annot[wh,2][1] wh.rf <- rf$chrom==as.character(chr) & rf$txStart>xmin & rf$txEnd<xmax par(mfrow=c(2,1)) plot(annot[wh,3],rat1[wh],type="l",xlab="",ylab="log2 Ratio",main=x,...) points(annot[...

...;- ellipse$call } else coord.ellipse <- NULL if (!is.null(ellipse.par)) { coord.ellipse.par <- ellipse.par$res npoint.ellipse.par <- ellipse.par$call } else coord.ellipse.par <- NULL if (is.null(xlim)) { xmin <- xmax <- 0 if (is.na(test.invisible[1])) xmin <- min(xmin, coord.ind[, 1]) if (is.na(test.invisible[1])) xmax <- max(xmax, coord.ind[, 1]) if (!is.null(res.mfa$ind.sup) & is.na(test.invisible[2]))...