Phoenix → R [General Sta­tis­tics]

posted by Helmut Homepage – Vienna, Austria, 2014-02-23 02:20  – Posting: # 12473
Views: 14,681

Hi Lucas,

export the stuff to R 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)


[image]

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

Cheers,
Helmut Schütz
[image]

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

Complete thread:

Activity
 Admin contact
20,250 posts in 4,262 threads, 1,398 registered users;
online 19 (0 registered, 19 guests [including 9 identified bots]).
Forum time (Europe/Vienna): 04:20 CET

No written law has ever been more binding than
unwritten custom supported by popular opinion.    Carrie Chapman Catt

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5