Michael, I really appreciate your help.
but I got the following error message when I wan trying to run the function
written by you:
Error in out[i, ] - apply(help[, c(grp1 + 1, grp2 + 5)], 2, sample, 1) :
number of items to replace is not a multiple of replacement length
I am not quite sure
Dear All,
I like to know how to simulate survival data using median (or mean) survival
time. Any help will be greatly appreciated.
Best wishes,
Kere
Kerenaftali Klein PhD| Biostatistician | Queensland Clinical Trials
Biostatistics Centre
The University of Queensland | School of Population
On 16 November 2010 16:10, wangwallace talentt...@gmail.com wrote:
Michael, I really appreciate your help.
but I got the following error message when I wan trying to run the function
written by you:
Error in out[i, ] - apply(help[, c(grp1 + 1, grp2 + 5)], 2, sample, 1) :
number of items
gauss--hermite weights the function f(x) to be integrated by exp(-nodes^2),
ie, it uses the 'trick' that exp(nodes^2) * exp(-nodes^2) = 1. so your sum
should be exp(qq$nodes^2) * exp(-qq$nodes^2) * f(qq$nodes^2) =
exp(qq$nodes^2) * qq$weights * f(qq$nodes). hence
sum(exp(qq$nodes^2) *
Hi Kere,
Step 1 is choose an appropriate distribution :) Do you have one in
mind ? Or are you interested in examining the effects of survival
times having the same mean but generated with alternative
distributions ?
One ready-rolled alternative is the SimSurv method in package prodlim.
To find
Hi,
The pdf function would not let me change the paper size and gives me
the following warning:
pdf(figure.pdf, width=6, height=10)
Warning message:
‘mode(width)’ and ‘mode(height)’ differ between new and previous
== NOT changing ‘width’ ‘height’
If I use the option paper = a4r,
Dear list,
I am new to this list and I am new to the world of R. Additionally I am not
very firm in statistics either but have to deal. So here is my problem:
I have a dataset (which I attach at the end of the post) with a binomial
response variable (alive or not) and three fixed factors
Update - sorry for the stupid question, let's say it's pretty late.
For those who may be as tired as I am and get the same warning, the
paper size should be given as an integer!
On 16 November 2010 04:17, Emmanuel Levy emmanuel.l...@gmail.com wrote:
Hi,
The pdf function would not let me
Hi,
who of you in this forum uses R (http://www.r-project.org/) with the
multicore, SNOW or CUDA packages, so for advanced calculations that need
more power than a workstation CPU? On which hardware do you compute these
scripts? At home/ at work or do you have data center access somewhere?
The
When using odfWeave on an OpenOffice input document, I can not open the output
document. I get the message
Format error discovered in the file in sub-document content.xml at 2,4047
(row,col)
Can anyone help me on this? (Apologies if this has been discussed before; I
have not been able to find
From: soren.hojsga...@agrsci.dk
To: r-h...@stat.math.ethz.ch
Date: Tue, 16 Nov 2010 11:32:06 +0100
Subject: [R] odfWeave - Format error discovered in the file in sub-document
content.xml at 2, 4047 (row, col)
When using odfWeave on an OpenOffice
Dear invaluable R-list,
my present problem is arranging/removing some panels in a lattice plot.
Please consider the following:
df.data -
cbind.data.frame(expand.grid(SUBJ=1:5,
TREAT=LETTERS[1:4],
REF=letters[1:4]
R-helpers,
I am hoping to find someone who uses both R and program Stata for GLMs.
I am a beginner R user, finding my own way through; learning code etc. at the
same time as learning the statistics I need to complete my project.
What I have is the code from Stata and am trying to reproduce
On 15/11/2010 10:53 PM, Steve Lianoglou wrote:
Hi,
I'm using R-2.12 on a linux 64bit machine.
When I run a chunk of code inside a foreach() %do% { ...} or %dopar%
{...} (with doMC backend) I keep getting a segfault. Running the
*same* code within lapply(something, function(x) ... ) doesn't
for a simple scatterplot:
plot(X ~ Y, type = 'p', col = 'red')
this produces red-edged circles, but I want to fill in the circles.
can this be done? I checked '? plot' already but couldn't find what I was
looking for.
cheers
--
View this message in context:
Greg,
Two thoughts:
1. It might be possible that 'vars' is a reserved word of sorts and if you
change the name of your vector RF might be happier
2. A way that works for me is to call importance as follows:
sel.imp - importance(sel.rf, class=NULL, scale=TRUE, type=NULL)
and then use the
Hi janesense,
You might find it useful to go through the tutorial listed here for help:
http://www.r-statistics.com/2010/04/repeated-measures-anova-with-r-tutorials/
Cheers,
Tal
Contact
Details:---
Contact me:
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
On 11/16/2010 04:53 AM, Steve Lianoglou wrote:
Hi,
I'm using R-2.12 on a linux 64bit machine.
When I run a chunk of code inside a foreach() %do% { ...} or %dopar%
{...} (with doMC backend) I keep getting a segfault. Running the
*same* code
Hey,
the parameter in the plot command you are looking for is pch.
For filled circles type:
plot(X ~ Y, type = 'p', col = 'red', pch=16)
The 16 determines the filled circles.
Cheers,
Mango
-
B.-Markus Schuller aka Mango
Sensory Ecology Group
DrCJones matthias.goddard at gmail.com writes:
for a simple scatterplot:
plot(X ~ Y, type = 'p', col = 'red')
this produces red-edged circles, but I want to fill in the circles.
plot(X ~ Y, type = 'p', col = 'red', pch=16)
Are you sure you want X~Y and not Y~X?
Columbine Caroline Waring caquilegia at hotmail.com writes:
I am hoping to find someone who uses both R and program Stata for GLMs.
[snip]
What I have is the code from Stata and am trying to reproduce the same
analysis in R - my program of choice.
. glm count md ms rf sg,
Fantastic thanks!
--
View this message in context:
http://r.789695.n4.nabble.com/scatterplot-with-filled-circles-tp3044690p3044772.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Are you sure you want X~Y and not Y~X?
You are right of course - sorry that was just a typo ;)
--
View this message in context:
http://r.789695.n4.nabble.com/scatterplot-with-filled-circles-tp3044690p3044776.html
Sent from the R help mailing list archive at Nabble.com.
Dear Mike,
Good point - thanks. The lines that caused the error mentioned above are simply:
=
x - 1:10
x
@
I could add that the document 'simple.odt' (which comes with odfWeave) causes
the same error - but at row=109, col=1577
sessionInfo()
R version 2.12.0 (2010-10-15)
Platform:
Try adding pch = 16 to your plot call.
HTH,
Dennis
On Tue, Nov 16, 2010 at 4:28 AM, DrCJones matthias.godd...@gmail.comwrote:
for a simple scatterplot:
plot(X ~ Y, type = 'p', col = 'red')
this produces red-edged circles, but I want to fill in the circles.
can this be done? I checked
Can i use readLines to extract only the linees with a specific word within?
If yes, how?
Tnx for help.
--
View this message in context:
http://r.789695.n4.nabble.com/Question-about-readLines-tp3044701p3044701.html
Sent from the R help mailing list archive at Nabble.com.
Hello everyone I'm a new R user and have some questions about lists
indexing, names and the lapply function.
first have a look to my data (names changed) it's a rainfall record for 67
stations,
str(data)
List of 67
$ e.somename.xx.txt :'data.frame': 456 obs. of 5
On Nov 16, 2010, at 7:37 AM, romzero wrote:
Can i use readLines to extract only the linees with a specific
word within?
If yes, how?
Just use:
lines - readLines( appropriate-arguments )
desired - lines[grep(word, lines]
Tnx for help.
--
View this message in context:
x - readLines('yourFile')
x - x[grepl('myWord', x)]
On Tue, Nov 16, 2010 at 7:37 AM, romzero romz...@yahoo.it wrote:
Can i use readLines to extract only the linees with a specific word within?
If yes, how?
Tnx for help.
--
View this message in context:
Dear R Users,
I'm looking for a package that allows to test hypothesis about a
homogeneity of odds ratio in k 2x2 tables. I know that Breslow-Day is
suitable but does anybody could me point out a package? I found diffR,
but as far as I see this package is for IRT theory.
Best,
Robert
Here's an article in Forbes about R in case anyone's interested.
http://blogs.forbes.com/smcnally/2010/11/10/names-you-need-to-know-in-2011-r-data-analysis-software/
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
On Nov 16, 2010, at 10:08 AM, Robert Ruser wrote:
Dear R Users,
I'm looking for a package that allows to test hypothesis about a
homogeneity of odds ratio in k 2x2 tables. I know that Breslow-Day is
suitable but does anybody could me point out a package? I found diffR,
but as far as I see this
Thanks for sharing the questions and responses!
Is it possible to appreciate how much the coefficients matter in one
or the other model?
Say, using Biau's example, using coxph, as.factor(grade2 ==
high)TRUE gives hazard ratio 1.27 (rounded).
As clinician I can grasp this HR as 27% relative
On 15.11.2010 18:35, francois fauteux wrote:
Hello,
Is it possible to remove those extra margins on the sample axis from
plot.dendrogram:
par(oma=c(0,0,0,0),mar=c(0,0,0,0))
ddr-as.dendrogram(hclust(dist(matrix(sample(1:1000,200),nrow=100
On Nov 16, 2010, at 9:18 AM, David Winsemius wrote:
On Nov 16, 2010, at 10:08 AM, Robert Ruser wrote:
Dear R Users,
I'm looking for a package that allows to test hypothesis about a
homogeneity of odds ratio in k 2x2 tables. I know that Breslow-Day is
suitable but does anybody could me
Thank you very much for this clarification,
Sincerely,
Fabon Dzogang.
2010/11/15 Uwe Ligges lig...@statistik.tu-dortmund.de:
On 03.11.2010 16:26, Fabon Dzogang wrote:
Hi,
I run R 2.10.1 under ubuntu 10.04 LTS (Lucid Lynx) and klaR version 0.6-4.
I compute a model over a 2 classes
Hi R experts,
I am looking for a simple error handling approach, whereby I could stop
function execution with a customized error message. For example:
for (i in 1:10) {
if (i == 5)
# I'd like to be able to stop right here with an error message I have
complete control over
}
The problem
Thanks, I'll try that.
Regarding your question Why are you using assignments to indicate the
return values of functions?,
I'd probably answer that it is a matlab vice... But I'll change my habits
soon enough!
Best regards,
Eduardo Horta
On Tue, Nov 16, 2010 at 12:27 PM, Erich Neuwirth
Thanks Mark, for sharing such a great article. Really feeling proud that I
am also part of that great community. BTW can anyone share the code on how
to draw the picture displayed there (blur surface plot)?
Thanks,
-Original Message-
From: r-help-boun...@r-project.org
Hi
Here is one way (If I understood what you did ask).
test-read.table(clipboard, header=T)
test
SubID CSE1 CSE2 CSE3 CSE4 WSE1 WSE2 WSE3 WSE4
1 165626224
2 264726623
3 355555545
4
On Nov 16, 2010, at 10:58 AM, Aleksey Naumov wrote:
Hi R experts,
I am looking for a simple error handling approach, whereby I could
stop
function execution with a customized error message. For example:
for (i in 1:10) {
if (i == 5){
#Then insert:
cat(Informative error message) ;
I have modified my code since asking my original question. The
classifier is now generated correctly (with a good, low error rate, as
expected). However, I am running into two issues:
1) I am getting an error at the prediction stage, I get only NA's when I
try to run data down the forest;
2) I
Il giorno mar, 16/11/2010 alle 10.58 -0500, Aleksey Naumov ha scritto:
for (i in 1:10) {
if (i == 5)
# I'd like to be able to stop right here with an error message I
have
complete control over
}
I'm far from expert, but perhaps this is near what you are looking for
for(i in
From: bogaso.christo...@gmail.com
To: marklee...@gmail.com; r-help@r-project.org
Date: Tue, 16 Nov 2010 21:48:10 +0530
Subject: Re: [R] R article in forbes
Thanks Mark, for sharing such a great article. Really feeling proud that I
am also part
Thank you both, these are very helpful hints.
Chris, could you please suggest how to modify what you sent to also show the
same labels as (horizontal) row names? I have not yet mastered the details
of R graphics...
Many thanks in advance,
Lara
On Mon, Nov 15, 2010 at 3:25 PM, Chris Stubben
David,
Thank you for your help. I agree, stop right here is ambiguous, let me
clarify.
I am at a terminal session, calling my top-level function f(). At any point
inside f() I'd like to be able to generate an error message, print it to the
terminal and stop execution of f() (but still remain in
Hi,
I want to merge 4 data frames with one column in common but I am
getting error message while using this function. Can any one help me out.
merge_all(Br,Ki,Lu,Pr,by=Genes)
Error: could not find function merge_all
I have installed the package: reshape but I still get this error
--
A bit embarrassed to post this seemingly trivial question, but I can't find
anything in the archive that's quite relevant:
a1=1
a2=2
obs=objects(pattern=glob2rx(a?))
I want to utilize 'obs' as a function argument to produce something like:
sum(a1,a2)
Obviously, sum(obs) doesn't work, but
yes, the data.frame is exactly the same as the one I posted earlier.
I was trying to see if the loop function works. And I got that message. here
below is the syntax I was trying to run, followed by the error message at
the end:
sampleX-function(X,nGrp1,nsamples){if(nGrp1=4)stop(can't sample
Hi I am a newbie in R.
I have data with dim of 20.
How to use lm if i want to do regression with the whole design matrix? My y
is the first column.
the left are xs.
Thanks a lot.
--
View this message in context:
Ottorino,
Thank you for your suggestion. I am not really looking to do any debugging,
but merely to issue an informative error message and stop execution of my
top-level function, called at an R terminal. See me response to David for
more details...
Thank you,
Aleksey
On Tue, Nov 16, 2010 at
On Nov 16, 2010, at 10:34 AM, Kevin Ummel wrote:
A bit embarrassed to post this seemingly trivial question, but I
can't find anything in the archive that's quite relevant:
a1=1
a2=2
obs=objects(pattern=glob2rx(a?))
I want to utilize 'obs' as a function argument to produce something
Hello
What's the differente betwen using plot and using print in order to
plot a graph?
For example in order to plot the result of a histogram.
cheers
--
View this message in context:
http://r.789695.n4.nabble.com/plot-vs-print-tp3045256p3045256.html
Sent from the R help mailing list archive
Dear David and Bogasso,
The image is from
http://socserv.mcmaster.ca/jfox/Courses/R-course/index.html, as you'd see
by clicking on the graph in the Forbes article, and it was indeed produced
with scatter3d() in the car package.
Best,
John
John Fox
Senator
Hi,
It would have been nice to see what you've tried.
Here is what I got (maybe not the easiest, but looks like it works)
a1 - 1
a2 - 2
obs - objects(pattern=glob2rx(a?))
sum(unlist(mget(obs, envir=.GlobalEnv)))
[1] 3
It also works with vectors:
a1 - 1:10
a2 - 11:20
sum(unlist(mget(obs,
I write about R every weekday at the Revolutions blog:
http://blog.revolutionanalytics.com
and every month I post a summary of articles from the previous month
of particular interest to readers of r-help.
In case you missed them, here are some articles related to R from the
month of October:
On Tue, Nov 16, 2010 at 10:30 AM, poko2000 quan.poko2...@gmail.com wrote:
Hi I am a newbie in R.
I have data with dim of 20.
How to use lm if i want to do regression with the whole design matrix? My y
is the first column.
the left are xs.
Thanks a lot.
Do you have the data stored in a
Tim,
Thanks for the suggestions. I tried them both, and while neither solved my
problem directly, they helped me to get to a solution.
What I've realized is that partialPlots, for some reason, always looks to
the Global environment when evaluating the x.var argument. So any
variables created
Hi,
It really depends, but almost always to plot a graph you would use
plot() (or some similar graphing function). print() is usually more
for displaying data in the console, not graphically. There are some
cases (e.g., inside a for loop), where you need to wrap the call to
plot (or similar
On Tue, Nov 16, 2010 at 11:20 AM, Mike Marchywka marchy...@hotmail.com wrote:
In many cases, the problem is not the R code but actually
getting the data you need to make real-life decisions or form opinions.
I sent a link to one of my rants to the OP on this topic
but it is a bit political so
On Tue, Nov 16, 2010 at 11:09 AM, skan juanp...@gmail.com wrote:
What's the differente betwen using plot and using print in order to
plot a graph?
For example in order to plot the result of a histogram.
Could you give an example?
It seems that you are referring to graphics functions in
Hi Abhijit,
[I've cc'd R-help to keep the discussion on the list]
On Tue, Nov 16, 2010 at 8:06 AM, Abhijit Dasgupta
adasgu...@araastat.com wrote:
Seth,
I was looking for something like this too. I've a question. If
you're reading the data from a connection, does R start reading the
next
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-do-lattice_002ftrellis-graphics-not-work_003f
On Nov 16, 2010, at 12:09 PM, skan wrote:
Hello
What's the differente betwen using plot and using print in
order to
plot a graph?
For example in order to plot the result of a histogram.
Giovanni,
Both matrices describing the points (A and B in my example) are the
same size, so the resulting matrix will always be square. Also, the
equation I'm using is essentially the following identity:
Var(A + B) = Var(A) + Var(B) + Cov(A, B) + Cov(B, A)
All the covariance matrices that
It is often more effective to send questions about lmer or glmer to
the r-sig-mixed-mod...@r-project.org mailing list, which I am cc:ing
on this response.
On Tue, Nov 16, 2010 at 3:25 AM, Annika annika.opperb...@jyu.fi wrote:
Dear list,
I am new to this list and I am new to the world of R.
On Tue, Nov 16, 2010 at 10:58 AM, Aleksey Naumov anau...@gmail.com wrote:
Hi R experts,
I am looking for a simple error handling approach, whereby I could stop
function execution with a customized error message. For example:
for (i in 1:10) {
if (i == 5)
# I'd like to be able to
Just increase the margins on the left side and add the rownames
x - cor(matrix(rnorm(600), 60, 100))
rownames(x)-paste(row, 1:100)
op-par(mar=c(1,5,1,1), xpd=TRUE)
image(t(x[nrow(x):1,]), axes=FALSE)
text(-0.01, seq(0,1,length=nrow(x) ), rownames(x), pos = 2, offset = 0,
cex = .7)
Another
Hi,
Did you also load the package (you'll need to load it every session
you want to use it)?
## load reshape package
library(reshape)
## now try
merge_all(Br, Ki, Lu, Pr, by=Genes)
If this does not resolve your problem, please run: sessionInfo()
at your console and report the output (it will
Thanks every one first.
I have a dataset like this
YearMonth Day HourDIR SPD (m/s) SKYCOVERTEMP
(C)RH(%)
20091 1 0 310 13.858243 -5 23
20091 1 1 330 14.305283
On Tue, Nov 16, 2010 at 9:40 AM, Jeff Bassett jbass...@cs.gmu.edu wrote:
Giovanni,
Both matrices describing the points (A and B in my example) are the
same size, so the resulting matrix will always be square. Also, the
equation I'm using is essentially the following identity:
Var(A + B) =
Try this:
aggregate(TEMP ~ Year + Month + Day, DF, mean)
On Tue, Nov 16, 2010 at 4:48 PM, facehappywy yiw...@gmail.com wrote:
Thanks every one first.
I have a dataset like this
YearMonth Day HourDIR SPD (m/s) SKYCOVER
TEMP (C)RH(%)
20091 1
Dear All,I would like to do multiple regression in R. I used: lm(y~X),
where y is a n by 1 vector, and X is a n by m matrix. I kept getting the error
message:Error in model.frame.default(formula = y ~ X, : invalid type (list)
for variable 'X'. However, when I used:
On Nov 16, 2010, at 1:04 PM, Kevin Ummel wrote:
Sorry, I shouldn't have used 'sum' as an example; I am looking for a
solution in the case of functions that do not result in simple
vectors or matrices.
The real-world example is something like the following using the
'raster' package,
Ravi Varadhan and I have been looking at UCMINF to try to identify why it gives
occasional
(but not reproducible) errors, seemingly on Windows only. There is some
suspicion that its
use of DBLEPR for finessing the Fortran WRITE() statements may be to blame.
While I can
find DBLEPR in Venables
Try this:
y - rnorm(100)
X - matrix(runif(100 * 10), ncol = 10)
lm(y ~ ., data = cbind.data.frame(y, X))
On Tue, Nov 16, 2010 at 5:07 PM, Tianchan Niu n...@isis.georgetown.eduwrote:
Dear All,I would like to do multiple regression in R. I used: lm(y~X),
where y is a n by 1 vector, and X is a
Tianchan -
Your X is not a matrix -- it's a dataframe. Probably the
simplest solution is to use
lm(y~as.matrix(X))
but you should also learn the difference between a data frame
and a matrix.
- Phil Spector
My problem is that I have a data set for every day of measurement in a
seperate file and I want to plot one parameter of the data for all the days
in one graph. I tried to use for loop but only the last data remains in the
program memory, I don`t know how to plot each day`s data continusly after
Sorry, I shouldn't have used 'sum' as an example; I am looking for a solution
in the case of functions that do not result in simple vectors or matrices.
The real-world example is something like the following using the 'raster'
package, though (I think) any command producing an S4 object could
Hi there:
I am a total beginner in R, and I have a simple question:
I have a table with thousands of lines that represent locations, and two
columns: latitude and longitude. I need to randomly sample 1000 lines. How
do I do it? I know the command sample, but it samples elements
independently, not
My objective is to look at differences in two species of fish from morphometric
measurements. My morphometric measurements are head length, eye diameter,
snout length, and measurements from tail to each fin. I want to use
discrimanant function analyis to determine if there are differences
Fabulicious! It worked!!!
One more question, in the following data frame as posted above:
SubIDCSE1 CSE2 CSE3 CSE4 WSE1 WSE2 WSE3 WSE4
1 6 5 6 2 6 22 4
2 6 4 7 2 6 6
Helllo,
I am having trouble with running the library Playwith in the R-2.12.0. running
under 32bit Windows XP. After calling the library the error message The
procedure entry point gdk_cairo_reset_clip could not be located in the dynamic
library libgdk-win32-2.0-0.dll. occurs and the R asks for
On Nov 16, 2010, at 12:13 PM, Mariana wrote:
Hi there:
I am a total beginner in R,
Time to read the Posting Guide. (Especially since Nabble appears to
deficient in providing an introductory document, at least judging from
the email behavior of some of its users.)
and I have a simple
Thanks, David. That does, indeed, work. It didn't occur to me that a list would
do the job as an argument.
Thanks for the fix!
kevin
On Nov 16, 2010, at 6:58 PM, David Winsemius wrote:
On Nov 16, 2010, at 1:04 PM, Kevin Ummel wrote:
Sorry, I shouldn't have used 'sum' as an example; I
Dear List,
this may be a Newbi question and may have been asked several times, but i am
too stupid to find the posts.
I have a plot of values against POSIXct time steps. If I want to add a second x
axis to the top margin of the plot, only numbers are at the tickmarks. Is there
a
Using David's examples:
do.call(stack, lapply(obs, get))
On Tue, Nov 16, 2010 at 4:04 PM, Kevin Ummel kevinum...@gmail.com wrote:
Sorry, I shouldn't have used 'sum' as an example; I am looking for a
solution in the case of functions that do not result in simple vectors or
matrices.
The
On second glance, while it works for 'stack', it doesn't appear to work for
'overlay':
overlay(sapply(obs,get))
Error in function (classes, fdef, mtable) :
unable to find an inherited method for function overlay, for signature
list, missing
At this point, this may be more of a 'raster'
Hi,
Can anyone suggest a clever way to compute a rolling weekly average of the
columns in a matrix? The column bit is straightforward use apply given a
function which does what you want on a column. With regard to a particular
column, the obvious way is to run a for loop indexing the last 7
Try this:
DF[sample(seq(nrow(DF)), 1000),]
Where DF is your data
On Tue, Nov 16, 2010 at 3:13 PM, Mariana mvale@gmail.com wrote:
Hi there:
I am a total beginner in R, and I have a simple question:
I have a table with thousands of lines that represent locations, and two
columns:
Hi Josh,
Thanks for your reply.
While I am loading the package, it says package plyr is
required and I tried to install plyr package but I am unable to do so
because it is giving an error message as shown below.
Please help me.
library(reshape)
Loading required
assuming your data takes the form of
locationlatitudelongitude
string num num
string2 num num
try:
sub - dat[sample.int(length(dat$location), 1000),]
--
Jonathan P. Daily
Technician - USGS Leetown Science Center
11649 Leetown Road
Tena koe Nasrin
Try points() instead of plot() in your second and subsequent calls to plot().
points() and lines() adds to the current plot by default. Of course you may
have difficulties with setting the x and y limits by that's another matter.
HTH
Peter Alspach
-Original
Hi Josh,
I thought of giving up and started writing code in excel using
VBA but then I saw ur message and gave a try in R. I got it, the location
from which I am calling the package plyr is not working out so called from
different CRANmirror location and then I was able load the
I am using \Sexpr to include a variable in a title of a Sweave document:
\documentclass[a4paper]{article}
echo=false, results=hide=
#mytitlevar - Stuff # case 1, everything is find
mytitlevar - Stuff_first # case 2, f is turned into sub-text
@
\title{MyTitle: \\ \Sexpr{mytitlevar} }
Greetings,
May you please suggest a package or function to use for fitting a GLMM
(generalized linear mixed model) with spatially correlated random effects?
Thank you,
Elijah DePalma
[[alternative HTML version deleted]]
__
A coefficient of -0.4 means that survival times are multiplied by
exp(-0.4), that is, people survival only 67% as long.
-thomas
On Wed, Nov 17, 2010 at 4:32 AM, Vincent Vinh-Hung anhx...@gmail.com wrote:
Thanks for sharing the questions and responses!
Is it possible to appreciate how
On Wed, 17 Nov 2010, Paolo Rossi wrote:
Hi,
Can anyone suggest a clever way to compute a rolling weekly average of the
columns in a matrix? The column bit is straightforward use apply given a
function which does what you want on a column. With regard to a particular
column, the obvious
Ben,
Thank you, it WAS a typo of sorts.
Officially I tried:
glm(count~md+ms+rf+sg+offset(log(Eff)),family=poisson,data=DepthHabGen)
glm(count~md+ms+rf+sg, offset=(log(Eff)),family=poisson,data=DepthHabGen)
(which of course are the same as eachother)
glm(count~md+ms+rf+sg,
On Nov 16, 2010, at 2:32 PM, Jannis wrote:
Dear List,
this may be a Newbi question and may have been asked several times,
but i am too stupid to find the posts.
I have a plot of values against POSIXct time steps. If I want to add
a second x axis to the top margin of the plot, only
On Nov 16, 2010, at 2:33 PM, Paolo Rossi wrote:
Hi,
Can anyone suggest a clever way to compute a rolling weekly average
of the
columns in a matrix? The column bit is straightforward use apply
given a
function which does what you want on a column. With regard to a
particular
column,
1 - 100 of 164 matches
Mail list logo