[R] Help with response CIs for lme

]] [2,], hsb.int[[1]][1,]+hsb.int[[1]][3,], hsb.int[[1]][1,]+hsb.int[[1]] [2,] + hsb.int[[1]][3,] + hsb.int[[1]][4,])) If so, is there an easier way to write it? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Re: [R] show unique combination of two factor

wanted. Thanks, Aimin Michael Dewey http://www.aghmed.fsnet.co.uk __ 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

Re: [R] random effect question and glm

-1,data = p5) m000 - glm(sc ~1+aa-1,data = p5) See above thanks, Aimin Yan Michael Dewey http://www.aghmed.fsnet.co.uk __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] means over factors in mlm terms

() for R, where mlm methods are still somewhat limited. -Michael Richard M. Heiberger wrote: ## Thank you, Michael, for the example. I believe you are looking ## for the natural extension to mlm of the proj function. ## Here is the lapply workaround from which that extension ## might

[R] Summary, was Re: Confidence interval for relative risk

At 14:43 10/11/2006, Michael Dewey wrote: After considerable help from list members and some digging of my own I have prepared a summary of the findings which I have posted (see link below). Broadly there were four suggestions 1 - Wald-type intervals, 2 - transforming the odds ratio confidence

Re: [R] Newbie problem ... Forest plot

parameter what happens, and how does that differ from what you wanted to happen? Michael Dewey http://www.aghmed.fsnet.co.uk __ 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

Re: [R] Newbie problem ... Forest plot

on the forest plot. Greetings... Peter Boolton MCA adcc.. Michael Dewey http://www.aghmed.fsnet.co.uk __ 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

Re: [R] how to create this design matrix?

[EMAIL PROTECTED] wrote: Dear Michael, This looks like a multivariate simple regression -- that is, 12 response variables, one predictor. If the data are in the matrix X, then lm(X[,1:12] ~ X[,13]) should do the trick. I hope this helps, John John Fox

[R] how to create this design matrix?

column of the 1st row of data; d1, 12 means the 12th column of the 1st row of data. But now since I have N samples(rows) of Y data; how do I conceptually do a 3D design matrix? On 11/15/06, John Fox [EMAIL PROTECTED] wrote: Dear Michael, -Original Message- From: Michael [mailto:[EMAIL

[R] Creating a table

A 2 1 1 4 B 1 00 1 colsum 31 15 so that the frequencies of df$loc are shown in the table for different time intervals. Thanks in advance for any hint, Michael

[R] how to compute p-value?

Hi all, I just want to understand how R computes p-value in a simple linear regression model? The reason is that in Matlab in the function which evaluate standard errors for multivariate normal regression, it just provide estimates and standard errors, without giving out p-value, It computes

[R] how to create this design matrix?

Hi all, I have a multiple-linear regression problem. There are 13 columns of data, the whole data matrix is: n x 13, where n is the number of samples. Now I want to regress EACH of the first 12 columns onto the 13th column, with 2-parameter linear model y_i = b0 + b1 * x_i, where i goes from

Re: [R] [R-SIG-Finance] looking for functions that can test/estimate CAPM, APT, Fama's factor model, etc.

are not very complete I guess.) Thanks a lot Brian and thanks everybody... On 11/13/06, Brian G. Peterson [EMAIL PROTECTED] wrote: On Sunday 12 November 2006 22:41, Michael wrote: I am also looking for interesting statistical experiments about testing and estimating CAPM, APT, Fama

Re: [R] Confidence interval for relative risk

did. Best, -- Wolfgang Viechtbauer Department of Methodology and Statistics University of Maastricht, The Netherlands http://www.wvbauer.com/ -Original Message- From: [EMAIL PROTECTED] [mailto:r-help- [EMAIL PROTECTED] On Behalf Of Michael Dewey Sent: Friday, November 10

Re: [R] Confidence interval for relative risk

At 12:35 12/11/2006, Peter Dalgaard wrote: Michael Dewey [EMAIL PROTECTED] writes: At 15:54 10/11/2006, Viechtbauer Wolfgang (STAT) wrote: Thanks for the suggestion Wolfgang, but whatever the original authors did that is not it. Did you ever say what result they got? -p

[R] looking for functions that can test/estimate CAMPM, APT, Fama's factor model, etc.

Hi all, I am also looking for interesting statistical experiments about testing and estimating CAPM, APT, Fama models, etc. using R using financial series data... please give me some pointers... I have been searching the R archives for the past a few hours and I vaguely got to know that there are

[R] Confidence interval for relative risk

unless I can help them with what they should have done. Is there another package I should have tried? Is there some other way of doing this? Am I doing something fundamentally wrong-headed? Michael Dewey http://www.aghmed.fsnet.co.uk __ R-help

Re: [R] Repeated Measures MANOVA in R

I have an in press paper on HE plots, http://www.math.yorku.ca/SCS/Papers/heplots.pdf that describes methods to visualize the dimensionality of effects in MLMs. The implementation is in SAS, but there's a link to a rudimentary R function in the paper. -Michael A. Bolu Ajiboye wrote: Can R do

[R] Help with plot(augPred()) plot(comparePred()) in nlme

_ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019+1-434-982

Re: [R] Reformat a data frame

similar things before along those lines. Michael Jerosch-Herold Thorsten Muehge [EMAIL PROTECTED] 11/07/06 1:10 AM Hello Experts, how do I reformat a data frame in the way described below: df1: ID desc resist thick temp 1 4711 100 5 20 2 4712 101 4 21 3 4711 99 3 19 4 4712 98 7 22 TO df2: id

Re: [R] have I an actual matrix?

, 558, 134, 53, 9) names(nums) - c(2002, 2002, 2002, 2002, 2002, 2003, 2003, 2003, 2003, 2004, 2004, 2004, 2004) barplot(nums) or require(gplots) barplot2(nums, plot.grid = TRUE, las = 1) # offers many nice options _ Professor Michael Kubovy University of Virginia

