Re: [R] install local packages

2006-03-23 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Hello all, I'm trying to install the local package under window system. Two ways I've tried: 1. using the menupackages install package(s) from local zip files My .zip file is mclust.zip. But it shows Errors which are: Error in

Re: [R] install local packages

2006-03-23 Thread Aleš Žiberna
1.: Did you try to extract the zip file manually? It seem that something is wrong with the zip file, not the procedure. Try downloading (compiling) it again. 2.: You should not use \ in your paths, you should either \\ of / instead. Best, Ales Ziberna [EMAIL PROTECTED] pravi: Hello all,

Re: [R] New to R

2006-03-23 Thread Petr Pikal
Wellcome. Hope you will enjoy it. Petr. PLEASE do read the posting guide! On 22 Mar 2006 at 15:24, [EMAIL PROTECTED] wrote: From: [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date sent: Wed, 22 Mar 2006 15:24:29 -0600 Subject:

Re: [R] sqlSave

2006-03-23 Thread Jacques VESLOT
Dear Tobias, I finally succeeded in exporting my dataframes from R into Access with this script: library(RODBC) canal - odbcConnectAccess(D:/Floristique.mdb) sqlSave(channel=canal, dat=flore, tablename=Floristique, rownames=F, safer=F, fast=F, varTypes=c(dates=Date)) odbcClose(canal)

Re: [R] NLME Covariates

2006-03-23 Thread Robin Martin
HLM question? Is there a minmum number of observations required for a category..I have individusals in work teams.I have incomplete data for all the teams ..sometimes I only have data for one person in a team.I assume that HLM can't work here! But what would be the mimimal.at the moment I

Re: [R] a question related to clustering

2006-03-23 Thread Berton Gunter
Sounds like you should go to the Bioconductor website http://www.bioconductor.org/ and get plugged into their extensive and professionally written software instead of shaking and baking your own. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA The business of the

Re: [R] install local packages

2006-03-23 Thread Peter Ehlers
Where did you get mclust.zip? CRAN shows mclust_2.1-11.zip. Why not just let R install mclust directly from the CRAN mirror at SFU? Use the menu: Packages Install package(s). As to your install.packages(), you need forward slashes, contriburl=, and you probably want a different destdir (or just

Re: [R] install local packages

2006-03-23 Thread Ulises M. Alvarez
Take a look at: R for Windows FAQ http://cran.us.r-project.org/bin/windows/base/rw-FAQ.html [EMAIL PROTECTED] wrote: Hello all, I'm trying to install the local package under window system. Two ways I've tried: 1. using the menupackages install package(s) from local zip files

Re: [R] install local packages

2006-03-23 Thread Don MacQueen
1. My best guess, and it is only a guess, is that the package wasn't built properly. That's assuming that when you used the menu, you navigated to the zip file. 2. Please review the documentation for the install.packages() function. When repos is NULL, your first argument is supposed to be a

[R] lme convergence

2006-03-23 Thread Chaudhary, Mohammad A.
Dear All: lme(sbp ~ cengirth, data = family, random= ~ 1 | familyid) converges but lme(sbp ~ cengirth, data = family, random= ~ cengirth | familyid) does not. I get the following message: Error in lme.formula(sbp ~ cengirth, data = family, random = ~cengirth | :

[R] about histograms

2006-03-23 Thread linda.s
x - rnorm(1000) hist(x) hist(x,nclass=100) The two histograms look different and check the help file but still confused... sorry if it is too simple. Linda __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] rnorm

2006-03-23 Thread linda.s
just curious, why put r in the rnorm for The Normal Distribution? Linda __ 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

[R] RMySQL's column limit

2006-03-23 Thread Mark Van De Vyver
Dear R-users, First, thank you to the developers for the very useful R-library RMySQL. While using this library a recieved an error message: RS-DBI driver: (could not run statement: Too many columns) The statement that generated the error was: dbWriteTable(dbcon, simdataseries, template,

[R] Estimation of skewness from quantiles of near-normal distribution

2006-03-23 Thread Leif Kirschenbaum
I have summary statistics from many sets (10,000's) of near-normal continuous data. From previously generated QQplots of these data I can visually see that most of them are normal with a few which are not normal. I have the raw data for a few (700) of these sets. I have applied several tests

[R] about read txt file

2006-03-23 Thread linda.s
I have txt file and use the following command to read it to R: z - read.table(c:/temp/test.txt) and I type z but the output is totally a mess; the test.txt is strutured with the first line as the variable names and when use excel to open it, all the vaules are well positioned. why in R it is

[R] new to R

2006-03-23 Thread linda.s
z - read.table(c:/temp/q.txt) Warning message: incomplete final line found by readTableHeader on 'c:/temp/q.txt' what does that mean? my q.txt is like: 1, 2, 3, 33, 44, 88, 23, 43, 69, __ R-help@stat.math.ethz.ch mailing list

[R] R 2.2.1: install packages

2006-03-23 Thread Hannu Kahra
Hi, when using a laptop and trying to install packages for R for Windows 2.2.1 I get the following error utils:::menuInstallPkgs() --- Please select a CRAN mirror for use in this session --- Warning: unable to access index for repository http://cran.dk.r-project.org/bin/windows/contrib/2.2

Re: [R] rnorm

2006-03-23 Thread Chuck Cleland
linda.s wrote: just curious, why put r in the rnorm for The Normal Distribution? From ?rnorm Value: 'dnorm' gives the density, 'pnorm' gives the distribution function, 'qnorm' gives the quantile function, and 'rnorm' generates random deviates. -- Chuck Cleland, Ph.D. NDRI,

[R] read txt?

2006-03-23 Thread linda.s
which kind of txt file can be read into R? can anyone give me an example? I used a txt file and the result is: Warning message: incomplete final line found by readTableHeader on 'c:/temp/q.txt' Linda __ R-help@stat.math.ethz.ch mailing list

Re: [R] about read txt file

2006-03-23 Thread Cuvelier Etienne
linda.s a écrit : I have txt file and use the following command to read it to R: z - read.table(c:/temp/test.txt) and I type z but the output is totally a mess; the test.txt is strutured with the first line as the variable names and when use excel to open it, all the vaules are well

Re: [R] about read txt file

2006-03-23 Thread linda.s
On 3/23/06, Cuvelier Etienne [EMAIL PROTECTED] wrote: linda.s a écrit : I have txt file and use the following command to read it to R: z - read.table(c:/temp/test.txt) and I type z but the output is totally a mess; the test.txt is strutured with the first line as the variable names

Re: [R] read txt?

2006-03-23 Thread Sumanta Basak
What kind of file are you trying to read in R? If it is comma delimited, then try, Z-read. table(c:/temp/q.txt, sep=,, header=TRUE) ##[header=TRUE if you have variable names in your file ## -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of linda.s

[R] Sampling in R

2006-03-23 Thread Noam Friedman
Hi all! I was wondering if you can help me with a little sampling issue I'm having in R. I have a database with 100 observations. I need to sample n=9 sample size, 200 times (i.e. get 200 samples of size 9). N (pop. size) =100 Each sample can't contain the same observation more than one

Re: [R] about histograms

2006-03-23 Thread Petr Pikal
Hi On 23 Mar 2006 at 0:02, linda.s wrote: Date sent: Thu, 23 Mar 2006 00:02:40 -0800 From: linda.s [EMAIL PROTECTED] To: [EMAIL PROTECTED] Copies to: R-help@stat.math.ethz.ch Subject:[R] about histograms x -

Re: [R] Sampling in R

2006-03-23 Thread Thomas Petzoldt
Does the following what you want: x - 1:100 s - matrix(0, nrow=200, ncol=9) for (i in 1:200) { s[i, ] - sample(x, 9) } m - rowMeans(s) hist(m) The default behavior of sample is without replacement. Thomas Noam Friedman wrote: Hi all! I was wondering if you can help me with a

Re: [R] Sampling in R

2006-03-23 Thread Sean Davis
On 3/23/06 5:52 AM, Noam Friedman [EMAIL PROTECTED] wrote: Hi all! I was wondering if you can help me with a little sampling issue I'm having in R. I have a database with 100 observations. I need to sample n=9 sample size, 200 times (i.e. get 200 samples of size 9). N (pop. size)

Re: [R] Sampling in R

2006-03-23 Thread pedro cañadilla
I prefer the function resample as it is given in the sample help, in this way: resample - function(x, size, ...) { if length(x = 1) {if (!missing(size) size == 0) x[FALSE] else x} else sample(x, size, ...) } Anyway, if an observation can not be twice in the sample, you have to use

[R] RGui: windows-record and command history

2006-03-23 Thread Thomas Steiner
a) How can I set the recording of all windows()-history forever to true? I want something like windows(record = TRUE) but not just for the window that opens then, but for all windows I will open ever. b) Scrolling up in RGui (windows 2000) to see past commands is nice, but: Is it possible to type

Re: [R] gray level values

2006-03-23 Thread Thomas Harte
# clear everything: rm(list=ls(all=TRUE)) detach() graphics.off() # make a test matrix: ( m- matrix((1:12)/12, 3, 4) ) library(gplots) # display the test matrix: image(m, col=colorpanel(12, darkblue, yellow, white), axes=FALSE) # here's the action part (no

[R] Best Solution: difference between 2 dates: IN MONTHS the way Mothers compute it

2006-03-23 Thread Smith, Phil
Folks: If you need to determine the difference between 2 dates in months, the best solution is described by David Reiner, below. Phil Smith CDC -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Sent: Fri 3/10/2006 3:52 PM To:

Re: [R] Sampling in R

2006-03-23 Thread Liaw, Andy
From: Sean Davis On 3/23/06 5:52 AM, Noam Friedman [EMAIL PROTECTED] wrote: Hi all! I was wondering if you can help me with a little sampling issue I'm having in R. I have a database with 100 observations. I need to sample n=9 sample size, 200 times (i.e. get 200 samples

[R] X11, fonts, R-2.0.1, R-2.2.1 and R-devel

2006-03-23 Thread Xavier Fernández i Marín
Hello, I am having some problems with the X11 display in a gentoo linux laptop with R compiled manually. (https://stat.ethz.ch/pipermail/r-help/2006-March/089701.html) Whether I can open the X11 device and use it when I am using 'ion' as a window manager, I can't open it using 'gnome', due to a

[R] Cross correlation in time series

2006-03-23 Thread Hufkens Koen
Hi list, I'm working on time series of (bio)physical data explaining (or not) the net ecosystem exchange of a system (+_ CO2 in versus CO2 out balance). I decomposed the time series of the various explaining variable according to scale (wavelet decomposition). With the coefficients I got from

[R] front- end problem while using nnet and tune.nnet

2006-03-23 Thread madhurima bhattacharjee
Dear R people, I am using tune.nnet from e1071 package to tune the parameters for nnet. I am using the following syntax: tuneclass -c(rep(1,46),rep(2,15)) tunennet -tune.nnet(x=traindata,y=tuneclass,size=c(50,75,100),decay=c(0,0.005,0.010),MaxNWts = 2) Here traindata is the training

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Duncan Murdoch
On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set the recording of all windows()-history forever to true? I want something like windows(record = TRUE) but not just for the window that opens then, but for all windows I will open ever. options(graphics.record=TRUE) will make that

[R] Create graphs from text console

2006-03-23 Thread Rainer M Krug
Hi I am using R 2.2.1 under Linux (SuSE 10) and would like to know if it is possible to create graphs, i.e. jpeg(filename=fn) try( coplot( mor$I_Morisita ~ mor$Year | mor$RunStr, show.given=FALSE) ) dev.off() from a text console? It gives me an error message on the jpeg() command:

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Martin Maechler
Duncan == Duncan Murdoch [EMAIL PROTECTED] on Thu, 23 Mar 2006 08:48:52 -0500 writes: Duncan On 3/23/2006 7:35 AM, Thomas Steiner wrote: b) Scrolling up in RGui (windows 2000) to see past commands is nice, but: Is it possible to type eg wi and the arrow up

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Philippe Grosjean
Duncan Murdoch wrote: On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set the recording of all windows()-history forever to true? I want something like windows(record = TRUE) but not just for the window that opens then, but for all windows I will open ever.

[R] Default lag.max in ACF

2006-03-23 Thread Rafal Stankiewicz
Hi, The default value for lag.max in ACF implementation is 10*log10(N) There several publications recommending setting lag.max to: - N/4 (Box and Jenkins, 1970; Chatfield, 1975; Anderson, 1976; Pankratz, 1983; Davis, 1986; etc.) - sqrt(N)+10 (Cryer, 1986) - 20=N=40 (Brockwell and Davis) Why

[R] clogit question

2006-03-23 Thread Dan Chan
Hi, I am playing with clogit(case~spontaneous+induced+strata(stratum),data=infert) from clogit help file. This line works. 1. But, why strata(stratum) doesn't have a coefficient like spontaneous and induced? 2. When I remove strata(stratum) from the command, this function seems to keep

[R] lme plot

2006-03-23 Thread Leonardo D Bacigalupe
Hi all, I have a questions regarding mixed effects models: 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? library(nlme)

Re: [R] RMySQL's column limit

2006-03-23 Thread David James
Hi, Mark Van De Vyver wrote: Dear R-users, First, thank you to the developers for the very useful R-library RMySQL. While using this library a recieved an error message: RS-DBI driver: (could not run statement: Too many columns) The statement that generated the error was:

[R] still unclear about the parameters of Moran's I

2006-03-23 Thread zhijie zhang
I have read the introductions on Moran of SPDEP package, but still unclear about the parameters of Moran's I, and can't calculate the Moran's I. For example,I have a dataset like the following(only an example): longitude latitudex 110.23 32.53 10 109.52 33.2120

[R] invalid variable type in model.frame within a function

2006-03-23 Thread Ingmar Visser
Dear expeRts, I came across the following error in using model.frame: # make a data.frame jet=data.frame(y=rnorm(10),x1=rnorm(10),x2=rnorm(10),rvar=rnorm(10)) # spec of formula mf1=y~x1+x2 # make the model.frame mf=model.frame(formula=mf1,data=jet,weights=rvar) Which gives the desired output:

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Gabor Grothendieck
On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set the recording of all windows()-history forever to true? I want something like windows(record = TRUE) but not just for the window that opens then, but for all windows I will open

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Peter Ehlers
Duncan Murdoch wrote: On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set the recording of all windows()-history forever to true? I want something like windows(record = TRUE) but not just for the window that opens then, but for all windows I will open ever.

[R] gam y-axis interpretation

2006-03-23 Thread Bliese, Paul D LTC USAMH
Sorry if this is an obvious question... I'm estimating a simple binomial generalized additive model using the gam function in the package mgcv. The model makes sense given my data, and the predicted values also make sense given what I know about the data. However, I'm having trouble

Re: [R] clogit question

2006-03-23 Thread Thomas Lumley
On Thu, 23 Mar 2006, Dan Chan wrote: Hi, I am playing with clogit(case~spontaneous+induced+strata(stratum),data=infert) from clogit help file. This line works. Yes, this is one of the nice features of R (that the examples work). 1. But, why strata(stratum) doesn't have a coefficient

Re: [R] clogit question

2006-03-23 Thread vito muggeo
Dear Dan, I think you need more (theorical) background here.. clogit() in package survival performs conditional logistic regression where you have several groups (the strata, the matched sets). There is an intercept for each stratum in the model, but you do not obtain them since estimation is

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Duncan Murdoch
On 3/23/2006 10:29 AM, Gabor Grothendieck wrote: On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set the recording of all windows()-history forever to true? I want something like windows(record = TRUE) but not just for the window

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Gabor Grothendieck
On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 10:29 AM, Gabor Grothendieck wrote: On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set the recording of all windows()-history forever to true? I want something

[R] nlme for groupedData with inner and outer factors

2006-03-23 Thread Dan Bebber
Hello, I am having trouble specifying a suitable nlme model. My data structure is described by gd - groupedData(ppath ~ lcut | exp, outer = ~ bait, inner = ~ weight, data = d) i.e. the response (ppath) of several subjects (sub) was measured at levels of a continuous variable (lcut). Subjects

Re: [R] generating multivariate autocorrelated time series

2006-03-23 Thread Paul Gilbert
If you can specify an ARMA or state-space model for the vector process, then you can do this with simulate in the dse bundle. Paul Gilbert Thomas Petzoldt wrote: Hello expeRts, for an application in hydrology I need to generate multivariate (log)normally distributed time series with given

[R] conservative robust estimation in (nonlinear) mixed models

2006-03-23 Thread dave fournier
Conservative robust estimation methods do not appear to be currently available in the standard mixed model methods for R, where by conservative robust estimation I mean methods which work almost as well as the methods based on assumptions of normality when the assumption of normality *IS*

Re: [R] R 2.2.1: install packages

2006-03-23 Thread Uwe Ligges
Hannu Kahra wrote: Hi, when using a laptop and trying to install packages for R for Windows 2.2.1 I get the following error utils:::menuInstallPkgs() --- Please select a CRAN mirror for use in this session --- Warning: unable to access index for repository

Re: [R] Create graphs from text console

2006-03-23 Thread Uwe Ligges
Rainer M Krug wrote: Hi I am using R 2.2.1 under Linux (SuSE 10) and would like to know if it is possible to create graphs, i.e. jpeg(filename=fn) try( coplot( mor$I_Morisita ~ mor$Year | mor$RunStr, show.given=FALSE) ) dev.off() from a text console? It gives me an error

[R] warning message using lmer()?

2006-03-23 Thread Mingyu Feng
Dear all, I use lmer to fit a mixed effect model.It give some warnings. What can I do about this? Here is the function and the warning message: model.growth.mcas5 - lmer(response ~ monthElapsed + (monthElapsed|studentID), +data= mcas5, family=binomial(link=logit), method='ML')

Re: [R] How to compare areas under ROC curves calculated with ROCR package

2006-03-23 Thread Frank Samuelson
The seROC routine you included is an very good approximation to the standard error of the Mann-Whitney-Wilcoxon/Area under the ROC curve statistic. It is derived from negative exponential models, but works very well in general (e.g. Hanley and McNeil, Diagnostic Radiology, 1982, v. 143, p. 29). A

[R] AIC mathematical artefact or computation problem ?

2006-03-23 Thread lionel humbert
Dear R user, I have made many logistic regression (glm function) with a second order polynomial formula on a data set containing 440 observation of 96 variables. I’ve made the plot of AIC versus the frequency (presence/observations) of each variable and I obtain a nearly perfect arch effect

[R] Intercepts in linear models.

2006-03-23 Thread Rolf Turner
A colleague asked me if there is a way to specify with a ***variable*** (say ``cflag'') whether there is an intercept in a linear model. She had in mind something like lm(y ~ x - cflag) where cflag could be 0 or 1; if it's 0 an intercept should be fitted, if it's 1 then no intercept.

Re: [R] invalid variable type in model.frame within a function

2006-03-23 Thread Gabor Grothendieck
Examine the source code of lm to determine the proper way of doing this. Following that gives: makemodelframe - function(formula,data,weights) { mf - match.call() mf[[1]] - as.name(model.frame) eval(mf, parent.frame()) } On 3/23/06, Ingmar Visser [EMAIL PROTECTED]

Re: [R] DIfference between weights options in lm GLm and gls.

2006-03-23 Thread Spencer Graves
In my tests, gls did NOT give the same answers as lm and glm, and I don't know why; perhaps someone else will enlighten us both. I got the same answers from lm and glm. Since you report different results, please supply a replicatable example. I tried the following:

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Martin Maechler
Martin == Martin Maechler [EMAIL PROTECTED] on Thu, 23 Mar 2006 15:04:23 +0100 writes: Duncan == Duncan Murdoch [EMAIL PROTECTED] on Thu, 23 Mar 2006 08:48:52 -0500 writes: Duncan On 3/23/2006 7:35 AM, Thomas Steiner wrote: Martin b) Scrolling up in RGui

Re: [R] Time-Series, multiple measurements, ANOVA model over time points, analysis advice

2006-03-23 Thread Spencer Graves
If this were my problem, I might start by considering each stimulus-response pair as a one observation, and I'd break the MEG into separate time series, each starting roughly 1 second before the stimulus and ending roughly 1 second after. If you've averaged many of these, I'm

Re: [R] invalid variable type in model.frame within a function

2006-03-23 Thread Thomas Lumley
On Thu, 23 Mar 2006, Ingmar Visser wrote: Dear expeRts, I came across the following error in using model.frame: # make a data.frame jet=data.frame(y=rnorm(10),x1=rnorm(10),x2=rnorm(10),rvar=rnorm(10)) # spec of formula mf1=y~x1+x2 # make the model.frame

Re: [R] Create graphs from text console

2006-03-23 Thread McGehee, Robert
Also check out the GDD package created by Simon Urbanek if bitmap does not fit your needs. On some systems, bitmap is slow or produces an inferior quality plot, in part because anti-aliasing is not supported (at least on my system). GDD, however, produces excellent anti-aliased graphs using the GD

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Duncan Murdoch
On 3/23/2006 10:34 AM, Peter Ehlers wrote: Duncan Murdoch wrote: On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set the recording of all windows()-history forever to true? I want something like windows(record = TRUE) but not just for the window that opens then, but for all windows I

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Duncan Murdoch
On 3/23/2006 10:46 AM, Gabor Grothendieck wrote: On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 10:29 AM, Gabor Grothendieck wrote: On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set the recording of all

Re: [R] Create graphs from text console

2006-03-23 Thread Liaw, Andy
From: Uwe Ligges Rainer M Krug wrote: Hi I am using R 2.2.1 under Linux (SuSE 10) and would like to know if it is possible to create graphs, i.e. jpeg(filename=fn) try( coplot( mor$I_Morisita ~ mor$Year | mor$RunStr, show.given=FALSE) ) dev.off() from a text

Re: [R] Estimation of skewness from quantiles of near-normal distribution

2006-03-23 Thread kumar zaman
This pertains to the first paragraph, you can use Dagostino test which is an omnibus test combining both skewness and kurtosis and has a high power, istead of only skewness of the data. Try ?dagoTest Ahmed Leif Kirschenbaum [EMAIL PROTECTED] wrote: I have summary statistics from many

[R] R code to generate Sweave input ? Sweave squared ?

2006-03-23 Thread Dirk Eddelbuettel
For some recurrent tasks, I would like to programmatically generate input for Sweave. While I could do that with many languages, in particular some starting with the letter P, I wouldn't mind advancing two more positions in the alphabet and use R itself to generate input for Sweave. In other

Re: [R] Intercepts in linear models.

2006-03-23 Thread Gabor Grothendieck
You can use update to modify a formula: fo - y ~ x fo - update(fo, . ~ . -1) willl remove the intercept. See ?update.formula On 3/23/06, Rolf Turner [EMAIL PROTECTED] wrote: A colleague asked me if there is a way to specify with a ***variable*** (say ``cflag'') whether there is an

[R] Writing a function to fit ALSOS models. problem with normalization?

2006-03-23 Thread David Hugh-Jones
Dear all, Below is my attempt at a function to fit Alternate Least Squares Optimal Scaling models, as described in Young (1981) Quantitative Analysis of Qualitative Data and Jacoby (1999) Levels of Measurement and Political Research: An Optimistic View. I would welcome any comments on coding

Re: [R] DIfference between weights options in lm GLm and gls.

2006-03-23 Thread Sundar Dorai-Raj
Hi, Spencer, For your call to gls you actually want: fit.gls.w - gls(y~x, data=DF, weights=varFixed(~1/w)) HTH, --sundar Spencer Graves wrote: In my tests, gls did NOT give the same answers as lm and glm, and I don't know why; perhaps someone else will enlighten us both. I got

[R] Accessing specific values of linear regression model

2006-03-23 Thread Kelly Vincent
I am just starting to learn R, and I'm struggling with a lot of the basic stuff. The online documentation is generally good, but I have been unable to find the answer for my problem this time. I am running linear regression on a particular set of data and plotting it. I've got the regression

Re: [R] DIfference between weights options in lm GLm and gls.

2006-03-23 Thread Spencer Graves
Hi, Sundar: Thanks, Sundar. That should have been obvious to me. However, I hadn't used varFixed before, and evidently I thought about it for only 1 ms instead of the required 2. With that change, I get the same answers for all three. Best Wishes, spencer Sundar

Re: [R] DIfference between weights options in lm GLm and gls.

2006-03-23 Thread Peter Dalgaard
Spencer Graves [EMAIL PROTECTED] writes: In my tests, gls did NOT give the same answers as lm and glm, and I don't know why; perhaps someone else will enlighten us both. The weights argument in gls (gnlslmenlme) specifies the variance, not the actual weight which is the reciprocal.

Re: [R] RMySQL's column limit

2006-03-23 Thread Mark Van De Vyver
My apologies, I'm using R v2.2.1 via the JGR Editor GUI, RMySQL v0.5.7 under windows XP with latest patches. MySQL is server version: 4.1.12 I'm pretty sure it's running on a linux box. The dimension of template is 2000, I know the error kicks in at 3000, but haven't iterated to find the exact

Re: [R] RGui: windows-record and command history

2006-03-23 Thread P Ehlers
Duncan Murdoch wrote: On 3/23/2006 10:34 AM, Peter Ehlers wrote: Duncan Murdoch wrote: On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set the recording of all windows()-history forever to true? I want something like windows(record = TRUE) but not just for the window that

[R] comparative density estimates

2006-03-23 Thread Michael Friendly
I have two series of events over time and I want to construct a graph of the relative frequency/density of these events that allows their distributions to be sensibly compared. The events are the milestones items in my project on milestones in the history of data visualization [1], and I want to

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Paul Murrell
Hi Duncan Murdoch wrote: On 3/23/2006 10:46 AM, Gabor Grothendieck wrote: On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 10:29 AM, Gabor Grothendieck wrote: On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 7:35 AM, Thomas Steiner wrote: a) How can I set

