留言 | /cn/about/

访客 *hillluck* @ 2010-05-21 08:53:20 写道：

访客 *th* @ 2010-05-25 04:55:30 写道：

访客 *nancyhui* @ 2010-05-27 12:41:03 写道：

访客 *kappler* @ 2010-06-04 14:47:59 写道：

访客 *GentleYang* @ 2010-06-21 07:21:47 写道：

访客 *SelinZhang* @ 2010-06-23 21:41:47 写道：

访客 *yixuan* @ 2010-07-07 19:49:56 写道：

访客 *王梓骅* @ 2010-07-08 01:42:58 写道：

访客 *nicole* @ 2010-07-19 11:40:42 写道：

访客 *王彦飞* @ 2010-07-20 21:56:10 写道：

访客 *lywzh* @ 2010-08-08 07:48:18 写道：

访客 *Lan Wangsen* @ 2010-08-23 05:13:30 写道：

quakes$Magnitude <- equal.count(quakes$mag, 4)
cloud(depth ～ lat * long | Magnitude, data = quakes, 
      zlim = rev(range(quakes$depth)),
      screen = list(z = 105, x = -70), panel.aspect = 0.75,
      xlab = "Longitude", ylab = "Latitude", zlab = "Depth")
cloud(depth ～ lat * long | Magnitude, data = quakes, 
      zlim = rev(range(quakes$depth)), panel.aspect = 0.75,
      screen = list(z = 80, x = -70), zoom = 0.7,
      scales = list(z = list(arrows = FALSE, distance = 2)), 
      xlab = "Longitude", ylab = "Latitude", 
      zlab = list("Depthn(km)", rot = 90))
p <-
    cloud(depth ～ long + lat, quakes, zlim = c(690, 30),
          pch = ".", cex = 1.5, zoom = 1,
          xlab = NULL, ylab = NULL, zlab = NULL,
          par.settings = list(axis.line = list(col = "transparent")),
          scales = list(draw = FALSE))
npanel <- 4
rotz <- seq(-30, 30, length = npanel)
roty <- c(3, 0)
update(p[rep(1, 2 * npanel)], 
       layout = c(2, npanel),
       panel = function(..., screen) {
           crow <- current.row()
           ccol <- current.column()
           panel.cloud(..., screen = list(z = rotz[crow], 
                                          x = -60,
                                          y = roty[ccol]))
       })
state.info <-
    data.frame(name = state.name,
               long = state.center$x,
               lat = state.center$y,
               area = state.x77[, "Area"],
               population = 1000 * state.x77[, "Population"])
state.info$density <- with(state.info, population / area)
cloud(density ～ long + lat, state.info,
      subset = !(name %in% c("Alaska", "Hawaii")),
      type = "h", lwd = 2, zlim = c(0, max(state.info$density)),
      scales = list(arrows = FALSE))
library("maps")
state.map <- map("state", plot=FALSE, fill = FALSE)
panel.3dmap <- function(..., rot.mat, distance, xlim, ylim, zlim,
                        xlim.scaled, ylim.scaled, zlim.scaled)
{
    scaled.val <- function(x, original, scaled) {
        scaled[1] + (x - original[1]) * diff(scaled) / diff(original)
    }
    m <- ltransform3dto3d(rbind(scaled.val(state.map$x, xlim, xlim.scaled),
                                scaled.val(state.map$y, ylim, ylim.scaled),
                                zlim.scaled[1]), 
                          rot.mat, distance)
    panel.lines(m[1,], m[2,], col = "grey76")
}
cloud(density ～ long + lat, state.info,
      subset = !(name %in% c("Alaska", "Hawaii")),
      panel.3d.cloud = function(...) {
          panel.3dmap(...)
          panel.3dscatter(...)
      }, 
      type = "h", scales = list(draw = FALSE), zoom = 1.1,
      xlim = state.map$range[1:2], ylim = state.map$range[3:4],
      xlab = NULL, ylab = NULL, zlab = NULL,
      aspect = c(diff(state.map$range[3:4]) / diff(state.map$range[1:2]), 0.3),
      panel.aspect = 0.75, lwd = 2, screen = list(z = 30, x = -60),
      par.settings = list(axis.line = list(col = "transparent"),
                          box.3d = list(col = "transparent", alpha = 0)))
data(Cars93, package = "MASS")
cor.Cars93 <- 
    cor(Cars93[, !sapply(Cars93, is.factor)], use = "pair")
data(Chem97, package = "mlmRev")
Chem97$gcd <-
    with(Chem97, 
         cut(gcsescore,
             breaks = quantile(gcsescore, ppoints(11, a = 1))))