Re: [R] Colour-coded Editor for R Code

Tinn-R From: [EMAIL PROTECTED] on behalf of Jon Minton Sent: Mon 06/11/2006 9:11 AM To: R-help@stat.math.ethz.ch Subject: [R] Colour-coded Editor for R Code Do any of you know any simple programming editors for R scripts which offer basic colour-coding and

Re: [R] struggling to plot subgroups

to min; returning Inf 4: no non-missing arguments to max; returning -Inf A second issue, I would also like to offset the second set of points along the x-axis so that the error bars will be visible. Any help is appreciated. Regards, Sumit Michael Dewey http://www.aghmed.fsnet.co.uk

Re: [R] Equivalent to Stata command

Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019+1-434-982-4751 Fax:+1-434

[R] CI -- fixed effects -- response -- mixed models

to you for guidance on how to proceed. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office

Re: [R] Matrix dimnames

Daniel Gatti wrote: How does one assign names to the columns of a matrix so that the columns can be accesses using the '$' operator? It seem that x$c1 below should return column 1. x = matrix(1:4,2) dimnames(x) = list(c(r1, r2), c(c1, c2)) x c1 c2 r1 1 3 r2 2 4 x$c1

Re: [R] Saving a function

Matt, patopatasfrias wrote: Hey, I'm a real novice... how do I go about saving a function that I created eg. If I created this function: Fun-function(x1,X2) {..}, how can I store it so that the next time I use R I can load it and not have to type it out again If you search the

Re: [R] no possible to load a package correctly dowloaded by the R prompt

('partitions') You could have said install.packages('untb', dependencies = TRUE), and this would have done the job for you. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102

[R] Help with unlist

Dear r-helpers, I have a list whose elements are str(durCut[[1]]) Ord.factor w/ 5 levels vLowlowmed..: 3 2 5 2 2 2 4 4 3 5 ... How do I unlist durCut into an ordered factor? str(unlist(durCut)) int [1:3024] 3 2 5 2 2 2 4 4 3 5 ... _ Professor Michael Kubovy

[R] Package development on Windows - latest guidance?

I would like to start building R packages under Windows XP. I have programming experience and a minimal but working knowlege of many Unix (-like) programming tools. The package functions (for now) will be from R source, not C or Fortran. I've installed Rtools, Perl, the MS hhc, and so on. I am

Re: [R] Package development on Windows - latest guidance?

Thanks to Gabor G., Duncan M., and Hong O. for helpful replies. I've made some progress, but have two questions. Can anyone explain *how* R CMD searches for latex? I have provided a batch file (shell script) and and alias (I use a Windows command shell with aliasing) that each provide the

Re: [R] Package development on Windows - latest guidance?

Thanks again to D. M. and G. G. for help. I have built a small package successfully, although latex is not working quite right. I might be able to get it working, but I'm unclear what the latex compiler is used for when building a package. Otherwise, I'll just install MikTeX and treat it as a