Re: [R] Intercepts in linear models.

2006-03-23 Thread Liaw, Andy
update() simply calls the model fitting function again using the new arguments, which means one redundant model fit. What I would do is construct the formula beforehand, e.g., myformula - if (cflag) y ~ x else y ~ x - 1 or something like that. Andy From: Gabor Grothendieck You can use

Re: [R] Accessing specific values of linear regression model

2006-03-23 Thread Liaw, Andy
You said you were fitting regression model, yet you used summary(aov(model))? It's rather unusual to use aov() for regression models... If you used lm() to fit the regression model, the following might help: y - rnorm(10); x - runif(10) fit - lm(y ~ x) anova(fit) Analysis of Variance Table

Re: [R] Intercepts in linear models.

2006-03-23 Thread Kjetil Brinchmann Halvorsen
Rolf Turner wrote: A colleague asked me if there is a way to specify with a ***variable*** (say ``cflag'') whether there is an intercept in a linear model. She had in mind something like lm(y ~ x - cflag) something like lm( if (cflag) y ~ x-0 else y ~ x, ... Kjetil where cflag

[R] kmeans Clustering

2006-03-23 Thread Mark Hempelmann
Dear WizaRds, My goal is to program the VS-KM algorithm by Brusco and Cradit 01 and I have come to a complete stop in my efforts. Maybe anybody is willing to follow my thoughts and offer some help. In a first step, I want to use a single variable for the partitioning process.

Re: [R] conservative robust estimation in (nonlinear) mixed models

2006-03-23 Thread Spencer Graves
I know of two fairly common models for robust methods. One is the contaminated normal that you mentioned. The other is Student's t. A normal plot of the data or of residuals will often indicate whether the assumption of normality is plausible or not; when the plot indicates

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Michael H. Prager
on 3/23/2006 2:42 PM Paul Murrell said the following: Hi Duncan Murdoch wrote: On 3/23/2006 10:46 AM, Gabor Grothendieck wrote: On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 10:29 AM, Gabor Grothendieck wrote: On 3/23/06, Duncan Murdoch

[R] PCA, Source analysis and Unmixing, environmental forensics

2006-03-23 Thread Mike Bock
I am using R for environmental forensics (determination of the sources and/or groupings in mixtures of organic chemicals in the field). The goal is to determine in there are groups of samples with similar/dissimilar compositions, and to assign samples to a potential source or a mixture of sources

Re: [R] comparative density estimates

2006-03-23 Thread Achim Zeileis
Michael, very nice and interesting plots! One alternative idea to compare the proportion of milestone items (that does not really answer the bandwith question) in Europe and North America might be a conditional density plot. After running your R source code, you could do: where -

[R] outer() function

2006-03-23 Thread Kyle LaMalfa
Greetings R-help community, I am relatively new to R, which may be why I am having trouble understanding this problem. I am trying to use outer() to generate a graphable surface of a function. If there is a better way to do this, I would appreciate the insight. Otherwise, could someone suggest

Re: [R] PCA, Source analysis and Unmixing, environmental forensics

2006-03-23 Thread Gabor Grothendieck
Review all the task views at: http://cran.r-project.org/src/contrib/Views/ On 3/23/06, Mike Bock [EMAIL PROTECTED] wrote: I am using R for environmental forensics (determination of the sources and/or groupings in mixtures of organic chemicals in the field). The goal is to determine in

Re: [R] RGui: windows-record and command history

2006-03-23 Thread Paul Murrell
Hi Michael H. Prager wrote: on 3/23/2006 2:42 PM Paul Murrell said the following: Hi Duncan Murdoch wrote: On 3/23/2006 10:46 AM, Gabor Grothendieck wrote: On 3/23/06, Duncan Murdoch [EMAIL PROTECTED] wrote: On 3/23/2006 10:29 AM, Gabor Grothendieck wrote:

Re: [R] conservative robust estimation in (nonlinear) mixed models

2006-03-23 Thread Berton Gunter
Ok, since Spencer has dived in,I'll go public (I made some prior private remarks to David because I didn't think they were worth wasting the list's bandwidth on. Heck, they may still not be...) My question: isn't the difficult issue which levels of the (co)variance hierarchy get longer tailed

Re: [R] outer() function

2006-03-23 Thread Thomas Lumley
On Thu, 23 Mar 2006, Kyle LaMalfa wrote: Greetings R-help community, I am relatively new to R, which may be why I am having trouble understanding this problem. I am trying to use outer() to generate a graphable surface of a function. If there is a better way to do this, I would appreciate

Re: [R] RMySQL's column limit

2006-03-23 Thread Jeffrey Horner
Mark Van De Vyver wrote: My apologies, I'm using R v2.2.1 via the JGR Editor GUI, RMySQL v0.5.7 under windows XP with latest patches. MySQL is server version: 4.1.12 I'm pretty sure it's running on a linux box. It turns out that this may be a MySQL limit:

Re: [R] outer() function

2006-03-23 Thread Gabor Grothendieck
This is a FAQ: http://cran.r-project.org/doc/manuals/R-FAQ.html#Why-does-outer_0028_0029-behave-strangely-with-my-function_003f On 3/23/06, Kyle LaMalfa [EMAIL PROTECTED] wrote: Greetings R-help community, I am relatively new to R, which may be why I am having trouble understanding this

Re: [R] new to R

2006-03-23 Thread Francisco J. Zagmutt
Hi Linda Did you already get a reply to your question? If not, try adding a new line at the end of the text (just hit enter after 69,the last number in your data and save the file). You also want to use the argument sep in read.table Since you have a comma at the end of each row you can

Re: [R] new to R

2006-03-23 Thread linda.s
On 3/23/06, Francisco J. Zagmutt [EMAIL PROTECTED] wrote: Hi Linda Did you already get a reply to your question? If not, try adding a new line at the end of the text (just hit enter after 69,the last number in your data and save the file). You also want to use the argument sep in read.table

[R] pairs function from graphics

2006-03-23 Thread Jungyeon Yoon
Hi, I have two questions about options for panel, upper.panel, lower.panel in pairs(). 1. I see examples using panel.smooth, panel.corr, panel.hist. Are there any plots that can be used here? 2. Can we have different plots in upper diagonal? For example, in a matrix of scatterplots, I want to

Re: [R] RMySQL's column limit

2006-03-23 Thread Mark Van De Vyver
Thanks for pointing that out. I've run the create table uploaded with that bug report and it fails for the same reason using MySQL 4.1.12 on linux and 5.0.18-nt So this is not a limitation of RMySQL. Regards Mark On 3/24/06, Jeffrey Horner [EMAIL PROTECTED] wrote: Mark Van De Vyver wrote: My

  1   2   >