Displaying 6 results from an estimated 6 matches similar to: "data manipulation"
2009 Dec 07
1
anova/factor
Dear Wiza[R]ds,
I have the following data in a data.frame. I need to do an anova with
multiple comparison but I don't know how to factor the groups for analysis.
There are 3 groups, 1,2 and 3 labelled in column 1. Help appreciated with
thanks in advance.
Group SI Sif SG Io I2 lol2
1 9.08e-05 9.08e+00 0.060842287 1.798556446 32.574500
2009 Nov 28
1
Plot fitted vs observed values
Dear Wiza[R]ds,
# I have the following experimentally observed data:
csdata <- data.frame(
time=c(0,1,3,9,20),
conc=c(638.697,395.69,199.00,141.58,112.16)
)
# weighting resp means response
wt.MM<- function(resp, time,A1,a1,A2,a2)
{
pred <- A1*exp(-a1*time)+A2*exp(-a2*time)
(resp - pred) / sqrt(pred)
}
# Fit using nls
cs.wt <- nls( ~ wt.MM(conc, time,A1,a1,A2,a2),
2009 Nov 29
1
Plotting observed vs. fitted values
Dear Wiza[R]ds,
I am very grateful to Duncan Murdoch for his assistance with this problem.
His help was invaluable. However, the problem has become a little more
complicated for me. Now, in each plot, I need to plot the observed and
fitted values of a supine and upright posture experiment. Here is what I
have and how far I got.
# tritiated (3H)-Norepinephrine(NE) disappearance from plasma
#
2009 Nov 18
4
Switch Help
Dear Rexperts,
Given,
aar <-function(command) {
switch(command,
{scrn = cat("scrn :Screening","\n")}
{dx = cat("dx :Diagnosis","\n")}
{df = cat("df :Don't Forget","\n")}
)
}
I want to be able to do:
aar("dx") # function does cat("dx :Diagnosis","\n")
2009 Nov 29
3
Plotting observed vs. Predicted values, change of symbols
Dear Wiz[R]ds,
I am deeply grateful for the help from Duncan Murdoch, Gray Calhoun, and
others. We are almost there. For whatever reason, I can't change the symbol
from a circle to a triangle in the upright posture plots. Any ideas? I have
included the problem in full.
# tritiated (3H)-Norepinephrine(NE) disappearance from plasma
# concentrations supine and upright
# supine
datasu <-
2010 Jan 06
1
boxplot help
Dear Rexperts,
I am trying to add a '+' identifying the mean in a boxplot using the
following
sizelist <- split(size, grp)
centers <- boxplot(sizelist, style.bxp = "att", medpch = "o",
ylab = "Prostate Volume (cm3)")
points(centers, unlist(lapply(sizelist, mean)), pch = "+")
But, I get error
Error in xy.coords(x, y) :