[R] glht for aov with Error() term

for repeated measurements ANOVAs? If so, how would I do that? Best regards, Michael Zehetleitner For example (senseless, but it makes the point): - snip - amod - aov(breaks ~ tension, data = warpbreaks) amod.err - aov(breaks ~ tension + Error(wool/tension), data

Re: [R] latent class model

. Karen [[alternative HTML version deleted]] Michael Dewey http://www.aghmed.fsnet.co.uk __ 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

[R] [R-pkgs] X2R: R-compatible output routines from C, ADMB, Fortran to R

. Michael H. Prager - [EMAIL PROTECTED] Andi Stephens Southeast Fisheries Science Center National Marine Fisheries Service, NOAA 101 Pivers Island Road Beaufort, North Carolina 28516 USA Jennifer L. Martin - [EMAIL PROTECTED] Northeast Fisheries Science Center National Marine Fisheries Service, NOAA

Re: [R] Error when naming rows of dataset

('latestWithNumber.txt', header=T) Error in row.names-.data.frame(`*tmp*`, value = c(1.1, 1.2, 1.3, : duplicate 'row.names' are not allowed Yongchuan Michael Dewey http://www.aghmed.fsnet.co.uk __ R-help@stat.math.ethz.ch mailing list https

Re: [R] Error when naming rows of dataset

0 2226.0825 8.7 2.8 11 12 0 2226.0825 8.7 2.9 12 13 0 2226.0825 8.7 2.1 13 14 0 2226.0825 8.7 From: Michael Dewey [EMAIL PROTECTED] Date: Wed 25/10/2006 6:38 PM SGT To: yongchuan [EMAIL PROTECTED

[R] plot(all.effects)

= coln, nrow = 2, ncol = 2, more = last[i]) } Error in plot.effect(x[[(i - 1) * cols + j]], row = i, col = j, nrow = rows, : formal argument row matched by multiple actual arguments _ Professor Michael Kubovy University of Virginia Department of Psychology

[R] setting options in a list

Hi, I have a question regarding setting options that use a list. There are none that I'm aware of in R-base, however, it can make sense for some custom situations. For example: options(myoptions=list(lwd=1,density=NULL,angle=45,col=green)) Getting values is simple: getOption(myoptions)$col

[R] Constructing predictions from HPDinterval() after lmer()

on 5 and 129 DF, p-value: 2.2e-16 rd.eff - all.effects(rd) rd.ci - data.frame(y = rd.eff[[1]]$fit, Lower = rd.eff[[1]] $lower, Upper = rd.eff[[1]]$upper) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Re: [R] plotting 95% confidence bands on a simple linear regression model from lm()

) ?summary.effect _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729

[R] Finding out about objects and classes

When R help simply states something like: Value: An object of class 'loess'. How do I find out more about that class? Shouldn't there be a link in the help file or something? ATB Mick __ R-help@stat.math.ethz.ch mailing list

Re: [R] glm cannot find valid starting values

Mauricéia | `. `'` CP: 126, CEP: 39401-089, Montes Claros - MG - Brasil | `- Fone: (38) 3229-8190 | [EMAIL PROTECTED] | ICQ#: 5692561 | LinuxUser#: 205366 Michael Dewey http://www.aghmed.fsnet.co.uk __ R-help@stat.math.ethz.ch mailing list https

[R] Hide line ends behind unfilled circles?

(2, 3, 3, 4, 6, 4, 5, 5, 6, 7, 7, 7, 8, 8)] y1 - yy[c(2, 3, 3, 4, 6, 4, 5, 5, 6, 7, 7, 7, 8, 8)] plot(yy ~ xx, pch = aa, cex = 3) segments(x0, y0, x1, y1) Can anyone suggest a way of insuring that the lines are hidden behind the unfilled circles? _ Professor Michael

[R] Compact presentation of multiple figures

), mar = c(1.1, 2.1, 1.1, 0.1)) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office

Re: [R] Compact presentation of multiple figures

=y4) plot(y ~ x, main=5, xlab=, ylab=, xaxt=s) plot(y ~ x, main=6, xlab=, ylab=) _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

Re: [R] glm cannot find valid starting values

valid starting values: please specify some Have you tried specifying some starting values as it asks? Try start=c(-3,1) for instance Michael Dewey http://www.aghmed.fsnet.co.uk __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] Log-scale in histogramm