ChemTab <- xtabs(～ score + gcd + gender, Chem97)
ChemTabDf <- as.data.frame.table(ChemTab)
env <- environmental
env$ozone <- env$ozone^(1/3)
env$Radiation <- equal.count(env$radiation, 4)
cloud(ozone ～ wind + temperature | Radiation, env)
splom(env[1:4])
fm1.env <- lm(ozone ～ radiation * temperature * wind, env)
fm2.env <- 
    loess(ozone ～ wind * temperature * radiation, env,
          span = 0.75, degree = 1)
fm3.env <- 
    loess(ozone ～ wind * temperature * radiation, env,
          parametric = c("radiation", "wind"), 
          span = 0.75, degree = 2)
library("locfit")
fm4.env <- locfit(ozone ～ wind * temperature * radiation, env)
w.mesh <- with(env, do.breaks(range(wind), 50))
t.mesh <- with(env, do.breaks(range(temperature), 50))
r.mesh <- with(env, do.breaks(range(radiation), 3))
grid <- 
    expand.grid(wind = w.mesh, 
                temperature = t.mesh,
                radiation = r.mesh)
grid[["fit.linear"]] <- predict(fm1.env, newdata = grid)
grid[["fit.loess.1"]] <- as.vector(predict(fm2.env, newdata = grid))
grid[["fit.loess.2"]] <- as.vector(predict(fm3.env, newdata = grid))
grid[["fit.locfit"]] <- predict(fm4.env, newdata = grid)
wireframe(fit.linear + fit.loess.1 + fit.loess.2 + fit.locfit ～ 
                                     wind * temperature | radiation, 
          grid, outer = TRUE, shade = TRUE, zlab = "")
levelplot(fit.linear + fit.loess.1 + fit.loess.2 + fit.locfit ～ 
                                     wind * temperature | radiation, 
          data = grid)
contourplot(fit.locfit ～ wind * temperature | radiation, 
            data = grid, aspect = 0.7, layout = c(1, 4),
            cuts = 15, label.style = "align")
levelplot(volcano)
contourplot(volcano, cuts = 20, label = FALSE)
wireframe(volcano, panel.aspect = 0.7, zoom = 1, lwd = 0.01)
## plot(levelplot(volcano), split = c(1, 1, 1, 3), more = TRUE)
## plot(contourplot(volcano, cuts = 20, label = FALSE), split = c(1, 2, 1, 3), more = TRUE)
## plot(wireframe(volcano, panel.aspect = 0.7, zoom = 1, lwd = 0.01), 
##      split = c(1, 3, 1, 3), more = FALSE)
data(Chem97, package = "mlmRev")
Chem97$gcd <-
    with(Chem97, 
         cut(gcsescore,
             breaks = quantile(gcsescore, ppoints(11, a = 1))))
ChemTab <- xtabs(～ score + gcd + gender, Chem97)
ChemTabDf <- as.data.frame.table(ChemTab)
data(Cars93, package = "MASS")
cor.Cars93 <- cor(Cars93[, !sapply(Cars93, is.factor)], use = "pair")
levelplot(cor.Cars93, 
          scales = list(x = list(rot = 90)))
ord <- order.dendrogram(as.dendrogram(hclust(dist(cor.Cars93))))
levelplot(cor.Cars93[ord, ord], at = do.breaks(c(-1.01, 1.01), 20),
          scales = list(x = list(rot = 90)))
tick.at <- pretty(range(sqrt(ChemTabDf$Freq)))
levelplot(sqrt(Freq) ～ score * gcd | gender, ChemTabDf, 
          shrink = c(0.7, 1),
          colorkey = 
          list(labels = list(at = tick.at, labels = tick.at^2)),
          aspect = "iso")
cloud(Freq ～ score * gcd | gender, data = ChemTabDf, type = "h",
      aspect = c(1.5, 0.75), panel.aspect = 0.75)
library("latticeExtra")
cloud(Freq ～ score * gcd | gender, data = ChemTabDf,
      screen = list(z = -40, x = -25), zoom = 1.1,
      col.facet = "grey", xbase = 0.6, ybase = 0.6,
      par.settings = list(box.3d = list(col = "transparent")),
      aspect = c(1.5, 0.75), panel.aspect = 0.75,
      panel.3d.cloud = panel.3dbars)
library("copula")
grid <-
    expand.grid(u = do.breaks(c(0.01, 0.99), 25),
                v = do.breaks(c(0.01, 0.99), 25))
grid$frank  <- with(grid, dcopula(frankCopula(2),    cbind(u, v)))
grid$gumbel <- with(grid, dcopula(gumbelCopula(1.2), cbind(u, v)))
grid$normal <- with(grid, dcopula(normalCopula(.4),  cbind(u, v)))
grid$t      <- with(grid, dcopula(tCopula(0.4),      cbind(u, v)))
wireframe(frank + gumbel + normal + t ～ u * v, grid, outer = TRUE, 
          zlab = "", screen = list(z = -30, x = -50), lwd = 0.01)
