Displaying 4 results from an estimated 4 matches for "ptemp".
Did you mean:
temp
2010 Jan 21
3
cross validation function translated from stata
...free_psa"
local prediction1 = "base"
local prediction2 = "full"
g `prediction1'=.
g `prediction2'=.
quietly g u = uniform()
sort `event' u
g set = mod(_n, 10) + 1
forvalues j=1(1)10{
quietly logit `event' `predictors1' if set~=`j'
quietly predict ptemp if set==`j'
quietly replace `prediction1' = ptemp if set==`j'
drop ptemp
quietly logit `event' `predictors2' if set~=`j'
quietly predict ptemp if set==`j'
quietly replace `prediction2' = ptemp if set==`j'
drop ptemp
}
tempfile dca`i'
quietly dca `event'...
2011 Apr 12
1
How to set the dimension of a matrix correctly?
...eq(-105.5,-93.5,by=0.07273)pgrid <- expand.grid(x=plon,y=plat)pdim <- dim(pgrid) # (10230,2)
#plot(pgrid, cex=0.5)
lat <- coords[,1] lon <- coords[,2]
ppt1 <- ppt[,1:xx[2]] # 1:528
pptpred <- matrix(0,ncol=xx[2],nrow=1)
############# Only test one period data##################
ptemp <- ppt1[,3] ll <- which(ptemp>0)
ppt2 <- matrix(0,nrow=length(ll),ncol=3) # (lon,lat,ptemp)
ppt2[,1] <- lat[ll] # y-axis ppt2[,2] <- lon[ll] # x-axis ppt2[,3] <- ptemp[ll] # ppt
pptd <- as.geodata(ppt2)
bin1 <- variog(pptd) # plot(bin1) # fig1
bin2 <-...
2010 Mar 11
2
about IRT simulation
..." part doesn't run. we don't know what is wrong
with it.
Thanks so much~~~
I <- 10
J <- 5
response <- matrix(0, 10, 5)
pij <- function(a,b,theta)
{
a <- rnorm(J, 0.8, 0.04)
a
b <- rnorm(J, 0, 1)
b
theta <- rnorm(I, 0,1)
theta
for( i in 1:I ) {
for( j in 1:J ) {
ptemp <- runif(1)
pij <- exp(a[j]*(theta[i]-b[j]))/(1+exp(a[j]*(theta[i]-b[j])))
response[i,j]<-ifelse(pij(b=b[j], a=a[j], theta[i]) < ptemp , 0 ,1)
}
}
}
response
helena
2011 Jun 01
0
Anyone have experience with kinship pedigree plot?
...t.pedigree.html
symbolsize controls symbolsize. Default=1.
width default=8. For a packed pedigree, the minimum width allowed in
the realignment of pedigrees.
density defines density used in the symbols. Takes up to 4 different values.
mar no comment(s)
my pedigree command is similar to:
plot(ptemp, width=5, symbolsize=1);
however, 'width', doesnt seem to actually do anything as far as I can
see (ie. the plot looks the same no matter what value I use). the
example command in the kinship doc uses:
plot(x, id=x$id, sex=x$sex, status=x$status,
affected=x$affected, cex=1, col=rep(1, le...