Martín Quiroga mquiroga at 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 |
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
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
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
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]
Sorry there was some garbage at the beginning. Should just
be:
On Sat, Mar 14, 2009 at 6:03 AM, Gabor Grothendieck
ggrothendi...@gmail.com 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
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
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
Add perl = TRUE
On Sat, Mar 14, 2009 at 5:42 AM, Daren Tan darenta...@gmail.com 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,
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
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]]
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()
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
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:
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
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 darenta...@gmail.com
To: r-help@r-project.org
Sent: Saturday, March 14, 2009 6:42 PM
Subject: [R] gsub
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
On Fri, Mar 13, 2009 at 2:39 PM, Martín Quiroga mquir...@ssdfe.com.ar 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
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
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
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
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
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
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
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.
--
baked Toast baked.toast at 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,
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 -q64
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
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).
--
David
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 best
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
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:
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 -
-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
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,
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:
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
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:
On Fri, Mar 13, 2009 at 7:58 AM, Ross Culloch ross.cull...@dur.ac.uk 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
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
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
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
check out the interaction.plot. This *may* be what you are looking
for.
?interaction.plot
On Mar 15, 8:14 am, David Kaplan dkap...@education.wisc.edu 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
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
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
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:
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
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 -
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 models
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
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=UTF8s=booksqid=1237069199sr=8-1
).
I add an region indicator (has numeric values
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
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 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
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
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
G'day Edna,
On Sat, 14 Mar 2009 22:52:38 -0500
Edna Bell edna.bel...@gmail.com 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
G'day Edna,
On Sat, 14 Mar 2009 23:00:12 -0500
Edna Bell edna.bel...@gmail.com 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
58 matches
Mail list logo