I've had the same problem and solved it by removing the cases with the new
levels - they need to be handled some other way, either by building a new
model or reassigning the factor level to one in the training set.
Nagu wrote:
>
> Hi,
>
> I get the following error when I try to predict the pr
Hello,
nowadays I am employed at my thesis with the link cauchit, nevertheless I
have not found theory it brings over of this one, if someone could help one
would be grateful for it very much.
You forgive for my English, I am from Chile.
Regards
[[alternative HTML version deleted]]
___
Wouldn't you agree that maybe the word 'normal' has something to do with the
concept of 'average' as it follows from a given CLT applying to partial sums,
S_n,
of a given sequence having identical means, mu, and variances, sig^2, that for
finite n, a we have the following weak approximation:
S_n/
>From what you've provided its nearly impossible to help you.
I think that I understand enough to provide a solution the problem
that I believe you are posing.
First some comments:
(1) You can't expect that passing a character string to the argument of
tapply that expects a function to work
I was wondering if someone might be able to tell me what formula R's
influence.measures function uses for determining whether the hat value
it computes is influential (i.e., the true/false value in the "hat"
column of the returned is.inf data frame). The reason I'm asking is
that its results disag
I have seen this sort of appearance on screen when plotting very
dense grids. When I use a different device, such as pdf(), no such
artifacts appear. This supports Uwe's suggestion.
-Don
At 7:38 PM +0100 3/8/08, Uwe Ligges wrote:
>Paul Winward wrote:
>> Hi,
>>
>> I have a plot with about 2 m
?lattice
take a look at the lattice package
On Sat, Mar 8, 2008 at 6:44 PM, questions? <[EMAIL PROTECTED]> wrote:
> I am using par(mfrow(4,2)...) to plot multiple graphs in one frame.
> However, the plots naturally partitioned into two sets, the first two
> rows and the 3rd and 4th row are sligh
At 03:54 PM 3/8/2008, David Winsemius wrote:
>"Robert A. LaBudde" <[EMAIL PROTECTED]> wrote in
>news:[EMAIL PROTECTED]:
>
> > Given a set of data with a number of variables plus a response, I'd
> > like to obtain a randomized subset of the rows such that the
> > marginal proportions of each variabl
I am using par(mfrow(4,2)...) to plot multiple graphs in one frame.
However, the plots naturally partitioned into two sets, the first two
rows and the 3rd and 4th row are slightly different,
I want them to be in one frame, but I also want to draw a line (or
something to make top two rows/bottom two
I am using par(mfrow(4,2)...) to plot multiple graphs in one frame.
However, the plots naturally partitioned into two sets, the first two
rows and the 3rd and 4th row are slightly different,
I want them to be in one frame, but I also want to draw a line (or
something to make top two rows/bottom two
Hi There,
I'm working on some cluster analyses on a large data-set using hclust with
Wards method and Manhattan (city block) distance measures. I've created
dendrograms to illustrate the clustering criteria, but would like to create
a plot to examine for the classic elbow criterion to use in det
That worked perfectly. Thank you!!!
Original message
>Date: Fri, 7 Mar 2008 18:13:24 -0300
>From: "Henrique Dallazuanna" <[EMAIL PROTECTED]>
>Subject: Re: [R] Combine two columns
>To: "[EMAIL PROTECTED]" <[EMAIL PROTECTED]>
>Cc: r-help@r-project.org
>
>Try:
>
>transform(x, AL1Al2 =
Martin Kaffanke wrote:
> Am Samstag, den 08.03.2008, 19:34 +0100 schrieb Uwe Ligges:
>> Martin Kaffanke wrote:
>>> Hi there!
>>>
>>> I have two little different data. One is a computer test on people, the
>>> other is a paper and pencil test. two boxplots show me that the data is
>>> almost the
"Robert A. LaBudde" <[EMAIL PROTECTED]> wrote in
news:[EMAIL PROTECTED]:
> Given a set of data with a number of variables plus a response, I'd
> like to obtain a randomized subset of the rows such that the
> marginal proportions of each variable are maintained closely in the
> subset to that of
Svyatkovskiy Alexey <[EMAIL PROTECTED]> wrote in
news:[EMAIL PROTECTED]:
> Hello,
>
> Tell me please, can I get a sourcse code for 'lda' from MASS?
> Is MASS package still present in the latest R-2.6.2 version (I
> couldn't find it)?
>
>
"lda" is in the VR (super-)package.
--
David Winsemi
?table
e.g., table(your.matrix[,1])
On Sat, Mar 8, 2008 at 3:15 PM, Donna Tucker <[EMAIL PROTECTED]> wrote:
>
> Hello,
> I would like to know how to count the number (cardinality) of a specific
> element in a single row of a matrix. At this time I have 30X3 matrix. The
> first column is the
Hello,
I would like to know how to count the number (cardinality) of a specific
element in a single row of a matrix. At this time I have 30X3 matrix. The
first column is the treatment number for each data point. I would like to know
how many of each treatments are in this matrix. i.e. I wa
Thanks to those who have replied to my original query. However, I'm
still confused on how obtain estimates, standard error and F-tests for
main effect and interaction contrasts which agree with the SAS code
with output appended below.
for example,
## Given the dataset (from Montgomery)
twoway <-
"Ng Stanley" <[EMAIL PROTECTED]> wrote::
> Seeking suggestions to compute uncentered (pearson correlation)
> correlation efficiently.
>
> corr from stats library works on x and y columns. dist from amap
> library works on x and y rows.
>
> My data layout is slightly different such that row(i)
See ?rollapply and ?rollmean in the zoo package
and ?filter and ?embed in the core of R.
On Sat, Mar 8, 2008 at 1:38 PM, Alberto Vieira Ferreira Monteiro
<[EMAIL PROTECTED]> wrote:
> This is probably very trivial - so I can't find an answer in the help files.
>
> When I have a series x (x[1], x[2]
On 08-Mar-08 18:04:54, Svyatkovskiy Alexey wrote:
> Hello,
> Tell me please, can I get a sourcse code for 'lda' from MASS?
> Is MASS package still present in the latest R-2.6.2 version (I couldn't
> find it)?
> Thankful in advance,
> Alexey.
MASS is part of VR (for "Venables and Ripley"), a "bundl
It is in the VR bundle (collection of packages), which is 'recommended',
so it is included in all binary distributions of R; if you have the source
code you find VR in R-2.6.2/src/library/Recommended
On Sat, 8 Mar 2008, Svyatkovskiy Alexey wrote:
> Hello,
>
> Tell me please, can I get a sourcse c
Am Samstag, den 08.03.2008, 19:34 +0100 schrieb Uwe Ligges:
>
> Martin Kaffanke wrote:
> > Hi there!
> >
> > I have two little different data. One is a computer test on people, the
> > other is a paper and pencil test. two boxplots show me that the data is
> > almost the same.
> >
> > So now
This is probably very trivial - so I can't find an answer in the help files.
When I have a series x (x[1], x[2], ... x[n]) and I want to construct a
new series y (y[1], y[2], ... y[m]) such that y's are either interpolations
of the x's (when m > n) or a weighted mean (when m < n), is there
any dir
Hello,
Tell me please, can I get a sourcse code for 'lda' from MASS?
Is MASS package still present in the latest R-2.6.2 version (I couldn't find
it)?
Thankful in advance,
Alexey.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/list
Paul Winward wrote:
> Hi,
>
> I have a plot with about 2 million points. It should be a fairly dense
> plot over the reals, but I'm seeing many narrow, white vertical lines
> uniformly spread in my plot where nothing was plotted. I'm using a
> period as the plot character. I've played aro
Ingmar Visser wrote:
> Hi,
> When running R CMD check I'm getting a number of errors that I don't
> quite follow and don't know where to start looking for an answer, any
> hints appreciated.
>
> R CMD check trunk
> * checking for working latex ... OK
> * using log directory '/Users/ivisser/D
Martin Kaffanke wrote:
> Hi there!
>
> I have two little different data. One is a computer test on people, the
> other is a paper and pencil test. two boxplots show me that the data is
> almost the same.
>
> So now I'd like to know if I could handle all data as one, by testing
> with ks.test:
JStainer wrote:
> an example from my R table will calculating the average LOOCV for two
> treatments ALL and AML
>
> table
> ALL AML
> 11.2 .3
> 2.87.3
> 31.1.5
> 41.2.7
> 53.21.2
> 61.11.1
> 7.90 .99
> 81.1.32
> 92.
Please do not cross post.
Note from yesterday:
https://stat.ethz.ch/pipermail/r-help/2008-March/156586.html
On Sat, Mar 8, 2008 at 12:40 PM, John McHenry <[EMAIL PROTECTED]> wrote:
> Hi WizaRds,
>
> (I've cross-posted to r-sig-finance because a lot of people there use 'zoo'
> objects).
>
> I'm
On Sat, 8 Mar 2008, Paul Smith wrote:
> Dear All,
>
> I am trying to install the rgl package on R 2.6.1 running on Fedora 8
> (Linux), but I am experiencing the following problems:
>
> -
> configure: error: X11 not found but required, configure aborted.
>
Hi WizaRds,
(I've cross-posted to r-sig-finance because a lot of people there use 'zoo'
objects).
I'm trying to plot 2 time series lined up with par(mfrow=c(2,1)), so that the
time index is the same for the two series. The data have monthly frequency.
An example will illustrate:
# data for t
Hi there!
I have two little different data. One is a computer test on people, the
other is a paper and pencil test. two boxplots show me that the data is
almost the same.
So now I'd like to know if I could handle all data as one, by testing
with ks.test:
> ks.test(el$angststoer, fl$angsts
Given a set of data with a number of variables plus a response, I'd
like to obtain a randomized subset of the rows such that the marginal
proportions of each variable are maintained closely in the subset to
that of the dataset, and possibly maintaining as well the two-factor
interaction margina
short answer: non-visible means that they are not exported from
a package with a namespace so that they are "not found" if you just
type the name. However, you can see these functions if you use
getAnywhere(fname)
longer answer: see the section on namespaces in the writing R
extensions
On Sat, 2008-03-08 at 08:07 -0600, Douglas Bates wrote:
> On Sat, Mar 8, 2008 at 2:57 AM, Alexandra Bremner
> <[EMAIL PROTECTED]> wrote:
> > I am attempting to model data with the following variables:
>
> > timepoint - n=48, monthly over 4 years
> > hospital - n=3
> > opsn1 - no of outcomes
Dale,
You might find it fruitful to look at the help pages for fit.contrast
() and estimble() functions in the gmodels package, and the contrast
() functions in the Hmisc package.
-Greg
On Mar 7, 2008, at 4:20PM , Thompson, David ((MNR)) wrote:
> Dale,
>
> Other than the first SAS contras
Thank you all very much!
Bob
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
> project.org] On Behalf Of Muenchen, Robert A (Bob)
> Sent: Saturday, March 08, 2008 10:49 AM
> To: [EMAIL PROTECTED]
> Subject: [R] Non-visible functions are asterisked
>
> Dear R-Help
It means that the package has not exported it so normally you would
not access it directly. For example, for summary.princomp you would
normally pass a princomp object to the summary generic.
If you did want to access summary.princomp directly, since it has not been
exported, you can use the form
G'day Bob,
On Sat, 8 Mar 2008 10:48:49 -0500
"Muenchen, Robert A (Bob)" <[EMAIL PROTECTED]> wrote:
> What does the non-visible mean?
Type at the command line
> summary.aov
and hit return. Note, leave out the `()'.
Now type at the command line
> summary.princomp
That should illustrate what
On 08/03/2008 10:48 AM, Muenchen, Robert A (Bob) wrote:
> Dear R-Helpers,
>
> I suspect I'm about to ask a FAQ, but I haven't been able to find an
> answer in the FAQ, AItR or an R Site Search. When I look at the methods
> of summary (below) it says, "Non-visible functions are asterisked". I
> loo
Dear R-Helpers,
I suspect I'm about to ask a FAQ, but I haven't been able to find an
answer in the FAQ, AItR or an R Site Search. When I look at the methods
of summary (below) it says, "Non-visible functions are asterisked". I
looked at the help file for summary.princomp, which did not comment on
thank you so much. I must have spend a few hours last night trying to figure
out.
thank you again,
claire
> Date: Sat, 8 Mar 2008 09:13:31 -0500> From: [EMAIL PROTECTED]> To: [EMAIL
> PROTECTED]> Subject: Re: [R] Errors in in my code> CC: r-help@r-project.org>
> > It says the file you are t
Hi,
I must have accidentally deleted my previous post. I am having a really
difficult time calculating the LOOCV (leave out cross validation).
table in excel
genes ALL AML p.value
1 1.2 .3 .01
2 .87.3 .03
3 1.1.5 .05
4 1.2
It says the file you are trying to access it not there. Do this in
two steps so you can see what is happening:
fileName <- paste(dataDirectory, annotation_HGU133Plus2, sep = "/")
cat( fileName)
... now use in it your read
My guess is that what you think is in "annotation_HGU133Plus2" is not
the
Christophe -- You'll want to coerce golubMerge to an ExpressionSet.
> library(Biobase)
> library(golubEsets)
> data(golubMerge)
> golubUpdated <- as(golubMerge, "ExpressionSet")
Warning messages:
1: In asMethod(object) :
missing or mis-shaped 'se.exprs' in original object; creating ExpressionSet
Here is a way of doing it. You check to see which ones meet the
criteria and then delete them.
> # create indices of columns to match
> first.col <- seq(1, to=ncol(SS) - 2)
> sec.col <- seq(3, to=ncol(SS))
>
> # determine rows that met criteria
> meet.crit <- apply(SS, 1, function(.row) any((.row
On Sat, Mar 8, 2008 at 2:57 AM, Alexandra Bremner
<[EMAIL PROTECTED]> wrote:
> I am attempting to model data with the following variables:
> timepoint - n=48, monthly over 4 years
> hospital - n=3
> opsn1 - no of outcomes
> total.patients
> skillmixpc - skill mix percentage
> nurse.hours.p
On Sat, Mar 8, 2008 at 7:00 AM, Keizer_71 <[EMAIL PROTECTED]> wrote:
> Hello,
> I am to run this R script but i keep getting this error.
> > expr<-exprs(golubMerge)
> Warning message:
> The exprSet class is deprecated, use ExpressionSet instead
> I tried to find information on the website
Horn's parallel analysis method of deciding how many factors to retain has
been implemented in R in the paran package. This update to paran adds an
option to produce a graph of the analysis along the lines of what Horn
presented in his seminal 1965 paper. It also makes a few minor improvements
Hello,
I am to run this R script but i keep getting this error.
> expr<-exprs(golubMerge)
Warning message:
The exprSet class is deprecated, use ExpressionSet instead
I tried to find information on the website but no luck. (exprSet...etc)
thank you.
--
View this message in context:
http://ww
Hello,
I am to run this R script but i keep getting this error.
> expr<-exprs(golubMerge)
Warning message:
The exprSet class is deprecated, use ExpressionSet instead
I tried to find information on the website but no luck. (exprSet...etc)
thank you.
--
View this message in context:
http://ww
Murray Richardson utoronto.ca> writes:
>
> Hello,
>
> I am looking for advice on a task I am trying to complete.
>
> I have a 4 column dataframe defining the start and end coordinates of
> line edges (from a CGAL alpha shapes function to define concave hulls
> from point clusters). I would
Uwe Ligges statistik.tu-dortmund.de> writes:
>
>
> Rainer M Krug wrote:
> > Hi
> >
> > I have a raster which I would like to plot in a greyscale instead of
> > colour. Is this possible, and how?
>
> Probably, and we may tell you how given you tell us what "raster" means:
>
> > ?raster
> No
JStainer wrote:
>
> Hi,
>
> I am trying to find out the best way to calculate the average LOOCV in R
> for several classifier for, KNN, centroid classification, DLDA and SVM.
>
> I have four types of diseases and 62 samples.
>
> Is there a R code available to do this?
>
>
>
--
View thi
an example from my R table will calculating the average LOOCV for two
treatments ALL and AML
table
ALL AML
11.2 .3
2.87.3
31.1.5
41.2.7
53.21.2
61.11.1
7.90 .99
81.1.32
92.1 1.2
JStainer wrote:
>
> Hi,
>
>
an example from my R table will calculating the average LOOCV for two
treatments ALL and AML
table
ALL AML
11.2 .3
2.87.3
31.1.5
41.2.7
53.21.2
61.11.1
7.90 .99
81.1.32
92.1 1.2
JStainer wrote:
>
> Hi,
>
>
Dear All,
I am trying to install the rgl package on R 2.6.1 running on Fedora 8
(Linux), but I am experiencing the following problems:
-
configure: error: X11 not found but required, configure aborted.
ERROR: configuration failed for package 'rgl'
** Removi
Hi,
I am trying to find out the best way to calculate the average LOOCV in R for
several classifier for, KNN, centroid classification, DLDA and SVM.
I have four types of diseases and 62 samples.
Is there a R code available to do this?
--
View this message in context:
http://www.nabble.com/
On Sat, 8 Mar 2008, Bob Green wrote:
> Hello,
>
> I have just installed v 2.6.2 on a new computer running Windows XP
> and tried to perform 'update packages' via the menu option on the R console.
>
> Any advice on the following problem is much appreciated.
>
> Bob
>
>
> Below are the warning and
On 3/5/08, Martin Kaffanke <[EMAIL PROTECTED]> wrote:
> Now I'd like to have it one sided, means only the left bottom side to be
> printed (the others are the same) and I'd like to have * where the
> p-value is lower than 0.05 and ** lower than 0.01.
Look here [1], at "Visualizing Correlations"
Hi,
When running R CMD check I'm getting a number of errors that I don't
quite follow and don't know where to start looking for an answer, any
hints appreciated.
R CMD check trunk
* checking for working latex ... OK
* using log directory '/Users/ivisser/Documents/projects/
depmixProject/depmi
Hello,
I have just installed v 2.6.2 on a new computer running Windows XP
and tried to perform 'update packages' via the menu option on the R console.
Any advice on the following problem is much appreciated.
Bob
Below are the warning and error messages received. A search of the
hard drive
Ruqiang Liang gmail.com> writes:
>
> I have a table like the following. I want to fit Cm to Vm like this:
> Cm ~Cl+Q1*b1*38.67*exp(-b1*(Vm-Vp1)*0.03867)/(1+exp(-b1*(Vm-Vp1)*0.03867))^2+
Q2*b2*38.67*exp(-b2*(Vm-Vp2)*0.03867)/(1+exp(-b2*(Vm-Vp2)*0.03867))^2
>
> I use nls, with start=list(Q1=2e-3,
I am attempting to model data with the following variables:
timepoint - n=48, monthly over 4 years
hospital - n=3
opsn1 - no of outcomes
total.patients
skillmixpc - skill mix percentage
nurse.hours.per.day
Aims
To determine if skillmix affects rate (i.e. no.of.outcomes/total.patien
Hi Everyone,
I am very new to R. When i run the code yesterday, it was working fine. I
was able to find gene annotation.
somehow, today, when i try to run it again, it is giving me errors message
that i don't understand. It says that it cannot open file. what file is it
looking for?
###
66 matches
Mail list logo
