Hi R,
Can we download Reuters (3000 Xtra) data from R? Does ODBC package help
me in this? Or otherwise, is there a way to extract daily closing prices
data of Reuters from R?
Thank you very much,
Shubha
This e-mail may contain confidential and/or privileged i...{{dropped:13}}
Joe Trubisz wrote:
Hi...
Is this possible in R?
I have 2-sets of data, that were collected simultaneously using
2-different data acquisition schemes.
The x-values are the same for both.
The y-values have different ranges (16.4-37.5 using one method,
557-634 using another).
In theory, if
CE.KA wrote:
Hi R users,
Is there a Package in R to
- set up a questionnaire?
- enter data?
Hi CE.KA,
I don't know, but I have written a general purpose questionnaire program
in Tcl-Tk that will administer the questionnaire and record the
responses. If you don't find anything, let me know.
John Sorkin wrote:
R 2.8.0
windows XP
I would like to divide the rows of data frame into five groups and then get the
mean of one column within the five groups. I have accomplished this using the
code below, but I hope there is an easier way, i.e. some function that I can
call
# create five
http://idisk.mac.com/jdeleeuw-Public/jacobi
This is paper/software for various techniques based on Jacobi plane
rotations.
There is R code for
-- classical cyclical Jacobi Eigen diagonalization
-- Jacobi-based SVD diagonalization
-- approximate simultaneous diagonalization of symmetric
Rounding can do no good because
round(8.8,1)-round(7.8,1)1
# still TRUE
round(8.8)-round(7.7)1
# FALSE
What you might do is compute a-b-1 and compare it to a very small number:
(8.8-7.8-1) 1e-10
# TRUE
K
On Wed, Dec 10, 2008 at 11:47 AM, emma jane [EMAIL PROTECTED] wrote:
Thanks Greg,
Hello R-users.
I want to know if there is in R the distribution and density function of the
Coxian. If there isn't, It will help if someone has written this functions.
Thanks
Borja
[[alternative HTML version deleted]]
__
R-help@r-project.org
Dear all,
I was trying to rotate a plot using R. I tried most of the examples offered
here in this forum, but for some reason it is not working for a Scatterplot
or in my case: plot (x,y). With a histogram I had no problems.
Is it not possible to rotate a simple plot?
Thanks a lot,
Michael
Dear Mr. Ripley,
indeed that it true. sfInit() currently have a bug on Windows depending
on the usage of the Linux tools and the broken Exceptionhandling. Too
bad I never tested it accordingly on Windows (as we do not have any
Windows machines in our institute).
snowfall 1.62 is in the pipe
I'm trying to use a lattice function within a function and have problems
passing the groups argument properly. Let's say I have a data frame
d - data.frame(x = rnorm(100), y = c(a, b))
and want to plot variable x in a densityplot, grouped by the variable y, then
I would do something like
I was trying to rotate a plot using R. I tried most of the examples offered
here in this forum, but for some reason it is not working for a Scatterplot
or in my case: plot (x,y). With a histogram I had no problems.
Is it not possible to rotate a simple plot?
I'm not 100% sure what exactly you
Hi,
in the first example, your feasible set is just one point (the one
that fulfills the 3 equations) and thus there is only this one point
which can maximize the objective function. In the second case, the
feasible set is a line. But the simplex algorithm tries to find an
optimizing value of the
Dear R gurus,
I have some climatic data for a region of the world. They are monthly averages
1950 -2000 of precipitation (12 months), minimum temperature (12 months),
maximum temperature (12 months). I have scaled them to 2 km x 2km cells, and
I have around 75,000 cells.
I need to feed them
Hello,
For entering data alone, you would not need a package. For simple
questionnaires, you could write a function.
It could go like this. For example you want to record people's names and
their ages:
# Sets up an empty database
database - c()
enter_data - function() {
show(Enter
Try the pseudo inverse:
m - rbind(c(1, 1, 1), c(1, 0, 1), c(0, 1, 0))
b - c(5, 2, 3)
library(MASS)
ginv(m) %*% b
On Thu, Dec 11, 2008 at 2:20 AM, Chris Line [EMAIL PROTECTED] wrote:
I have a set of linear equations and would like to find any feasible
solution. A simplex solution works in
You can have look to
*S. Dray*. On the number of principal components: A test of
dimensionality based on measurements of similarity between matrices.
/Computational Statistics and Data Analysis/, 52:2228-2237, 2008.
which is implemented in the testdim function of the ade4 package.
Cheers.
On Wed, 10 Dec 2008 18:29:43 +0100, Peter Dalgaard [EMAIL PROTECTED] wrote:
You (as many before you) have overlooked the ave() function, which can
replace the ordering as well the do.call(c,tapply())
Majority of questions on this list concern data manipulation. Many are repetitive.
Hi all,
Is there a simple function already implemented for getting the ISO
weeks of a Date object?
I couldn't find one, and so wrote my own function to do it, but would
appreciate a pointer to the default way. If a function is not yet
implemented, could the code below be of interest to submit to
Dear List,
I have two dataframes with overlapping colnames and want to merge them.
Actually, what I want is more similar to rbind, but the dataframes
differ in their columns. Here are the examples:
df1 - data.frame(A = c(1,2), B = c(m,f), C = c(at home, away))
df2 - data.frame(A = c(2), C =
Dear Stefan,
Why not use merge() if you want to merge two datasets? ;-)
df1 - data.frame(A = c(1,2), B = c(m,f), C = c(at home, away))
df2 - data.frame(A = c(2), C = c(at home))
merge(df1, df2, all = TRUE)
HTH,
Thierry
Dear List,
I am writing a function in R with the facility to store models for later
use in scoring.
It would be very useful if I could include in the name of the file
stored the name of
the model object being stored, this name being chosen by the user in the
function
call. A simple function to
strftime(x, %V)
E.g.
strftime(as.POSIXlt(Sys.Date()), %V)
is 50, and you might want as.numeric() on it.
Note that this is OS-dependent, and AFAIR Windows does not have it.
On Thu, 11 Dec 2008, Gustaf Rydevik wrote:
Hi all,
Is there a simple function already implemented for getting the ISO
A slightly simpler version is
format(Sys.Date(), %V)
On Thu, 11 Dec 2008, Prof Brian Ripley wrote:
strftime(x, %V)
E.g.
strftime(as.POSIXlt(Sys.Date()), %V)
is 50, and you might want as.numeric() on it.
Note that this is OS-dependent, and AFAIR Windows does not have it.
On Thu, 11 Dec
format(d, %U) and format(d, %W) give week numbers using
different conventions. See ?strptime
On Thu, Dec 11, 2008 at 7:43 AM, Gustaf Rydevik
[EMAIL PROTECTED] wrote:
Hi all,
Is there a simple function already implemented for getting the ISO
weeks of a Date object?
I couldn't find one, and
Dear Stefan,
You could use rbind.fill() from reshape package :
install.packages(reshape)
library(reshape)
df1 - data.frame(A = c(1,2), B = c(m,f), C = c(at home, away))
df2 - data.frame(A = c(2), C = c(at home))
rbind.fill(df1, df2)
AB C
1 1m at home
2 2faway
3 2 NA
re pseudo inverse
On the point of generalised inverses - GINV is usually taken to mean the
moore-penrose pseudo inverse - this is the least squares projection.
There are others - e.g. the Drazin inverse which amounts to diagonalisation
- of course this inverse may not be available in R.
Gerard
Dear List Members,
this question is not directly related to R but follows up a comment in
the documentation of 'cov'. In the case of missing values in the input
matrix 'x' there are various options specified via the parameter
'use'. One of them being 'pairwise.complete.obs'. In the 'Details'
Not an R package, but EpiData is free software designed to to exactly this.
It's a wonderful piece of software. Define fields, add annotations, provide
defaults, provide allowable values or ranges, calculate one field based on
entry to another, conditional skipping from one question to another
Hi,
Good idea, what do you say we try and write a page on this in the R
wiki?
I started the topic:
http://wiki.r-project.org/rwiki/doku.php?id=guides:overview-data-manip
Once the content is there, it wouldn't be much of an effort to create
a reference-card format if required.
Best
Hi,
I've estimated a simple kernel density of a univariate variable with
density(), but after I would like to find out the CDF at specific
values.
How can I do it?
thanks for your help, with it I am very close to finish my first
little bit more serious work in R,
Viktor
On Thu, 11 Dec 2008 14:28:31 +0100, Viktor Nagy wrote:
VN Hi,
VN
VN I've estimated a simple kernel density of a univariate variable with
VN density(), but after I would like to find out the CDF at specific
VN values.
VN How can I do it?
VN
Answer 1.
Use approfun to interpolate the outcome from
Hi Kenn,
Well I think your use or round isn't optimal solution. If you using
round(x,1)-round(x,1) you create 2 problems
First: error propagation because you make 2 round.
Second: you don't using guard digits approach.
The optimal use of round is using in last calculation:
Look this
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
You (as many before you) have overlooked the ave() function, which can
replace the ordering as well the do.call(c,tapply())
Majority of questions on this list concern data manipulation. Many are
repetitive. Overlooking like that will always happen unless some
comprehensive data
On Thu, 11 Dec 2008 15:11:06 +0100, Adelchi Azzalini wrote:
AA
AA In practical terms you can still use density(), as indicated above, but
AA selecting a suitably smaller bandwith compared to the one used for
AA density estimation.
PS.
of course this is numerically most inefficient...
Adelchi
If you're intending to create a model using PCs as predictors, select
the PCs based on whether they contribute significanctly to the model
fit.
In chemometrics (multivariate stats in chemistry, among other things),
if we're expecting 3 or 4 PC's to be useful in a principal component
regression,
Peter Dalgaard wrote:
Wacek Kusnierczyk wrote:
Peter Dalgaard wrote:
Johannes Huesing wrote:
Stavros Macrakis [EMAIL PROTECTED] [Wed, Dec 10, 2008 at
04:59:25AM CET]:
So I conclude that what is really meant by R semantics are based on
Scheme
semantics is R has functions as first-class
dear Hadley and Bert,
thank you very much for your suggestions. I asked one question and I
learned 2 things:
1. Hadley,
library(plyr)
ddply(data, .(V1), colwise(cl))
that is exactly what I was searching for.
2. Bert,
?tapply says that the first argument is an **atomic** vector. A
factor is
If however I wanted to call the function densityplot within a function and
pass the groups argument as an argument of that function, how would I have
to proceed? It is not as straightforward as
f - function(data, groupvar) {
densityplot(~ x, data, groups = groupvar)
}
probably
Dear Prof. Ripley,
Thank you for your quick response.
(A)
link-sqrt is a name and not accepted. link=sqrt is a literal character
string, and is.
I am not entirely sure whether I understand that statement but this is what I
found out. If I specify family=gaussian(link=sqrt), the glm()
Hi,
I agree -- and my examples using round were meant as bad and dangerous
examples. Using round at the last step is better and may solve the problem,
but in your example ...
round(8.8-7.8,1)==1
[1] TRUE
... you have to know in advance how many decimal places can possibly make a
difference (is
Sys.setlocale(,C)
x1 - as.character(date()) # I use date to record the time, and save
it to sqlite database, to it converted to character
x1_2 - strptime(x1, %a %b %d %H:%M:%S %Y)
x2 - as.character(date())
x2_2 - strptime(x2, %a %b %d %H:%M:%S %Y)
X-c(x1_2,x2_2)
order(X) ## I want to get the
Hello!
I am pretty new to R statistics and trying to figure out how to show
the barplot and the single data points in one single plot.
I am using the tplot code from the following site for showing the
single datapoints:
http://biostat.mc.vanderbilt.edu/twiki/pub/Main/TatsukiRcode/
You might take a look at the transformBy function in the doBy package
For example,
new.df=transformBy(~group,data=my.df, new=y/max(y))
David Freedman
baptiste auguie-2 wrote:
Dear list,
I have a data.frame with x, y values and a 3-level factor group,
say. I want to create a new
Hi everyone,
I want to write a wrapper function that uses the hist() function. Now I want
to allow the hist breaks argument as optional in my function. So when my
function contains the breaks argument I want the hist() function to use it,
if not, I want the hist() function to use its default for
Hi Shubha
I have created an extension DLL for downloading time series data from Reuters.
You can download it from here:
http://www.theresearchkitchen.com/blog/archives/287
There is also a short manual available at the same location:
Hi Shubha,
I'm replying offlist because although I don't have an answer to your
question, i thought i would point you in the direction of Rmetrics
(http://www.rmetrics.org/rmetricsPackages.htm) as there might be
something in there that may be helpful if you don't get a more direct
answer from
On Thu, 11 Dec 2008 15:19:03 +0100, hadley wickham [EMAIL PROTECTED] wrote:
You (as many before you) have overlooked the ave() function, which can
replace the ordering as well the do.call(c,tapply())
Majority of questions on this list concern data manipulation. Many are
repetitive.
Dear Alessia,
I am very new to R and wanted to know if there is a package that, given very
long nucleotide sequences, searches and identifies short (7-10nt) motifs.. I
would like to look for enrichment of certain motifs in genomic sequences.
I tried using MEME (not an R package, I know),
On 12/11/2008 11:01 AM, Mark Heckmann wrote:
Hi everyone,
I want to write a wrapper function that uses the hist() function. Now I want
to allow the hist breaks argument as optional in my function. So when my
function contains the breaks argument I want the hist() function to use it,
if not, I
Hi all,
Just on this question :
can I assume any R internal defined function can be used to describe the
link (e.g. = arctan) so long as its increasing and monotone?
How might abs work for example - (except at 0)?
And/or finally, can I define any old function in R called myfun and use
Thank you
ryancw wrote:
Not an R package, but EpiData is free software designed to to exactly
this. It's a wonderful piece of software. Define fields, add annotations,
provide defaults, provide allowable values or ranges, calculate one field
based on entry to another, conditional
Sorry for the last post. I didn't use the latest version of R. It
works under Linux as well for R-2.8.0 patch.
Best
On Thu, Dec 11, 2008 at 11:34 PM, ronggui [EMAIL PROTECTED] wrote:
Sys.setlocale(,C)
x1 - as.character(date()) # I use date to record the time, and save
it to sqlite database,
do you mean something like the following:
f - function (x) {
deparse(substitute(x))
}
x - 5
y - 6
z - 7
f(x)
f(y)
f(z)
I hope it helps.
Best,
Dimitris
Philip Whittall wrote:
Dear List,
I am writing a function in R with the facility to store models for later
use in scoring.
It would
Please do read the help page: fortune(WTFM) applies.
On Thu, 11 Dec 2008, Gerard M. Keogh wrote:
Hi all,
Just on this question :
can I assume any R internal defined function can be used to describe the
link (e.g. = arctan) so long as its increasing and monotone?
How might abs work for
A quick question really:
I have a database with extremely long integer IDs (eg
588848900971299297), which is too big for R to cope with internally
(it appears to store as a double), and when I do any frequency tables
erroneous results appear. Does anyone know of a package that extends
OK, this should be trivial but I'm not finding it. I want to compress
the test,
if (i==7 | i==10 | i==30 | i==50) {}
into something like
if (i in c(7,10,30,50)) {}
so I can build a excludes vector
excludes - c(7,10,30,50)
and test
if (i in excludes) {}
However, I'm not finding a clue on
Hi,
We're using stats4 for a logistic regression. The code is
chdreg.logit2 - glm(chd ~ age + sex, family = binomial)
summary(chdreg.logit2)
oddsratios - coef(chdreg.logit2)
exp(oddsratios)
# Calculate model predicted values
pred - predict(chdreg.logit2,type=response)
The glm part runs fine,
have a look at merge(), e.g.,
df1 - data.frame(A = c(1,2), B = c(m,f), C = c(at home, away))
df2 - data.frame(A = c(2), C = c(at home))
merge(df1, df2, all = TRUE, sort = FALSE)
I hope it helps.
Best,
Dimitris
Stefan Uhmann wrote:
Dear List,
I have two dataframes with overlapping
Have you looked at the R interface to gmp?
http://cran.r-project.org/web/packages/gmp/index.html
Rory Winston
RBS Global Banking Markets
Office: +44 20 7085 4476
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Aaron Robotham
Sent: 11 December 2008
On Thu, Dec 11, 2008 at 2:10 PM, Prof Brian Ripley
[EMAIL PROTECTED] wrote:
A slightly simpler version is
format(Sys.Date(), %V)
On Thu, 11 Dec 2008, Prof Brian Ripley wrote:
strftime(x, %V)
E.g.
strftime(as.POSIXlt(Sys.Date()), %V)
is 50, and you might want as.numeric() on it.
Take a look at ?any.
On Thu, Dec 11, 2008 at 3:11 PM, David B. Thompson, Ph.D., P.E.,
D.WRE, CFM [EMAIL PROTECTED] wrote:
OK, this should be trivial but I'm not finding it. I want to compress the
test,
if (i==7 | i==10 | i==30 | i==50) {}
into something like
if (i in c(7,10,30,50)) {}
Hi,
This has been asked before but not sufficiently answered from what I
could find. How do you create combinations
with repetitions (multisets) in R?
If I have
set - array(1:3)
And I want to choose all combinations of picking 2 numbers, I want to
get a print out like
[,1] [,2]
[1,]1
hi all,
I want to do a plot and put the legend on the left of y axis this is my
code:
x-seq(1980,2005,1)
plot(x,tfa_ita,type=l,col=1,xlim=c(1979,2005),ylim=c(0.2,1.7),xlab=,ylab=,main=Totale
Attivita` Finanziarie)
lines(x,tfa_spa,type=l,col=2)
lines(x,tfa_aus,type=l,col=3)
A quick question really:
I have a database with extremely long integer IDs (eg
588848900971299297), which is too big for R to cope with internally
(it appears to store as a double), and when I do any frequency tables
erroneous results appear. Does anyone know of a package that extends
OK, this should be trivial but I'm not finding it. I want to compress
the test,
if (i==7 | i==10 | i==30 | i==50) {}
into something like
if (i in c(7,10,30,50)) {}
so I can build a excludes vector
excludes - c(7,10,30,50)
and test
if (i in excludes) {}
Works for me.
Hi,
It is generally not the case that the best PC set, say, the top k PCs (where
k p, p being the number of predcitors) contain the best predictor subset
in linear regression. Hadi and Ling (Amer Stat, 1998) show that it is even
possible to have an extreme situation where the first (p-1) PCs
Hi,
Perhaps you can use expand.grid and then remove the mirror combinations,
values - 1:3
tmp - expand.grid(values, values)
unique.combs - tmp[tmp[, 1]=tmp[, 2], ]
unique.combs[do.call(order, unique.combs), ] # reorder if you wish
Var1 Var2
111
412
713
522
8
i %in% c(7,10,30,50)
On Thu, Dec 11, 2008 at 12:11 PM, David B. Thompson, Ph.D., P.E., D.WRE, CFM
drdbthomp...@gmail.com wrote:
OK, this should be trivial but I'm not finding it. I want to compress the
test,
if (i==7 | i==10 | i==30 | i==50) {}
into something like
if (i in c(7,10,30,50))
If they are IDs, you presumably don't need to perform arithmetic on them, so
why not store them as strings? If you're reading them with read.table, see
the colClasses parameter. I am not sure how to do this in RODBC;
as.isthere (as in read.table) does not affect columns that look like
numbers --
Dear R users --
I think this question was asked before but there was no reply to it.
I would appreciate any suggestion any of you might have. I am
interested in plotting several implicit functions (F(x,y,z)=0) on
the same fig. Is there anyone who has an example code of how to do
this?
Thank
Dear R-user
I need a function to approximate a complex integration. My function is:
aprox2=function(s,x,rate){
dexp(x,rate)*exp(-s*x)
}
where argument s is a complex number. I can't use the integrate function
because it's only used with numeric arguments
Does anyone know some function to
Dear Reuben,
On Thu, Dec 11, 2008 at 12:53 PM, baptiste auguie ba...@exeter.ac.uk wrote:
Hi,
Perhaps you can use expand.grid and then remove the mirror combinations,
values - 1:3
tmp - expand.grid(values, values)
unique.combs - tmp[tmp[, 1]=tmp[, 2], ]
unique.combs[do.call(order,
On Thu, 11 Dec 2008, Reuben Cummings wrote:
Hi,
This has been asked before but not sufficiently answered from what I
could find. How do you create combinations
with repetitions (multisets) in R?
If I have
set - array(1:3)
Why wrap 1:3 in array() ??
And I want to choose all combinations
Hi,
I often have the problem of combining data sets of different lengths
together.
Simple example: I have data frame a, with two columns C1 and C2
and another data frame b with only one column V1.
Data frame b is much bigger than a, but C1 of a has the same
levels as V1 of b. (so in other
On Wed, Dec 10, 2008 at 6:39 PM, Bert Gunter gunter.ber...@gene.com wrote:
...?tapply says that the first argument is an **atomic** vector. A factor
is not an atomic vector. So tapply interprets it as such by looking only at
its representation, which is as integer values.
What is the
Hi,
How can I check if a certain ... argument has been passed on to my
user-defined function or not?
foo - function(data, ...)
{
### here I want to check whether xlab was passed with the ... arguments
### or if the ... arguments did not contain an xlab argument
}
I tried missing(xlab) ,
Try this:
foo - function(data, ...)
{
### here I want to check whether xlab was passed with the ... arguments
### or if the ... arguments did not contain an xlab argument
args - list(...)
return(ifelse(xlab %in% names(args), Exists, Missing))
}
On Thu, Dec 11, 2008 at 4:22 PM, Mark Heckmann
See
?match.call
and note the expand.dots arg.
HTH,
Chuck
On Thu, 11 Dec 2008, Mark Heckmann wrote:
Hi,
How can I check if a certain ... argument has been passed on to my
user-defined function or not?
foo - function(data, ...)
{
### here I want to check whether xlab was passed
Gabor Grothendieck ggrothendieck at gmail.com writes:
format(d, %U) and format(d, %W) give week numbers using
different conventions. See ?strptime
Gabor,
the results of format(aDate, W) appear to be incorrect anyway, see:
format(as.Date(2008-01-01), %W) #- 00
There is never a
According to the definition in ?strptime (which is not the same as the
ISO definition):
format(x, %W) returns
Week of the year as decimal number (00–53) using Monday as the first
day of week (and typically with the first Monday of the year as day 1
of week 1). The UK convention.
The first day
Replies inline below.
Best regards,
-- Bert
___
From: macra...@gmail.com [mailto:macra...@gmail.com] On Behalf Of Stavros
Macrakis
Sent: Thursday, December 11, 2008 10:53 AM
To: Bert Gunter
Cc: Patrizio Frederic; r-help@r-project.org
Subject: Re: [R] repeated
I'm new to R so forgive me if this seems like a simple question:
So I have table where the row titles are string variables. When I plot the
data with rows along the x-axis, the data is ordered alphabetically as
opposed to the order of the table.
How can I preserve the row order of the table in
It would be easier to answer your question if we knew what your
data look like, what R commands you've tried, and what result
you want.
One possibility: plot the data against 1:nrow(yourdata), and add
the row names as labels.
Sarah
On Thu, Dec 11, 2008 at 2:35 PM, qroberts lvaic...@bu.edu
for a simple example:
x - list()
x[[a]] - list(a=c(1,2,3),b=c(3,4,5))
x[[b]] - list(a=c(6,7,8),b=c(9,10,11))
lapply(x,sum)
this fails w/
Error in FUN(X[[1L]], ...) : invalid 'type' (list) of argument
Just wondering if I have overlooked something obvious.
one can also do:
Gabor Grothendieck ggrothendieck at gmail.com writes:
According to the definition in ?strptime (which is not the same as the
ISO definition):
format(x, %W) returns
Week of the year as decimal number (00–53) using Monday as the first
day of week (and typically with the first Monday of
This will recursively lapply and may or may not be
what you are looking for:
rapply(x, sum)
On Thu, Dec 11, 2008 at 2:59 PM, Whit Armstrong
armstrong.w...@gmail.com wrote:
for a simple example:
x - list()
x[[a]] - list(a=c(1,2,3),b=c(3,4,5))
x[[b]] - list(a=c(6,7,8),b=c(9,10,11))
Perhaps you mean is that the definition ought be otherwise but
at least according to one standard the definition is correct:
http://www.opengroup.org/onlinepubs/009695399/functions/strptime.html
On Thu, Dec 11, 2008 at 3:01 PM, Hans W. Borchers hwborch...@gmail.com wrote:
Gabor Grothendieck
Hello,
I have a file with two dependent variables (three and five) and one
independent variable. I do i.mod - lm(cbind(three, five) ~ species,
data=i.txt) and get the following output:
Coefficients:
three five
(Intercept) 9.949 9.586
species -1.166 -1.156
I do a
On Thu, Dec 11, 2008 at 2:55 AM, Thomas Zumbrunn t.zumbr...@unibas.ch wrote:
I'm trying to use a lattice function within a function and have problems
passing the groups argument properly. Let's say I have a data frame
d - data.frame(x = rnorm(100), y = c(a, b))
and want to plot variable x in
On Thu, 11 Dec 2008, Whit Armstrong wrote:
for a simple example:
x - list()
x[[a]] - list(a=c(1,2,3),b=c(3,4,5))
x[[b]] - list(a=c(6,7,8),b=c(9,10,11))
lapply(x,sum)
this fails w/
Error in FUN(X[[1L]], ...) : invalid 'type' (list) of argument
Just wondering if I have overlooked something
Bert,
Thanks for your reply. I suspect we agree more than you might think
Comments inline below. I've snipped out parts.
-s
On Thu, Dec 11, 2008 at 2:45 PM, Bert Gunter gunter.ber...@gene.com wrote:
Rationale? -- you'll have to ask the developers
Hmm. It would be
Joe Trubisz wrote:
Hi...
Is this possible in R?
I have 2-sets of data, that were collected simultaneously using
2-different data acquisition schemes.
The x-values are the same for both.
The y-values have different ranges (16.4-37.5 using one method, 557-634
using another).
In theory, if
Dear all R users,
I am going to use R to process some of my physiological data about eye.
The problem is the recording machine does not sample in a reliably
constant rate: the time intervals between data sampled can vary from
9msec to ~120msec, while most around in the 15-30msec range.
The below
Thanks, Gabor and Prof. Ripley.
Sorry for the oversight.
I grepped the lapply help for recursive prior to sending my question.
why does it appear as *r*ecursive in the help file? or is that just
a formating problem on my machine?
-Whit
On Thu, Dec 11, 2008 at 3:13 PM, Prof Brian Ripley
Is there an equivalent to MPlus's Full Information Maximum Likelihood (FIML)
missing data estimator for R? If so, is there a way to take covariance
structures produced by such a package and perform multiple regression with
these?
If you are unfamiliar with Mplus' FIML below is a link to their
replies inline below.
Bert Gunter wrote:
Replies inline below.
[bert (?)]...?tapply says that the first argument is an **atomic**
vector. A
factor is not an atomic vector. So tapply interprets it as such by looking
only at its representation, which is as integer values.
On Thu, 11 Dec 2008, Whit Armstrong wrote:
Thanks, Gabor and Prof. Ripley.
Sorry for the oversight.
I grepped the lapply help for recursive prior to sending my question.
why does it appear as *r*ecursive in the help file? or is that just
a formating problem on my machine?
It is marked as
Hi all,
Is there a function to extract row names from a data frame based on
row names from another data frame? I can write a loop function to do
this, but this may be inefficient in terms of processing.
thanks for any information,
Wade
__
Hi, I am simulating 2-dimensional data using the RandomFields library and
the gaussRF function therein. While this is done with the code below, I
would like the landscape to be continuous or smooth at the edges. That is, I
would like the upper edge to smoothly connect to the lower edge AND the
1 - 100 of 144 matches
Mail list logo