Log y-axis on histograms are conceptually wrong, but aren't a bad idea either. It is conceptually safer to show this using density. Consider an exponential distribution, which could look better with a log y-axis: x - rexp(1,.1) xd - density(x,from=0) par(mfrow=c(2,1)) plot(xd)

[R] Object attributes in R

Hi, I have questions about object attributes, and how they are handled when subsetted. My examples will use: tm - (1:10)/10 ds - (1:10)^2 attr(tm,units) - sec attr(ds,units) - cm dat - data.frame(tm=tm,ds=ds) attr(dat,id) - test1 When a primitive class object (numeric, character, etc.) is

Re: [R] Incompatability of lme4 and Matrix packages

this newer version. Mark Lyman Hi Mark, I had the same problem (on a Mac), which I solved by install.packages('Matrix') install.packages('lme4') On Oct 9, 2006, at 11:29 AM, Mark Lyman wrote: _ Professor Michael Kubovy University of Virginia Department

[R] Tinn-R functionality for Mac?

Dear readers: I like very much how one can hightlight R-code sections in Tinn-R and send them to R for execution. Is there an editor for the Mac that similarly allows you to highlight code and send it to the R-console? Thank you in advance for recommendations! Michael Jerosch-Herold

Re: [R] Installing Lindsey's packages

Hi Ken, My message was not well-focused. Here's the real problem: mk% R CMD INSTALL rmutil.tgz ERROR: cannot extract package from 'rmutil.tgz' On Oct 8, 2006, at 5:23 AM, Ken Knoblauch wrote: Hi Michael, What works for me is to retrieve the tgz from the garbage and to compile it from

Re: [R] Simulate p-value in lme4

Unfortunately the latter fails (even if I replace grps with n.grps): Error in slot(model2, frame) : object model2 not found In any event, I would be eager to hear more discussion of the pros and cons of these approaches. _ Professor Michael Kubovy University of Virginia

[R] latex and anova.lme problem

_ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab:B019+1-434-982-4751

[R] 'weaver' package problem

2.3-8 1.3-2 1.0-10 0.65-3 nlme lme4 Matrixlattice MASS JGR iplots JavaGD rJava 3.1-77 0.9975-1 0.9975-2 0.14-9 7.2-29 1.4-11 1.0-40.3-5 0.4-10 _ Professor Michael Kubovy University

[R] Installing Lindsey's packages

, locale: C)? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982

[R] glm and plot.effects

0.2.20.1-1 2.29 0.9975-0 0.9975-1 1.4-11 iplots JavaGD rJava MASSlattice 1.0-40.3-5 0.4-10 7.2-29 0.14-9 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville

[R] Arrowheads at line-ends

Hi, I'm quite new to R and I'm trying to do a simple plot with two lines, one upper and one lower. The problem is that the lines are reaching limiting values at x=0 for the upper line and x=1 for the lower. Thus, they are undefined at that values. I would like to indicate that with

Re: [R] error from pmvnorm

? Thanks a lot. Yonghai Li Michael Dewey http://www.aghmed.fsnet.co.uk __ 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

Re: [R] R 2.3.1 and SPSS 14.0

read.spss should work as it has. I haven't had any problems importing data files created in SPSS 14.0. Problems can arise when importing data files created by the SPSS Data Entry product (currently at version 4.0), and this is covered in the R Data Import/Export manual. Michael

[R] Documentation patch for 'match' and 'palette'

Here is a patch to improve documentation for finding useful, yet newish, functions: 'findInterval' and 'colorRamp'. I think that it is worthwhile to mention these in the 'seealso' section of the similar 'match' and 'palette' documents. I had difficulty finding these functions at first, as they

[R] Loading rgdal library

