search for: x_b

...ning was to long. Thanks in advance, Anne KERVAHU rm(list=ls()) nbr=10 #list of data x_a<-c(0,0.5,1,1.5,2,3,4,6,8,15,20,40,60,80,120) y_a<-c(5.4771213,5.0791812,4.8450980,4.3010300,4,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.361...

...tas: I am trying to plot a marginal distribution on the plane behind a persp() plot. My existing code is: library(MASS) X <- mvrnorm(1000,mu=c(0,0),Sigma=matrix(c(1,0,0,1),2)) X.kde <- kde2d(X[,1],X[,2],n=25) # X.kde is list: $x 1*n, $y 1*n, $z n*n persp(X.kde,phi=30,theta=60,xlab="x_b",ylab="x_a",zlab="f") ->res Any suggestions are very appreciated. Thank you, Colin [[alternative HTML version deleted]]

...the sum of squares. The model is number[2]-(x/number[1])^number[3] Number [2] and number [1] change as the 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, th...

...n help me. It’s simple but I suppose I am missing something. 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...

...ary outcome on x1, but the data format is now different. How can I actually feed R with mydatanew 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:...

Dear R developers, I am a little disappointed that my bug report only made it to the wishlist, with the argument: Well, it does not say it has. Only relevant to prediction intervals. predict.lm does calculate prediction intervals for linear models from weighted regression, so they should be correct, right? As far as I can see they are bound to be wrong in almost all cases, if no weights