wireframe(frank + gumbel + normal + t ～ u * v, grid, outer = TRUE, 
          zlab = "", screen = list(z = -30, x = -50), 
          scales = list(z = list(log = TRUE)), lwd = 0.01)
kx <- function(u, v)
    cos(u) * (r + cos(u/2) * sin(t*v) - sin(u/2) * sin(2*t*v))
ky <- function(u, v) 
    sin(u) * (r + cos(u/2) * sin(t*v) - sin(u/2) * sin(2*t*v))
kz <- function(u, v) 
    sin(u/2) * sin(t*v) + cos(u/2) * sin(t*v)
n <- 50
u <- seq(0.3, 1.25, length = n) * 2 * pi
v <- seq(0, 1, length = n) * 2 * pi
um <- matrix(u, length(u), length(u))
vm <- matrix(v, length(v), length(v), byrow = TRUE)
r <- 2
t <- 1
wireframe(kz(um, vm) ～ kx(um, vm) + ky(um, vm), shade = TRUE,
          screen = list(z = 170, x = -60),
          alpha = 0.75, panel.aspect = 0.6, aspect = c(1, 0.4))
data(USAge.df, package = "latticeExtra")
str(USAge.df)
library("RColorBrewer")
brewer.div <-
    colorRampPalette(brewer.pal(11, "Spectral"),
                     interpolate = "spline")
levelplot(Population ～ Year * Age | Sex, data = USAge.df,
          cuts = 199, col.regions = brewer.div(200),
          aspect = "iso")

访客 *Lan Wangsen* @ 2010-08-23 05:23:12 写道：

访客 *Lan Wangsen* @ 2010-08-23 09:43:42 写道：

访客 *Lan Wangsen* @ 2010-08-23 09:47:36 写道：

访客 *Lan Wangsen* @ 2010-08-23 09:53:10 写道：

访客 *兰旺森* @ 2010-08-24 13:27:54 写道：

访客 *兰旺森* @ 2010-08-24 13:37:15 写道：

访客 *王少斌* @ 2010-10-09 18:18:51 写道：

访客 *lg21c* @ 2010-10-10 05:36:51 写道：

访客 *lg21c* @ 2010-10-10 06:44:45 写道：

访客 *lg21c* @ 2010-10-10 07:29:21 写道：

访客 *lg21c* @ 2010-10-10 07:42:02 写道：

访客 *lg21c* @ 2010-10-10 07:51:05 写道：

访客 *苏菲的世界* @ 2010-10-10 17:12:10 写道：

访客 *Cc* @ 2010-10-11 16:03:44 写道：

访客 *lg21c* @ 2010-10-12 12:14:12 写道：

访客 *lg21c* @ 2010-10-12 12:45:18 写道：

访客 *lg21c* @ 2010-10-13 14:21:07 写道：

访客 *Lower-Control-Line* @ 2010-10-15 18:40:54 写道：

访客 *wangshaobin* @ 2010-10-21 07:26:06 写道：

访客 *fangjie* @ 2010-11-18 08:14:34 写道：

访客 *huzhao* @ 2010-12-03 09:56:35 写道：

访客 *橙子* @ 2011-01-04 16:11:38 写道：

访客 *橙子* @ 2011-01-06 13:04:13 写道：

访客 *尚磊* @ 2011-01-26 10:20:21 写道：

访客 *刘兵* @ 2011-02-02 21:55:32 写道：

访客 *yixuan* @ 2011-02-03 13:22:47 写道：

访客 *尚磊* @ 2011-02-22 21:08:06 写道：

访客 *尚磊* @ 2011-02-24 19:23:49 写道：

访客 *李毅* @ 2011-04-12 13:58:12 写道：

访客 *Hao Zhang* @ 2011-04-14 01:25:37 写道：

访客 *任丽烨* @ 2011-04-22 06:12:19 写道：

访客 *张静* @ 2011-05-09 16:35:15 写道：

访客 *张静* @ 2011-05-09 16:40:45 写道：

访客 *室光* @ 2011-05-10 20:07:39 写道：

访客 *FUgui* @ 2011-05-13 19:42:33 写道：

model volatility;

{
mu ～ dnorm(0,0.1);
phistar ～ dbeta(20,1.5);
itau2 ～ dgamma(2.5,0.025);
k～dchisqr(8)I(2,50);

theta0 ～ dnorm(mu,itau2);

thmean[1] <- mu + phi*(theta0-mu);
theta[1] ～ dnorm(thmean[1],itau2);

for (t in 2:n) { thmean[t] <- mu + phi*(theta[t-1]-mu);
                     theta[t] ～ dnorm(thmean[t],itau2);
                   }
for (t in 1:n) { yisigma2[t] <- 1/exp(theta[t]);
                     y[t] ～ dt(0,yisigma2[t] ,k);
                   }


phi <- 2*phistar-1;
tau <- sqrt(1/itau2);

}

