2017年5月2日星期二

linear vs nolinear


A regression model is called nonlinear, if the derivatives of the model with respect to the model parameters depends on one or more parameters. This definition is essential to distinguish nonlinear from curvilinear regression. A regression model is not necessarily nonlinear if the graphed regression trend is curved. A polynomial model such as y = b0 + b1x + b2x2 + e appears curved when y is plotted against x. It is, however, not a nonlinear model. To see this, take derivatives of y with respect to the parameters b0b1, and b2: dy/db0 = 1, dy/db1 = x, dy/db2 = x2 None of these derivatives depends on a model parameter, the model is linear. In contrast, consider the log-logistic model y = d + (a - d)/(1 + exp{b log(x/g)}) + e Take derivatives with respect to d, for example: dy/dd = 1 - 1/(1 + exp{b log(x/g)}). The derivative involves other parameters, hence the model is nonlinear

.http://www.ats.ucla.edu/stat/sas/library/SASNLin_os.htm

Like Father Like Son


2017年4月27日星期四

The idea of utility -a very interesting explanation

https://stats.stackexchange.com/a/276376/61705

"The idea of utility is very simple. It converts the numerical outcome into the utility to you. Say, you get 10 apples or 5 oranges, or some combination of two. How would you compare apples to oranges? Easy, I offer you to take either of two choices, and you pick one. So, theory goes that there's some kind of a function 
U(x1,x2,,xn) such that when you feed it a basket of goods 1,2,,n with quantities xi of each good, this function returns a scalar value, which reflects the utility of this basket to you. You pick the one with a greater utility.
So, if I feed it U(10,0) or U(0,5), where the good 1 is apples, and the good 2 is oranges, it'll return two different numbers, the greater one will indicate which basket has more utility to you. YOu can combine goods in the same basket U(1,1) - one apple and one orange.
So, in a nutshell, the utility function is a concept that allows us to compare apples to oranges. The theory states that you pcik the baskets based on their utility, which is calculated by plugging the basket descriptions into this utility function."

2017年4月18日星期二

A function to calculate Asymptotic Relative Efficiency (ARE) of signed-rank Wilcoxon and contaminated normals

#sigma1=1 and sigma2=3

ARE_W_t<-function(epsilon){

effi<-(((1-epsilon)^2/(2*sqrt(pi))+epsilon^2/(6*sqrt(pi))+epsilon*(1-epsilon)/(2*sqrt(pi)))*sqrt(12))^2*(1+epsilon*8)

return (effi)

}

2017年4月16日星期日

A function for Walsh averages

library(Rfit)
x<-c(-4,-2,5,9)
x
walsh(x)

mywalsh<-function (x)
{
n <- length(x)
w <- vector(n * (n + 1)/2, mode = "numeric")
ind <- 0
for (i in 1:n) {
for (j in i:n) {
ind <- ind + 1
w[ind] <- 0.5 * (x[i] + x[j])
}
}
return(w)
}

mywalsh(x)