I am trying to load the rgdal library in Mac OS X 10.4.7 (pismo g3 500 mhz). I already loaded sp and compiled rgdal successfully, such that rgdal shows up in the package manager GUI. But when I try to load rgdal, I get the following error message: Error in dyn.load(x, as.logical(local),

Re: [R] Best use of LaTeX listings package for pretty printing R code

- 4 \end{lstlisting} That was \lstinline|x - 22| \lil{q - 'cat'}. \end{document} _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

[R] Fatal error: unable to restore saved data in .RData --- no package called 'nlme'

I am using Windows XP and R 2.3.1. During my lastest session I updated my packages and now when I try to start R 2.3.1 I get the following error message: Fatal error: unable to restore saved data in .RData in an error window and Error in loadNamespace(name): there is no package called 'nlme' in

Re: [R] Variable im Data Frame Namen

to answer, but I find that when I am tempted to do this it is almost always better to make a list (in your case of data frames) and then access the elements of the list Example: week - 28 test(week) - data.frame(a,b,s,c); test28 Michael Dewey http://www.aghmed.fsnet.co.uk

[R] [R-pkgs] cairoDevice 1.2 on CRAN

, Michael Lawrence ___ R-packages mailing list R-packages@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-packages __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Call for Abstracts: Tools for Intelligent Data Analysis (GfKl 07)

will be invited to submit a paper for a book published by Springer (Studies in Classification, Data Analysis, and Knowledge Organization). Deadline for submission of abstracts: Fr. November 10, 2006 Session Chairs: Michael Hahsler [EMAIL PROTECTED] Kurt Hornik [EMAIL PROTECTED] For information

[R] justRMA

I am trying to preprocess a large dataset of affymetrix data. Creating an affybatch is not possible with the computer I am running it on, so I have used the justRMA command to run RMA. I have read the affy document describing the justRMA command and the help documentation but I am unclear as

Re: [R] prediction interval for new value

incentive to develop some facility with the concept and code by working a few out first. HTH Regards, Michael Grant --- Peter Dalgaard [EMAIL PROTECTED] wrote: Sachin J [EMAIL PROTECTED] writes: Google search gave me this: http://ewr.cee.vt.edu/environmental/teach/smprimer

[R] Access Rows in a Data Frame by Row Name

I have created a data frame using the read.table command. I want to be able to access the rows by the row name, or a vector of row names. I know that you can access columns by using the data.frame.name$col.name. Is there a way to access row names in a similar manner? [[alternative

Re: [R] Reading fixed column format

read every dataset. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434

Re: [R] Command equivalent of rgui File, Save to File?

Duncan Murdoch wrote: On 9/11/2006 3:57 PM, Michael Prager wrote: R 2.3.1 on Windows XP Professional. I am writing some scripts to generate examples. The Rgui menu item File, Save to File is helpful. Is there perhaps an equivalent R function that can be incorporated into a script

[R] Command equivalent of rgui File, Save to File?

R 2.3.1 on Windows XP Professional. I am writing some scripts to generate examples. The Rgui menu item File, Save to File is helpful. Is there perhaps an equivalent R function that can be incorporated into a script? Mike Prager Southeast Fisheries Science Center, NOAA Beaufort, North Carolina

Re: [R] Successive Graphs

[EMAIL PROTECTED] wrote: Hello! I have written an R script on a Windows platform where I calculate eight result matrices I plot using matplot. I would like to display the resulting plots successively, rather than simultaneously, and I was wondering if anyone could point me in the right

[R] plot image matrix with row/col labels

, 35 Ouvriers, 36 Journaliers, 37 Domestiques, 38 Chevaux, 39 Chiens, 40 Moralite ) dimnames(mat) - list(rnames, cnames) # how to plot the image matrix with row/col names??? # show the image matrix plot(imagematrix(mat)) -- Michael Friendly Email: friendly AT yorku DOT ca

Re: [R] Consulta sobre bases en R

. [[alternative HTML version deleted]] Michael Dewey http://www.aghmed.fsnet.co.uk __ 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

Re: [R] fitting a gaussian to some x,y data

for my confusing and perhaps incorrect usage of terminology. Cheers, Michael __ 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

Re: [R] fitting a gaussian to some x,y data

Whoops, I forgot to add, many thanks to all who replied publicly and privately. I am very appreciative of all comments and suggestions. Mark Leeds was especially kind in terms of clarifying why my description of the situation was so confusing. Cheers, Michael On Aug 26, 2006, at 2:51 PM

[R] fitting a gaussian to some x,y data

to want the data in some other format or something because it is complaining about singular gradient. I'm sure this isn't hard and the answer must be out there somewhere but I can't find it. I appreciate any assistance. Cheers, Michael filepath - system.file(data, infile , package=datasets) mm

Re: [R] fitting a gaussian to some x,y data

it with a gaussian. Does that make sense? Thanks! Michael On Aug 25, 2006, at 2:04 PM, MARK LEEDS wrote: hi : i'm not clear on what you mean but someone else might be ? if you say ( x,y), then it sounds like you are talking about a bivariate normal distribution. to fit a regular one dimensional

Re: [R] negatively skewed data; reflecting

(Lsoc) + 1). But if you do a Google search on 'ladder powers negative skew' you'll find another answer. The idea of the ladder of powers is due to J W Tukey. _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

Re: [R] HH and Rcmdr.HH packages available

.tar and then installed it with the command sudo R CMD INSTALL Desktop/HH_1.4.tar.gz _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall

[R] Rgraphviz fails to load

RuuidJGR JavaGD rJava MASS lattice 1.10.6 1.10.01.4-40.3-40.4-5 7.2-27.1 0.13-10 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400

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

-hamburg.de/instkrim/kriminologie/Mitarbeiter/Enzmann/Enzmann.html Michael Dewey http://www.aghmed.fsnet.co.uk __ 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

Re: [R] Vector Join

At 13:23 13/08/2006, Michael Zatorsky wrote: Hi, I'm working on producing a simple cumulative frequency distribution. Thanks to the help of the good people on this list I now have four vectors that I'd like to join/relate into a table. e.g. v1 - myHistogram$breaks # classes v2

[R] Vector Join

class in a format I can write out to a text file as: v1v2v3v4 v1v2v3v4 etc... Any advice will be appreciated. Regards Michael. Coming soon: Celebrity Survivor - 11 celebrities, 25 days, unlimited

[R] jpeg() and JGR

_ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab

[R] R2HTML: request for an extended example

Dear R-helpers, If you know of an extended example of the use of R2HTML, in which the various constructs are present in one place, could you please point me to it or send it to me? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS

Re: [R] R2HTML: request for an extended example

On Aug 9, 2006, at 9:23 AM, Tobias Verbeke wrote: Michael Kubovy wrote: Dear R-helpers, If you know of an extended example of the use of R2HTML, in which the various constructs are present in one place, could you please point me to it or send it to me? There was an article in the R

[R] Capabilities and JGR

solution there. Version 2.3.1 (2006-06-01) powerpc-apple-darwin8.6.0 JGR 1.4.2 _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick Road

[R] R2HTML incomplete echo

(x)? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office:B011+1-434-982-4729 Lab

[R] Frequency Distribution

only downloaded R about 2 hours ago, apologies if this is obviously documented somewhere I missed.) Regards Michael. Send instant messages to your online friends http://au.messenger.yahoo.com __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch

