Displaying 20 results from an estimated 4000 matches similar to: "wish list: generalized apply"
2007 Jan 19
4
Newbie question: Statistical functions (e.g., mean, sd) in a "transform" statement?
Greetings listeRs -
Given a data frame such as
times
time1 time2 time3 time4
1 70.408543 48.92378 7.399605 95.93050
2 17.231940 27.48530 82.962916 10.20619
3 20.279220 10.33575 66.209290 30.71846
4 NA 53.31993 12.398237 35.65782
5 9.295965 NA 48.929201 NA
6 63.966518 42.16304 1.777342 NA
one can use "transform" to
2007 Dec 16
4
improving a bar graph
Hello,
Below is the code for a basic bar graph. I was seeking advice
regarding the following:
(a) For each time period there are values from 16 people. How I can
change the colour value so that each person has a different colour,
which recurs across each of the three graphs/tie epriods?
(b) I have seen much more sophisticated examples using lattice (e.g
each person has a separate
2004 Jul 04
2
Random intercept model with time-dependent covariates, results different from SAS
Dear list-members
I am new to R and a statistics beginner. I really like the ease with which I can
extract and manipulate data in R, and would like to use it primarily. I've
been learning by checking analyses that have already been run in SAS.
In an experiment with Y being a response variable, and group a 2-level
between-subject factor, and time a 5-level within-subject factor. 2
2011 Oct 06
2
Titles changing when a plot is redrawn
I ran into a problem with titles on graphs. I wanted a graph with
multiple subplots, with each having a title that involved both
a Greek letter and an identifier for each graph. Below is a
simplified version of code to do this. The graph appears fine,
with the first graph having "i=1" in the title, and the second
graph having "i=2" in the title. However, when I resize the
2012 Nov 30
1
help on "stacking" matrices up
Dear All,
#I have the following code
Dose<-1000
Tinf <-0.5
INTERVAL <-8
TIME8 <-matrix(c((0*INTERVAL):(1*INTERVAL)))
TIME7 <-matrix(c((0*INTERVAL):(2*INTERVAL)))
TIME6 <-matrix(c((0*INTERVAL):(3*INTERVAL)))
TIME5 <-matrix(c((0*INTERVAL):(4*INTERVAL)))
TIME4 <-matrix(c((0*INTERVAL):(5*INTERVAL)))
TIME3 <-matrix(c((0*INTERVAL):(6*INTERVAL)))
TIME2
2010 May 20
1
Strange behaviour when using diff with POSIXt and POSIXlt objects
Dear list,
I´m calculating time differences between series of time stamps and I noticed
something odd:
If I do this...
> time1=strptime("2009 05 31 22 57 00",format="%Y %m %d %H %M")
> time2=strptime("2009 05 31 23 07 00",format="%Y %m %d %H %M")
>
> diff(c(time1,time2),units="mins")
Time difference of 10 mins
.. I get the correct
2016 Dec 08
2
require(..., quietly=TRUE) does not suppress warning
Hi,
The `quietly` argument of `require` is documented as follows:
quietly: a logical. If ?TRUE?, no message confirming package
attaching is printed, and most often, no errors/warnings are
printed if package attaching fails.
However:
> require(foo, quietly=TRUE)
Warning message:
In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, :
2012 Jun 24
2
Win 64 package build - ERROR: loading failed for 'x64'
I have developed an R package that works under Win32, but when I attempt to build it on Win64,
I get ERROR: loading failed for 'x64'
More precisely, I developed and tested the package under Win32 and it works. But when I move
to a 64 bit Windows 7 (Home Premium) system, and attempt to build both 32 bit and 64 bit packages,
the 32 bit package seems to build, but the 64 bit build
2012 Apr 15
2
xyplot type="l"
Probably a stupidly simple question, but I wouldn't know how to google it:
xyplot(neuro ~ time | UserID, data=data_sub)
creates a proper plot.
However, if I add
type = "l"
the lines do not go first through time1, then time2, then time3 etc but in
about 50% of all subjects the lines go through points seemingly random
(e.g. from 1 to 4 to 2 to 5 to 3).
The lines always start at time
2005 Nov 15
1
Repeates Measures MANOVA for Time*Treatment Interactions
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
2016 Dec 08
0
wish list: generalized apply
> On Dec 8, 2016, at 12:09 PM, John P. Nolan <jpnolan at american.edu> wrote:
>
> Dear All,
>
> I regularly want to "apply" some function to an array in a way that the arguments to the user function depend on the index on which the apply is working. A simple example is:
>
> A <- array( runif(160), dim=c(5,4,8) )
> x <- matrix( runif(32), nrow=4,
2008 Jun 06
1
calling a C function with a struct
I am trying to call a precompiled C function that uses a struct as one of
it's arguments.
I could write a wrapper function in C, but I was hoping there is some way
to
pack fields into an array of type raw that could be passed directly to the
function.
Here is some more detail. The C struct is simple, but has mixed types:
struct STRUCT1 {
long type;
long nx;
double
2016 Dec 08
1
require(..., quietly=TRUE) does not suppress warning
Well, I'm getting a warning (not an error) when the package doesn't exist.
I interpreted "most often" to mean that suppressing warnings/errors is why you'd most often use this argument, as most packages don't emit startup messages.
And technically there isn't a problem with attaching the package, since we don't even try to attach packages that don't exist.
2011 May 01
1
Mean/SD of Each Position in Table
I have 100+ .csv files which have the basic format:
> test
X Substance1 Substance2 Substance3 Substance4 Substance5
1 Time1 10 0 0 0 0
2 Time2 9 5 0 0 0
3 Time3 8 10 1 0 0
4 Time4 7 20 2 1 0
5 Time5
2010 Dec 06
3
0.5 != integrate(dnorm,0,20000) = 0
Hello:
The example "integrate(dnorm,0,20000)" says it "fails on many
systems". I just got 0 from it, when I should have gotten either an
error or something close to 0.5. I got this with R 2.12.0 under both
Windows Vista_x64 and Linux (Fedora 13); see the results from Windows
below. I thought you might want to know.
Thanks for all your work in creating
2010 Jun 17
2
constrOptim( ): conflict between help page and code
There is a contradiction between what the help page says and what constrOptim actually
does with the constraints. The issue is what happens on the boundary.
The help page says
The feasible region is defined by ?ui %*% theta - ci >= 0?,
but the R code for constrOptim reads
if (any(ui %*% theta - ci <= 0))
stop("initial value not feasible")
The following example
2005 Aug 10
2
Treatment-response analysis along time
Dear R people,
I wonder if you could give me a hand with some of my data. I have a very
typical analysis in biology, however it is difficult for me to find the
right way to analyse. I had a group of animals, I gave them a treatment, and
I measure a variable along time -one??s per day- along 5 days,for
example(fake data):
Animals Time1 Time2 Time3 Time4
1 1 5 3
2012 Mar 28
3
Connect lines in a dot plot on a subject-by-subject basis
I am trying to plot where data points from a give subject are connected by a line. Each subject is represented by a single row of data. Each subject can have
up to five observations. The first five columns of mydata give the time of observation, columns 6-10 give the values at each time point. Some subjects have
all data, some are missing values.
The code I wrote to draw the plot is listed below.
2007 May 13
2
Some questions on repeated measures (M)ANOVA & mixed models with lme4
Dear R Masters,
I'm an anesthesiology resident trying to make his way through basic
statistics. Recently I have been confronted with longitudinal data in
a treatment vs. control analysis. My dataframe is in the form of:
subj | group | baseline | time | outcome (long)
or
subj | group | baseline | time1 |...| time6 | (wide)
The measured variable is a continuous one. The null hypothesis in
2010 Jul 16
1
Nested if help
Hello,
I am trying to find a direct way to write a nested if of sorts to find data
for a specific time range for a specific day (across a range of days) and
have exhausted my abilities with the manuals I have at hand. I have a good
deal of data of this approximate form:
day time price
1 1am 5
1 2am 7
1 3am 9
1 4am 12
2 1am 5
2 2am 7
2