search for: x_c

...he data changes but for all the set of data number[3] must be identical. I have 3 set of data (x1,y1), (x2,y2), (x3,y3). x_a<-c(0,0.5,1,1.5,2,3,4,6) y_a<-c(5.4,5,4.84,4.3,4,2,1.56,1.3) x_b<-c(0,1,2,3,4,5,6,7,8,9,10,11,12) y_b<-c(5.34,4.98,4.66,4.06,3,3.4,2.7,2.9,2.6,2.6,1.9,1.3,1.4) x_c<-c(0,3,6,8,10,12,14,16,18,20,24,26,28,30) y_c<-c(5.5,5.1,4.3,4,3.7,3.2,3.04,2.778,2.56,2.25,1.78,1.44,1.32,1.2) x<-c(x_a,x_b,x_c) y<-c(y_a,y_b,y_c) long<-c(0,8,21,35) Hence, the sum of squares is: Sce= sum( sum((y- number[4]-(x/number[1])^number[7])^2)+ sum((y- number[5]-(x/n...

...le parameters DET = xnew*xnew - xnew*Mz + Cov_xy; Center = cbind(Mxz*(Myy-xnew)-Myz*Mxy , Myz*(Mxx-xnew)-Mxz*Mxy)/DET/2; Par = cbind(Center+centroid , sqrt(Center[2]*Center[2]+Mz+2*xnew)); return(Par) } #EXAMPLE library(plotrix) # Create a Circle of radius=10, centre=5,5 R = 10; x_c = 5; y_c = 5; thetas = seq(0,pi,(pi/64)); xs = x_c + R*cos(thetas) ys = y_c + R*sin(thetas) # Now add some random noise mult = 0.5; xs = xs+mult*rnorm(rnorm(xs)); ys = ys+mult*rnorm(rnorm(ys)); plot(xs,ys,pch=19,cex=0.5,col="red",xlim=c(-10,20),ylim=c(-10,20),asp=1) # real circle draw.cir...

...I generate a data matrix X, say of order n*p, where n represents independent row-vectors and p correlated col vectors. Let the row representation be X = (X’_1, . . ., X’_n)’. I generate the differences of these vectors as D_{ab} = X_a – X_b, with different indices a and b, and similarly D_{cd} = X_c – X_d, and make a quadratic form D’_{ab}D_{cd} = A_{abcd}, say. Please note that, here a is unequal b, c is unequal d, a is unequal c, b is unequal d. By slightly shuffling the same set of vectors, I generate two other similar quadratic forms, say A_{acbd} and A_{adcb}, where the indices are agai...

Hi, I have some doubt about how qualitative factors are coded in R. For instance, I consider a response y, a quantitative factor x and a qualitative factor m at 3 levels, generated as follow : y_c(6,4,2.3,5,3.5,4,1.,8.5,4.3,5.6,2.3,4.1,2.5,8.4,7.4) x_c(3,1,3,1,2,1,4,5,1,3,4,2,5,4,3) m_gl(3,5) lm(y~x+m) Coefficients: (Intercept) x m2 m3 3.96364 0.09818 0.44145 0.62291 In literature, 2 usual implicit coding process are suggested : m1=0 or m1=-m2-m3. Does R use one of these process ? (I've...

Hello. > x_c(1,2,5,7) > plot(x,type="h",lwd=10) produces for me a plot where the tops of the high-density lines are rounded. Is there a way to get the high-density lines so that their tops are flat? Thanks, Michael -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

...2.30103,1.5563025,1.30103,1 ,1.6434527,0.60206,0.60206,0.30103,0,0.4771213) x_b<-c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14) y_b<-c(5.3424227,4.9867717,4.6627578,4.0606978,3.07182,3.4771213,2.7993405,2 .9395193,2.69897,2.6127839,1.9777236,1.3222193,1.4149733,0.7781513,1.322219) x_c<-c(0,3,6,10,12,14,16,18,20,24,26,28,30,34) y_c<-c(5.5185139,5.1461280,4.3617278,3.771513,3.20412,3.0413927,2.7781513,2. 5682017,2.255272,1.7823917,1.447158,1.3222193,1.2787536,0.69897) #number of kinetics nb=3 #complete data set x<-c(x_a,x_b,x_c) y<-c...

...the sample data to a formula that I supply it. Reading the documentation for loess - it seems that I can do this by supplying a formula.. but I have not had much luck. Perhaps I am reading the documentation incorrectly. Here is my current loess line code: mg.lo = loess(y ~ x, mg, span=0.25) mg$x_c <- predict(mg.lo) which produces a curve like the attached file loess_curve.png. and I would like to have a curve "something like" the arccot in the attached file arccotangent.png (from wikipedia), but all x values are positive - meaning just the best fit of the sample data to that t...

...anew to perform again a logistic regression on x1 only? (2) This is a bit more conceptual. Let us say that you have a set of products A,B,C,D,E,F. Each product has a list of features: x_A for product A, x_B for B etc... Each customer has its own set of parameters (age, sex, income etc..) I call x_cust. Finally, the customer is confronted with two products (e.g. A and D; combinations may vary, I call each combination of two products a scenario) and asked which one he would like to buy. Bottom line: your data are in the format 1 x_A x_cust 0 x_D x_cust meaning that a certain customer chos...

...aster.ms QA Contact: xorg-team at lists.x.org Created attachment 107505 --> https://bugs.freedesktop.org/attachment.cgi?id=107505&action=edit Xorg log Here you can see a video of it and also more important information is in the youtube description: https://www.youtube.com/watch?v=X_c-26mB304 Like the text says, there seems to be some kind of 2D drawing action in Xorg that takes maybe a thousand times longer with nouveau than usual. But this drawing is rarely used, i.e. mostly only on window creation, so the delay doesn't get as apparent in most cases you run into (I suppos...

...= min(x), u = max(x), n = as.integer(n))) Cpretty2 <- function (x, n = 5) { ##-- return what .C("pretty... doe n <- length(p <- pretty(x,n)) c(l=p[1], u=p[n], n=n-1) } Cpretty(pi) all(Cpretty1(pi) == Cpretty2(pi)) all(Cpretty1(c(0,pi)) == Cpretty2(c(0,pi))) for(i in 1:100) { x_c(runif(1),pi) for(k in -14:14) { ok <- all(Cpretty1(x*10^k) == Cpretty2(x*10^k)) if(!ok) cat("NOT ok: x=",deparse(x)," k=",k,"\n") } } Gpretty <- function (x, n = 5) { plot.new(); op <- par(); on.exit(par(op)) par(lab=c(n, par("lab")...