Displaying 4 results from an estimated 4 matches for "rsm2".
Did you mean:
rsm
2011 Apr 01
5
Predicción de valor máximo en superficie de respuesta, con paquete rsm
...áximo predicho para la
variable de respuesta en los valores de x1 y x2 estimados.
-------------------------------------------------------------
library("rsm")
ChemReact
CR <- coded.data(ChemReact, x1 ~ (Time - 85)/5, x2 ~ (Temp - 175)/5)
CR
# Realizando un ajuste de segundo orden:
CR.rsm2 <- rsm(Yield ~ Block + SO(x1, x2), data = CR)
CR.rsm2
summary(CR.rsm2)
# Graficando:
contour(CR.rsm2, ~ x1 + x2, at = summary(CR.rsm2)$canonical$xs, image =
TRUE, img.col = terrain.colors(50))
xs <- canonical(CR.rsm2)$xs
points (xs[1], xs[2], pch=2, col="red")
points(CR[1:7,], col=...
2004 Jun 02
1
factor -> numeric
...n types in the same variable
(another variable deifnes the observation type). Some
of these observation types are factors, so R naturally
classifies the entire variable as a factor. I want to
select a subset and convert the values to numeric
values, but it isn't working; as shown below:
> rsm2<-subset(rsm,PlantStrataId==2,select=c(Value))
> rsm2$Value
[1] 70 30
15 50
[26] 2
[51]
6...
2002 Apr 25
1
An unexpected exception has been detected in native code outside the VM
...library(ctest) print(Rvals) file <- "zz" dev <- "png" print(dev) draw <- "pairs_corr_scatter" if (draw == "pairs_corr_scatter") { source("/export/home/users/ruser/rsc/pcscatter01.R") NODE = rnorm(77) RSM2 = rnorm(77) CADT = rnorm(77) VTS = c(NODE - c(rnorm(12), rep(0, 65))) xRS = c(NODE - c(rnorm(27), rep(0, 50))) reda <- data.frame(NODE = NODE, RSM2 = RSM2, CADT = CADT, VTS = VTS, xRS = xRS) pcscatter01(file = file, dev = dev, reda = reda, title = &...
2006 Sep 19
0
How to interpret these results from a simple gamma-frailty model
...e
rs<- age*beta[1]+tc*beta[2]+hdlc*beta[3]+sbp*beta[4]+diab*beta[5]+smok*beta[6]
pevent <- exp(-10*lambda)^(exp(rs)) # Original 10-year risk
rsM1<- (cf[i,1,1]*age+cf[i,1,2]*tc+cf[i,1,3]*hdlc+cf[i,1,4]*sbp+cf[i,1,5]*diab+cf[i,1,6]*smok)
peventM1 <- exp(-10*lambda)^(exp(rsM1))
rsM2<- (cf[i,1,1]*age+cf[i,1,2]*tc+cf[i,1,3]*hdlc+cf[i,1,4]*sbp+cf[i,1,5]*diab+cf[i,1,6]*smok)
peventM2 <- exp(-10*lambda)^(estv*exp(rsM1))
# Proportion of more accurate predictions
pred <- sum(abs(pevent-peventM2) < abs(pevent-peventM1))/n
print(paste("Theta = ",l," prop...