Folks,
I often build up R objects starting from NULL and then repeatedly
using rbind() or cbind(). This yields code like:
a <- NULL
for () {
onerow <- craft one more row
a <- rbind(a, onerow)
}
This works because rbind() and cbind() are forgiving when presented
with a NULL arg: they act like
G'day Edna,
On Sat, 14 Mar 2009 23:00:12 -0500
Edna Bell wrote:
> Anyway, my question now is: what determines if a function is a
> builtin vs. a closure, please?
The help page of typeof states that:
The possible values are listed in the structure 'TypeTable' in
'src/main/util.c'. Current valu
G'day Edna,
On Sat, 14 Mar 2009 22:52:38 -0500
Edna Bell wrote:
> Dear R Gurus:
Well, I guess I can answer nevertheless. :)
> How do I find the functions which are primitives, please?
?is.primitive
Thus, the following code would give you all the primitive functions in
package base:
R> pos <
Dear R Gurus:
I'm working slowly through "R Programming for Bioinformatics", which
is really interesting!
Anyway, my question now is: what determines if a function is a
builtin vs. a closure, please?
For instance:
> typeof(sqrt)
[1] "builtin"
> typeof(mean)
[1] "closure"
>
Thanks,
Edna Bell
_
Dear R Gurus:
How do I find the functions which are primitives, please?
Thanks,
Edna Bell
__
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
Dear R-helpers,
I have produced several models, named model1, model2, model3, etc...
I would like to extract several elements from each model's object, e.g. at
minimum the estimates, SEs, and P values of each model's intercept and
slopes, model R-squared, and AIC...
...and then produce a new obj
Exactly what type of cleaning do you want to do on them? Can you read
in the data a block at a time (e.g., 1M records), clean them up and
then write them back out? You would have the choice of putting them
back as a text file or possibly storing them using 'filehash'. I have
used that technique
Dear R helpers:
I am a newbie to R and have a question related to cleaning large data frames
in R.
So far, I have been using SAS for data cleaning because my data sets are
relatively large (handling multiple files, each could be as large as 5-10
G).
I am not a fan of SAS at all and am eager to mo
Hi everyone,
I'm working on a multi-level model based on one of the "radon" example from
Gelman & Hill's book (
http://www.amazon.com/Analysis-Regression-Multilevel-Hierarchical-Models/dp/052168689X/ref=pd_bbs_1?ie=UTF8&s=books&qid=1237069199&sr=8-1
).
I add an "region" indicator (has numeric val
Dear hyena,
Actually, looking at this a bit more closely, the first models dedicate 6
parameters to the correlational and variational structure of the three
variables that you mention -- 3 variances and 3 covariances; the second
model also dedicates 6 parameters -- 3 factor loadings and 3 error va
Dear hyena,
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On
> Behalf Of hyena
> Sent: March-14-09 5:07 PM
> To: r-h...@stat.math.ethz.ch
> Subject: [R] SEM model testing with identical goodness of fits
>
> HI,
>
>I am testing several
Use reshape in the base of R or melt/cast in the reshape package.
Using the former:
Lines <- "ID Year V1
1 19801
1 19812
1 19826
1 19834
2 19805
2 19815
2 19825
2 19836"
# DF <- read.table("myfile.dat", header = TRUE)
DF <- read.table(textConnect
On Mar 14, 2009, at 5:46 PM, David Winsemius wrote:
On Mar 14, 2009, at 5:14 PM, David Kaplan wrote:
Greetings all,
I have two questions. I have a data set that is arranged in the
example below. I wish to obtain a plot of the performance of each
ID over Year on v1. It's not clear how
On 14/03/2009 12:02 PM, Pedro Mardones wrote:
Dear all;
Does anyone know how to add grid lines to a persp plot? I've tried
using lines(trans3d..) but the lines of course are superimposed into
the actual 3d surface and what I need is something like the plot shown
in the following link:
http://ther
On 12/03/2009 6:51 AM, Duncan Murdoch wrote:
On 11/03/2009 11:35 PM, Oliver wrote:
hi, all
I am looking at R package RGL to draw a colored mesh/surface plot like
this one (from matlab).
http://www.mathworks.com/access/helpdesk/help/techdoc/visualize/cbar.gif
The key features I am looking for i
On Mar 14, 2009, at 5:14 PM, David Kaplan wrote:
Greetings all,
I have two questions. I have a data set that is arranged in the
example below. I wish to obtain a plot of the performance of each
ID over Year on v1. It's not clear how I set this up?
ID Year V1 1 19801 1 1981
check out the interaction.plot. This *may* be what you are looking
for.
?interaction.plot
On Mar 15, 8:14 am, David Kaplan wrote:
> Greetings all,
>
> I have two questions. I have a data set that is arranged in the example
> below. I wish to obtain a plot of the performance of each ID over Ye
Greetings all,
I have two questions. I have a data set that is arranged in the example
below. I wish to obtain a plot of the performance of each ID over Year
on v1. It's not clear how I set this up?
ID Year V1
1 19801
1 19812
1 19826
1 19834
2 19805
2
HI,
I am testing several models about three latent constructs that
measure risk attitudes.
Two models with different structure obtained identical of fit measures
from chisqure to BIC.
Model1 assumes three factors are correlated with each other and model
two assumes a higher order factor exi
I don't have a solution, but just noticed that levelplot has a
colorkey=FALSE mode that might solve the problem. The filled.contour
help page suggests to me that you would need to alter the code of
filled.contour to get what you want, but it appears to be an .Internal
function.
On Mar 13
On Fri, Mar 13, 2009 at 7:58 AM, Ross Culloch wrote:
>
> Hi Kingsford,
>
> Thanks yet again for your help! I have tried this, and once again i have
> failed! I have put the code that i've used below (i'm sure you'll note some
> bad practice)
It doesn't appear you changed the factor to a character
On Mar 14, 2009, at 12:02 PM, Pedro Mardones wrote:
Dear all;
Does anyone know how to add grid lines to a persp plot? I've tried
using lines(trans3d..) but the lines of course are superimposed into
the actual 3d surface and what I need is something like the plot shown
in the following link:
htt
Thanks Joe. This is good. But can you please tell me how to set up the
variable as either fixed effect or random effect?
*as per the documentation, rhierBinLogit model is:*
P = exp(Xij BETAi) / (1 + exp(Xij BETAi)
where, BETAi = Z Delta[i,]
for keyword i and week j.
Here all X are modeled as r
Dear Colleagues:
I need to make a graphic that uses the Nimbus rather than Helvetica font
family so that the font can be embedded in the encapsulated postscript
file. This is to satisfy a requirement from a journal for electronic
submission of figures. I do the following:
postscript('lowes
Ideas... not a solution. Plot the grid within your ranges using
something along the lines, literally and figuratively, based on the
second example of persp's help pages. For the z=8 grid lines on that
example you could use:
for (ix in seq(-10,10, by=5)) lines (trans3d(x=ix, y=seq(-10,10,
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Hi All,
I found an interesting thread discussing R and other packages here:
http://anyall.org/blog/2009/02/comparison-of-data-analysis-packages-r-matlab-scipy-excel-sas-spss-stata/
Plenty of well-reasoned comments.
I thought it may be informative
library(chron)
Lines <- " ID Target Actual Date Time
1 2030 0 -57.0 12/20/08 17:03:00
2 2030 90 90.0 12/20/08 18:41:00
3 2030 45 43.8 12/21/08 14:36:00
4 2030 0 -23.8 12/21/08 19:30:00
5 2030 90 90.2 12/21/08 21:48:00
6 2030 45 48.6 12/22/08 13:02:00"
DF
Dear all;
Does anyone know how to add grid lines to a persp plot? I've tried
using lines(trans3d..) but the lines of course are superimposed into
the actual 3d surface and what I need is something like the plot shown
in the following link:
http://thermal.gg.utah.edu/tutorials/matlab/matlab_tutorial
I am working with time-date data to plot temporal activity.
I would like to create a "template" to repeat the date and time for
all potential data sets that may have differing time zone origins.
Reading the > ?strptime help section I note that time zone is a
variable to add.
Looking at the
with the example below, i think the following solves your problem:
text = ',,,apple,,orange,lemon,strawberry'
gsub('^,+|,+(?=,)|,+$', '', text, perl=TRUE)
vQ
Tom wrote:
> sub(",$","",gsub("^,*|,*$|(,)*", "\\1",
> ",,,apple,,orange,lemon,strawberry"))
>
> It may not be the
Read this from the as.POSIXct help page more carefully:
"Usage
as.POSIXct(x, tz = "", ...)"
"tz
A timezone specification to be used for the conversion, if one is
required. System-specific (see time zones), but "" is the current
timezone, and "GMT" is UTC (Universal Time, Coordinated)."
--
On Mar 14, 2009, at 8:52 AM, Aimin Yan wrote:
I have a data set like this:
> head(FormatedData)
ID Target Actual Date Time
1 2030 0 -57.0 12/20/08 17:03:00
2 2030 90 90.0 12/20/08 18:41:00
3 2030 45 43.8 12/21/08 14:36:00
4 2030 0 -23.8 12/21/08 19:30:00
5 20
Search the R mailing list archive for other reports of success on
AIX. I compiled a 64-bit version of R-2.8.0 successfully with the
following options:
export OBJECT_MODE=64
export CC="xlc_r -q64"
export CFLAGS="-O -qstrict"
export CXX="xlC_r -q64"
export CXXFLAGS="-O -qstrict"
export F77="xlf_r -q
baked Toast gmail.com> writes:
> Now the problem is that when I
> want to import this data to R with the RODBC package (the used code is
> displayed on the bottom of this message) it converts this marker to a
> number and cuts of the starting zero's.
Use parameter as.is in odbcQuery, which is b
I cannot. (And you should know from reading the Posting Guide that
homework questions are discouraged on R-help, anyway. And you might
also note that informative subject lines are requested.)
I can suggest that you look at the code and examples in the fExtremes
and extRemes packages.
--
D
It's not entirely clear from the help pages, but reading Thereau and
Gramsch' "Modeling Survival Data" it appears to me that the basehaz
function in the survival package returns the estimated baseline
_cumulative_ hazard function H(t|mean covariates). Such a function
will be monotonic upwa
Hello,
The following book can now be ordered from Amazon.com and Springer.com.
Mixed effects models and extensions in ecology with R. (2009).
Zuur, AF, Ieno, EN, Walker, N, Saveliev, AA, and Smith, GM.
All R code and data sets can be downloaded from:
http://www.highstat.com/books.htm
Kind r
I am working with time-date data to plot temporal activity.
I would like to create a "template" to repeat the date and time for
all potential data sets that may have differing time zone origins.
Reading the > ?strptime help section I note that time zone is a
variable to add.
Looking at the
> if you don't mind about the efficiency of your program, just use ?
> sapply on a range of wavelengths. I'm not particularly keen on the
> implementation of the FORTRAN program you used compared to the one I
> suggested, it seems like a very crude convolution.
>
Thank you very much for your sugge
On 14 Mar 2009, at 13:08, Jinsong Zhao wrote:
Hi,
For a good discussion of the link between colour and spectra I would
suggest,
http://www.fourmilab.ch/documents/specrend/
which provides an open-source C code to perform the conversion you
ask
for. I asked for some advice on how to wrap a
> Hi,
>
> For a good discussion of the link between colour and spectra I would
> suggest,
>
> http://www.fourmilab.ch/documents/specrend/
>
> which provides an open-source C code to perform the conversion you ask
> for. I asked for some advice on how to wrap a R function around this
> code last wee
On Fri, Mar 13, 2009 at 2:39 PM, Martín Quiroga wrote:
> Hi everyone! I am a biologist from Argentina and have to solve this problem.
> I have an insect population obtained from 10 different nests and need to
> know its sex ratio. But as I cannot ensure insects independence I need to
> run a model
I have a data set like this:
> head(FormatedData)
ID Target Actual Date Time
1 2030 0 -57.0 12/20/08 17:03:00
2 2030 90 90.0 12/20/08 18:41:00
3 2030 45 43.8 12/21/08 14:36:00
4 2030 0 -23.8 12/21/08 19:30:00
5 2030 90 90.2 12/21/08 21:48:00
6 2030 45
sub(",$","",gsub("^,*|,*$|(,)*", "\\1",
",,,apple,,orange,lemon,strawberry"))
It may not be the best solution, but it was my first thought
Tom
- Original Message -
From: "Daren Tan"
To:
Sent: Saturday, March 14, 2009 6:42 PM
Subject: [R] gsub and regex to tidy comma-limited
Hi,
For a good discussion of the link between colour and spectra I would
suggest,
http://www.fourmilab.ch/documents/specrend/
which provides an open-source C code to perform the conversion you ask
for. I asked for some advice on how to wrap a R function around this
code last week but sad
Hi there,
I try to plot visible light spectrum (380nm~780nm) with color
corresponding to the specific wavelength. However, I don't find a
function that could do this.
Another question, it's possible to plot a color space chromaticity
diagram like this:
http://upload.wikimedia.org/wikipedia/common
Dear all,
I use the RODBC package to import data from a Microsoft Access
database. When importing the different variables I use a marker, a
unique code which consists of a series of numbers starting with
several (not a fixed number) of zero's. Now the problem is that when I
want to import this dat
OK, the underlying problem is that the projection code (i.e. the mapproj package plus some
parts of the maps package) does not clip its output to the specified limits.
Let me explain:
1) the maps databases are made up of line segments which are combined to form
polygons
2) when you call map() w
hi
could you help to write program in R to solve
1-write function to time MLE of extreme Value type 1 distribution
2 find the Extreme Value 1 by using MOM, OLS,PWM,LS\
_
[[alternative HTML version deleted]]
__
Hi Jim, that works, which is great!!! Alexis also gave me an answer that
works if you wanted to look at that post, but that was based on adding
colours to existing Sets from BrewerColor. However, your answer is the first
i've got where you can actuall pick your entire palette, it's something that
Add perl = TRUE
On Sat, Mar 14, 2009 at 5:42 AM, Daren Tan wrote:
> I am cleaning up comma-limited values, so that only one comma
> separates each value. Using the example below, as much as I try with
> regex, I can't remove the last comma. I hope to have a one-liner
> solution, if possible.
>
>
Ross Culloch wrote:
Hi Kingsford,
Thanks yet again for your help! I have tried this, and once again i have
failed! I have put the code that i've used below (i'm sure you'll note some
bad practice) if that is any use to help explain where i'm going wrong, it
seems to run fine and feeds back just
zoo makes no assumptions regarding times and
dates other than they are ordered and have certain
methods. It doesn't even know the difference between
a time and a date nor does it know about any time
or date classes. (The exception is interface routines
such as read.zoo.)
If the dates uniquely spe
Sorry there was some garbage at the beginning. Should just
be:
On Sat, Mar 14, 2009 at 6:03 AM, Gabor Grothendieck
wrote:
> If its truly 1-1, i.e. they have the same number of comonents and
> the i-th component of the first series corresponds to the i-th
> component of the second seriesm, then j
I am cleaning up comma-limited values, so that only one comma
separates each value. Using the example below, as much as I try with
regex, I can't remove the last comma. I hope to have a one-liner
solution, if possible.
gsub("^,*|,*$|(,)*", "\\1", ",,,apple,,orange,lemon,strawberry")
[1] "a
johnhj wrote:
Hii Jholtman,
I will make a graph of the median values and not to print to the console.
I tried to plot with:
plot(V3 ~ grup, data = median) ??
but I get an error message.
I also tried to print the output of median<-with(x, tapply(V3, grup,
median)) to a text file with the X a
Jeff:
As a followup to this question -- I have a pair of tables that I want to do
a 1 to 1 join on, but the date field contains the full time, down to the
second, to base the join on (e.g. in a given day, there are going to MANY
observations, but not at the exact same time). I might be missing
so
Hello ,
I am hoping for some advice regarding obtaining the values for the
hazard function in a cox regression that I have undertaken. I have a
model in the following form, analysed with the package survival (v.
2.34-1) and a log-log plot obtained using Design (v. 2.1-2).
For two variables,
Martín Quiroga ssdfe.com.ar> writes:
> The response variable has a binomial distribution (males or females).
> Ive been reading for a while and found the MASS and lmer packages that will
> allow me to do such a thing with my data. I found the script I should write
> lmer(y~fixed+(time | random)
59 matches
Mail list logo