search for: x_a

Displaying 5 results from an estimated 5 matches for "x_a".

Did you mean: _a
2013 Oct 15
1
plotting a marginal distribution on the plane behind a persp() plot
...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]]
2005 Apr 01
1
optim problem, nls regression
Hi, I try to fit a non linear regression by minimising the sum of 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&l...
2013 Feb 18
3
Generating QFs from same sample
...ity can 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 quad...
2005 Apr 08
0
TR: The results of your email commands
...scrib but I don't know where. so if you can help me, it will be great.... I have try an other program where y and x are directly calculate in sce function but the time of running 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.47712...
2010 Oct 31
1
Questions about Probit Analysis
...egression of the binary 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 t...