Re: [R] Newbie Question: Using R with PHP and MySQL

MASFERFC Team 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 the d

[R] Lost in substitute: nlsList and nlme with dynamic formula

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 = nlsList(conc~SSbiexp(time,A1,lrc

Re: [R] for plots

jim holtman 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 li

Re: [R] R graphics display window

Felipe Carrillo 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 = 12) Dieter ___

Re: [R] lmer coefficient distributions and p values

Daniel Lakeland 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). The res

Re: [R] to combine bwplot + srt option?

KOITA Lassana - STAC/ACE 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 cannot, since bwplo

Re: [R] Import of Access data via RODBC changes column name ("NO" to "Expr1014") and the content of t he column

Maciej Hoffman-Wecker 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 <- odbcConnectAccess("./te

Re: [R] Using sunflowerplot to add points in a xyplot panel

Ronaldo Reis Junior 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 func

Re: [R] Problem with

Daniel Stahl 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 __ R-help@stat

Re: [R] p value statistic for lm object

Arjun Bhandari 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 <- c(4.17,5.5

Re: [R] bayesm - question about 'rscaleUsage function'

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) > out=rscaleUsage(Data=surveydat,Mcmc=

Re: [R] multiple nls - next fit even after convergence problem

GOUACHE David 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 going th

Re: [R] Invert Likert-Scale Values

Alexis Delevett 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 linear reg

Re: [R] [R-pkgs] new package plotAndPlayGTK

Felix Andrews 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

Re: [R] Fitting exponential curve to data points

Andrew Clegg 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) comments on

Re: [R] custom point shapes

baptiste Auguié 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 a

[R] What is parameter "display" (windows/grDevices)

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

[R] Print grid/ggplot to a metafile

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

Re: [R] eps in odfWeave

Dieter Menne 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,

Re: [R] eps in odfWeave

Kuhn, Max 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... suggestions before

[R] eps in odfWeave

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" #pl

Re: [R] Setting up a blank table with column names in the hard drive

Yuchen Luo 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 the results o

Re: [R] latex of ftable (Hmisc?)

Chuck Cleland optonline.net> writes: > How about a Sweave approach? Something like this in the *.Rnw file: > > \documentclass[letter]{article} ... > > ftable(symptom ~ treatment + sex) > > > > \end{document} > > Then Sweave() the *.Rnw file to produce a *.tex file in the working > d

[R] latex of ftable (Hmisc?)

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 <- f

Re: [R] Linear Mixed Models with nlme, more than one random effect

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 if you currently want to use lme4. Dieter

[R] User Defined R Functions in Sweave/Latex

M. Jankowski 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: > basicplot

Re: [R] How to draw several ROC curves on a common graph

Tirthadeep 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 ___

Re: [R] Stacked barchart color

hadley wickham gmail.com> writes: > > On 6/12/07, Dieter Menne 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 exce

[R] Stacked barchart color

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. barcha

Re: [R] how to ignore error messages?

Jianping Jin 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 paramet

Re: [R] Standard errors of the predicted values from a lme (or lmer)

Fränzi Korner 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, say, fixed fac

Re: [R] Standard errors of the predicted values from a lme (or lmer)-object

Fränzi Korner 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 residuals(lmeresu

Re: [R] Function to Sort and test AIC for mixed model lme?

Ken Nussear 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 dropterm or

Re: [R] RColorBrewer Package Color index

W. Katy 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 f

Re: [R] Multivariate Outlier Detection - Robust

Fiona Sammut 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 R. I also

Re: [R] Need an working examples

Aurel Razvan Duica 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 really simple ex

Re: [R] ANOVA results in R conflicting with results in ot her software packages

Simon Blomberg 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

Re: [R] SweaveInput and absolute paths

Kevin R. Coombes 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 pat

Re: [R] Multiple trellis plots without intervening space

David Lindelof 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 bit agai

Re: [R] xyplot source file only work with echo=TRUE

, but xyplot dont work. Why it dont 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.general&followup=83924> See:

Re: [R] Regarding Vista

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 z

Re: [R] Wikibooks

ppe Grosjean informed me: "Your page is way too long and is a rather crude copy and paste from the long thread in the mailing list." I disagree. Why do you have a "large guides" section? And taking into account the amount of work I put into reformatting the transcript, I

Re: [R] Regarding Vista

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 installation of chm-files

Re: [R] equivalent datatypes

Duncan Murdoch 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 https://stat.ethz.ch/mailman/listin

Re: [R] tcltk crashing R after the (ab)use of tkwait

Alberto Monteiro 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 cli

Re: [R] logistic regression on contingency table

es the response can also 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 mailin

[R] LD50 contrasts with lmer/lme4

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 50

Re: [R] Google, hard disc drives and R

Tim Churches 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 wh

[R] fit.contrast (package gmodels) and lmer

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) library(lm

Re: [R] Suppresing default text in pairs.lmList() in package = nlme

Michael Kubovy 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 sources

Re: [R] Near function?

Bart Joosen 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) should

Re: [R] "lme" in R and Splus-7

yyan liu 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 in th

Re: [R] coefficients of each local polynomial from loess() or locfit()

Liu, Delong (NIH/CIT) [C] 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 chang

Re: [R] coefficients of each local polynomial from loess() or locfit()

Liu, Delong (NIH/CIT) [C] 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 posti

Re: [R] coefficients of each local polynomial from loess() or locfit()

Liu, Delong (NIH/CIT) [C] 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 o

Re: [R] Help with filled.contour()

Michael Kubovy 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 > * >

Re: [R] Help with filled.contour()

Michael Kubovy 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 smooth;

Re: [R] zero random effect sizes with binomial lmer

Daniel Ezra Johnson 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 Subject and

Re: [R] Types of objects?

Raymond Balise 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 https://stat.ethz.ch/mailman/listinfo/r-hel

Re: [R] Crosstab from sql dump

Alan Gibson 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

Re: [R] how to do muliple comparisons in linear mixed models

liu, jcheng 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 __ R-help@stat.math.ethz.ch mai

Re: [R] exporting a table to latex

Rainer M Krug 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 scientifi

Re: [R] solving non linear system of eq

Paulino Perez Rodriguez 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 https://stat.ethz.

Re: [R] Switch and integer

Knut Krueger 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)) > switch (

Re: [R] What training algorithm does nnet package use?

Wee-Jin Goh 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 said,

Re: [R] error in format.df when using dec

Rainer M Krug 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 th

Re: [R] Hmisc, latex cdec not as expected - again

Henrik Frederik Thomsen 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 upgr

Re: [R] Using a deriv function in nlme

It might 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/p

Re: [R] Rank transformation and the linear mixed model

Bruno L. Giordano 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 https://stat.ethz.ch/

Re: [R] Prediction intervals for predict.gls

Spencer Jones 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 __

Re: [R] Compare linear regressios for significant differenc es of the slopes

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 ..

Re: [R] Lmer, heteroscedasticity and permutation, need help please

while lmer 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 comme

[R] toggleDoc: clutter reduction mainly for S4 documentation files

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. Version

Re: [R] write data to pdf

Franco Mendolia 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 anyway better to thi

Re: [R] Adding locfit confidence intervals in trelis xyplot

Deepayan Sarkar gmail.com> writes: > > On 10/18/06, juan f poyatos 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 function that does. S

Re: [R] Nested source() errors

Pierce, Ken 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 seems to

Re: [R] creating bins for a plot

Jeffrey Stratford 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 - 1.0) a

