Displaying 20 results from an estimated 500 matches similar to: "Importing data from text file with mixed format"
2008 Apr 28
2
F values from a Repeated Measures aov
Hi Folks,
I have repeated measures for data on association time (under 2
acoustic condtions) in male and female frogs as they grow to adulthood
(6 timepoints). Thus, two within-subject variables (Acoustic
Condition: 2 levels, Timepoint: 6 levels) and one between-subject
variable (Sex:male or female).
I am pretty sure my distributions depart from normality but I would
first like to simply run a
2009 Sep 02
2
Average over data sets
Hello,
I have a number of files output1.dat, output2.dat, ... , output20.dat,
each of which monitors several variables over a fixed number of
timepoints. From this I want to create a data frame which contains the
mean value between all files, for each timepoint and each variable.
The code below works, but it seems like I should be able to do the
second part without a for loop. I played
2010 Sep 16
1
ANOVA - more sophisticated contrasts
dear list,
i am using a multifactorial design with two treatments (factor A: drugs,
three levels; factor B: theraphy, two levels) and a time factor (three
levels, different timepoint). hypothetically, i measured the same subjects
for all treatements and timepoints, so its a repeated measurement design.
now i ran an anova in R and also some Tukey post-hoc tests using glht. but
what i am actually
2010 Sep 19
1
boyplots nearly identical but still highly significant effect?
dear list,
i am running a within-design ANOVA with 4 factors (4,4,2 and 2 levels each).
the last one is a time factor comprising two different treatment timepoints.
i fit a mixed-effects model using lme and apply the anova function to the
outcome. according to this analysis, there are highly significant main
effect on the first and the time factor. i then checked the boxplots for the
two 4-level
2009 Oct 08
3
how to convert list into a vector
Hi ,
I want to convert a my list:
> myList
[[1]]
[1] 1 2 3 4 5
[[2]]
[1] 6 7 8
[[3]]
[1] 9 10 11
into something like c(1,2,3,4,5,6,7,8,9,10,11)
I realized that this is possible by
> c(do.call("cbind",myList))
But only when list[[1]] through list[[3] have equal length of vectors within
them
Any ideas?
Thank you for time,
-Daniel
--
View this message in context:
2007 Nov 20
1
Vectorization/Speed Problem
Hi,
I cannot find a 'vectorized' solution to this 'for loop' kind of problem.
Do you see a vectorized, fast-running solution?
Objective:
Take the value of X at each timepoint and calculate the corresponding value
of Y. Leading 0's and all 1's for X are assigned to Y; otherwise Y is
incremented by the number of 0's adjacent to the last 1. The frequency and
2013 Nov 08
1
how to derive true surface area from `computeContour3d' (misc3d package) -- follow up
regarding my previous mail for this topic, I have in the meantime
identified my misconception.
actually, `computeContour3d' returns the vertices just fine in the correct
coordinate frame.
the misconception was caused basically by assuming that the `level'
argument was a fractional
threshold relative to the maximum of the array. so I believed that the
rendered cube actually
is
2010 Jul 08
0
0 exit status on packages update try
I am not entirely sure what is going on. I am sure that I am missing
something. maptools won't detach... Thanks for all of your help.
R 2.11.1
OS Ubuntu 10.04
R --vanilla
install.packages("sp", dep=TRUE)
Loading required package: reshape
Loading required package: plyr
Loading required package: grid
Loading required package: proto
This is vegan 1.17-3
Loading required package:
2010 Jun 10
1
do faster ANOVAS
Dear all R users,
I want to realize 800 000 ANOVAS and to store Sum of Squares of the effects. Here is an extract of my table data
Product attribute subject rep t1 t2 t3 … t101
P1 A1 S1 R1 1 0 0 … 1
I want to realize 1 ANOVA per timepoint and per attribute, there are 101 timepoints and 8 attributes so I want to realize 808 ANOVAS. This will be an ANOVA with two factors :
Here is one example:
2013 Sep 05
2
Problem installing "sp" package
Hi,
I get an error message (see below) when trying to install the sp package
in R version 2.15.2 running under Ubuntu 12.04 LTS.
I have run sudo apt-get build-dep r-base, but to no help.
Any help will be greatly appreciated.
Terje
> install.packages("sp")
Installing package(s) into ?/home/terjej/R/x86_64-pc-linux-gnu-library/2.15?
(as ?lib? is unspecified)
trying URL
2008 Mar 08
1
analysing mixed effects/poisson/correlated data
I am attempting to model data with the following variables:
timepoint - n=48, monthly over 4 years
hospital - n=3
opsn1 - no of outcomes
total.patients
skillmixpc - skill mix percentage
nurse.hours.per.day
Aims
To determine if skillmix affects rate (i.e. no.of.outcomes/total.patients).
To determine if nurse.hours.per.day affects rate.
To determine if rates vary between
2008 Feb 05
1
Extracting level-1 variance from lmer()
All,
How does one extract the level-1 variance from a model fit via lmer()?
In the code below the level-2 variance component may be obtained via
subscripting, but what about the level-1 variance, viz., the 3.215072 term?
(actually this term squared) Didn't see anything in the archives on this.
Cheers,
David
> fm <- lmer( dv ~ time.num*drug + (1 | Patient.new), data=dat.new )
2008 Sep 11
1
plot of all.effects object
All,
I'm trying to plot an all.effects() object, as shown in the help for
all.effects and also Crawley's R book (p.178, 2007). The data has a repeated
measures structure, but I'm using all.effects for the simple lm() fit here.
Below is a reproducible example that yields the error message.
fm.ex = lm(dv ~ time.num*drug*X, data = dat.new)
fm.effects = all.effects(fm.ex, xlevels =
2007 Jul 09
2
ANOVA: Does a Between-Subjects Factor belong in the Error Term?
I am executing a Repeated Measures Analysis of Variance with 1 DV (LOCOMOTOR
RESPONSE), 2 Within-Subjects Factors (AGE, ACOUSTIC CONDITION), and 1
Between-Subjects Factor (SEX).
Does anyone know whether the between-subjects factor (SEX) belongs in the
Error Term of the aov or not? And if it does belong, where in the Error Term
does it go? The 3 possible scenarios are listed below:
e.g.,
1.
2005 Feb 07
1
Incorrect behavior for ordering timepoints in "reshape" (PR#7669)
Full_Name: Dav Clark
Version: 2.0.1
OS: OS X 10.3
Submission from: (NULL) (128.122.87.35)
When the timepoints that reshape uses (in direction="long") are negative or
fractional, the time label is assigned incorrectly. It is easier to give an
example than to describe the problem abstractly:
Assume you have a data.frame header with values related to peri-stimulus time
like this:
2004 May 05
1
Granule Pos of start of page...
OK... i've come across a problem trying to get the granule pos of the start
of the page... it's not so crucial with single stream ogg files... but now
that i have theora+vorbis in a file, i'm finding that when i seek to a
position, i have no way to determine the relative offsets of the different
streams at the new seek point and hence the av is out of sync.
So given a page, is it
2013 Apr 01
1
Help Please, ggplot2
library(ggplot2)
a<- read.table("data", header=T)
b = na.omit(a)
ggplot(data=b) + geom_line(aes(x=timepoint, y=value,group=sample, colour= factor(sample))) +? geom_point(aes(x=timepoint, y=value, group=s
ample)) + facet_wrap(~bio, scales = "free",ncol = 5) +theme_bw() + opts(legend.direction = "horizontal",??? legend.position = "top",????
2009 Jun 08
1
Looking for easy way to normalize data by groups
Hi,
i do have a dataframe representing data from a repeated experiment. PID
is a subject identifier, Time are timepoints in an experiment which was
repeated twice. For each subject and all three timepoints there are 2
sets of four values.
df <- data.frame(PID = c(rep("A", 12), rep("B", 12), rep("C", 12)),
Time = rep(c(0, 0, 0, 0, 30, 30, 30,
2008 May 28
1
heatmap-changing column or row names
Dear R Community,
I am trying to create an heatmap for the following set of data:
##example of data matrix
o4
V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18
green 27 28 29 29 28 28 26 25 25 23 23 22 22 21 21 22 22 22
yellow 6 8 8 7 7 7 6 6 6 7 7 7 6 6 6 6 6 5
red 15 15 15 15 15 15 14 13 12 11 12 10 9 8 7 6 8 9
pink
2012 Feb 20
1
Reporting Kaplan-Meier / Cox-Proportional Hazard Standard Error, km.coxph.plot, survfit.object
What is the best way to report the standard error when publishing
Kaplan-Meier plots? In my field (Vascular Surgery), practitioners
loosely refer to the "10% error" cutoff as the point at which to stop
drawing the KM curve. I am interpreting this as the *standard error
of the cumulative hazard*, although I'm having a difficult time
finding some guidelines about this (perhaps I am