DeaR
I am trying to use a dynamically create formula with nlsList and nlme, but I
cannot get the environment of the string-generated formal to work similarly
to the manually entered one.
Any idea?
Dieter
#-
library(nlme)
# Pinheiro/Bates p 280
fm1Indom.lis =
MASFERFC Team masferfc at gmail.com writes:
I'm currently running MySQL 5 queries with PHP 5.2.3 and returning results
to the end-user in web-page tables of ca. 50 rows by 5 columns. I'd like to
(more or less) simultaneously return to the browser a couple of canned
charts and graphs based on
Felipe Carrillo mazatlanmexico at yahoo.com writes:
How can the R graphics window be customized
programmatically? Either minimized,maximized or change
the size of the default that ships with R.
Assuming Windows as OS (which you forgot to mention)
windows(width = 2, height = 2, pointsize =
jim holtman jholtman at gmail.com writes:
Turn 'Recording on for the plots.
windows(record=TRUE)
or select from the GUI.
I mean if the object has several plots how can I
get those?
gam.object - gam(y ~ s(x,6) + z,data=gam.data)
plot(gam.object,se=TRUE)
Or use something like
KOITA Lassana - STAC/ACE lassana.koita at aviation-civile.gouv.fr writes:
Could someone help me to combine bwplot and srt option (exemple srt = 45
degree or srt 90 degree)? My graphic contains 146 boxplots, I would like
to label all of them. As you know, labels are not readable.
You
Daniel Lakeland dlakelan at street-artists.org writes:
We have used the lmer package to fit various models for the various
experiments that she has done (random effects from multiple
measurements for each animal or each trial, and fixed effects from
developmental stage, and genotype etc).
Daniel Stahl daniel_stahl at operamail.com writes:
I would like to do a large number of e.g. 1000 paired ttest using the
by-function.
Technical stuff aside: do you really want to do 1000 paired t-tests?
Followed by a *** pick-nic?
Mmm...
Dieter
Ronaldo Reis Junior chrysopa at gmail.com writes:
Hi,
I use panel.points to add points to a xyplot graphic. But I like to use the
sunflowerplot to plot my points because this is very superimposed. It is
possible to use this? I try but it dont work directly. It may be need to put
this
Maciej Hoffman-Wecker Maciej.Hoffman-Wecker at bioskin.de writes:
...
But I succeeded in developing a tiny example, that reproduces the bug
(wherever it is).
I generated a small Access data base test2.mdb with one table Tab1 and
following columns:
.
library(RODBC)
.con -
Arjun Bhandari arb.em at adia.ae writes:
I conduct a univariate regression with lm function. I would like to get
the p value for the regression. Is there a method that would enable me to
extract the p value into a variable.
# From lm docs (Sorry, Annette, for misusing your data)
ctl -
paulandpen at optusnet.com.au writes:
I get this part of the code, it all makes sense.
##
if(nchar(Sys.getenv(LONG_TEST)) != 0) {R=1000} else {R=5}
{
data(customerSat)
surveydat = list(k=10,x=as.matrix(customerSat))
Mcmc1 = list(R=R)
set.seed(66)
Alexis Delevett adelevet at yahoo.com writes:
I am using R to process some community survey data. Several item responses are
recorded via a 7-point
Likert-Scale. As I have coded the responses, 1 represents high agreement, and
7 high disagreement. This
of course impacts the coefficients in a
GOUACHE David D.GOUACHE at arvalisinstitutduvegetal.fr writes:
...
This is my original function :
comp.fit.2-function(tab)
{
fit.log-nls(surf.vert.tot ~ 100/(1+exp(((log(81))/a)*(sum.T.levee-b))),
start=list(
...
I've thought of just storing the 3 model results in a list, and then
Felix Andrews felix at nfrac.org writes:
Package plotAndPlayGTK provides a graphical user interface for R
.
This package is
based on RGtk2, and so requires the GTK+ libraries.
Warning for Windows users: When I installed GTK+ from
Andrew Clegg andrew.clegg at gmail.com writes:
... If I want to demonstrate that a non-linear curve fits
better than an exponential, what's the best measure for that? Given
that neither of nls() or optim() provide R-squared.
To supplement Karl's comment, try Douglas Bates' (author of nls)
baptiste Auguié ba208 at exeter.ac.uk writes:
I'm new to R, but a search through the list didn't quite solve this
problem: I want to draw a few ellipses (or any custom shape for that
matter) at given locations in a graph. I know how to plot points from
my data, set the point type pch to
In the docs of window/grDevices:
All these devices are implemented as windows devices, the _display_
parameter selects which is actually used.
What is the display parameter?
Dieter Menne
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch
Dear UseRs called Hadley, or Paul,
I am trying to print an edited ggplot2/grid graphics to a metafile. With the
commented line below it works, but when I edit the plot by uncommenting the
line, it fails, because it's illegal to have 2 graphics in a metafile. It
works with pdf, but even then I get
Dear Weavers,
Does someone have an example of using eps or any other vector graphics with
odfWeave? It tried the example below (and commented variants) with
simple.odt in the examples directory, and got an error.
Dieter
#---
library(odfWeave)
plotInfo - getImageDefs()
plotInfo$type = eps
Kuhn, Max Max.Kuhn at pfizer.com writes:
Thanks, Max.
What is the version of odfWeave? If you use those specifications in the
current version (0.5.9):
Yes, version 0.5.9 (it was further below in my message, I forgot to put the
odfWeave flag before it.
I had tried your horizontal...
Dieter Menne dieter.menne at menne-biomed.de writes:
library(odfWeave)
plotInfo - getImageDefs()
plotInfo$type = eps
plotInfo$plotWidth = 4
plotInfo$plotHeight = 4
plotInfo$device = postscript
plotInfo$args = list( horizontal = FALSE,onefile = FALSE, paper = special)
setImageDefs
Dear latexRs,
I tried to make a latex printout of a simple categorial ftable. It should
look like the output of print.ftable. Any ideas how to get the syntax of
summary.formula right. Or some alternative? As far I see, xtable does not
have method for ftable.
Dieter
library(Hmisc)
n=500
sex -
Chuck Cleland ccleland at optonline.net writes:
How about a Sweave approach? Something like this in the *.Rnw file:
\documentclass[letter]{article}
...
ftable(symptom ~ treatment + sex)
at
\end{document}
Then Sweave() the *.Rnw file to produce a *.tex file in the working
Yuchen Luo realityrandom at gmail.com writes:
This should be a very common operation and I believe there should be a nice
way in R to handle it. I couldn't find it in the manual or by searching on
line. I am wondering if I could ask for some help in this community.
I am trying to record
M. Jankowski mjankowski at gmail.com writes:
My Sweave code is intended to make lots of plots and create a *.pdf
document. Sweave is perfect for this. The only problem is that I find
myself using the same R code, within my Sweave input file over an
over.
This is the R function:
joris.dewolf at cropdesign.com writes:
Have a look at Rnew volume 5/1 (http://cran.r-project.org/doc/Rnews/) where
Doug Bates explains this nicely. Consider using lme4 for your purpose.
But be aware of http://finzi.psych.upenn.edu/R/Rhelp02a/archive/76742.html
if you currently want to use
Tirthadeep sabya231 at gmail.com writes:
library(ROCR)
plot(roc1)
plot(roc2)
gives two plots on two different graph. Now i want to merge on a single
graph.
plot(roc1)
points(roc2)
Use parameter add=TRUE in all but the first call
Dieter
Dear Latticer,
I want to give individual colors to all elements in a simple stacked
barchart. I know why the example below does not work (and it is a excellent
default), but is there any workaround for this?
Dieter
# This only colors red and green, but I want blue and gray for Peatland.
hadley wickham h.wickham at gmail.com writes:
On 6/12/07, Dieter Menne dieter.menne at menne-biomed.de wrote:
Dear Latticer,
I want to give individual colors to all elements in a simple stacked
barchart. I know why the example below does not work (and it is a excellent
default
Jianping Jin jjin at email.unc.edu writes:
...
I wrote a code to iterate a non-linear fit with a set of data. The entire
process didn't implemented to the end because an error message, singular
gradient. I knew that some sub-sets (columns) do not fit my formula well
and may result in
Fränzi Korner fraenzi.korner at oikostat.ch writes:
sorry for not being more specific. I would like to use R to get a prediction
(with standard error) of the response in a mixed model at selected values of
the fixed-effects factors. Hence, in a mixed model, say, for response body
size with,
Fränzi Korner fraenzi.korner at oikostat.ch writes:
how do I obtain standard errors of the predicted values from a lme (or
lmer)-object?
Not totally clear what you mean. intervals(lmeresult) gives the confidence
intervals for the coefficients. Otherwise, you can do some calculations with
Ken Nussear knussear at mac.com writes:
I'm running a series of mixed models using lme, and I wonder if there
is a way to sort them by AIC prior to testing using anova
(lme1,lme2,lme3,lme7) other than by hand.
You can try something like the following. However, also consider using
W. Katy linlicc at hotmail.com writes:
Can anyone kindly help me with the interpretation of colors in the
RColorBrewer Package (brewer.pal) ?
For example, in brewer.pal command, when picking up the RdBu palette,
does red means larger or smaller value?
Conventionally, the redder colors
Fiona Sammut fiona.sammut at um.edu.mt writes:
Whilst reading a book on robust statistics, in particular a section on
multivariate outlier detection, I came across the functions outmve,
outmgv and outpro. I tried the help.search command to try and find
documentation on these 3 commands in
Aurel Razvan Duica AurelRazvan.Duica at amdocs.com writes:
Can you please provide me with an working example of clustering in R
(not very complicated) and with an example of survey in R.
We would appreciate this very much. We need them as soon as possible.
Sooner is not possible: there are
David Lindelof lindelof at ieee.org writes:
I'm trying to create a 10x7 rectangular grid of individual plots, each
of which is a contour plot that shows a slice from a three-dimensional
matrix.
The problem is that you construct individual plots, and then try to glue these
together. This is a
Kevin R. Coombes krc at mdacc.tmc.edu writes:
Is there a way to turn off the automatic inclusion of ./ at the
beginning of a path specified in an \SweaveInput{} instruction?
Giving a full file path name works for me (Windows, R 2.4.1)
\SweaveInput{C:/tmp/MyTitle.rnw}
Note the if the
Simon Blomberg s.blomberg1 at uq.edu.au writes:
.
I won't discuss the dangers of types of sums of squares and different
contrast codings. That would be tempting the wrath of the gods. See
section 7.18 in the R FAQ. John Fox's Companion book also has a brief
discussion (p. 140).
And
work without echo=T? It is possible
to write a source that dont need echo=T on source command?
Most likely another variant of
http://post.gmane.org/post.php?group=gmane.comp.lang.r.generalfollowup=83924
See: Why do lattice/trellis graphics not work?
Dieter Menne
gyadav at ccilindia.co.in writes:
I am facing the same problem in my case. R 2.4.1 have installed
successfully, but when i try to install the packages from a local zip
file. It gives the following error message.
+++
utils:::menuInstallLocal()
Error in
Alberto Monteiro albmont at centroin.com.br writes:
library(tcltk)
tt - tktoplevel()
done - tclVar(0)
but - tkbutton(tt, text=OK, command=function() tclvalue(done) - 1)
tkpack(but)
tkwait.variable(done)
works as fine as long as I click the OK. However, if I close
the window (by clicking
Duncan Murdoch murdoch at stats.uwo.ca writes:
On 3/29/2007 6:45 AM, Schmitt, Corinna wrote:
- cell array
I don't know what a cell is in C.
This is MathLabish.
Dieter
__
R-help@stat.math.ethz.ch mailing list
christian.ritter at shell.com writes:
The Vista issue is not innocent as it threatens the life of R within large
corporations. So any posts on how R
runs under Vista and what has to be done to make it work and what cannot be
done etc will be very useful.
Main problem seems to be the
be specified
as a factor (when the first level denotes failure and all others success)
or as a two-column matrix with the columns giving the numbers of successes
and failures.
Dieter Menne
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch
Dear R-list,
I have a data set from 20 pigs, each of which is tested at crossed 9 doses
(logdose -4:4) and 3 skin treatment substances when exposed to a standard
polluted environment. So there are 27 patches on each pig. The response is
irritation=yes/no.
I want to determine equally effective
Tim Churches tchur at optushome.com.au writes:
A recent paper from Google Labs, interesting in many respects, not the
least the exclusive use of R for data analysis and graphics (alas not
cited in the approved manner):
http://labs.google.com/papers/disk_failures.pdf
...
For all of you
Friends of lme,
following the documentation of fit.contrast in package gmodels, it should be
possible to use it with lmer. At least with sim.lmer=TRUE and lmer(not 2),
because the latter has no mcmcsamp yet.
I could not get this to work. UseR error?
--
library(gmodels)
Michael Kubovy kubovy at virginia.edu writes:
I would like to suppress the text 'Scatter Plot Matrix' that appears
under the plot. Could someone please suggest how?
Must be as special Virginia Brand of lmList. That Test does not turn up in my
output, and the only place I found it in the
Bart Joosen bartjoosen at hotmail.com writes:
Hi,
I have an integer which is extracted from a dataframe, which is sorted by
another column of the dataframe.
Now I would like to remove some elements of the integer, which are near to
others by their value. For example:
integer: c(1,20,2,21)
yyan liu zhliur at yahoo.com writes:
I used the function lme in R and Splus-7. With the same dataset and same
argument for the function, I got
quite different estimation results from these two software. Anyone has this
experience before?
From the FAQ:
R by default uses treatment contrasts
Liu, Delong (NIH/CIT) [C] liud2 at mail.nih.gov writes:
Instead, I am more interested in obtaining the
coefficient estimates of local polynomial from loess(). Is it
straightforward to modify loess() so that the coefficient estimates can
be put into the return list of loess()?
No need to
Liu, Delong (NIH/CIT) [C] liud2 at mail.nih.gov writes:
I want to extract estimated coeffiicents of each local polynomial at
given x from loess(), locfit(), or KernSmooth(). Can some experts
provide me with suggestions? Thanks.
Try
cars.lo - loess(dist ~ speed, cars)
str(cars.lo)
List
Liu, Delong (NIH/CIT) [C] liud2 at mail.nih.gov writes:
I want to extract estimated coeffiicents of each local polynomial at given x
from loess(), locfit(), or
KernSmooth(). Can some experts provide me with suggestions? Thanks.
Before you start on your own, also note Brian Ripleys recent
Michael Kubovy kubovy at virginia.edu writes:
I tried and it gave a strange result. See
http://people.virginia.edu/~mk9y/mySite/twoGaussian.R
and
http://people.virginia.edu/~mk9y/mySite/twoGaussian.pdf
*
Session Info
*
Raymond Balise iteachtyping at gmail.com writes:
Is there a list of object types used by R? Is there a function that I can
use to tell the type of an object?
typeof()
Dieter
__
R-help@stat.math.ethz.ch mailing list
Daniel Ezra Johnson johnson4 at babel.ling.upenn.edu writes:
I am fitting models to the responses to a questionnaire that has
seven yes/no questions (Item). For each combination of Subject and
Item, the variable Response is coded as 0 or 1.
I want to include random effects for both
Michael Kubovy kubovy at virginia.edu writes:
The following plot is a first approximation to what I need:
.. see below
***
It has four drawbacks. I would like:
(1) to remove the white grid;
(2) to remove the white contours;
(3) its appearance to be
Alan Gibson alan.gibson at gmail.com writes:
Im looking for a simple function to produce a crosstab from a dumped
sql query result.
xtabs in stats is a simple solution, but package reshape is much more flexible
and comes with a good introduction.pdf.
Dieter
liu, jcheng liujcheng at gmail.com writes:
i want to compare the several main effects in a linear model. i wonder
how to do the multiple comparisons for less confident intervals.
Package multcomp by Torsten Hothorn.
Dieter
__
Rainer M Krug RKrug at sun.ac.za writes:
I am using the latex() command from the Hmisc package to export table
to latex. For formating of the values, I use the format() function. But
I don't manage to get the format of the values right: I don't want
scientific notation, but if I say
Knut Krueger Knut-krueger at einthal.de writes:
Switch is working with character , with integer, but not in the third
example
...
# not working
test = c(3,9,3,9,3,9,8,9,8,9,8,9,3,9,3,9,5,9,3,9,1,2,7,9,3,9,1,1,
2,2,3,9,2,1,2,5, 9,8,9,1,2)
count1 - 0
for (i in 1:length(test))
Paulino Perez Rodriguez perpdgo at colpos.mx writes:
How can I solve a non linear system of equations using R?
Have your tried to enter nonlinear equations into http://search.r-project.org?
Dieter
__
R-help@stat.math.ethz.ch mailing list
Wee-Jin Goh wjgoh at brookes.ac.uk writes:
Just to add to this, I also need to know what language is the nnet
package written in? Is it in pure R or is it a wrapper for a C
library.
As usual, you can download the full source to find out what you want, but it's a
bit hidden. Simply
Rainer M Krug RKrug at sun.ac.za writes:
When using format.df with dec, the formating does not work (see code
below).
format.df(x, dec=2)
a b m 1 m 2
1 1.120 3 5 7
2 2.230 4 6 8
This is the same bug introduced in the latest version
Henrik Frederik Thomsen HFTH at DCE.AU.DK writes:
I have encountered at problem when using the Hmisc package to generate
latex-tabels.
I can't seem to control the number of digits in the latex-file. I have
used the Design and Hmisc packages for a while without this problem, but
after
Rainer M Krug wrote:
I have (8 measures * 96 groups) = 768 datasets for which I did linear
regressions using lm().
Now I want to compare the slopes for each of the 8 measures in each of
the 96 groups. As I understand , I can not use
anova(lm1, ..., lm8)
as the lm1 ... lm8 are
Spencer Jones ssj1364 at gmail.com writes:
I am using R 2.3.0, Is there a way to get prediction intervals using
predict.gls?
A r-project site search (alway recommended) would have given you
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/42932.html
Dieter
Bruno L. Giordano bruno.giordano at music.mcgill.ca writes:
I am looking for references about mixed models built on rank transformed
data.
To complement Frank H: Try lrm in package Design.
Dieter
__
R-help@stat.math.ethz.ch mailing list
by a similar problem as in
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/4555.html
Dieter Menne
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
is under development an does not support
this argument YET.
Dieter Menne
__
R-help@stat.math.ethz.ch 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 provide commented
The following perl-script is available from
http://www.menne-biomed.de/download/toggleDoc.zip
After it is tested on other systems (currently only Windows),
I plan to upload it to
http://cran.mirroring.de/other-software.html
toggleDoc.pl: Clutter reduction for S4-R documentation files.
Jeffrey Stratford stratja at auburn.edu writes:
I'm trying to plot the ratio of used versus unused bird houses
(coded 1 or 0) versus a continuous environmental gradient (proportion of
urban cover [purban2]) that I would like to convert into bins (0 -
0.25, 0.26 - 0.5, 0.51 - 0.75, 0.76 -
Pierce, Ken ken.pierce at oregonstate.edu writes:
Does anyone know of any issues with nesting source() calls within
multiple scripts? I have at least one script which always finds errors
when I source it but runs fine when run on its own. It containd source()
calls to other scripts and it
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
On 10/18/06, juan f poyatos jpoyatos at cnio.es wrote:
Dear all,
I am trying to include confidence intervals in a xyplot.
...
Well, panel.locfit doesn't have any options to draw confidence bands,
so you'll have to write a panel
Franco Mendolia franco.mendolia at gmx.de writes:
Is there a possibility in R to save data in pdf-format?
I do not want to save a plot but some lines of simple text.
I remember that about 4 years ago there was a thread on the subject, but I
could not find it any more. But, maybe, it's
to be a problem with SQLServer. The following example with
Access works correctly for me, even if I close the connection in between.
(Win2k, R 2.4.0)
Dieter Menne
-
library(RODBC)
channel - odbcConnectAccess(db1.mdb)
df - data.frame(T=1, S=10)
sqlSave(channel, df, test, rownames=F)
#odbcClose
Jenny Stadt jennystadt at yahoo.ca writes:
if(length(real.d)=30 length(real.b)=30
beta1*beta2*theta1*theta20 )
{ r - 1; corr - 1; }
real.d and real.b are two vectors, beta1,beta2,theta1,and theta2 are
constants. The error occurred like this:
Error in
if (length(real.d) = 30
quotes.
Creating the text there was not much fun, waiting times are simply
unacceptable and the Wiki only give an empty page when syntax errors (for
example from quotes) are detected.
Dieter Menne
__
R-help@stat.math.ethz.ch mailing list
https
Thomas Wutzler twutz at bgc-jena.mpg.de writes:
without grouped data (athough I recommend using grouped data after
reading the book), the call can be formulated as:
Well, this is probably the only point where I disagree with Douglas Bates. I
found using grouped data introduces a level of
Marc Schwartz (via MN mschwartz at mn.rr.com writes:
If the content of this particular transcript is likely to be static,
consider an alternative of making it available in a PDF document that is
linked on that page.
Then others can perhaps contribute by providing other relevant content
in
Frank E Harrell Jr f.harrell at vanderbilt.edu writes:
As far as I know, no one has implemented a latex method for anova.lme
objects.
Well, I have, based on your latex/Hmisc, but it's a bit to my private taste, so
I won't upload it to CRAN. Also has latex.glm, .lme, .summary.lme.
Jue.Wang2 at sanofi-aventis.com writes:
Does any one know why wilcox.exact sometimes doesn't agree with wilcox_test on
the estimate of the
difference of medians in two levels ?
Example removed
See:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/85893.html
Dieter
Gamer, Matthias gamer at uni-mainz.de writes:
Specifically, I have data from a psychometric function relating the
frequency a subject's binary response (stimulus present / absent) to the
strength of a physical stimulus. Such data is often modeled using a
cumulative gaussian function.
AgusSusanto gusanto at gmail.com writes:
I obtained error messages when I run these commands in UNIX, but I
obtained correct result when I run these command in WINDOWS. Can
somebody point out the problem and give the solution. Thanks.
dt-read.table(file=Fall.dat)
dim(dt)
[1] 1941
Anupam Tyagi AnupTyagi at yahoo.com writes:
Is there something in R that will display both observed values and their
influence on calculated statistics?
In general : yes. You may give the example in the lm documentation a first
try, and then rephrase your question.
Dieter
files to mark those parts that should only be included in the full
documentation.
(For fear of being rip...ed into pieces, going off for vacation now).
Dieter Menne
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
Right. the grid call shouldn't be necessary, axis.line controls the
panel borders. And tck can be a vector, to get rid of the ugly bumps
on top:
xyplot(x ~ x | g, type = l, lwd = lwd,
scales = list(tck = c(-1, 0)),
carlos riveira carlos.riveira at yahoo.com writes:
I am working with Probit regression (I cannot switch to logit) can anybody
help me in finding out how to obtain with R Finney's fiducial confidence
intervals for the levels of the predictor (Dose) needed to
produce a proportion of 50% of
nls not converging for zero-noise cases
Setzer.Woodrow at epamail.epa.gov writes:
No doubt Doug Bates would gladly accept patches ... .
The zero-noise case is irrlevant in practice, but quite often I have uttered
/(!! (vituperation filter on) when nls did not converge with real data. The
Yours truly dieter.menne at menne-biomed.de writes:
...
Recently, a colleague fitted gastric emptying
curves using GraphPad, with 100% success, and
nls failed for one third of these. When we
checked GraphPads output more closely, some of
the coefficients looked like 2.1 with a confidence
Earl F. Glynn efg at stowers-institute.org writes:
It's not clear to me why this problem cannot be fixed somehow. You
You might try optim instead of nls, which always (well, as far I used it)
converges. However, resulting coefficients may be totally off, and you should
use profiling to check
Andrew Robinson A.Robinson at ms.unimelb.edu.au writes:
I have up to 12 measures of a protein for each of 6 patients, taken
every two or three days. The pattern of the protein looks periodic,
but the height of the peaks is highly variable. I'm testing for
periodicity using a Monte Carlo
Ajay Narottam Shah ajayshah at mayin.org writes:
Consider this situation:
x1 - runif(100); x2 - runif(100); y - 2 + 3*x1 - 4*x2 + rnorm(100)
m1 - summary(lm(y ~ x1))
m2 - summary(lm(y ~ x2))
m3 - summary(lm(y ~ x1 + x2))
Now you have estimated 3 different competing models, and
(maxiter=15, tol=0.01), # nothing makes
sense here
+trace=TRUE)
This toy problem is exactly what the warning is for:
Warning
Do not use nls on artificial zero-residual data.
Add some noise and try again.
Dieter Menne
__
R-help
Dirk Enzmann dirk.enzmann at uni-hamburg.de writes:
I receive R-help messages in digest form that makes it difficult to
answer a post. I noticed that my answer is not added to the thread
(instead, a new thread is started) although I use the same subject line
(starting with Re: ) as the
John Kane jrkrideau at yahoo.ca writes:
I have a dataframe where I would like to order first
by variable, year, and then within that variable by
month.
http://www.ats.ucla.edu/STAT/r/faq/sort.htm
and
http://wiki.r-project.org/rwiki/doku.php?id=misc:rtips-status
Dieter
Joseph LeBouton lebouton at msu.edu writes:
Can anyone help me understand why an lm model summary would return an
r.squared of ~0.18 with an intercept term, and an r.squared of ~0.98
without the intercept? The fit is NOT that much better, according to
plot.lm: residuals are similar
Xin jasonshi510 at hotmail.com writes:
I try to list all of combination among (a,b,c,d,e,f,g,h,i,j,k)
Check function combinations in gtools.
Dieter
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Simone Vincenzi simone.vincenzi at nemo.unipr.it writes:
Dear R-list,
I have a statistical problem with the comparison of two short time-series of
density data in an ecological framework. I have to compare two short time
series (5 years, one value for each year) of species density data
1 - 100 of 287 matches
Mail list logo