Bioequivalence and Bioavailability Forum

Phoenix → R [General Sta­tis­tics]

posted by Helmut  – Vienna, Austria, 2014-02-23 03:20 (3708 d 18:58 ago) – Posting: # 12473
Views: 20,125

Hi Lucas,

export the stuff to and some quick/dirty code:
In the Phoenix-project navigate to Stud Res intra > Output Data > Result right-click > Export
Filename: StudResIntra
Filetype: Text (Tab delimied) (*.txt)
Save
The default location is the desktop (which conveniently is also the default import site in R). Start R and paste this code into the console:

res <- read.table("StudResIntra.txt", header=TRUE, sep="\t")
attach(res)
ResLim <- c(-max(abs(StudResIntra))*1.1, max(abs(StudResIntra))*1.1)
PreLim <- c(min(abs(Predicted))*1.1, max(abs(Predicted))*1.1)
SubLim <- c(min(abs(Subject))*1.1, max(abs(Subject))*1.1)
TR <- subset(res, Sequence == "TR")
RT <- subset(res, Sequence == "RT")
RL <- function() {
  abline(h=c(0, -2, 2, -3, 3), lty=c(1, rep(2, 4)),
    col=c(rep("black", 2), rep("red", 2)))
}
lab <- "Studentized residual"
split.screen(c(3, 2))
screen(1)
plot(0, 0, type="n", xlim=PreLim, ylim=ResLim, xlab="Predicted", ylab=lab,
  main="Scatterplot 1", cex.main=0.9, las=1)
RL()
points(TR$Predicted, TR$StudResIntra, pch=23, col="red", bg="#ffdddd")
points(RT$Predicted, RT$StudResIntra, pch=22, col="darkgreen", bg="#ddffdd")
screen(2)
plot(0, 0, type="n", xlim=SubLim, ylim=ResLim, xlab="Subject", ylab=lab,
  main="Scatterplot 2", cex.main=0.9, las=1)
RL()
points(TR$Subject, TR$StudResIntra, pch=23, col="red", bg="#ffdddd")
points(RT$Subject, RT$StudResIntra, pch=22, col="darkgreen", bg="#ddffdd")
screen(3)
HisLim <- c(0, max(hist(StudResIntra, breaks="FD", plot=FALSE)$density)*1.1)
hist(StudResIntra, breaks="FD", freq=FALSE, xlim=ResLim, ylim=HisLim, xlab=lab,
  ylab="Frequency", main="Histogram", cex.main=0.9, las=1, border="blue",
  col="#ddddff")
curve(dnorm(x, mean=mean(StudResIntra), sd=sd(StudResIntra)), add=TRUE)
box()
screen(4)
qqnorm(StudResIntra, xlim=ResLim, ylim=ResLim, xlab="Theoretical quantiles",
  ylab=lab, main="Normal Q-Q plot", cex.main=0.9, pch=21, col="blue", bg="#ddddff",
  las=1, plot.it=TRUE)
qqline(StudResIntra, datax=F, distribution=qnorm)
screen(5)
boxplot(StudResIntra, xlim=c(0.5, 2.5), ylim=ResLim, ylab=lab, main="Boxplot",
  cex.main=0.9, las=1, at=1, border="blue", col="#ddddff")
RL()
points(x=jitter(rep(2, length(StudResIntra)), factor=4), StudResIntra,
  pch=21, col="blue", bg="#ddddff")
close.screen(all = TRUE)
detach(res)



Good luck! If you have a Connect license you even don’t have to leave PHX.

—
Dif-tor heh smusma 🖖🏼 Довге життя Україна! 
Helmut Schütz


The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

• Qualitative analysis of the ANOVA residuals Lucas 2014-02-19 14:27 
  • Qualitative analysis of the ANOVA residuals Helmut 2014-02-19 15:22
    • Qualitative analysis of the ANOVA residuals Lucas 2014-02-19 15:58
      • Qualitative analysis of the ANOVA residuals Helmut 2014-02-19 16:25
        • Qualitative analysis of the ANOVA resid­uals Lucas 2014-02-19 17:20
        • Qualitative analysis of the ANOVA resid­u­als Lucas 2014-02-20 19:07
          • First hints for Phoenix Helmut 2014-02-20 20:11
          • Phoenix: resid­u­als Helmut 2014-02-21 19:25
            • Phoenix: resid­u­als BRN 2014-03-04 15:15
              • Phoenix: resid­u­als Helmut 2014-03-04 16:52
            • Phoenix: resid­u­als Rosy 2018-10-15 07:03
              • Software validation Helmut 2018-10-15 11:43
                • Software validation Rosy 2018-10-18 00:38
                  • Avoid raw residuals Helmut 2018-10-18 17:14
                    • Avoid raw residuals Rosy 2018-10-18 18:05
          • Phoenix → RHelmut 2014-02-23 02:20
            • Phoenix → R Lucas 2014-02-25 19:28
  • Qualitative analysis of the ANOVA residuals ElMaestro 2014-02-19 15:32
    • Qualitative analysis of the ANOVA residuals Helmut 2014-02-19 16:06
    • Qualitative analysis of the ANOVA residuals Lucas 2014-02-19 16:13
      • Qualitative analysis of the ANOVA resid­u­als ElMaestro 2014-02-19 20:40