Hi!
Replying to my own post: Thanks Duncan for your hint to try the newest
version! I did not think of this in the first place as I'm using the
newest Ubuntu feisty fawn release which got released a couple of days
ago. Nevertheless, the r-cran-hdf5 package is outdated and the newest
one available
Hi,
i have 656 attributes ind INTERVALL_VAR 119 in GROUP
and this morning i'm little confused why the inner loop hang if it
arrive 656th column.
My Task is a t-test and correlation with all columns in INTERVALL_VAR
for all attributes in GROUP.
many thanks regards,
christian
for( k in
sorry,
I don't undersatnd what happens
Annee_O;Id_Essai;Id_Rep;Id_Geno;Id_Cult;Lib_Geno;St_Cult;Id_Par;X_Par;Y_Par;Id_Cara;Surf_O;Val_O;Ori_O;Stade_O;Date_O;Id_Bloc;Id_TrT1
2004;1006003;1;55094;1012988;XF 338/1;;1;1;1;137;;9.4;P;;09/09/2004;1;0
2004;1006003;1;55094;1012988;XF
Petr == Petr Klasterecky [EMAIL PROTECTED]
on Mon, 23 Apr 2007 14:27:03 +0200 writes:
Petr Hi,
Petr you seem to have mixed 2 different things:
Petr 1) changing a working directory - see ?setwd, ?getwd
Petr However, this will NOT load another .Rdata file.
Petr 2) loading
Just a little caveat concerning the use of the clipboard to transfer data from
excel to R:
If I remember right, at least up to Excel 2003, copying data from a worksheet
(Edit/Copy) to the windows clipboard, copies the data in the number format in
which they are seen in Excel. That is, if you
Hi,
i understand my problem , because i overwrite my result's from previous
loop's again and agian, really stupid :-]
regards, christian
Hi,
i have 656 attributes ind INTERVALL_VAR 119 in GROUP
and this morning i'm little confused why the inner loop hang if it
arrive 656th column.
My
You will do better to use e.g. bitmap() and bitmap the output at
600 dpi, if that is what you are aiming at. The default line width
will be 1/96.
You haven't told us what OS 'my PC' is running, but people who say that
usually mean Windows. On the Windows' png device the default line width
elyakhlifi mustapha wrote:
sorry,
I don't undersatnd what happens
Annee_O;Id_Essai;Id_Rep;Id_Geno;Id_Cult;Lib_Geno;St_Cult;Id_Par;X_Par;Y_Par;Id_Cara;Surf_O;Val_O;Ori_O;Stade_O;Date_O;Id_Bloc;Id_TrT1
2004;1006003;1;55094;1012988;XF 338/1;;1;1;1;137;;9.4;P;;09/09/2004;1;0
Dear R-experts,
I imported successful a workspace from Matlab. The information is stored
in the variable data. If I use the command length(data) I get back a
number which corresponds to the number of variables which were imported.
Here it is 82. Now I only want to know the names of the imported
I have a question regarding the samr package.
For a 2 class unpaired problem, with sample 1 of size N1 and sample 2
of size N2, samr computes at most (N1+N2)! permutations of the two
samples (if the user-supplied parameter nperms allows it). However,
there are only (N1+N2)/(N1!*N2!) DISTINCT
Hallo Peter,
thank you. It is exact what I wanted. Now I can modify my program.
Corinna
Von: Peter Konings [mailto:[EMAIL PROTECTED]
Gesendet: Dienstag, 24. April 2007 10:55
An: Schmitt, Corinna
Betreff: Re: [R] information extraction
Hi Corinna,
Dear useRs,
The seqinR package is a library of utilities to retrieve and analyse
biological sequences.
A new version of seqinR, seqinR 1.0-7, has been released on CRAN.
Here is a summary of changes:
o A new *experimental* function extractseqs() to download
sequences thru zlib compressed
Dears helpers,
I have a question not specially about R, but about statistical modelling. This
is the problem, I want to conduct a regression analysis to explain the causes
of students fail in an exam. I have two variables the score obtain at the exam
and the categorical variable coding 0 in
I wanna display some data which there are subsets of a dataframe called don
but there are errors like this
L=as.numeric(levels(factor(don$Id_Cara)))
for(i in L){
+ donC(i)=subset(don, Id_Cara == i, select = c( Id_TrT1, Id_Geno, Id_Rep,
Val_O))
+ donC(i)
+ }
Erreur dans donC(i) = subset(don,
Hi
you shall:
use appropriate subject e.g. how to subset data
use [ ] brackets
initialize object donC before using it for assignment
and maybe try to look to some docummentation how to manipulate R objects
(Paul Johnsons R tips are easily found by Google and they help me a lot
during my first
Are you trying to create a 'list' of the subsets? If so, try:
L=as.numeric(levels(factor(don$Id_Cara)))
donC - list()
for(i in L){
donC[[i]] - subset(don, Id_Cara == i, select = c( Id_TrT1,
Id_Geno, Id_Rep, Val_O))
print(donC[[i]])
}
Or better using the 'names' of the factors
L -
hello
when I entered following command, I got NA values for some catagories.
tapply(slp_jeo2$slp,slp_jeo2$jeo,mean )
999 Ca Cka DCy Jh JKi Kk
14.06665 NA 14.60445 NA NA NA NA NA
KTa KTac Ku
I've rolled up R-2.5.0.tar.gz a short while ago. This is a development
release which contains a number of new features. In particular:
- Object name completion by integration of package 'rcompletion' by
Deepayan Sarkar
- New recommended package 'codetools' by Luke Tierney
- New
Ahmet,
Try the following syntax: tapply(slp_jeo2$slp,slp_jeo2$jeo,mean, na.rm =
T)
Steve
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of temiz
Sent: Tuesday, April 24, 2007 1:22 PM
To: r-help@stat.math.ethz.ch
Subject: [R] problem in tapply command
probably you're looking for
tapply(slp_jeo2$slp, slp_jeo2$jeo, mean, na.rm = TRUE)
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax:
Yannan Jiang wrote:
Hi there,
I am trying to run simulations using R with linear mixed model (lme). There
are two factors in my fixed effect model, educ (treatment and control) and
mth (visit 1, 2, and 3). What I want to obtain is the estimated treatment
difference (treatment -
Thank, Jim to look to my post . I dont explain pass/fail by the exam score but
I want to explain both of them by a set of explanotory variables.
y_{i}=a_{0}+\sum_{i}a_{i}x_{i}
or P(z_{i}=k)=F(a_{0}+\sum_{i}a_{i}x_{i}) z_{i}=1(c_{k}yc_{k+1}
I want to know if there is a
Mustapha,
donC is not a function if it is a list write donC[[i]] or donC[i] for indexing
Justin BEM
Elève Ingénieur Statisticien Economiste
BP 294 Yaoundé.
Tél (00237)9597295.
- Message d'origine
De : elyakhlifi mustapha [EMAIL PROTECTED]
À : R-help@stat.math.ethz.ch
Envoyé le :
hello,
If I wanna export data.frame from R to html have you got some ideas to do this?
___
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch
Use the package R2HTML.
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22
Ohttp://maps.google.com/maps?f=qhl=enq=Curitiba,+Brazillayer=ie=UTF8z=18ll=-25.448315,-49.276916spn=0.002054,0.005407t=kom=1
On 4/24/07, elyakhlifi mustapha [EMAIL PROTECTED] wrote:
hello,
If I
Dear all -
I imported (on a Mac) a big table with 2000 lines:
mydata - read.table(file(/Users/didiw/Desktop/R/all.txt), header = TRUE)
mydata[1:15,]
location Spezies Spec E.MPa. Phi No Trial
1LC PJ 13.27 7.51 1 1
2LC PJ 14.24 6.68 1 1
Dear R-Experts,
I just imported a workspace from Matlab. I know that I can get the names of the
imported variables with names(). It works. The variable ca consists of
several elements. I want to get the names of the elements to handle my output
better. But names(ca) doesn't work. Why? I did
I'm a new user of R-packages.
I'm in need of information with regard to Stochastic Kernel.
I would like to know if I can use R to compute stochastic Kernels according
to Quah methodology.
Best wishes,
Oriana
-
-
Hi
look at aggregate
something like
mymean-aggregate(mydata[,4:5], list( location, Spezies), mean)
mysd-aggregate(mydata[,4:5], list( location, Spezies), sd)
and then cbind(appropriate columns of resulting data frames)
Petr Pikal
[EMAIL PROTECTED]
[EMAIL PROTECTED] napsal dne
Hello,
my problem is that I don't know how long must be an array of double while
calling C from R.
R-Code:
array - c(1,1,1)
save - .C ( Ctest , a = array )
C-Code:
void Ctest ( double *array )
{ ...
array = (double*) realloc ( array , new_number * sizeof(double) ) ;
...
}
The length of
Henrique Dallazuanna wrote:
Use the package R2HTML.
Or xtable.
Jeff
--
http://biostat.mc.vanderbilt.edu/JeffreyHorner
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Is it so much effort to use the search on r-project.org?
elyakhlifi mustapha wrote:
hello,
If I wanna export data.frame from R to html have you got some ideas to do
this?
___
[[alternative HTML
Hi,
I would like to know if someone could clarify a question regarding
license for redistribution of raw data in the datasets package. If
possible, I would like to reuse those data in a package for scipy
(http://www.scipy.org), which is under a BSD license. The package
datasets itself is
ok,
I have problems with write function
F
Id_TrT1 Id_Geno Id_Rep Val_O
30 55094 185
90 55096 187
15 0 55098 192
21 0 55079 176
27 0 55095 192
33 0 55099 198
39 0 55092 1
Sven == Sven Knüppel [EMAIL PROTECTED]
on Tue, 24 Apr 2007 13:53:09 +0200 writes:
Sven Hello,
Sven my problem is that I don't know how long must be an array of double
while calling C from R.
Sven R-Code:
array - c(1,1,1)
save - .C ( Ctest , a = array )
Sven
Hallo,
just try it like a data frame, as I mailed you last time. Here an example:
dat - data.frame(Class=I(Id_TrT1), Levels=I(1), Values=I(2))
new.info - c(Class=Id_Geno, Levels=7 , Values=64208 64209 64210 64211
64212 64213 64214)
dat - rbind(dat, new.info)
dat
new.info - c(Class=
Hi,
I write several xyplot graphics on a source file. When I try to use
source(graphics.R) the source don't work, but if I use
source(graphics.R,echo=T) it work. Generally some commands work without
echo=T, but xyplot dont work. Why it dont work without echo=T? It is possible
to write a
[EMAIL PROTECTED] wrote:
I have data in a two dimensional table. each row of the data adds
upto 100 ( hence they are percentages ). it can be interpreted as
like this A - I are the matches and P - X are the players. Thus
Player P scored 20% of the runs during this season in Match C, 60% in
Hallo,
it might be a problem with the underscore in the names. Look up the syntax
conventions. Here is a example for the write function in my programs perhaps it
might help you:
helpLetters=
cat(\nThe workspace is now imported!, helpLetters ,\n\n)
cat(The names of the variables are the same as
Hi
had you consider to look at help page? How do you expect your written file
to be named?
From help page
write is a wrapper for cat, which gives further details on the format used
cat is useful for producing output in user-defined functions. It converts
its arguments to character strings,
I mis-spoke. It seems I had two collections of functions in the same
directory. One by Colin Goodall, and one by David Scott (I have no
record of where he is/was located). It is the *latter* collection
that does all its work from within Fortran.
I'll have another look at what Colin Goodall
Read FAQ 7.22
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-do-lattice_002ftrellis-graphics-not-work_003f
--sundar
Ronaldo Reis Junior said the following on 4/24/2007 6:38 AM:
Hi,
I write several xyplot graphics on a source file. When I try to use
source(graphics.R) the source don't
Ronaldo Reis Junior chrysopa at gmail.com writes:
I write several xyplot graphics on a source file. When I try to use
source(graphics.R) the source don't work, but if I use
source(graphics.R,echo=T) it work. Generally some commands work without
echo=T, but xyplot dont work. Why it dont
Quoting [EMAIL PROTECTED]:
Hi everyone,
I have a question about limmaGUI.I've just started to use the program
for microarray analysis.My problem is after loading data into Limma
(with GAl file and RNA targets file),I'm not able to create an M box
plot and everytime I'm trying to do that,
Dear r-helpers,
I have been bitten by a cryptic comment in the help page for xYplot().
Here is some code from Madeline Bauer, one of the authors, which is
essentailly the code on the help page for xYplot():
# This example uses the summarize function in Hmisc to
# compute the median and outer
Dear all,
I need to compute the ACF (autocorrel) of an AR6 process, given the values
of its parameters (w1,w2,w3,w4,w5,w6).
First, I notice that there is an error as soon as the sum of the wi equals 1
:
Error in drop(.Call(La_dgesv, a, as.matrix(b), tol, PACKAGE = base)) :
system is
Hi,
I've been trying here to install Rmpi on an SGI IA-64 machine with 64
processors, running SuSE Linux Enterprise Server 9, R 2.4.0 and
lam-mpi 7.1.3. While I've read of similar problems on this list, I
think I've got an entirely new set of error messages to contribute
(see below). I'm not sure
You do need to specify a stationary ARMA process for it to have an acf:
your example is not valid since
Mod(polyroot(c(1, -w)))
[1] 0.6726057 1.6256859 1.6256859 1.7448484 1.7448484 18.4778223
Have you actually looked at the reference on the help page? It is there
to be help-ful.
Like
On Tue, 24 Apr 2007, Hendrik Fuß wrote:
Hi,
I've been trying here to install Rmpi on an SGI IA-64 machine with 64
processors, running SuSE Linux Enterprise Server 9, R 2.4.0 and
lam-mpi 7.1.3. While I've read of similar problems on this list, I
think I've got an entirely new set of error
Hi folks,
Is there a function to show the size of an R object? eg. in Kbytes?
Couple months ago Bendix Carstensen posted this marvelous little function
lls(), which shows all objects in the current workspace by mode, class and
'size'. This is a wonderful enhancement to the build-in ls()
Dear R-users,
I have recently discovered the languageR package and its aovlmer.fnc
which offers a very slick way to test fixed effects hypotheses in lmer
models using mcmc.
However, the R help page for aovlmer.fnc states that the lmer object
must be based a continuous response variable. I
Hi Nicolas --
I think your code is assuming that all nodes have access to the same
set of variables. One solution is to write in a more completely
'functional' style
parNullDistribIntersection - function(g1, g2, perm=1000, cl=cl) {
n1 = nodes(g1)
parSapply(cl,
Aaaah, wonderful. Thanks Ben. lls appears below. Should these functions be a
permanent part of R? Just a thought.
H.
lls
function (pos = 1, pat = ) {
dimx - function(dd) if (is.null(dim(dd)))
length(dd)
else dim(dd)
lll - ls(pos = pos, pat = pat)
cat(formatC(mode, 1,
Hi,
Is 'object.size()' the function your are looking for ?
Best,
François
-Message d'origine-
De : [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] De la part de Horace Tso
Envoyé : mardi 24 avril 2007 17:59
À : r-help@stat.math.ethz.ch
Objet : [R] Size of an object in workspace
Hi folks,
Try the following and look at what they return:
str(ca)
dimnames(ca)
-- Tony Plate
[EMAIL PROTECTED] wrote:
Dear R-Experts,
I just imported a workspace from Matlab. I know that I can get the names of
the imported variables with names(). It works. The variable ca consists of
several
Hi,
I am looking for documentation, reference guides, etc. that explain the
output of functions... For example using cor.test(, method=pearson)
with Pearson's corr coeff the output is:
Pearson's product-moment correlation
data: a and b
t = 0.2878, df = 14, p-value = 0.
alternative
I draw dotplot using the following code:
sd.dotplot-dotplot(data.47.nmr$om_sd ~ as.factor(data.47.nmr$position)
|as.factor(data.47.nmr$pr), data = data.47.nmr,layout=c(1,1),
xlab=Position, xlim=range(data.47.nmr$position),ylab=Sd, main=Changes
of omega angle in different positions,
scales =
See ll() in R.oo (that is two L:s), e.g.
ll()
member data.class dimension objectSize
1 author character 1120
2 myfunc function NULL512
3 x matrix c(3,11)248
4 y array c(5,7,1)264
ll() is quite flexible so you can create
I have been playing around with sparse matrices in the Matrix
package, in particularly with the Cholesky factorization of matrices
of class dsCMatrix. And BTW, what a fantastic package.
My problem is that I have to carry out repeated Cholesky
factorization of a spares symmetric matrices, say
Hi,
I searched the archives and did not find a good solution to that.
assume I have 10 sets and I want to have the common character elements of them.
how could i do that?
--
Weiwei Shi, Ph.D
Research Scientist
GeneGO, Inc.
Did you always know?
No, I did not. But I believed...
---Matrix III
Hi,
If you look at the documentation for the function you are interested
in, in this case ?cor.test, it will generally give you an explanation
of the return values (often brief, and not too helpful if you aren't
already familiar with the test), but also one or more references
that you can turn to
On 24/04/07, Sarah Goslee [EMAIL PROTECTED] wrote:
Hi,
If you look at the documentation for the function you are interested
in, in this case ?cor.test, it will generally give you an explanation
of the return values (often brief, and not too helpful if you aren't
already familiar with the
Hi,
I am trying to avid the somewhat costly startup overhead of launching
a separate R executable for each client request on Linux.
My current architecture is such that My Java client explicitly calls
R in batch mode and passes it certain parameters. The initital
startup takes almost 10 seconds
Dear R-help,
This is about randomForest's handling of NA and NaNs in test set data.
Currently, if the test set data contains an NA or NaN then
predict.randomForest will skip that row in the output.
I would like to change that behavior to outputting an NA.
Can this be done with flags to
On Tue, Apr 24, 2007 at 02:15:43PM -0400, Markus Loecher wrote:
Hi,
I am trying to avid the somewhat costly startup overhead of launching
a separate R executable for each client request on Linux.
My current architecture is such that My Java client explicitly calls
R in batch mode and passes
Hi Markus,
Take a look at RSOAP (http://rsoap.sf.net). It was designed for
handling concurrent client connections to R, and minimizes the per-
connection startup time by pre-starting R and forking of processes as
requests come in. Each client can maintain a stateful connection, if
assume t2 is a list of size 11 and each element is a vector of characters.
the following codes can get what I wanted but I assume there might be
a one-line code for that:
t3 - t2[[1]]
for ( i in 2:11){
t3 - intersect(t2[[i]], t3)
}
or there is no such apply?
On 4/24/07, Weiwei Shi
RWebServices
http://wiki.fhcrc.org/caBioc/
offers a more structured approach to this -- map R functions and data
classes to their Java representation, expose Java as a web service,
service requests using persistent R workers.
Martin
Markus Loecher [EMAIL PROTECTED] writes:
Hi,
I am trying
On Tue, 24 Apr 2007, Weiwei Shi wrote:
Hi,
I searched the archives and did not find a good solution to that.
assume I have 10 sets and I want to have the common character elements of
them.
how could i do that?
list.of.sets - lapply(1:10,function(x) sample(letters,20)) # for example
Dear Weiwei Shi,
How about using recursion?
intersection - function(x, y, ...){
+ if (missing(...)) intersect(x, y)
+ else intersect(x, intersection(y, ...))
+ }
a - letters[1:4]
b - letters[2:5]
c - letters[3:6]
d - letters[4:7]
e - letters[5:8]
intersection(a, b)
[1] b c d
I don't think there's that sort of apply-reduce function in R, but for
this problem, the last line below happens to be a one-liner:
set.seed(1)
x - lapply(1:10, function(i) sample(letters, 20))
table(unlist(x))
a b c d e f g h i j k l m n o p q r s t u v w x
y
On 4/24/07, Aimin Yan [EMAIL PROTECTED] wrote:
I draw dotplot using the following code:
sd.dotplot-dotplot(data.47.nmr$om_sd ~ as.factor(data.47.nmr$position)
|as.factor(data.47.nmr$pr), data = data.47.nmr,layout=c(1,1),
xlab=Position, xlim=range(data.47.nmr$position),ylab=Sd, main=Changes
On 4/24/07, Weiwei Shi [EMAIL PROTECTED] wrote:
assume t2 is a list of size 11 and each element is a vector of characters.
the following codes can get what I wanted but I assume there might be
a one-line code for that:
t3 - t2[[1]]
for ( i in 2:11){
t3 - intersect(t2[[i]], t3)
}
I had a similar solution by using frequency but having more codes :(
I also like the recursive idea : I initially tried to use rapply,
however, which can only take one-arg function.
thanks, everyone.
-w
On 4/24/07, hadley wickham [EMAIL PROTECTED] wrote:
On 4/24/07, Tony Plate [EMAIL
On 4/24/07, Tony Plate [EMAIL PROTECTED] wrote:
I don't think there's that sort of apply-reduce function in R, but for
this problem, the last line below happens to be a one-liner:
Only if you have character data though:
x - lapply(1:10, function(i) sample(20, 15))
Hi,
I have my files and R session on a portable hard drive to be able to
work across computers since I have some extremely large files. I have
just started doing this and generally it works as I expect. But after a
while of using my R session, certain types of calls to open connections
on the
you could try something like the following:
t2 - lapply(1:11, function(i) c(a, sample(letters[1:5],
sample(10, 1), TRUE), b))
unq.vals - unique(unlist(t2))
ind - rowSums(sapply(t2, %in%, x = unq.vals)) == length(t2)
unq.vals[ind]
I hope it helps.
Best,
Dimitris
Dimitris
hadley wickham wrote:
On 4/24/07, Weiwei Shi [EMAIL PROTECTED] wrote:
assume t2 is a list of size 11 and each element is a vector of characters.
the following codes can get what I wanted but I assume there might be
a one-line code for that:
t3 - t2[[1]]
for ( i in 2:11){
t3 -
Hi,
My name is Katie and I was wondering if you could help me with my problem. I
am trying to write a function in R that computes the statistics (mean,
standard error, confidence intervals) for stratified samples. I am not that
familiar with R and I am having dificulties setting this function up.
Have you seen Pinheiro and Bates (2000) that lays out the nlme
package? It is very helpful.
Hank
On Apr 14, 2007, at 8:02 PM, Kyle. wrote:
You probably can do this with lme function, but I don't know that for
sure. aov (included in the stats package), with a call to the
Error function how I
On 4/24/07, Gardar Johannesson [EMAIL PROTECTED] wrote:
I have been playing around with sparse matrices in the Matrix
package, in particularly with the Cholesky factorization of matrices
of class dsCMatrix. And BTW, what a fantastic package.
My problem is that I have to carry out repeated
have u seen ?merge ?
- Original Message
From: Weiwei Shi [EMAIL PROTECTED]
To: R Help R-help@stat.math.ethz.ch
Sent: Tuesday, April 24, 2007 10:55:51 PM
Subject: [R] intersect more than two sets
Hi,
I searched the archives and did not find a good solution to that.
assume I have 10
Hi everybody,
I work with data with following pattern
comm
Date Value
1 4/10/2007 361.2
2 4/11/2007 370.1
3 4/12/2007 357.2
4 4/13/2007 362.3
5 4/16/2007 363.5
6 4/17/2007 368.7
7 4/18/2007
Dear all,
I am a little bit puzzled by the way round() works.
Consider the following code
a-123456.3678
round(a,digits=10)
[1] 123456.4
I would expect the outcome to be something like 123456.3678 or
123456.368, instead the computer gives me 123456.4 no matter how large
the digits are.
Can
library(gdata)
ll()
On 4/24/07, Horace Tso [EMAIL PROTECTED] wrote:
Hi folks,
Is there a function to show the size of an R object? eg. in Kbytes?
Couple months ago Bendix Carstensen posted this marvelous little function
lls(), which shows all objects in the current workspace by mode,
Hello,
I have a Problem to make Log-Returns of the dataset EuStockMarkets.
Is there any function which could calculate it for me?
data(EuStockMarkets)
Thanks!
Alin Soare
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch
diff(log(x))
And this should be made into a function.
H.
Soare Marcian-Alin [EMAIL PROTECTED] 4/24/2007 2:17 PM
Hello,
I have a Problem to make Log-Returns of the dataset EuStockMarkets.
Is there any function which could calculate it for me?
data(EuStockMarkets)
Thanks!
Alin Soare
Dear Peter, Hadley, et al.,
Just to bring it back around, here's a recursive Fold():
a - letters[1:4]
b - letters[2:5]
c - letters[3:6]
d - letters[4:7]
e - letters[5:8]
Fold - function(f, x, y, ...){
+ if (missing(...)) f(x, y)
+ else f(x, Fold(f, y, ...))
+ }
Fold(intersect,
'print' by default only prints 7 digits. See below:
round(a, digits=10)
[1] 123456.4
print(round(a, digits=10), digits=10)
[1] 123456.3678
On 4/24/07, Robert Duval [EMAIL PROTECTED] wrote:
Dear all,
I am a little bit puzzled by the way round() works.
Consider the following code
On Tue, 24 Apr 2007, Tomas Mikoviny wrote:
Hi everybody,
I work with data with following pattern
comm
Date Value
1 4/10/2007 361.2
2 4/11/2007 370.1
3 4/12/2007 357.2
4 4/13/2007 362.3
5 4/16/2007 363.5
6
see
?getReturns in fMultivar package (part of Rmetrics).
the function has an option between continuous(log returns) or discrete
calculation.
I used it with timeSeries objects. works well.
good luck.
A.
- Original Message -
From: Soare Marcian-Alin [EMAIL PROTECTED]
To:
I am really struggling with this question.
Three students take part of an experiment.
student smoking non-smokingcancer
110.5 7.5 6.5
2 9.5 6.5 8.4
3 8.5 7.2 5.5
the proper
Hi Keti
Before reinventing the wheel from scratch you might want to take a
look at the survey package
http://faculty.washington.edu/tlumley/survey/
best
robert
On 4/24/07, Keti Cuko [EMAIL PROTECTED] wrote:
Hi,
My name is Katie and I was wondering if you could help me with my problem. I
am
Hello,
I want to simulate 100 values of the ARMA Process with this function:
x[i] = 0.5 * x[i-1] + 0.2 * x[i-2] + x[i] + 0.9 * x[i-1] + 0.2 * x[i-2] +
0.3 * x[i-3]
which possibilities do I have?
Alin Soare
[[alternative HTML version deleted]]
Hello
Blind toy_Car toy_truck toy_boat
1 6.3 7.5 5.4
2 3.4 8.1 6.1
3 2.2 4.4 5.1
How do we calculate the F-statistic in R.
Any help is really appreciated.
--
View this message in context:
Dear List,
Trying to explain my situation as simply as possible
for me:
I am running a series of iteration on coxph model on
simulated data (newly generated data on each iteration
to run under coxph; in my example below- sim.fr is the
generated data). However, sometimes i get warning
messages
Hello all,
I know this is a pretty easy question but I can't find it in S poetry or
R help.
How can I make a negative number positive. Such as
-5 to be +5
I tried +(-5), but that didn't work.
So no, I don't mean taking a -5^2 just to get a positive number.
This is in a function so it's not just
On Tue, 2007-04-24 at 16:46 -0700, H. Paul Benton wrote:
Hello all,
I know this is a pretty easy question but I can't find it in S poetry or
R help.
How can I make a negative number positive. Such as
-5 to be +5
I tried +(-5), but that didn't work.
So no, I don't mean taking a -5^2
?try
Wrap each iteration in a try() call
Also ?tryCatch if you want to get fancy -- and can understand the rather
arcane docs.
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Mohammad Ehsanul
Karim
Sent:
Thanks
Research Technician
Mass Spectrometry
o The
/
o Scripps
\
o Research
/
o Institute
Rajarshi Guha wrote:
On Tue, 2007-04-24 at 16:46 -0700, H. Paul Benton wrote:
Hello all,
I know this is a pretty easy question but I can't find it in S poetry or
R help.
How can I
1 - 100 of 112 matches
Mail list logo
