Displaying 8 results from an estimated 8 matches for "yint".
Did you mean:
int
2010 Apr 19
3
nls for piecewise linear regression not converging to least square
...copied from real data.
y1<-c(2.37700445, 1.76209775, 0.09795576, 2.21834963, 6.62262243,
15.70471269, 21.92956392, 36.39401717, 32.43620195, 44.77442277)
x1<-c(24.6, 28.9, 33.2, 37.6, 42.0, 46.4, 50.9, 55.3, 59.8, 64.3)
dat <- data.frame(x1,y1)
nlmod <- nls(y1 ~ ifelse(x1 < xint+(yint/slp), yint, yint +
(x1-(xint+(yint/slp)))*slp),
data=dat, control=list(minFactor=1e-5,maxiter=500,warnOnly=T),
start=list(xint=39.27464924, yint=0.09795576, slp=2.15061064),
na.action=na.omit, trace=T)
##plotting the function
plot(dat$x1,dat$y1)
segments(x0=0, x...
2010 Apr 19
2
nls minimum factor error
...ties?
##Here's a sample dataset and code.
y2<-c(2.404529, 1.625661, 1.013981, 3.810921, 10.023745, 10.990817,
10.740636, 11.246827,17.022761, 21.430386)
x2<-c(25.0, 29.3, 33.8, 38.3, 42.8, 47.2, 51.6, 55.8, 60.4, 64.9)
dat <- data.frame(x2,y2)
nlmod <- nls(y2 ~ ifelse(x2 < xint+(yint/slp), yint, yint +
(x2-(xint+(yint/slp)))*slp),
data=dat, control=list(minFactor=1e-5,maxiter=500,warnOnly=T),
start=list(xint=40.49782, yint=1.013981, slp=0.8547828),
na.action=na.omit, trace=T)
##plotting the function
plot(dat$x2,dat$y2)
segments(x0=0, x1=coef...
2017 May 31
2
stats::line() does not produce correct Tukey line when n mod 6 is 2 or 3
...double *w, /* work and output: resid. & fitted */
+ int *indx, /* to get thirds of points independently of repeated values */
+ /* all the above of length */ int n,
+ double coef[2])
{
int i, j, k;
double xb, x1, x2, xt, yt, yb, tmp1, tmp2;
double slope, yint;
for(i = 0 ; i < n ; i++) {
- z[i] = x[i];
- w[i] = y[i];
+ z[i] = x[i];
+ w[i] = y[i];
+ indx[i] = i;
}
- R_rsort(z, n);/* z = ordered abscissae */
+ rsort_with_index(z, indx, n);/* z = ordered abscissae */
- tmp1 = z[il(n, 1./6.)];
- tmp2 = z[iu(n...
2002 Jan 25
2
selecting clusters of points
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
2017 May 31
4
stats::line() does not produce correct Tukey line when n mod 6 is 2 or 3
Seriously, if a method gives a wrong result, it's wrong. line() does NOT
implement the algorithm of Tukey, even not after the patch. We're not
discussing Excel here, are we?
The method of Tukey is rather clear, and it is NOT using the default
quantile definition from the quantile function. Actually, it doesn't even
use quantiles to define the groups. It just says that the groups
2006 Aug 02
0
Trying to use segmented in a function
...n: this package seems to be the one that most people use for
this, as far as I can see.
Since I want to run some simulations, I'm trying to write functions
that use segmented, and it's driving me mad. Here's a simple example:
simdata<-function
(Ns=200,Xmean=20,Xsd=5,SdYerr=0.5,Yint=0,threshold=20,slopebelow=0.5,slo
peabove=1)
{
Xs<-rnorm(Ns,Xmean,Xsd)
Yerr<-rnorm(Ns,0,SdYerr)
D<-ifelse(Xs<=threshold,0,1)
XminusX0<-Xs-threshold
Ys<-Yint+slopebelow*Xs+slopeabove*XminusX0*D+Yerr
plot(Xs,Ys)
linmod<-lm(Ys~Xs)
segment<-segmented(linmod,Z=Xs,psi=...
2011 Jul 18
1
Spam (Usuario R)
...@uv.es)
----------------------------------------------------------------------
Message: 1
Date: Fri, 15 Jul 2011 04:20:47 -0700 (PDT)
r-help-es@r-project.org, r-help-request@r-project.org,
mabel.rad@nielsen.com
Subject: Re: [R-es] Hawai''i,
Message-ID:
<1310728847.93356.yint-ygo-j2me@web161508.mail.bf1.yahoo.com>
Content-Type: text/plain; charset=utf-8
hello! http://gmidistribuidora.com.br/indexwx5.php?yluckyid=37du6
------------------------------
Message: 2
Date: Fri, 15 Jul 2011 13:45:35 +0200
From: Usuario R <r.user.spain@gmail.com>
To: r-help-es <r...
2011 Nov 21
0
[LLVMdev] How to get ELF section virtual starting address from MCSymbolRefExpr?
...dy 'help' to
llvmdev-request at cs.uiuc.edu
You can reach the person managing the list at
llvmdev-owner at cs.uiuc.edu
When replying, please edit your Subject line so it is more specific
than "Re: Contents of LLVMdev digest..."
Today's Topics:
1. llvm_anyint_ty clarification (Alberto Magni)
2. Re: How to make Polly ignore some non-affine memory accesses
(Bernhard Reutner-Fischer)
3. Re: How to make Polly ignore some non-affine memory accesses
(Tobias Grosser)
4. Re: How to make Polly ignore some non-affine memory...