search for: lamda1

Displaying 6 results from an estimated 6 matches for "lamda1".

Did you mean: lamda
2009 Jul 30
1
lmer() and "$ operator is invalid for atomic vectors"
...=TRUE) attach(data) new<-data.frame(lamda, age.o_1, Item.l2_2,Subject) new Subject <- factor(Subject) Age <- factor(age.o_1,levels=c("0","1"),labels=c("young","old")) Item <- factor(Item.l2_2,levels=c("1","2"),labels=c("lamda1","lamda2")) new2<-data.frame(lamda, Subject, Age, Item) new2 new3 <- na.omit(new2) new3 new4=as.data.frame(new3) > (lmer.lamda <- lmer(lamda ~ Age*Item + ( 1 | Subject), lamda, data=new4)) Error in family$family : $ operator is invalid for atomic vectors Here are th...
2023 Aug 25
2
Query on finding root
Sir, I want to solve the equation Q(u)=mean, where Q(u) represents the quantile function. Here my Q(u)=(c*u^lamda1)/((1-u)^lamda2), which is the quantile function of Davies (Power-pareto) distribution. Hence I want to solve , *(c*u^lamda1)/((1-u)^lamda2)=28353.7....(Eq.1)* where lamda1=0.03399381, lamda2=0.1074444 and c=26104.50. When I used the package 'Davies' and solved Eq 1, I got the answer u=0.39...
2023 Aug 27
1
Query on finding root
On Fri, 25 Aug 2023 22:17:05 +0530 ASHLIN VARKEY <ashlinvarkey at gmail.com> wrote: > Sir, Please note that r-help is a mailing list, not a knight! ?? > I want to solve the equation Q(u)=mean, where Q(u) represents the > quantile function. Here my Q(u)=(c*u^lamda1)/((1-u)^lamda2), which is > the quantile function of Davies (Power-pareto) distribution. Hence I > want to solve , *(c*u^lamda1)/((1-u)^lamda2)=28353.7....(Eq.1)* > where lamda1=0.03399381, lamda2=0.1074444 and c=26104.50. When I used > the package 'Davies' and solved Eq 1, I g...
2006 Jul 22
1
Why the contrain does not work for selecting a particular range of data?
.... I used two way to do this but the results are different. 1.Way one I use the data for x2=1 and run the program. It works for me. Tthe program is described as below: function (parameters,y1,x11) { p<-parameters[1] alpha1<-parameters[2] beta1<-parameters[3] delta1<-parameters[4] lamda1<-parameters[5] mu<-alpha1*((x11)^beta1)*exp(-delta1*(x11^lamda1)) ifelse(y1>0|x11>0, L<-lgamma(y1+p)+p*(log(p)-log(mu+p))+y1*(log(mu)-log(mu+p))-lfactorial(y1)-lgamma(p) ,Inf) L } This is working for me. 2 Way two: I select the data whose x2=1 in the whole range of dat...
2008 Jan 06
2
how to get residuals in factanal
In R factanal output, I can't find a function to give me residuals e. I mannually got it by using x -lamda1*f1 -lamda2*f2 - ... -lamdan*fn, but the e I got are not uncorrelated with all the f's. What did I do wrong? Please help. Yijun ____________________________________________________________________________________ Be a better friend, newshound, and
2009 Jul 02
2
constrained optimisation in R.
i want to estimate parameters with maximum likelihood method with contraints (contant numbers). for example sum(Ai)=0 and sum(Bi)=0 i have done it without the constraints but i realised that i have to use the contraints. Without constraints(just a part-not complete): skellamreg_LL=function(parameters,z,design) { n=length(z); mu=parameters[1]; H=parameters[2]; Apar=parameters[3:10];