I'd recommend having a look at the lrm function in the Design package.
Looking at the examples should show you how you can product this type
of plot using some of the other components of this package.
On 26/11/06, Aimin Yan [EMAIL PROTECTED] wrote:
I am trying to fit a logistic regression model
Thorsten Muehge wrote:
Hello,
I would like to plot the following matrix:
Wk=x achsis.
Para 1 = left y-axis as a barplot
para 2 right y-axis as a normal scatter plat.
I could not find such a solution in any of my documentation.
Can someone help me?
Thanks a lot
Thorsten
WkPara
I believe this is specific to the MacOS (and Windows) GUI, which you did
not mention you were using.
See ?flush.console for the solution. Your subject line is misleading: it
is the console that is delaying showing the result.
On Sat, 25 Nov 2006, Kevin Middleton wrote:
As part of a larger
Spencer,
I tried the mixed effects approach you suggest using the random effects
module of
AD Model Builder: (http://www.otter-rsch.ca/admbre/admbre.html). What are
94 unbounded parameters in Schnute et al (1998), now become realizations
of a Gaussian random variable, with the corresponding
Ricardo Rodríguez - Your XEN ICT Team wrote:
Hi all,
I've found quite usefull colored-grid created by image() but I'm facing a
doubt I am not able to solve.
Given the following data rectangle...
1 2 3 4 5 6 7 8 9 10 11 12 13 14
1 12 22 0 7 2 1 0 2 0 2 6
Hello!
I have a data set with 8 columns and in about 5000 rows. What I want to
do is to generate samples of this data set.
Samples of a special size, as example 200.
What is the easiest way to do this? No special things are needed, only
the random selection of 200 rows of the data set.
?sample should tell you what you need to know.
On 26/11/06, Alexander Geisler [EMAIL PROTECTED] wrote:
Hello!
I have a data set with 8 columns and in about 5000 rows. What I want to
do is to generate samples of this data set.
Samples of a special size, as example 200.
What is the easiest
Hello!
Douglas Bates wrote:
...
Not quite. There is now a capability in lmer to specify starting
estimates for the relative precision matrices which means that you can
use the estimates from one fit as starting estimates for another fit.
It looks like
fm1 - lmer(...)
fm2 - lmer(y ~
jim holtman[EMAIL PROTECTED] 24/11/2006 23:21
If you data is a matrix, then try:
image(1:5, 1:14, data.rect)
text(row(data.rect), col(data.rect), data.rect)
Thanks, Jim. It works great and perfectly fill my requirements. I better
understand now how text() works.
By the way, there is a line
Jim Lemon[EMAIL PROTECTED] 26/11/2006 12:05
Hi Ricardo,
This might be what you want (say your data frame is called my.df):
library(plotrix)
color2D.matplot(my.df,c(1,0),c(0,0),c(0,1))
text(rep(0.5:13.5,5),rep(seq(4.5,0.5,by=-1),14),
unlist(my.df),col=white)
and in fact it looks so neat that I
On 11/23/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
I'm trying to teach myself R, and by the way, re-learning statistics using
Crawley's Statistics: an introduction using R.
I've reached the regression chapter, and when it deals with non-linear
regresion using the nls library I face the
Hi,
I have a list of 20 elements, each of them of variable length and with a
structure like this:
lasker[[1]][1:10,]
Var1 Freq
1 1988-023
2 1988-031
3 1988-041
4 1988-052
5 1988-063
6 1988-071
7 1988-081
8 1988-091
9 1989-031
10 1989-04
If DF is the data frame in your post then turn it into a zoo
object, convert that to class ts and then back to zoo replacing
NA's with zero:
library(zoo)
z - zoo(DF$Freq, as.yearmon(as.Date(paste(DF$Var1, 01, sep = -
zz - as.zoo(as.ts(z))
zz[is.na(zz)] - 0
Suggest you read the
Gabor Grothendieck escribió:
If DF is the data frame in your post then turn it into a zoo
object, convert that to class ts and then back to zoo replacing
NA's with zero:
library(zoo)
z - zoo(DF$Freq, as.yearmon(as.Date(paste(DF$Var1, 01, sep = -
zz - as.zoo(as.ts(z))
On 11/26/06, Gregor Gorjanc [EMAIL PROTECTED] wrote:
Hello!
Douglas Bates wrote:
...
Not quite. There is now a capability in lmer to specify starting
estimates for the relative precision matrices which means that you can
use the estimates from one fit as starting estimates for
Hello Andrew Robinson and R-List
Thanks, Andrew, but this does not work. Puttting zero weights on
structural zeros, one's elsewhere in glm() does not deliver the
appropriate expected cell counts loglm() provides as one would expect.
If you run the code provided below, you'll see loglm() delivers
Hello Andrew,
I'm not sure how this is a problem. You can multiply the fitted
values by the weights again, if you wish, or ignore the structural
zeros altogether. The other predicted values all seem to me to be the
same.
I think that is a feature. Sometimes the goal for glm with structural
Hi all,
I just compiled package dna v. 0.2 from source and am having some
difficulty with the window argument of the dotmatrix function.
Some test code; say I'm simulating some nucleotide sequences:
sequence1-trunc(runif(200,0,4))
sequence2-trunc(runif(200,0,4))
I embed a nice 100
On Nov 26, 2006, at 5:23 PM, jim holtman wrote:
The phrase at the bottom (what is the problem that you are trying to
solve) is one that I have used for the last 20 years at work and I
have it as part of my signature line since a lot of people know me
from that phrase. It is one of the first
Hi-
I'm wrestling with some of my data apparently not being called into a GLM or
an LM. I'm looking at factors affecting fish annual catch rates (ie. CPUE)
over 30 years. Two of the factors I'm using are sea surface temperature and
sea surface temperature anomaly. A small sample of my data is
split - sample(2,nrow(dataframe),replace=T,prob=c(0.04,0.96))
dataframe[split==1,] # 200
dataframe[split==2,] # 4800
regards, christian
?sample should tell you what you need to know.
On 26/11/06, Alexander Geisler [EMAIL PROTECTED] wrote:
Hello!
I have a data set with 8 columns and
Douglas Bates wrote:
snip
Don't you find it somewhat disingenuous that you publish a comparison
between the AD Model Builder software that you sell and R - a
comparison that shows a tremendous advantage for your software - and
then you write I am not proficient in R?
I think there is
Sorry. broken Engrish.
I tried to give comment on Sweave by force.
rawcode=T is necessary for an option(If you need comment).
try
wget http://r.nakama.ne.jp/CommentOnSweave/SweaveProfile
R_PROFILE=SweaveProfile R CMD Sweave hoge.Rnw
--
EI-JI Nakama [EMAIL PROTECTED]
Hi,
Can someone please offer a procedure for going from CIs produced by
intervals.lme() to fixed-effects response CIs.
Here's a simple example:
library(mlmRev)
library(nlme)
hsb.lme - lme(mAch ~ minrty * sector, random = ~ 1 | cses, Hsb82)
(intervals(hsb.lme))
(hsb.new - data.frame
You will need a numeric axis, not a factor axis, and control of both
the at= and labels= of the axis. Here is an example using your
data and the lattice xyplot function.
try$visit.position - try$visit
levels(try$visit.position)
levels(try$visit.position)[1] - -20
levels(try$visit.position)
Thanks to Edzer and Roger,
I can now plot with increased font sizes. However, jpeg still does not
reproduce these, nor does it show up in high quality. What I would like
to do is produce some highresolution jpegs.
Any help would be appreciated
Thanx
Herry
R2.4 on Mandriva 10.2 linux.
Dr
[EMAIL PROTECTED] wrote:
Thanks to Edzer and Roger,
I can now plot with increased font sizes. However, jpeg still does not
reproduce these, nor does it show up in high quality. What I would like
to do is produce some highresolution jpegs.
Please specify the commands that you used,
Dears,
I am writing sice having question about your R package copula.
I am using this package to fit my data with some of the 22 copulas mentioned
in Nelsen 1999.
However, when trying random number generator I have got some difficulties.
For example when I have familie number A13 where cupola
28 matches
Mail list logo