Re: [R] sqlSave, fast=F option, bug?

eems 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, &q

Re: [R] if statement error

Jenny Stadt yahoo.ca> writes: > if(length(real.d)>=30 && length(real.b)>=30 && beta1*beta2*theta1*theta2>0 ) > > { 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) >

Re: [R] Transcript of "Conservative ANOVA tables"

Marc Schwartz (via MN 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 that s

Re: [R] syntax of nlme

Thomas Wutzler 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 indirection

[R] Transcript of "Conservative ANOVA tables"

ength 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

Re: [R] latex and anova.lme problem

Frank E Harrell Jr 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. http://www.menne-b

Re: [R] Estimate in Wilcox_test vs wilcox.exact

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

Re: [R] "which" command

AgusSusanto 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] 1

Re: [R] Fitting a cumulative gaussian

Gamer, Matthias 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. Well, m

Re: [R] plotting

Anupam Tyagi 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

[R] S4 documentation - Was: "exact Wilcoxon signed rank test "

econdary functions as "protected", I propose a keyword "protected" in the .rd 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 __

Re: [R] xyplot tick marks and line thickness

Deepayan Sarkar 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)), >par.settings =

Re: [R] Finney's fiducial confidence intervals of LD50

carlos riveira 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 responses(L

Re: [R] nls convergence problem

Yours truly 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 > interval i

Re: [R] nls convergence problem

nls not converging for zero-noise cases 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 dreaded "min

Re: [R] Presentation of multiple models in one table using xtable

Ajay Narottam Shah 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,

Re: [R] A model for possibly periodic data with varying amp litude [repost, much edited]

Andrew Robinson 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 simulation en

Re: [R] nls convergence problem

Earl F. Glynn 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 t

Re: [R] nls convergence problem

ls(Y ~ a*exp(b*X)+c, > +start=list(a=6,b=-0.5,c=1), > +control=nls.control(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 a

Re: [R] How to order or sort a data.frame

John Kane 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 _

Re: [R] How to reply to a thread if receiving R-help mails in digest form

Dirk Enzmann 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 original

Re: [R] failed to load gplots

Veronika Bókony 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) : > packa

Re: [R] Help with short time series

Simone Vincenzi 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 (it is the >

  1   2   3   >