date:20050330

[R] Recall() and sapply()

2005-03-30 Thread Robin Hankin

Hi.
I'm having difficulty following the advice given in help(Recall).  
Consider the two
following toy functions:

f1 - function(n){
  if(length(n)1){return(sapply(n,f1))}
  matrix(n,n,n)
}
f2 - function(n){
  if(length(n)1){return(sapply(n,Recall))}
  matrix(n,n,n)
}
f1() works as desired (that is, f(1:3), say, gives me a three element 
list whose i-th element
is an i-by-i matrix whose elements are all i).

But f2() doesn't.
How do I modify either function to use Recall()?  What exactly is 
Recall() calling here?


--
Robin Hankin
Uncertainty Analyst
Southampton Oceanography Centre
European Way, Southampton SO14 3ZH, UK
 tel  023-8059-7743
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] From FAQ 7.21 to a command like apply(sapply(list(f1,f2,f3),is.na),2,sum)

2005-03-30 Thread Heinz Tuechler

At 16:13 29.03.2005 -0800, Thomas Lumley wrote:
On Wed, 30 Mar 2005, Heinz Tuechler wrote:

 Dear all,

 Last December there was a thread regarding the famous FAQ 7.21 How can I
 turn a string into a variable? and asking what people want to do with
 these strings.
 My, certainly trivial application would be as follows:
 Assume I have a data.frame containing besides others also the columns f1,
 f2, ..., fn and I want to create a command like:
 apply(sapply(list(f1,f2,f3),is.na),2,sum)
 or
 summary(cbind(f1,f2,f3))

 Can I start from paste('f',1:3,sep='') to arrive at the abovementioned
 command?

No parse,as.name or other complications needed. It's all just indexing. 
Suppose your data frame is called dd

fs-paste('f',1:3,sep='')
apply(sapply(dd[,fs],is.na),2,sum)
summary(dd[,fs])

   -thomas


Thank you, Thomas, for your answer. I was curious if there was a simple way
to do this without referring to the data.frame, so that the resulting
command would correspond in its effect exactly to the abovementioned examples.
It's not urgent, but I will try further.

Many thanks

Heinz

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] From FAQ 7.21 to a command likeapply(sapply(list(f1,f2,f3),is.na),2,sum)

2005-03-30 Thread Heinz Tuechler

At 09:32 30.03.2005 +0200, Heinz Tuechler wrote:
At 16:13 29.03.2005 -0800, Thomas Lumley wrote:
On Wed, 30 Mar 2005, Heinz Tuechler wrote:

 Dear all,

 Last December there was a thread regarding the famous FAQ 7.21 How can I
 turn a string into a variable? and asking what people want to do with
 these strings.
 My, certainly trivial application would be as follows:
 Assume I have a data.frame containing besides others also the columns f1,
 f2, ..., fn and I want to create a command like:
 apply(sapply(list(f1,f2,f3),is.na),2,sum)
 or
 summary(cbind(f1,f2,f3))

 Can I start from paste('f',1:3,sep='') to arrive at the abovementioned
 command?

No parse,as.name or other complications needed. It's all just indexing. 
Suppose your data frame is called dd

fs-paste('f',1:3,sep='')
apply(sapply(dd[,fs],is.na),2,sum)
summary(dd[,fs])

  -thomas


Thank you, Thomas, for your answer. I was curious if there was a simple way
to do this without referring to the data.frame, so that the resulting
command would correspond in its effect exactly to the abovementioned
examples.
It's not urgent, but I will try further.

Many thanks

Heinz

Continuation:
Maybe not an elegant solution, but it seems to work:
apply(sapply(eval(parse(text=paste('list(',paste('f',1:3,sep='',
collapse=','),')'))) ,is.na),2,sum)
What I missed in my earlier attempts was collapse=','.

Heinz

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] locfit

2005-03-30 Thread . .

Hallo all,
can somebody tell me if is possible use the locfit function like a 
semiparametric model. I need to  use two explicative variables and one of them 
is a dummy variable.
Thanks Michele


-

 Giochi, Rubrica Scaricalo ora! 
[[alternative HTML version deleted]]

__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] Habituation model

2005-03-30 Thread Laurent Fanchon

Dear all,

I am looking at habituation of dogs trotting on a treadmill.
I record the ground reaction force and I analyze it with several discrete variables (maximum, minimum,...)
For each variable, I get between 40 and 50 data per sample.

I record data at time 1 min, 2 min, and 4 min a day, and I have 4 days of measurement (one day a week).
That means I have 12 samples : Day1_Min1, Day1_Min2, Day1_Min4, Day2_Min1,...

Furthermore, I have 28 dogs to analyze.

My questions is : when can I consider the data stabilized?
The problem is I am studying habituation with several sessions.
It seems logical that values from each first minute could be very dissimilar to others.
It is not a fully linear training system.
I have already done some ANOVAs with Day and Min factors. I found a significant effect.
It is quite logical because the dog is learning. The question is then when does it stop learning? or more precisely when is it trained enough to be analyzed?

I could do comparisons among all samples with Student test, but it is surely a simple approach.

I can presuppose the maximal allowed variability for each variable : in the region of 5%.

I am really new to both R and stats so if these questions are very simple and I am just missing something, suggestions about good texts or examples on R would be great.
I am generating data with Scilab, I have a single matrix corresponding to each dog. But I can change it if needed.

Any help would be greatly appreciated
Thanks
Laurent Fanchon
DVM, MS
Ecole Nationale Veterinaire d'Alfort
France
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

62 matches

Mail list logo