访客 *FUgui* @ 2011-05-13 19:48:03 写道：

访客 *FUgui* @ 2011-05-13 19:51:50 写道：

访客 *FUgui* @ 2011-05-13 20:47:38 写道：

访客 *Sid* @ 2011-06-10 10:51:39 写道：

访客 *isicic* @ 2011-06-23 05:22:07 写道：

访客 *ma* @ 2011-07-16 01:31:48 写道：

访客 *程贺* @ 2011-07-28 09:31:30 写道：

访客 *驴子* @ 2011-07-29 20:05:34 写道：

访客 *驴子* @ 2011-08-02 16:50:34 写道：

访客 *william_ou* @ 2011-08-09 14:52:05 写道：

访客 *huangshujia* @ 2011-08-13 14:09:10 写道：

访客 *cici93* @ 2011-08-19 22:17:00 写道：

访客 *chentong* @ 2011-08-21 05:14:17 写道：

访客 *WU Jun* @ 2011-09-15 07:26:21 写道：

访客 *lg21c* @ 2011-10-04 05:06:20 写道：

访客 *合乎* @ 2011-10-24 16:03:59 写道：

访客 *ypchen* @ 2011-11-16 10:32:49 写道：

访客 *gouweichao* @ 2011-11-20 19:54:00 写道：

访客 *Guo MengBiao* @ 2011-12-03 16:09:15 写道：

访客 *Cc* @ 2011-12-06 05:43:26 写道：

访客 *letranger* @ 2011-12-29 22:35:07 写道：

model IRT
{
#modelling
for (i in 1: nstud) {
for (j in 1: nqs) {
p[i,j]<-c[j]+(d[j]-c[j])*(exp(1.7*a[j]*(theta[i]-b[j]))/(1+exp(1.7*a[j]*(theta[i]-b[j]))))
r[i,j]～dbern(p[i,j])
}
theta[i]～dnorm(0,1)
}
#Priors
for (k in 1:nqs) {
a[k]～dlnorm(0, 8 )
b[k]～dnorm(0,  .25)
c[k]～dbeta(5, 17)
d[k]～dbeta(17, 5)
}
}

#initial value
list(a=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),
b=c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0),
c=c(0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1),
d=c(0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98,0.98))

#data
list(nstud = 489, nqs = 83, r=structure(.Data=c(0,0,0,1,0,0,…1,1),.Dim = c(489,83)))

访客 *Tlexander* @ 2012-01-15 18:03:51 写道：

访客 *yixuan* @ 2012-02-12 02:57:36 写道：

访客 *Cloud* @ 2012-02-21 12:36:26 写道：

访客 *Dengyishuo* @ 2012-02-22 12:58:56 写道：

访客 *tanglei* @ 2012-02-23 02:11:44 写道：

访客 *360512239* @ 2012-02-27 13:18:58 写道：

访客 *jing* @ 2012-03-03 13:29:58 写道：

访客 *Janechung* @ 2012-03-09 08:24:30 写道：

访客 *tanglei* @ 2012-03-21 05:26:12 写道：

访客 *李青龙* @ 2012-03-22 03:12:37 写道：

访客 *平平* @ 2012-03-29 16:10:21 写道：

访客 *培德 章* @ 2012-03-31 01:23:24 写道：

访客 *Zhuwenjie* @ 2012-04-06 03:19:07 写道：

访客 *Journeyxk* @ 2012-04-08 13:49:46 写道：

访客 *adophper* @ 2012-05-06 12:00:27 写道：

访客 *hzqim* @ 2012-05-08 11:28:07 写道：

访客 *Christopher Meng* @ 2012-06-18 03:49:03 写道：

访客 *可可* @ 2012-06-23 08:57:30 写道：

访客 *小熊猫* @ 2012-07-04 18:41:11 写道：

访客 *Terra_li* @ 2012-08-02 21:06:27 写道：

访客 *Huangxin * @ 2012-08-03 18:10:49 写道：

访客 *Shzlg5416* @ 2012-08-09 10:44:56 写道：

访客 *Shzlg5416* @ 2012-08-09 10:47:27 写道：

访客 *290692407* @ 2012-08-17 17:12:12 写道：

访客 *刘冬梅* @ 2012-08-22 20:02:56 写道：

访客 *Junyu* @ 2012-08-29 00:33:38 写道：

访客 *devbian* @ 2012-09-07 03:36:57 写道：

访客 *Junyu* @ 2012-09-10 17:10:01 写道：

访客 *394263788* @ 2012-09-13 04:00:26 写道：

访客 *pala* @ 2012-10-02 06:14:06 写道：

访客 *Guest* @ 2012-10-12 03:37:37 写道：

访客 *杨雪蓉* @ 2012-10-13 23:11:05 写道：