Hello,
I ran AIC for some competing models I created. I get df and an AIC score
from the AIC procedure. Can I use the models with the lowest AIC scores from
this procedure to choose my 'best' models? If not, what else do I need to do
(and know) and how can I do it in R to chose the 'best' models?
There is a new CRAN Task View for Time Series Analysis covering 68
packages on CRAN that provide time series functionality and data.
Check it out at
http://cran.r-project.org/web/views/TimeSeries.html
Base R ships with a lot of functionality useful for time series, in
particular in the stats pac
rajesh j wrote:
Hi,
I'd like to plot vertical lines on my plot using the locator.But since
locator returns two values I end up getting 2 vertical lines for every
click.Is there some way only the x coordinate from locator can be returned?
locator(1)$x
--
Dr Richard Rowe
Zoology & Tropical
Hi,
I'd like to plot vertical lines on my plot using the locator.But since
locator returns two values I end up getting 2 vertical lines for every
click.Is there some way only the x coordinate from locator can be returned?
--
Rajesh.J
"
Thanks to the help of Gregor Gorjanc I could solve the problem:
1: With the command
sudo chmod a+x /usr/local/bin/Sweave.sh
I had only
ls -l /usr/local/bin/Sweave.sh
-rwxr-x--x 1 root root 29997 2008-10-11 15:42 /usr/local/bin/Sweave.sh
So I added sudo chmod o+r /usr/local/bin/Sweave.sh
On Mon, Oct 13, 2008 at 11:47 PM, Frank E Harrell Jr
<[EMAIL PROTECTED]> wrote:
> Gabor Grothendieck wrote:
>>
>> On Mon, Oct 13, 2008 at 11:21 PM, Frank E Harrell Jr
>> <[EMAIL PROTECTED]> wrote:
>>>
>>> [EMAIL PROTECTED] wrote:
I recall a concept of Snout: sensitivity that is high enou
Gabor Grothendieck wrote:
On Mon, Oct 13, 2008 at 11:21 PM, Frank E Harrell Jr
<[EMAIL PROTECTED]> wrote:
[EMAIL PROTECTED] wrote:
I recall a concept of Snout: sensitivity that is high enough to
essentially rule out the presence of disease. And Spin: specificity that
is high enough to essent
Assuming system() does not do what you, you will need to explain to us
what you are trying to do -- please see the posting guide and references.
Note that 'is there a package to' questions are impossible to answer, as
no one know about all R packages.
On Mon, 13 Oct 2008, jcs wrote:
Is ther
Greg Snow wrote:
Actually, somewhat counterintuitively, ppv tends to me more affected by
specificity and npv by sensitivity. You can see this from the function
SensSpec.demo in the TeachingDemos package (also see the corresponding example
on the help page for tkexamp (same package)).
I don't
On Mon, Oct 13, 2008 at 11:21 PM, Frank E Harrell Jr
<[EMAIL PROTECTED]> wrote:
> [EMAIL PROTECTED] wrote:
>>
>> I recall a concept of Snout: sensitivity that is high enough to
>> essentially rule out the presence of disease. And Spin: specificity that
>> is high enough to essentially rule in th
John Sorkin wrote:
Of course Prof Baer is correct the positive predictive value (PPV) and the
negative predictive values (NPV) serve the function of providing conditional
post-test probabilities
PPV: Post-test probability of disease given a positive test
NPV: Post-test probability of no disease
If you don't have a strong reason to force the intercept to 0, then leave it in
(default) and just ignore it in the summary. Most of the time the intercept is
a nusiance parameter, it needs to be estimated to make the other estimates
meaningful, but tests on it are not of direct interest.
[EMAIL PROTECTED] wrote:
I recall a concept of Snout: sensitivity that is high enough to essentially
rule out the presence of disease. And Spin: specificity that is high enough
to essentially rule in the presence of disease.
So perhaps the below is backwards? The higher the sensitivity, th
Greg,
I don't think the intercept is 0.
My response variable is number of individuals within habitat. My independent
variables are descritpors of the habitat, I am building models to see which
are more 'important' than others. I will build a bunch and then use AIC to
pick the 'best' model. I was a
Is there a package that allows for ssh connections or scp in R? Does anyone
know how to do this?
--
View this message in context:
http://www.nabble.com/SSH-in-R-tp19962051p19962051.html
Sent from the R help mailing list archive at Nabble.com.
__
R-h
Robert W. Baer, Ph.D. wrote:
- Original Message - From: "Frank E Harrell Jr"
<[EMAIL PROTECTED]>
To: "John Sorkin" <[EMAIL PROTECTED]>
Cc: ; <[EMAIL PROTECTED]>;
<[EMAIL PROTECTED]>
Sent: Monday, October 13, 2008 2:09 PM
Subject: Re: [R] Fw: Logistic regresion - Interpreting (SENS) a
Look at the China Sea graphic in the vignette to the gridbase
package:
http://cran.r-project.org/web/packages/gridBase/vignettes/gridBase.pdf
On Mon, Oct 13, 2008 at 5:19 PM, Richard E. Glor
<[EMAIL PROTECTED]> wrote:
> Anybody know how to rotate shapes generated with the 'points' function? I'm
Try:
> mydata2 <- with(mydata, data.frame(Var1 = rep(Var1, Count),
+ Response=rep(Response, Count) )
+ )
> mydata2
Var1 Response
1A1
2A2
3A2
4B6
5B6
6B6
7A7
8A7
9A7
Hope this helps,
___
What question are you trying to answer? That is what you need to decide first.
In the no-intercept model you can test if the slope is non-zero, but you
probably already know the answer without needing data or a test (unless most
the y values are around 0 as well (both sides) and the science sa
For data such as
mydata <- read.table(textConnection("Var1 Response Count
A 1 1
A 2 2
B 6 3
A 7 3"),header=TRUE)
closeAllConnections()
is there a command to transform this data so that there is a single
observation per row?
A 1
A 2
A 2
B 6
B 6
B 6
A 7
A 7
A 7
Thanks,
Juliet
__
Greg,
I will not be using categorical variables. That sample code I provided was
from lm help. All of my variables are real numbers. The no intercept model,
testing all means equal to zero, seems less informative then testing whether
the means are equal to each other, in the intercept included mode
The rimage package has functions for reading in and plotting jpeg files that
you could use for displaying the photograph. If you then can find 2 points in
the image (not on the same horizontal or vertical line) for which you know the
coordinates in the coordinate system that you want to plot in
One clarification:
If groups is a numeric variable, then dropping the intercept forces it to 0 and
therefore fits a line that goes through the origin, but if group is categorical
(which seems likely with a name like group), then while it forces the intercept
to 0, the dummy variable encoding st
On Mon, 13 Oct 2008, Richard E. Glor wrote:
Anybody know how to rotate shapes generated with the 'points' function? I'm
trying to place points around a radial diagram such that the y-axis of
individual shapes are oriented with the radii of the circle rather than the
y-axis of the larger plot
I don't think that you can do this with points (though I would be happy to be
proven wrong), but you could look at the my.symbols function in the
TeachingDemos package. You need to create a function to plot the symbol that
you want with the rotation that you want, but then you can use it to pla
Actually, somewhat counterintuitively, ppv tends to me more affected by
specificity and npv by sensitivity. You can see this from the function
SensSpec.demo in the TeachingDemos package (also see the corresponding example
on the help page for tkexamp (same package)).
I don't think that Frank i
thanks Hadley.
On Mon, Oct 13, 2008 at 6:32 PM, hadley wickham <[EMAIL PROTECTED]> wrote:
> Hi Stephen,
>
> It will definitely be in the next release because I just wrote the
> code last week. I'm aiming for another release early November.
>
> Hadley
>
> On Mon, Oct 13, 2008 at 2:07 PM, stephen s
Dear R users,
How to change lattice panel label/text from the automatically generated
label (based on the conditioning) to our own set of label?
for example:
someStuff <- data.frame(area = rep(c("SOUTH", "NORTH", "EAST", "WEST"), each
= 25),
group = rep(c("A","B","C","D"), each =
I recall a concept of Snout: sensitivity that is high enough to essentially
rule out the presence of disease. And Spin: specificity that is high enough
to essentially rule in the presence of disease.
So perhaps the below is backwards? The higher the sensitivity, the greater the
NPV? And th
On Mon, Oct 13, 2008 at 3:38 PM, erwann rogard <[EMAIL PROTECTED]> wrote:
> hello,
>
> i need help with:
>
> data$f1<-as.factor(data$f1)
> data$f2<-as.factor(data$f2)
> s3<-equal.count(data$s2,number=3)
> densityplot(~y| f1 + f2 + s3,
> data=mydata
> )
>
> this produces 3 plots, *successively*, o
Indeed, however as I stated in my prior Email, the cases of a 0 or 1
prevalence are degenerative and are of little practical importance. And
as noted in my EMail message, I was talking about values of PPV and NPV
as a function of sensitivity and specificity when the prevalence is
fixed.
John
John
hello,
i need help with:
data$f1<-as.factor(data$f1)
data$f2<-as.factor(data$f2)
s3<-equal.count(data$s2,number=3)
densityplot(~y| f1 + f2 + s3,
data=mydata
)
this produces 3 plots, *successively*, one for each value of s3. i was
hoping it would produce one plot. is trellis limited to 2 condi
Hi Stephen,
It will definitely be in the next release because I just wrote the
code last week. I'm aiming for another release early November.
Hadley
On Mon, Oct 13, 2008 at 2:07 PM, stephen sefick <[EMAIL PROTECTED]> wrote:
> I would like to be able to do the xyplot in ggplot below. I read in
On Mon, 13 Oct 2008, Frank E Harrell Jr wrote:
useR wrote:
Hi R helpers,
One rather statistical question?
What would be the best startegy to shortlist thousands of continous
variables automaticaly using R
as the preparation for logistic regression modleing!
Thanks
The easiest approac
John Sorkin wrote:
Of course Prof Baer is correct the positive predictive value (PPV)
and the negative predictive values (NPV) serve the function of
providing conditional post-test probabilities PPV: Post-test
probability of disease given a positive test NPV: Post-test
probability of no disease g
useR wrote:
Hi R helpers,
One rather statistical question?
What would be the best startegy to shortlist thousands of continous
variables automaticaly using R
as the preparation for logistic regression modleing!
Thanks
The easiest approach is to use a random number generator.
Frank
--
Michael Just wrote:
Great,
Thanks,
Michael
On Mon, Oct 13, 2008 at 3:56 PM, Rolf Turner <[EMAIL PROTECTED]>wrote:
On 14/10/2008, at 9:42 AM, Michael Just wrote:
What is the difference when including or not including the intercept when
using lm()?
x.noint <- lm(weight ~ group - 1))# omittin
Of course Prof Baer is correct the positive predictive value (PPV) and the
negative predictive values (NPV) serve the function of providing conditional
post-test probabilities
PPV: Post-test probability of disease given a positive test
NPV: Post-test probability of no disease given a negative tes
One thing to try is to download Sweave.bat from
http://batchfiles.googlecode.com
and place it in the same directory as the Rnw file (or anywhere on your
path) and then from Windows console:
Sweave foo.Rnw
If MiKTeX is in a standard location Sweave.bat will find it and it
will locate R itself from
Hi Marko,
this may be helpful:
http://www.ingentaconnect.com/content/bpl/rssb/2008/0070/0001/art5;jsessionid=an2la3spa0n5h.alexandra?format=print
Happy modeling!
Stephan
useR schrieb:
Hi R helpers,
One rather statistical question?
What would be the best startegy to shortlist th
Sorry, I forgot to include a reproducible example on my last e-mail but here it
is: Since the file is large to be included here:
The path to the foo.Rnw examples is:
www.stat.umn.edu/~charlie/Sweave/foo.Rnw
and is suppossed to produce a pdf like this one:
http://www.stat.umn.edu/~charlie/Sweave/f
You should note that the author of the lrm function (at least the one in the
Design package, I don't know of others) is also one of the most vocal opponents
of stepwise regression methods. Using stepwise with lrm() is kind of like
borrowing someone's "down with violence" sign to hit them over t
Currently I find that if I call stl() repeatedly I can use the weights array
that is part of the stil output to detect outliers. I also find that if I
repeatedly call stl() (replacing the outliers after each call) that the
"remainder" portion of the stil output gets reduced. I am calling it like
An alternative is to use txtStart (and other functions mentioned in the same
help page) from the TeachingDemos package. This does the sinking, but can also
include the commands as well as allow you to insert comments. The etxt
variants allow you to postprocess the whole transcript into a postsc
Anybody know how to rotate shapes generated with the 'points'
function? I'm trying to place points around a radial diagram such
that the y-axis of individual shapes are oriented with the radii of
the circle rather than the y-axis of the larger plot area. Perhaps
something analogous to the
On Monday 13 October 2008, Waichler, Scott R wrote:
> I would like to use a map or aerial photo as a background to plotting
> solid lines and text, and semi-transparent color contours, in base and
> lattice graphics. Plot coordinates need to be consistent with the
> georeferenced background. For
Hi R helpers,
One rather statistical question?
What would be the best startegy to shortlist thousands of continous
variables automaticaly using R
as the preparation for logistic regression modleing!
Thanks
__
R-help@r-project.org mailing list
ht
I would like to use a map or aerial photo as a background to plotting
solid lines and text, and semi-transparent color contours, in base and
lattice graphics. Plot coordinates need to be consistent with the
georeferenced background. For example, a color contour plot would have
an gray-toned aeri
Great,
Thanks,
Michael
On Mon, Oct 13, 2008 at 3:56 PM, Rolf Turner <[EMAIL PROTECTED]>wrote:
>
> On 14/10/2008, at 9:42 AM, Michael Just wrote:
>
> What is the difference when including or not including the intercept when
>> using lm()?
>>
>> x.noint <- lm(weight ~ group - 1))# omitting interce
On 14/10/2008, at 9:42 AM, Michael Just wrote:
What is the difference when including or not including the
intercept when
using lm()?
x.noint <- lm(weight ~ group - 1))# omitting intercept
x <- lm(weight ~ group))
This has nothing to do with forcing the intercept to 0, correct?
On the cont
What is the difference when including or not including the intercept when
using lm()?
x.noint <- lm(weight ~ group - 1))# omitting intercept
x <- lm(weight ~ group))
This has nothing to do with forcing the intercept to 0, correct?
Thank you kindly,
Michael
[[alternative HTML version del
- Original Message -
From: "Frank E Harrell Jr" <[EMAIL PROTECTED]>
To: "John Sorkin" <[EMAIL PROTECTED]>
Cc: ; <[EMAIL PROTECTED]>;
<[EMAIL PROTECTED]>
Sent: Monday, October 13, 2008 2:09 PM
Subject: Re: [R] Fw: Logistic regresion - Interpreting (SENS) and (SPEC)
John Sorkin wrote
Hello,
I have the data set of 1 + 49 variables. One of them is binary other
are continous.
I would like to be able to fit the model with all 49 variables and
then run stepwise model selction.
I'd appriciate some code snippets...
__
R-help@r-project.
Thanks for the swift response.
cat it is.
Cheers,
Michael
On Mon, Oct 13, 2008 at 3:14 PM, Jorge Ivan Velez
<[EMAIL PROTECTED]>wrote:
>
> Dear Michael,
> You can use cat() as follows:
>
> sink("test.txt")
> cat('Here goes your text','\n','\n') # "\n" writes white spaces
> summary(x)
> sink()
>
Dear Michael,
You can use cat() as follows:
sink("test.txt")
cat('Here goes your text','\n','\n') # "\n" writes white spaces
summary(x)
sink()
See ?cat for more information.
HTH,
Jorge
On Mon, Oct 13, 2008 at 4:02 PM, Michael Just <[EMAIL PROTECTED]> wrote:
> Hello,
> How can I add notes (i
On 13-Oct-08 20:02:20, Michael Just wrote:
> Hello,
> How can I add notes (i.e. text) to a sink output?
>
> sink("test.txt")
>#This text will describe the test
> summary(x)
> sink()
>
> How can I add that text above to the sink output?
> Thanks,
> Michael
Anything on the lines of:
sink("test.
On 14/10/2008, at 9:02 AM, Michael Just wrote:
Hello,
How can I add notes (i.e. text) to a sink output?
sink("test.txt")
#This text will describe the test
summary(x)
sink()
How can I add that text above to the sink output?
?cat
##
Hello,
How can I add notes (i.e. text) to a sink output?
sink("test.txt")
#This text will describe the test
summary(x)
sink()
How can I add that text above to the sink output?
Thanks,
Michael
[[alternative HTML version deleted]]
__
R-help@r-p
Hi Matthieu,
> Does anybody have experience with Sweave run from Kile? I'm trying to
> make it run but have problems and don't know if the instructions are
> false or I do something wrong (my knowledge in bash and shell is too low
> to understand it)...
...
It would help if you stated that you us
R-help:
I have a variable ("ID_list") containing about 1800 unique numbers, and a
143066x29 data frame. One of the columns ("ID") in my data frame contains a
list of ids, many of which appear more than once. I'd like to find the
subset of my data frame for which "ID" matches one of the numbers
Dear Eva;
I shouldn't have sent my unhelpful reply to the entire list, since it
is now glaringly obvious that I did not carefully read your original
question. You are outside my experience, since I have not used lme4,
but I wonder if questions about over-dispersion shouldn't be handled
b
Liviu: Thanks for the links, I'll check them out. On a different note, have you
used MiKTEX at all? I have downloaded it but I don't know how to make it work.
Sweave and Stangle seem to work fine but when I use "texi2dvi" it crashes.
library(tools)
Sweave("C:/Program Files/R/R-2.7.2/bin/foo.Rn
John Sorkin wrote:
Frank,
Perhaps I was not clear in my previous Email message. Sensitivity and specificity do tell us about the quality of a test in that given two tests the one with higher sensitivity will be better at identifying subjects who have a disease in a pool who have a disease, and th
I would like to be able to do the xyplot in ggplot below. I read in
the archive that Hadley was working on this for the next release, and
I can not find the documentation (Aug. 23rd).
River.Mile <- c(215 ,202, 198, 190, 185, 179, 148, 119, 61)
Cu <- rnorm(9)
Fe <- rnorm(9)
Mg <- rnorm(9)
Ti <- rn
On Oct 13, 2008, at 5:36 AM, Dieter Menne wrote:
Prof Brian Ripley stats.ox.ac.uk> writes:
Yes: DF[is.na(match(row.names(DF), exclude_me)), ]
Assuming everything is possible in R: would it be possible to make
the below
work without breaking existing code?
a <- data.frame(x=1:10)
rowna
Frank,
Perhaps I was not clear in my previous Email message. Sensitivity and
specificity do tell us about the quality of a test in that given two tests the
one with higher sensitivity will be better at identifying subjects who have a
disease in a pool who have a disease, and the more sensitive t
John Sorkin wrote:
Jumping into a thread can be like jumping into a den of lions but here goes . .
.
Sensitivity and specificity are not designed to determine the quality of a fit (i.e. if your model is good), but rather are characteristics of a test. A test that has high sensitivity will proper
I don't know if I understand (small example with R command wouuld help),
but, assuming your data.frame is called 'df'
subset(df, ID %in% ID_list)
Question, is ID_list a "list" or a vector, and are they really "numbers"
or "factors"?
Kyle. wrote:
R-help:
I have a variable ("ID_list") contai
Do you have an example. I am not sure what you mean.
On Mon, Oct 13, 2008 at 9:48 AM, guria <[EMAIL PROTECTED]> wrote:
>
> How Graphs in R with leveling of point can be done?
> Please help.
> --
> View this message in context:
> http://www.nabble.com/graphs-in-R-tp19955281p19955281.html
> Sent f
Use the task scheduler in Windows and have a batch files executed.
On Mon, Oct 13, 2008 at 11:44 AM, Tony Breyal
<[EMAIL PROTECTED]> wrote:
> Dear R-Help,
>
> Is it possible to set R up to run a particular script at specific
> times of the day? trivial example: If the time is now 8:59:55am and I
>
(Sorry for the duplicate posting, first posting did not contain text in the
body of the message.)
Greetings,
I am attempting to plot the model fits of a generalized additive model using
the gam package (gam version 1.0, R v. 2.6.2). The gam object is created
without any apparent problem, but w
Dear R users,
I'm trying to estimate the matrix of regression parameters.
I need to do it numerically, so I used optim and nls.
I got the initial parameter estimates from least squares, and input them into
those functions.
But when I run the optim function, it stops in 30 seconds and shows
'con
Dear useRs,
I am pleased to announce the availability of the new package 'StatMatch'
(version 0.4)
http://cran.at.r-project.org/web/packages/StatMatch/index.html
'StatMatch' contains some functions to perform Statistical Matching.
Statistical Matching methods aim at integrate two samples, referr
R-help:
I have a variable ("ID_list") containing about 1800 unique numbers, and a
143066x29 data frame. One of the columns ("ID") in my data frame contains a
list of ids, many of which appear more than once. I'd like to find the
subset of my data frame for which "ID" matches one of the numbers i
Dear R-Help,
Is it possible to set R up to run a particular script at specific
times of the day? trivial example: If the time is now 8:59:55am and I
wish to run a function at 9am, I do the following:
my.function <- function(x) {
p1 <- proc.time()
Sys.sleep(x)
print('Hello R-Help!')
Dear Michal,
I had the same problem of you in installing quadprog packages ?
Did you resolve it?
Can you help me?
the error is:
* Installing *source* package 'quadprog' ...
** libs
gfortran -fpic -g -O2 -c aind.f -o aind.o
In order to use gfortran please type either:
source /usr/local/free/gfor
__
R-help@r-project.org 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, minimal, self-contained, reproducible code.
How Graphs in R with leveling of point can be done?
Please help.
--
View this message in context:
http://www.nabble.com/graphs-in-R-tp19955281p19955281.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
Dear All,
I have a query relating to the use of the ‘predict’ and ‘gamm’ functions. I am
dealing with large (approx. 5000) sets of presence/absence data, which I am
trying to model as a function of different of environmental covariates.
Ideally my models should include individual and colony as
Try this:
split(iris, factor(a, levels = 1:3))
On Mon, Oct 13, 2008 at 1:06 PM, Jay <[EMAIL PROTECTED]> wrote:
> Hello,
>
> I'll use part of the iris dataset for an example of what I want to
> do.
>
>> data(iris)
>> iris<-iris[1:10,1:4]
>> iris
> Sepal.Length Sepal.Width Petal.Length Petal.Widt
Thanks so much.
On Oct 13, 1:14 pm, "Henrique Dallazuanna" <[EMAIL PROTECTED]> wrote:
> Try this:
>
> a<-factor(c(3, 3, 3, 2, 3, 2, 2, 3, 2, 3), levels = 1:3)
> split(iris, a)
>
> lapply(split(iris, a), dim)
>
>
>
> On Mon, Oct 13, 2008 at 2:06 PM, Jay <[EMAIL PROTECTED]> wrote:
> > Hello,
>
> > I
On 10/13/08, Michael Lawrence <[EMAIL PROTECTED]> wrote:
>
>
> On Sun, Oct 12, 2008 at 4:50 PM, Dirk Eddelbuettel <[EMAIL PROTECTED]> wrote:
> >
> >
> > On 12 October 2008 at 12:53, cls59 wrote:
> > | On a related note... does anyone know good resources for binding a C++
> > | program to the R libr
Try this:
a<-factor(c(3, 3, 3, 2, 3, 2, 2, 3, 2, 3), levels = 1:3)
split(iris, a)
lapply(split(iris, a), dim)
On Mon, Oct 13, 2008 at 2:06 PM, Jay <[EMAIL PROTECTED]> wrote:
> Hello,
>
> I'll use part of the iris dataset for an example of what I want to
> do.
>
> > data(iris)
> > iris<-iris[1:
Hello,
I'll use part of the iris dataset for an example of what I want to
do.
> data(iris)
> iris<-iris[1:10,1:4]
> iris
Sepal.Length Sepal.Width Petal.Length Petal.Width
1 5.1 3.5 1.4 0.2
2 4.9 3.0 1.4 0.2
3 4.7
You could look at the xspline function. It approximates b-splines or Bezier
curves given control points and shape parameters. It can either plot the
spline or return a bunch of point on the curve for comparison to other values.
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data
Jumping into a thread can be like jumping into a den of lions but here goes . .
.
Sensitivity and specificity are not designed to determine the quality of a fit
(i.e. if your model is good), but rather are characteristics of a test. A test
that has high sensitivity will properly identify a large
On Mon, 2008-10-13 at 16:28 +0200, [EMAIL PROTECTED] wrote:
> > If you used daisy, is there a problem with converting the resulting
> > object to a full dissimilarity matrix and extracting the relevant
> > row/column you need for the target site?
>
> Well, the lost of efficiantcy is huge. I need t
Maithili Shiva wrote:
Dear Mr Peter Dalgaard and Mr Dieter Menne,
I sincerely thank you for helping me out with my problem. The thing is taht I
already have calculated SENS = Gg / (Gg + Bg) = 89.97%
and SPEC = Bb / (Bb + Gb) = 74.38%.
Now I have values of SENS and SPEC, which are absolute in n
and see page 8 of:
http://www.r-project.org/doc/Rnews/Rnews_2003-3.pdf
R as a Simulation Platform in Ecological Modelling
Thomas Petzoldt
"Conclusions
The examples described so far, plus the experience
with R as data analysis environment for measure-
ment and simulation data, allows to conclude
Dieter Menne wrote:
Maithili Shiva yahoo.com> writes:
I havd main sample of 42500 clentes and
based on their status as regards to defaulted / non - defaulted, I have
genereted the probability of default.
I have a hold out sample of 5000 clients. I have calculated (1) No of
correctly classif
There is the zoomplot function in the TeachingDemos package that allows you to
zoom in/out the current plot. But it is a bit of a kludge. The better option
is probably to just set the xlim and ylim arguments in a new plot command. You
can use the locator function as one way to find the coordi
All,
I'm getting the same error message as that discussed in a previous post (Feb
3, 2006). The reply to that post was to insure that the ld98 program was in
the system path (as also suggested in the help on ldBands). I have done
this but this does not change the result. Any advice much apprec
and http://cran.at.r-project.org/web/packages/deSolve/
2008/10/13 megha patnaik <[EMAIL PROTECTED]>
> See http://cran.r-project.org/web/packages/odesolve/index.html
>
> 2008/10/13 <[EMAIL PROTECTED]>
>
>
>>
>> Dear R-users,
>>
>> I am trying to perform some simulations from a model defined by ord
See http://cran.r-project.org/web/packages/odesolve/index.html
2008/10/13 <[EMAIL PROTECTED]>
>
>
> Dear R-users,
>
> I am trying to perform some simulations from a model defined by ordinary
> differential equations. I would be grateful if someone could indicate some
> functions/packages/examples
Dear R-users,
I am trying to perform some simulations from a model defined by ordinary
differential equations. I would be grateful if someone could indicate some
functions/packages/examples I could look at.
Thank you in advance.
Sebastien
__
R-help@r
I am trying to fit a discrete function to my dataset using nls().
fit<-nls(T2~form(SOA,t1weight,t2weight,d1weight),
start=list(t1weight=1,t2weight=1,d1weight=1), data=data1, trace=TRUE)
The problem is that my function ("form") includes a discrete function
and in that function I used the varia
__
R-help@r-project.org 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, minimal, self-contained, reproducible code.
__
R-help@r-project.org 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, minimal, self-contained, reproducible code.
sungard.com> writes:
> There are two good papers that illustrate how to compare classifiers
> using Sensitivity and Specificity and their extensions (e.g., likelihood
> ratios, young index, KL distance, etc).
>
> See:
> 1) Biggerstaff, Brad, 2000, "Comparing diagnostic tests: a simple
> graphi
Hi Maithili,
There are two good papers that illustrate how to compare classifiers
using Sensitivity and Specificity and their extensions (e.g., likelihood
ratios, young index, KL distance, etc).
See:
1) Biggerstaff, Brad, 2000, "Comparing diagnostic tests: a simple
graphic using likelihood ratio
1 - 100 of 127 matches
Mail list logo