Re: [R] Frequency Distribution

in mind in the package that you thought may help? All I'm looking to do is ask it to give me frequencies and cumulative frequencies for the whole dataset, using intervale widths of 100 or 1000 (in much the same way the data would have to have been binned before producing a histogram. Regards Michael

[R] Beamer and Sweave

problem or an Sweave problem? Suggestions? _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA 22904-4400 Parcels:Room 102Gilmer Hall McCormick RoadCharlottesville, VA 22903 Office

Re: [R] Beamer and Sweave

, Michael Kubovy [EMAIL PROTECTED] wrote: Dear R-helpers, Here is a minimal .Rnw file which shows that builds do not work in frames that contain chunks of verbatim code: \documentclass[]{beamer} \author{} \date{} \title{Title} \begin{document} \frame[containsverbatim]{\frametitle{Here

[R] bullseye or polar display of circular data

contexts you can define a circular attribute for your data. Are there plot routines for such circular data? thank you! Michael Jerosch-Herold __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] read.spss 'error reading system-file header'

the .sav file in question was created with the SPSS Data Entry product. If that is the case, then it is covered by section 3.1 of the R Data Import/Export document. Michael __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman

Re: [R] Great R documentation

aim and why I should use this one tell me where to go for more information but great vignettes, websites or book are also of interest. What is your favourite bit of R documentation, and why? Thanks, Hadley Michael Dewey http://www.aghmed.fsnet.co.uk

[R] Functions ,Optim, Dataframe

I think I need to clarify a little further on my original question. I have the following two rows of data: mydat-data.frame(d1=c(3,5),d2=c(6,10),p1=c(.55,.05),p2=c(.85,.35)) mydat d1 d2 p1 p2 1 3 6 0.55 0.85 2 5 10 0.05 0.35 I need to optimize the following function using optim for each row

[R] Functions ,Optim, Dataframe

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 I want to use these two steps on a dataframe:

<    1   2   3   4   5   6   7   8   9   10   >