Veronika Bókony vbokony at enternet.hu writes:
Hi everyone!
I get the following message anyhow I try to load the package gplots:
Error in lazyLoadDBfetch(key, datafile, compressed, envhook) :
ReadItem: unknown type 241
Error in library(pkg, character.only = TRUE) :
Daniel Coleman coleman.daniel at gene.com writes:
When I run gnls I get the error:
Error in nls(y ~ cbind(1, 1/(1 + exp((xmid - x)/exp(lscal, data = xy, :
step factor 0.000488281 reduced below 'minFactor' of 0.000976563
My first thought was to decrease minFactor but
Brian Lunergan ff809 at ncf.ca writes:
Evening all:
I'm taking a little time to experiment with R, Sweave, and Miktex/LaTex but
I've run up against some problems and -well- I hope that there are some on
the list who might have some suggestions. This will be kind of wordy as I
will
Antje niederlein-rstat at yahoo.de writes:
I'm not sure if I'm in the right place with my question...
I'm running R on Windows and wrote a function and saved it as .R file.
It looks like this:
bmi - function(weight, height) {
bmi - weight / height^2
bmi
}
If I want to use this
René J.V. Bertin rjvbertin at gmail.com writes:
I have some data in a data.frame where for each of a number of
subjects, I have scores for all of a number of symptoms.
Subjects are subdivided in a number of groups, which have unequal sizes.
I'd like to plot between-group differences in
Michael Papenfus mmpapenf at wisc.edu writes:
I have defined the following function:
fr-function(x) {
u-x[1]
v-x[2]
sqrt(sum((plnorm(c(3,6),u,v)-c(.55,.85))^2))
}
which I then solve using optim
y-optim(c(1,1),fr)
y$par
[1] 1.0029771 0.7610545
This works fine.
Now
Peter Eiger Peter.Eiger at gmx.net writes:
I've got a problem with RODBC and saving (sqlSave) of a dataframe in Access.
R 2.0.1 is running on windows XP.
When executing the examples in R help for the USArrests data set sqlSAve
works fine, but running
sqlSave() for a dataframe Adat
Linda Lei llei at bccrc.ca writes:
When dealing with the categorical covariates (for example 3 groups), it
will come out different results if we add the command factor in front
of the categorical covariate or not
The catch is here: the covariate is not a categorial, but a number. In that
Dear R-Listeners,
as the Sweave faq says:
http://www.ci.tuwien.ac.at/~leisch/Sweave/FAQ.html
creating several figures from one figure chunk does not work, and for
standard graphics, a workaround is given. Now I have a multipage trellis
plot with an a-priori unknown number of pages, and I don't
Erin Hodgess hodgess at gator.dt.uh.edu writes:
Has anyone put R on a web server any time, recently, please?
(Red Hat Linux)
The University of Montana put a version up in 2003, but
I was wondering if anyone had done so, please?
Also, where would I find information on such an
Douglas Bates bates at stat.wisc.edu writes:
.
I encourage users of lmer who wish to determine the precision of the
estimates of the variance components to create a Markov chain Monte
Carlo sample of the parameters and evaluate the HPDintervals.
sm1 - mcmcsamp(fm1, 5)
jenny tan jennytimp at hotmail.com writes:
May I know how does one superimpose the 95% confidence and prediction
intervals on the linear regression line of a scatterplot?
Check
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/8380.html
Dieter
__
Li, Hua HUL at stowers-institute.org writes:
I am using extract.lme.cov to extract the covariance matrix of lme. But
the results are not expected.
b - lme(travel~1,Rail,~1|Rail)
The default correlation for lme is no correlation within groups.
extract.lme.cov(b,Rail)
The part of
zhijie zhang epistat at gmail.com writes:
The default behavior in R when performing a regression model with missing
values is to exclude any case that contains a
missing value? How could i set the bahavior that R deal with missing values?
Difficult to say, since there is no such thing as a
Maria Wang mariayzw at hotmail.com writes:
There is a dataset with Latitude Longitude information and corresponding
value shaped like this:
Lati Longi Value
18.0001 -159.6667 123
18.0023 -159.6665 124
I
Mag. Ferri Leberl ferri.leberl at gmx.at writes:
I want to export data to LaTeX. As I want to employ the data as freely as
possible I want to avoid the xtable-command and instead generate some List
like
\MyOwnPrettyCommand{Adam}{Auer}{17}
\MyOwnPrettyCommand{Bertram}{Bauer}{14}
Xu, Xiuli (NIH/NHLBI) [E] xux3 at nhlbi.nih.gov writes:
I would really appreciate it if someone can give suggestions on how to
do spectra fitting in R using ordinary least square fitting and
non-negativity constraints. The lm() function works well for ordinary
least square fitting, but how
Matthias Braeunig mb.atelier at web.de writes:
How can I set the ylim in individual panels of a multiple plot.ts?
The panels are all scaled to fully fit the series inside. But I need
specific scales, and colors
Here is an example:
plot.ts(cbind(1:10,1:10/10),ylim=c(0,3)) #
flash johny flash.johny at gmail.com writes:
Hi,
I am working on project to build a model for analysing Mortgage-backed
securities.
I am doing this as a part of my summer project.
Try
http://cran.r-project.org/src/contrib/Views/
Dieter
__
Uwe Ligges ligges at statistik.uni-dortmund.de writes:
Anyway, these days Tinn-R seems to be the Windows way to go. It
already supports R, LaTeX and SWeave.
A somewhat underused feature in Tinn-R is background color. I use a very light
(5%) gray for the Rcode = sections in rnw-file, giving a
Erin Hodgess hodgess at gator.dt.uh.edu writes:
I have used the latex function from the Hmisc
package and it is just great!
However, I have a new question regarding that function:
is there an option for summary(lm(y~x)), please? There are
options for different types of objects, but I
Angel Roman aroman6 at gmail.com writes:
I've got a problem with a PHP script, in which I call the system function
(to call another processes). The call is :
system(R --slave --vanilla path/to/source/source.R);
...
but if I change the lines to:
png(file=/path/to/file/file.png);
Stephan Lindner lindners at umich.edu writes:
The problem is to create new variables from a data frame which
contains both individual and group variables, such as mean age for an
household. My data frame:
df
hhid h.age
1 1001002023
2 1001002023
...
where hhid is the
Peter Dalgaard p.dalgaard at biostat.ku.dk writes:
?ave
Finally, someone remembered it!
I feel ashamed. May I suggest to add a link under by and/or tapply ?
Dieter
__
R-help@stat.math.ethz.ch mailing list
Marc Bernard bernarduse1 at yahoo.fr writes:
I am wondering on how to use the abline.xyplot with xyplot such that I will
have different vertical lines for
each panel. More sepcifically, suppose that the xyplot generates 4 panels
defined by the combination of
two binary variables: X_1 and
Marc Bernard bernarduse1 at yahoo.fr writes:
I am looking for a function to extract, from an nlme object, the estimated
variance-covariance matrix of the random effects.
VarCorr
D.Menne
__
R-help@stat.math.ethz.ch mailing list
Michael Conklin michael.conklin at markettools.com writes:
I have searched through the help files and I have been unsuccessful in
solving this problem.
I am trying to create a small wrapper function that will go around a
call to a plot function and create a windows metafile in the
mahendra mahfood mahendra025 at yahoo.com writes:
i cant get the character output of R in delphi.
i have got the graphic of R in delphi. i wanna get the character output in
R for my delphi application. anybody knows???
the example source please?
I have written a little wrapper for Delphi
Lusk, Jeffrey J jjlusk at purdue.edu writes:
using this approach, but the output for the fixed effects doesn't report
a p-value or the degrees of freedom (unlike the examples listed in
Faraway's book, which I tried and got the same incomplete output). Any
idea how I can get the complete
Yves Magliulo ym at climpact.com writes:
* http://www.itp.phys.ethz.ch/econophysics/R/index.html
R-financial
Not that I want to discourage you, but be warned, on the download page it says;
Download ...
IMPORTANT WARNING: NOTE, THAT THE RMETRICS PACKAGES ARE NOT COMPLETE, IN
Dear R-Core,
after switching to 2.3.0, all my trusted do.call constructs that worked in
2.2 and earlier fail. I noted that changes were introduced to do.call, but I
could not find out how these relate to my problem.
The following example works in 2.2 and earlier, but fails because rownames
are
Farrel Buchinsky fbuchins at wpahs.org writes:
So how do you reply to someone's posting such that the reply lands up in the
original thread. I subscribe to the list with a daily single e-mail. I guess
if I got every message in real-time it would be easier since then you simply
reply to the
Anamika Chaudhuri canamika16 at yahoo.com writes:
I was trying to calculate an age adjusted rate for which we have the
formula as
sum(weight*rate)/sum(weight). I am using the framingham dataset called fram
here.I cant figure out how I
could change the logical statement I have for the
Leon amstat2006 at gmail.com writes:
Hi, All
I am using is.integer function to examine whether an object is an integer or
not, but I get such results,
x-3
is.integer(x)
[1] FALSE
x - 3
typeof(3)
[1] double
This may not be a wise decision in hindsight, but probably was made to avoid
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
I am using the Deepayan's Sweave trick to set graphics parameters for all
graphs:
ltheme = canonical.theme(color=TRUE)
sup = trellis.par.get(superpose.line)
ltheme$superpose.line$col = c('black',red,blue,#e3,green,
gray)
zhongmiao wang zhongmiao at gmail.com writes:
I am studying the effect of schools on student achievement growth over
time. School effect is random. The effects of schools in prior years
are assumed to be persistent till the current year. Thus, the total
school effect in the second year is
zhongmiao wang zhongmiao at gmail.com writes:
I keep getting the error message Error in lme.formula(fixed = Score ~
factor(time) - 1, data = ldata, random = list(dumid = mat), :
iteration limit reached without convergence (9)
Is there a way to set up starting values, especially the
I am using the Deepayan's Sweave trick to set graphics parameters for all
graphs:
ltheme = canonical.theme(color=TRUE)
sup = trellis.par.get(superpose.line)
ltheme$superpose.line$col = c('black',red,blue,#e3,green,
gray)
Works perfectly, there is only a minor nuissance that
Fredrik Karlsson dargosch at gmail.com writes:
I am using Sweave to generate figures and tables for my thesis. Now, I
have run into the problem that all of a sudden, Sweave produces not
only an \includegraphics command, but a
\include command as well. This breakes LaTeX processing.
Check
andrea valle valle at di.unito.it writes:
I'm using scatterplot3d. I want to add dotted lines parallel to the 3
axes, in order to create a sort of cage in the 3D box (so to represent
subdomains in my dataspace).
I saw that points3d can be used to draw points or lines to the existing
Brian Quinif bquinif at gmail.com writes:
I used the round() function to round some estimates to 3 decimal
places. I then sent put the rounded estimates in a matrix and used
latex() to make a LaTeX table from them. However, in my table, there
are estimtes which only have 2 decimal places.
marilyn ford m.ford at griffith.edu.au writes:
I have just started using R and am stumped over something:
I am using the lme4 package, using lmer (and putting it in to fit0.lmer
) and then I do an anova,
but I am not getting the full table. This is what I get:
fit0.lmer -
Katherine A Grieve grieve at u.washington.edu writes:
I am wondering how to obtain SE estimates for fixed effects from a nonlinear
mixed effects model?
It's a bit of fiddling, but there is an example coming close to what you want
on page 373 of Pinheiro/Bates, and in file
skim033 at student.ucr.edu writes:
Hi, I am writing a function that includes 'sum' function
such as:
f-function(x){
c-c(-1,0,1)
f-sum(c+x)
}
expecting f to be -1+x+0+x+1+x=3x. But I found out that f is
sum(x). So, f is always a scalar, which means that f(c(0,1))
is not a vector as
Uwe Ligges ligges at statistik.uni-dortmund.de writes:
You get the user coordinates of the plotting region by
par(usr)
Now simply make the legend right of that plotting region, e.g. with
x corrdinates at
par(usr)[2] + epsilon
and y coordinates at
mean(par(usr)[3:4])
I always
Matt Goff goff at nawwal.org writes:
The problem is that boxplot is displaying groups that are empty in the
plot.
Call factor() again on the groups, which will drop levels. You can do that in an
extra line, or on-the-fly:
data-data.frame(values=c(1:25), groups=rep(c(A,B,C,D,E),
Andreas Svensson andreas.svensson at bio.ntnu.no writes:
I had a suspicion that you can't have the lme4 package loaded when using
lme (from the nlme package), and lo! I get the full summary of lme only
if lme4 is NOT loaded.
Yes, currently the two don't coexist well, so better make sure
G. Alex Janevski galexski at umich.edu writes:
points(y~x, pch=*, col=black)
lm(y~x)
fm=lm(y~x)
abline(fm, col=red)
This works. The problem arises in that I would like to run my simulation
multiple times, to plot the data points together on the same plot, and
more importantly the
Brian Quinif bquinif at gmail.com writes:
What I need to do is make tables with the std. errors beneath the
estimates in parentheses (standard econ style). How can I make a
dataframe with that format?
Best prepare the output in R as strings; there are alternatives in LaTeX, but
it can be
CG Pettersson cg.pettersson at evp.slu.se writes:
Fitting and reducing a first model for grain yield went smoothly. When I
wanted to look at the fixed effects with estimable() I got an error
message claiming that I was using the wrong variable names, estimable
wanted the variable names in
Gabor Grothendieck ggrothendieck at gmail.com writes:
odbcConnectionDbase(c:\\documents and settings\\egc\\desktop\\test.dbf)
or use simple forward slashes, which are easier to read and work even under
Windows. (Just for the record, looks it was not the problem)
Dieter
To create a more end-user readable table captions for modeld, I would like
to get rid of the I(...) construct in formulae (what's the hell does the
I(..) mean in the contrast table)
Example:
effect ~ I(sqrt(nitro))*treat + I(nitro^2)
should giv
effect ~ sqrt(nitro)*treat + nitro^2
In know,
Sylvain Clément sylvain.clement at univ-lille3.fr writes:
I search the archive for finding a simple
solution for using TukeyHSD with a multistratum
aov result (a repeated emasure anova). The Question have
been asked but I've found no clear answer.
res-aov(y~Fa*Fb+Error(Subject/(Fa*Fb)) )
Gabor Grothendiek wrote:
read.table(textConnection(sub(pat, '\\1 \\2', patid)), as.is = TRUE)
V1 V2
1 ALAN 334
2 AzD 44
3 NA
Looks like you are suffering from a creativity push. I like that one, with no
success I had tried the variant
sub(pat, c(\\1, \\2') patid))
before.
I can't get the PERL subexpression translated to R. Following, for example,
B. Ripley's
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/58984.html
I am using sub, but it looks like an ugly substitute. Assume I want to
extract the first alpha part and the first numeric part, but only if they
are
Dieter Menne dieter.menne at menne-biomed.de writes:
patid=c(ALAN334,AzD44,44AZD)
txt =sub(([[:alpha:]]+)([[:digit:]])+,\\1,patid)
num =sub(([[:alpha:]]+)([[:digit:]])+,\\2,patid)
Sorry, a ) was at the wrong place. Here the corrected version.
patid=c(ALAN334,AzD44,44AZD)
txt
Gabor Grothendieck ggrothendieck at gmail.com writes:
In the third case there is no match so there are no
substitutions. Handle it separately:
pat = ^([[:alpha:]]+)([[:digit:]]+)
result - cbind(txt = sub(pat, \\1, patid), num = sub(pat, \\2, patid))
result[regexpr(pat, paid) 0,] - NA
Gabor Grothendieck ggrothendieck at gmail.com writes:
We could use sapply to reduce it slightly:
result - sapply(sprintf(\\%d, 1:2), sub, pattern = pat, x = patid)
result[regexpr(pat, patid) 0,] - NA
Looks like we should make a generalized wrapper for it.
Dieter
Leonardo D Bacigalupe L.Bacigalupe at sheffield.ac.uk writes:
I'm trying to plot e.g. fitted vs residuals for each level of the
random effects, and i'm getting the same error.
I guess this might be a problem of the graphic capabilities of R.
Is there any way to obtain those plots?
Peter Dalgaard p.dalgaard at biostat.ku.dk writes:
alist - lapply(as.list(allV), as.name)
names(alist) - allV
# next line
alist - c(as.list(as.name(data.frame)), alist)
mode(alist) - call
data - eval(alist, sys.parent(1))
Uwe Ligges ligges at statistik.uni-dortmund.de writes:
I want to attach a fixed legend to mosaicplot, e.g. -3..3 independently of
the data in the graphics. This must be hidden somehow in shading/legend, but
I could not get the right handle to it.
Example:
plot(1:10)
Hello, friends of mosaicplot,
I want to attach a fixed legend to mosaicplot, e.g. -3..3 independently of
the data in the graphics. This must be hidden somehow in shading/legend, but
I could not get the right handle to it.
Dieter
__
Peter Tait petertait at sympatico.ca writes:
Is it possible to remove the NA level from a factor with in
read.table()? If not what is the most efficient way of doing this?
inrange-read.table(C://...,header=T,sep=|,colClasses=c( id=factor))
See parameters na.strings in read.table.
Michael Kubovy kubovy at virginia.edu writes:
When I issue the command
latex(anova(raw1.lmer0, raw1.lmer, raw1.lmerI), file = 'raw1.tex',
rownamesTexCmd = c('baR', 'addit', 'multip'), longtable = F, dcolumn
= T, booktabs = T, t able.env = F, colheads = NULL, colnamesTexCmd = c
Michael Kubovy kubovy at virginia.edu writes:
Unfortunately this doesn't work.
rowname = NULL
suppresses row names altogether. I want to substitute different row
names
This example works for me. Take into account that rownamesTexCmd adds
the leading \, so you cannot use it
Rick Bilonick rab45+ at pitt.edu writes:
When I try to run the example from Variogram with an lme object, I get
an error (although summary works):
library(nlme)
fm1 - lme(weight ~ Time * Diet, BodyWeight, ~ Time | Rat)
Variogram(fm1, form = ~ Time | Rat, nint = 10, robust =
BORGULYA Gábor borgulya at gyer2.sote.hu writes:
I couldn't figure and couldn't google out how to make construct a pair of
\Sexpr s or a LaTeX macro that would include
5\cdot 10^{-12}
into the LaTeX output istead of
5e-12 .
a =1.3452e-12
asp =
Dear R-Listers,
(well... called Depayan)
The standard example
barchart(yield ~ variety | site, data = barley,
groups = year, layout = c(1,6), stack = TRUE,
auto.key = list(points = FALSE, rectangles = TRUE, space = right),
ylab = Barley Yield (bushels/acre),
scales = list(x =
Felipe Martínez-Pastor Felipe.Martinez at uclm.es writes:
I am analyzing some data obtained after measuring some parameters at
different times in samples obtained from many subjects. The model is
quite simple: aov(parameter ~ Time + Error(Subject/Time))
Now I want to make a pairwise
voodooochild at gmx.de writes:
i want to plot a count process in the following way, but i don't know
how i can do that. i have data for example x-(0,2,6,2,8,4,.) and
dates y which is a vector of weekly dates for example
(01/01/06, 08/01/06, 15/01/06, 22/01/06, ), now i want to plot
Carsten Steinhoff carsten.steinhoff at gmx.de writes:
I want to find the parameters mu and sigma that minimize the following
function.
It's important, that mu and sigma are strictly positive.
-
optimiere = function(fmean,smean,d,x,mu,sigma)
{
merk = c()
for (i in
Mahdi Osman m_osm at gmx.net writes:
Hi list,
I am fitting microarray data (intensity) model using the lme package in R
environment. I have 5 fixed variables in the model. One of the fixed
variables is genes. I am trying to get p-values for different genes. But I
am getting only one
Denis Chabot chabotd at globetrotter.net writes:
Quantreg package
But I'd like to extract the predicted values. The help for
predict.qss1 indicates this:
..
predict.qss1(object, newdata, ...)
and states that newdata is a data frame describing the observations
at which
Original question...
is it possible to pass additional parameters to the model function that are
known and groupwise constant with nlsList? I could not find something like a
keep this fixed option in the documentation and the code (my fault...?)
Spencer Graves spencer.graves at pdf.com
into groups and use
globals to pass the constant parameters, but it is ugly code and won't work
when an over-all nlme is needed.
Dieter Menne
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Aleš Žiberna ales.ziberna at gmail.com writes:
?par
Look for arguments
mfrow and mfcol
Original Post
[mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of R Granell,
Medicine
Is there a way of obtaining more than one histogram in the graph-window so
you can edit all the plots
Michael comtech.usa at gmail.com writes:
plot(testError, col=red)
lines(testVar, col=black)
Only one plot (the red one) appear on the Window, the black line did not
appear...what's wrong?
We can only guess, because you did not supply the data. But I am quite sure,
that the lines-data
Paul Johnson pauljohn32 at gmail.com writes:
echo=F,quiet=T,print=F=
library(Design, verbose=F)
at
But in the output there is still one page of banner information, which
starts like this:
Loading required package: Hmisc
Hmisc library by Frank E Harrell Jr
Type library(help= Hmisc
Laura Quinn laura at env.leeds.ac.uk writes:
I was wondering whether there is a way to select random samples from a
data matrix, retaining the indexing for the rows and columns?
Is this you are looking for? Dieter
# create data
n = 10
world = array(rnorm(n*n*n),c(n,n,n))
nsamples=12
#
Dear Latticers,
I want to position two trellis plots of different forms
side-by-side. The plot types are slightly different,
aspect=1 required, but panels should look the same in
both plots. Current workaround uses a guessed factor.
Any way to improve this?
Dieter
library(lattice)
n1 = 20
Thanks, Deepayan
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
If that's the only reason, I would suggest rbind-ing them and then use
scales = list(x = free)
I probably will go for this.
If you want the first two columns to have the same x-limits, you can specify
xlim =
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
Documentation said that using panel.width with non-standard aspect (which
must
be 1) leads to undefined effects.
Yes, because panel.width and panel.height together define the aspect,
and will override whatever the 'aspect' argument
Dear RWeavers,
to speed up report generation with Sweave, I am looking for an option
skipifexists for figures. A named figure chunk should not be executed if
the file to be produced already exists.
I know that this could be done with hooks, but I want to avoid to re-invent
this wheel.
Dieter
John Logsdon j.logsdon at quantex-research.com writes:
I am trying to extract the solution from a simple lme calculation.
y-c(-1.118,-.5,.5,1.118,10)
gp-factor(c(rep('one',4),'two'))
res-lme(y~1,rand=~1|gp)
...
Random effects:
Formula: ~1 | gp
(Intercept) Residual
N.W.Galwey nw.galwey at ukonline.co.uk writes:
The syntax below works, and gives the expected results, in R version 2.0.0,
provided that I have loaded packages Matrix, latticeExtra and lme4.
However, I cannot get it to work in version 2.2.1. Can anyone tell me how
to fit this model in the
Martin Maechler maechler at stat.math.ethz.ch writes:
One thing -- very useful for you -- I forgot to add:
You can easily quickly revert the other cbind/rbind
activation by using
methods:::bind_activation(FALSE)
so you don't need to unload lme4 or Matrix, and you can
reactivate
A larg program which worked with lme4/R about a year ago failed when I
re-run it today. I reproduced the problem with the program below.
-- When lme4 is not loaded, the program runs ok and fast enough
-- When lme4 is loaded (but never used), the do.call fails
with infinite recursion after 60
Peter Dalgaard p.dalgaard at biostat.ku.dk writes:
A larg program which worked with lme4/R about a year ago failed when I
re-run it today. I reproduced the problem with the program below.
-- When lme4 is loaded (but never used), the do.call fails
with infinite recursion after 60
Dieter Menne dieter.menne at menne-biomed.de writes:
Dear Rweavers,
When generating reports with Sweave, I would like to quote some results in
the abstract (Something like The treatment effect is 10 mbar, see page
100).
Currently, I used verbatimwrite and friends to write to multiple
Dear Rweavers,
When generating reports with Sweave, I would like to quote some results in
the abstract (Something like The treatment effect is 10 mbar, see page
100).
Currently, I used verbatimwrite and friends to write to multiple files to be
included, but I wonder is there is a more elegant
Ken Beath kbeath at efs.mq.edu.au writes:
I'm maximising a reasonably complex function using nlme (version
3.1-65, have also tried 3.1-66) and am having trouble with fixed
parameter estimates slightly away from the maximum of the log
likelihood. I have profiled the log likelihood and
David Whiting david.whiting at ncl.ac.uk writes:
The nearest I got to was to be able to colour row. Take a look at Table
10 here:
http://biostat.mc.vanderbilt.edu/twiki/pub/Main/StatReport/latexFineControl.pdf
After some private emailing with David and downloading the latest version of
Zev Ross zev at zevross.com writes:
Using glmmPQL you can extract the full table of estimates, SE, p-values
etc using as an example:
mymodel-glmmPQL(mymodel here)
summary(mymodel)[[18]]
How can I pull this table out of a lmer object in lme4?
Looks like it is not so easy, the code
Dear latex/R-Sweavers,
Using the codel below, I can color text in individual cells for latex
output.
Is there a similar way to get a background shading? My attempts failed
because I did not get the closing brace at the right place with Hmisc/latex.
library(Hmisc)
x -
Michael Kubovy kubovy at virginia.edu writes:
I run
e7.lmeb3 - lme(baLO ~ bar + factor( delta), data = e7, random = ~ 1
| sub, method = ML)
... cut
how can I get the CIs for the fixed effects in the 20 cells of the
bar * delta design?
A typo, maybe? Your design is
Leo Gürtler leog at anicca-vijja.de writes:
hlm - groupedData(laut ~ design | grpzugeh, data = imp.not.I)
the grouped data object is located at and can be downloaded:
www.anicca-vijja.de/lg/hlm_example.Rdata
...
2) What is the cause of the convergence error which seems to
Ann Hess hess at stat.colostate.edu writes:
Here is my model statement:
Model-lm(log2PM~P+T+P*T)
where P has 16 levels, T(treatment) has 12 levels and I am interested in
looking at different treatment comparisons.
With so many levels, you should do something against fishing for error
Using Hmisc, I get the following output from format.df (which will feed to
latex)
Min Q1 Med Mean Q3 Max
Age 30 32.5 36 37.50 41.00 49
Height 174 175.0 178 179.00 181.00 188
Weight 68 74.0 78 76.67 79.75 83
Automatic decimal points adjustments does not
Dieter Menne dieter.menne at menne-biomed.de writes:
Using Hmisc, I get the following output from format.df (which will feed to
latex)
Min Q1 Med Mean Q3 Max
Age 30 32.5 36 37.50 41.00 49
Height 174 175.0 178 179.00 181.00 188
Weight 68 74.0 78
Mahdi Osman m_osm at gmx.net writes:
How can I stop my xemacs from lauching R whenever it starts? I want to use
the usual M-x R instead.
This is something I have to edit in the init file.
Remove
(start-R)
at the end of the init file (assuming you use the John-Fox Version).
Dieter
101 - 200 of 287 matches
Mail list logo
