similar to: how to deal with NaN in "boot" object

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

Hi there, I am a Chinese R user. I hope to display Chinese character in a plot, and than save it in PostScript format. I have read the article titled "Non-Standard Fonts in PostScript and PDF Graphics", especially the section about CJK fonts. I also tried the code: > pdf("chinese.pdf", width=3, height=1) > grid.text("\u4F60\u597D", y=2/3,

Hi there, According to ?par, 'font' is an integer which specifies which font to use for text, that 1 corresponds to plain text (the default), 2 to bold face, 3 to italic and 4 to bold italic. When I test Chinese character in pdf(), I found that 1 to bold face, 2 to italic, 3 to bold italic, 4 to symbol. and I don't find how to set plain text. In the following code, the font to use

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 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

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,

It looks like an encoding problem. It works fine for me with R encoding set to UTF-8 Here is part of my sessionInfo() results [1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C [3] LC_TIME=en_CA.UTF-8 LC_COLLATE=en_CA.UTF-8 [5] LC_MONETARY=en_CA.UTF-8 LC_MESSAGES=en_CA.UTF-8 I would suggest issuing the command sessionInfo() and seeing what your encoding is. On Tue, 20 Oct 2020 at 08:22,

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

Hi there, I don't know why the following codes are return different results. > ifelse(3 > 2, 1:3, length(1:3)) [1] 1 > if (3 > 2) 1:3 else length(1:3) [1] 1 2 3 Any hints? Best, 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

Hi there, Why the same string is displayed in different form? > abc[,1] [1] "?land" "Afghanistan" > abc name 1 <c5>land 2 Afghanistan And more... > dput(abc, "aa.txt") > dget("aa.txt") name 1 <c5>land 2 Afghanistan > dget("aa.txt")[,1] [1] "<c5>land"

Hi there, I wrote a function that wraps MASS::boxcox as: bc <- function(vec) { lam <- boxcox(lm(vec ~ 1)) lam <- lam$x[which.max(lam$y)] (vec^lam - 1)/lam } When I invoke it as: > x <- runif(20) > bc(x) Error in eval(predvars, data, env) : object 'vec' not found I have googled, and rewrote the above function as: bc <- function(vec) { dat <<-

Hi there, I wrote a function that wraps MASS::boxcox as: bc <- function(vec) { lam <- boxcox(lm(vec ~ 1)) lam <- lam$x[which.max(lam$y)] (vec^lam - 1)/lam } When I invoke it as: > x <- runif(20) > bc(x) Error in eval(predvars, data, env) : object 'vec' not found I have googled, and rewrote the above function as: bc <- function(vec) { dat <<-

Hi Jinsong, Maybe "srt" will do what you want. As you noted non-orthogonal rotations may not work on some devices. You may be able to download mirror fonts from places like fontmirror.com, but it more complicated. Jim On Thu, Mar 30, 2023 at 7:25?PM Jinsong Zhao <jszhao at yeah.net> wrote: > > Hi there, > > > I just noticed that "crt" is a graphical

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

Also see the "arctext" function in the plotrix package. Jim On Fri, Mar 31, 2023 at 11:12?AM Jim Lemon <drjimlemon at gmail.com> wrote: > > Hi Jinsong, > Maybe "srt" will do what you want. As you noted non-orthogonal > rotations may not work on some devices. You may be able to download > mirror fonts from places like fontmirror.com, but it more

ifelse is vectorized. On Wed, Dec 13, 2017 at 7:31 AM, Jinsong Zhao <jszhao at yeah.net> wrote: > Hi there, > > I don't know why the following codes are return different results. > > > ifelse(3 > 2, 1:3, length(1:3)) > [1] 1 > > if (3 > 2) 1:3 else length(1:3) > [1] 1 2 3 > > Any hints? > > Best, > Jinsong > >

Hi there, I tried to export the names of country to a csv file with write.csv(). In the resulted file, ?land was coverted to <c5>land. Is there any way could prevent this happening? Thanks! > abc [1] "?land" > write.table(abc, file = "") "x" "1" "<c5>land" Best, Jinsong

hi there, I find the following codes produce strange output. plot(1:10, xlab = expression(NO[3]^-~(mg/L))) you will notice that the unit, mg/L is in superscript format. That means that "~" is not for space. However, the help page of plotmath does not mention this behavior. Best, Jinsong [[alternative HTML version deleted]]

Hi there, I just noticed that "crt" is a graphical parameter, which is stated as follows in the help page of "par": A numerical value specifying (in degrees) how single characters should be rotated. It is unwise to expect values other than multiples of 90 to work... However, I did not find any code example to demo the usage of "crt", and even know how to present a