similar to: How to code the factor for ANOVA

Hi there, The following data is obtained from a long-term experiments. > mydata <- read.table(textConnection(" + y year Trt + 9.37 1993 A + 8.21 1995 A + 8.11 1999 A + 7.22 2007 A + 7.81 2010 A + 10.85 1993 B + 12.83 1995 B + 13.21 1999 B + 13.70 2007 B + 15.15 2010 B + 5.69 1993 C + 5.76 1995 C + 6.39 1999

Hi there, Is there a function that can do a Fisher's LSD multiple comparisons in a two-way ANOVA? I hope to get a result similar with TukeyHSD(). Especially, I hope to know the significance of comparisons between the interactions of two factors. In the following example: x <- c(76, 84, 78, 80, 82, 70, 62, 72, 71, 69, 72, 74, 66, 74, 68, 66, 69, 72, 72, 78, 74, 71, 73, 67, 86, 67, 72,

Hi all, I analysed my data with lme and after that I spent a lot of time for mean separation of treatments (post hoc). But still I couldn’t make through it. This is my data set and R scripts I tried. replication fertilizer variety plot height 1 level1 var1 1504 52 1 level1 var3 1506 59 1 level1 var4 1509 54 1 level1 var2 1510 48 2 level1 var1 2604 47 2 level1 var4 2606 51 2 level1 var3

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,

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 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 R folk, First off I want to thank those of you who responded with comments for my R quick and dirty stats tutorial. They've been quite helpful, and I'm in the process of revising them. When it comes to repeated measures MANOVA, I'm in a bit of a bind, however. I'm beginning to see that all of the documentation is written for psychologists, who have a slightly

Hi! How to calculate the correct SE of mean in a nested or spliplot anova? Nested example: --------------------- m <- aov(Glycogen~Treatment+Error(Treatment/Rat/Liver)) > m Call: aov(formula = Glycogen ~ Treatment + Error(Treatment/Rat/Liver)) Grand Mean: 142.2222 Stratum 1: Treatment Terms: Treatment Sum of Squares 1557.556 Deg. of Freedom 2 Estimated

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

Sorry for bothering you - but I have a little problem. I`m completely new in R and trying to use lmer. That works, but I m not sure if my code is right, so I would just like to get an opinion from the experts. I have a Growthrate of plants (RGR Fixed Factors: Fertilizer, Status (alien or native plants) and their interaction Random Factors: Age is a covariable; than I have plantfamily and

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

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

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

Dear R Community, Below is a problem with a simple ggplot2 graph. The code returns the error message below. Error: stat_count() must not be used with a y aesthetic. My code is below and the data is attached as a ?text? file. # Graph of the probabilities library(digest) library(DT) datatable(probability) str(probability) probability$Fertilizer <- as.factor(probability$Fertilizer)

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 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 James, If you want to specify the y-values, you need to use stat="identity" as below: ggplot(probability, aes(x=Fertilizer, y=prob)) + geom_bar(stat="identity", aes(fill=Treatment)) best, huzefa On Tue, Apr 12, 2016 at 1:02 PM, James Henson <jfhenson1 at gmail.com> wrote: > Dear R Community, > > Below is a problem with a simple ggplot2 graph. The code

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