search for: hzau

...n CJK font, and it's not what I hope to get. I hope the Latin character is displayed in Helvetica. Any suggestions? Thanks in advance! Regards, Jinsong -- Jinsong Zhao, Ph.D. College of Resources and Environment Huazhong Agricultural University Wuhan 430070, P.R. China E-mail: jszhao at mail.hzau.edu.cn

..."\u4F60\u597D", family="GB1", font = 4) ## symbol text(5,5, "is 'hello' in Chinese", font = 4) # bold italic dev.off() -- Jinsong Zhao, Ph.D. College of Resources and Environment Huazhong Agricultural University Wuhan 430070, P.R. China E-mail: jszhao at mail.hzau.edu.cn

...- num = 1000 + (cidfontIndex - 1)*5 + 1 + face; + num = 1000 + (cidfontIndex - 1)*5 + face; } else { cidfontfamily = NULL; Regards, Jinsong -- Jinsong Zhao, Ph.D. College of Resources and Environment Huazhong Agricultural University Wuhan 430070, P.R. China E-mail: jszhao at mail.hzau.edu.cn

Hello, I plot a 3D plot/surface using cloud function in lattice package. Several days ago, I post a question about how to remove the framebox out the 3D plot/surface. However, I think I don't make me understood clearly. I just hope to romove the outer rectangle or box around the box. Thanks again. Regards, Jinsong Zhao

Hi, Is there a package that could do response surface analysis equivalent to SAS RSREG procedure? Thanks! Regards, Jinsong

Hi, I hope to use R to perform some hypothesis testing, such as t.test() and var.test(). However, I don't find a function that could do test on means with variance known (i.e., u test or z test in some textbook), and a function that could do test on variances of normal distribution (i.e. chi-squared test). Thanks in advance for any hints. Best wishes, Jinsong

Hi R users, Is is possible for me to use the try function with boot? I would to do the bootstraping with a nonlinear model(it works well when R < 1000). But it does not work very well (when R is large) thus I try to use "try" to resolve. I put the try function in two cases: case1: put the try in front of the boot > c1.try<-try(boot(c1data, statistic = c1.fun,

Hi there, I hope to draw a plot like this: http://www.sg-chem.net/swizard/Ru-bqdi-spectra.gif is it possible to draw it using R? thanks for any suggestions. regards, Jinsong

Dear all, I have a very simple set of data and I would like to analyze them with probit analysis. dose event trial 0.0 3 15 1.1 4 15 1.3 4 15 2.0 3 15 2.2 5 15 2.8 4 15 3.7 5 15 3.9 9 15 4.4 8 15 4.8 11 15 5.9 12 15 6.8 13 15 The dose should be transformed with log10(). I use glm(y ~ log10(dose), family=binomial(link=probit)) to do probit analysis, however, I have to exclude the

Dear all, I have encountered a problem when perform stepwise regression. The dataset have more 9 independent variables, but 7 observation. In R, before performing stepwise, a lm object should be given. fm <- lm(y ~ X1 + X2 + X3 + X11 + X22 + X33 + X12 + X13 + X23) However, summary(fm) will give: Residual standard error: NaN on 0 degrees of freedom Multiple R-Squared: 1, Adjusted

Hello list. I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without modifications. How did I try it? Created a (non-root) build environment (not a mock ) Installed the kernel.scr.rpm and did a rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee prep-out.log The build failed at the end: Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL Checking

Hi all, The package crossdes could contruct a complete sets of mutually orthogonal latin squares. The construction works for prime powers only. I hope to know whether there is a way to construct a mutually orthogonal Lation square for 10 or other numbers that could not be prime powers. Thanks for any suggestions. Best wishes, Jinsong Zhao

Hi there, I use boot() to do bootstrap simulation on my statistic, I get a boot object, named "obj". For certain reasons, there are some "NaN" in obj$t. Now, I hope to get confidence interval using boot.ci(), it give the following error: Error in if (const(t, min(1e-08, mean(t)/1e+06))) { : missing value where TRUE/FALSE needed Now, what I could do is I have to edit the

Hi there, I hope to use a string as an input in my function, however, I don't want to input the quotation mark. Is it possible? Thanks in advance. Regards, Jinsong

Hi I have a data frame, including x1, x2, x3, and y. I use lm() to fit second-order linear model, like the following: ft <- lm(y ~ x1 + x2 + x3 + I(x1 * x1) + I(x1 * x2) + I(x1 * x3) + I(x2 * x2) + I(x2 * x3) + I(x3 * x3), mydata) if the independent variable number is large, the formula will be very long. Is there a easy way to write formula like the above one? I have read the R

Hello, I hope to plot my data in a dot diagram, something like the output of the last example of layout function in graphics package. The only difference is that I need the rectangle in histogram is replaced by dot. Thanks in advance. Best wishes, Jinsong Zhao

Hi there, I hope to build a model Y ~ X1 + X2 + X3 + X4 with X1 has two instrumental variable A and B, and X2 has one instrumental variable A. I have searched the R site and mailling list, and known that the tsls() from sem package and ivreg() from AER package can deal with instrumental variable regression, however, I don't know how to formula the model. Any suggestion will be really

Hi there, When I do bootstrap on a maximum likelihood estimation, I try the following code, however, I get error: Error in minuslogl(alpha = 0, beta = 0) : object "x" not found It seems that mle() only get data from workspace, other than the boot.fun(). My question is how to pass the data to mle() in my case. I really appreciated to any suggestions. Best wishes, Jinsong

Hi, When I run the following code, r <- c(3,4,4,3,5,4,5,9,8,11,12,13) n <- rep(15,12) x <- c(0, 1.1, 1.3, 2.0, 2.2, 2.8, 3.7, 3.9, 4.4, 4.8, 5.9, 6.8) x <- log10(x) fr <- function(c, alpha, beta) { P <- c + (1-c) * pnorm(alpha + beta * x) P <- pmax(pmin(P,1),0) -(sum(log(choose(n,r))) + sum(r * log(P)) + sum((n -r)* log(1-P))) } fit <- mle((fr), start = list(c

Hi there, I just installed R 2.11.1 on my PC, which runs a Windows XP Home. The installation is successful, however, when I double click on the R icon, I get the following error message: R for Windows GUI front-end has encountered a problem and needs to close. We are sorry for the inconvenience. Error signature is: AppName: rgui.exe AppVer: 2.111.52157.0 ModName: msvcrt.dll ModVer: