[EMAIL PROTECTED] wrote:
> Hello all,
>
> I'm trying to install the local package under window system. Two ways I've
> tried:
>
> 1. using the menu"packages" "install package(s) from local zip files...".
>
>My .zip file is "mclust.zip". But it shows Errors which are:
>
>"Err
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 al
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:
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(can
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 have
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 statis
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
j
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 menu"packages" "install package(s) from local zip files..."
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
c
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 |
:
i
> 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
PLEAS
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
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, over
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
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 tot
> 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
https://stat.ethz.ch/ma
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
War
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.
N
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
https://stat.
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
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
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 lind
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 tim
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 <-
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 w
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
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
replace=
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 typ
# 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 pa
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: Smith
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 (
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
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 the
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
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
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:
Erro
> "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
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.
>
>
> options(gra
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)
Wh
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 ru
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)
a
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:
>
> d
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
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:
>
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
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.
>
>
> options(gr
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 interp
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 coef
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
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
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"
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
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
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* satisfie
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
> http://cran.dk.
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
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'
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
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 wi
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.
Th
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]>
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:
> "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>
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 guessin
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
> mf=model.frame(formula=m
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
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
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 reco
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
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
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 w
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
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
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
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 li
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
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 (&gnls&lme&nlme) specifies the variance,
not the actual weight which is the r
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 exac
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
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
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 Steine
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 up
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 Ta
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
>
> whe
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.
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
proble
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:
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
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 <- factor(c(rep(
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
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
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
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 dist
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 appre
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:
http://bugs.mysql.com/b
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
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 eit
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 re
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
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:
>
1 - 100 of 120 matches
Mail list logo