Displaying 2 results from an estimated 2 matches for "x1bk".
Did you mean:
x1b
2005 Aug 08
1
bug found in predict.locfit in locfit package (PR#8057)
...043032, 0.449906,-0.951917,-0.077104,
-0.356833,-0.286042, 0.065750, 0.159677,-0.075792)
x2=c(-0.3050, 1.0125, 0.2050, 0.1025, 0.9550,
0.6975, 1.5550, 0.0225, 0.2575, 0.3725,
2.0075, 2.1275, 0.7200, 0.2950, 0.2875,
-0.2800,-0.6050, 0.2125,-0.5525,-1.7850)
ndat=length(y)
ybk=y
x1bk=x1
x2bk=x2
######## Joint probability function of y, x1 and x2
# fit joint probability function
fityxv=locfit(~y+x1+x2,alpha=1,deg=1)
fyxv=predict(fityxv,where="data")
######## Marginal distribution of gxv
# fit marginal distribution of y
fitxv=locfit(~x1+x2,alpha=0.5,deg=1)
gxv=predict(...
2005 Oct 05
0
bug found in predict.locfit in locfit package ( PR#8057)
...-0.356833,-0.286042, 0.065750, 0.159677,-0.075792)
> x2=c(-0.3050, 1.0125, 0.2050, 0.1025, 0.9550,
> 0.6975, 1.5550, 0.0225, 0.2575, 0.3725,
> 2.0075, 2.1275, 0.7200, 0.2950, 0.2875,
> -0.2800,-0.6050, 0.2125,-0.5525,-1.7850)
>
> ndat=length(y)
> ybk=y
> x1bk=x1
> x2bk=x2
> ######## Joint probability function of y, x1 and x2
> # fit joint probability function
> fityxv=locfit(~y+x1+x2,alpha=1,deg=1)
> fyxv=predict(fityxv,where="data")
>
> ######## Marginal distribution of gxv
> # fit marginal distribution of y
> fitx...