Maithili Shiva maithili_shiva at 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 classified
Alex Karner aakarner at ucdavis.edu writes:
I realize these limitations. However, I know that my actual dataset is
reasonably well behaved in the range I want to predict, and I'm not using
the predicted values for any further analysis, only for schematic purposes
in the plot.
I'm still
Hi,
I have a really large graph and would like to zoom in on portions of the
graph and post them as blocks below the graph.Is there an add on package to
do this?
--
Rajesh.J
I skate to where the puck is going to be, not where it has
On Sun, 12 Oct 2008, [EMAIL PROTECTED] wrote:
On 12-Oct-08 22:09:46, jim holtman wrote:
Seem to work fine on my R 2.7.2 version of Windows:
png(file=myplot.png, bg=transparent, units='cm',
width=12,height=15, res=200)
plot(1:10)
rect(1, 5, 3, 7, col=white)
dev.off()
Did you check the
On Sun, Oct 12, 2008 at 1:39 AM, Felipe Carrillo
[EMAIL PROTECTED] wrote:
I am working on a publication and I have heard about LaTEX but I haven't
actually tried to learn about it until today. I've found a few
There are two more packages that might be of interest:
RReportGenerator [1] and
On Sun, 12 Oct 2008, David Winsemius wrote:
Two follow-up questions:
A) I get an error message when using Harrell's describe() function on one of
my variable, telling me that sum() is not meaningful for a difftime object.
Why should sum() not be meaningful for a collection of interval
Dieter Menne wrote:
Maithili Shiva maithili_shiva at 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
Peter Dalgaard p.dalgaard at biostat.ku.dk writes:
But if he has a hold out sample, isn't he already cross-validating??
I wonder if you're answering the right question there. Could he just
be looking for Sp=Gg/(Gg+Bg), Se=Bb/(Gb+Bb)? (If I got the notation
right.)
You are right. My
I have already used aidsEst() function in the micEcon package for the
estimation of elasticities.
But with the LES, I plan to estimate the minimal consumption
level...which is not possible with the AIDS model, do I?
Arne Henningsen wrote:
Hi Marie!
On Friday 10 October 2008 12:40:23,
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 nature. My question
Thanks for your reply Mark,
but no, using predict on the new data.frame does not help here.
I had first thought that the probelm was due the explanatory variable
(age) and the offset one (date) being very similar (highly correlated, I am
trying to tease their effect apart, and hoped offset
On Mon, 13 Oct 2008, Peter Dalgaard wrote:
Dieter Menne wrote:
Maithili Shiva maithili_shiva at 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
If you want to remove the N, then you can work with the indices:
x
[1] NA B NA B B NA N A B B A NA A N N N
A B B A
# if you want the indices of the non-N, then
(indx - which(is.na(x) | x != N))
[1] 1 2 3 4 5 6 8 9 10 11 12 13 17 18 19 20
x[indx]
[1] NA B NA B B NA A B B A NA
Is this what you want:
x.df
timestamp id
1 2008-05-27 22:57:00 763830873067
2 2008-05-27 23:00:00 763830873067
3 2008-05-27 23:01:00 763830873067
4 2008-05-27 23:01:00 763830873067
5 2008-06-05 11:34:00 763830873067
6 2008-05-29 23:08:00 765253440317
7 2008-05-29
Not sure exactly what you are trying to do since you did not provide
commented, minimal, self-contained, reproducible code. Let me take a
guess in that you also have to test for NAs:
x - sample(c(N, A, B, NA), 20, TRUE)
x
[1] A A B NA N NA NA B B N N N B A NA A
B NA A NA
x != N
[1]
Hi the list,
I need to compute Gower distance between a specific individual and all
the other individual.
The function DAISY from package cluster compute all the pairwise
dissimilarities of a population. If the population is N individuals,
that is arround N^2 distances to compute.
I need
Hello everybody,
I would like to read a SAS Data data1.sas7bdat in R! Is this possible?
Thank you a lot in advance ;),
Stefo
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Stefo Ratino wrote:
Hello everybody,
I would like to read a SAS Data data1.sas7bdat in R! Is this possible?
Only if you have access to SAS. The file format is proprietary and noone
has bothered to decipher it.
Thank you a lot in advance ;),
Stefo
library(foreign)
Rdata=read.ssd(Z:/MyFolder,data1,sascmd = C:/Program Files/SAS/SAS
9.1/sas.exe)
Rdata
Y D1 D2 D3
1 100 1 0 0
2 101 1 0 0
3 105 1 0 0
4 200 0 1 0
5 201 0 1 0
6 205 0 1 0
7 300 0 0 1
8 301 0 0 1
9 305 0 0 1
where 'data1' is the SAS datafile to be
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 library?
RCpp, at http://rcpp.r-forge.r-project.org, formerly known as
Hello
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)...
I discovered recently Sweave and wanted to run it
I use Sweave on a Mac from my LaTeX editor and it requires a similar shell
script to typeset the document. Mine uses R it's self (instead of sh) and
looks like:
#!/usr/bin/Rscript
args - commandArgs(T)
fname - strsplit(args[1],'\\.')[[1]][1]
Sweave(paste(fname,'Rnw',sep='.'))
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 to compute the distance
several time on data base that count 1000 or
The functions `gam' in packages `mgcv' or `gam', will do this, as will the
ssanova* functions in `gss' (there are numerous alternatives)
On Friday 10 October 2008 23:32, Verschuere Benjamin wrote:
Hi,
I was wondering if there is a function in R which performs bivariate non
parametric
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
Pedro.Rodriguez at 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
Try this:
x$Timestamp - as.POSIXct(strptime(as.character(x$Timestamp), %d/%m/%y
%H:%M))
x$time - as.numeric(x$Timestamp as.POSIXct(strptime(07-08-2008 00:00,
%d-%m-%Y %H:%M)))
with(x, tapply(Timestamp, user_id, function(x)diff(range(x), units=secs),
simplify = F))
On Mon, Oct 13, 2008 at 7:57
__
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.
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
__
Hi All,
I have a data frame which has columns comprised mainly of NAs. I know
there are functions na.pass and na.omit etc which can be used in these
situations however I can't them to work in this case. I have a function
which returns the data according to some rule i.e. removal of N in this
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 I could
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 ordinary
Hi,
To simplify, suppose I have 2 observations each day for three days. I
would like to define the correlation structure of these 6 observations
as follows: the correlation of 2 observations on the same day is, say,
alpha, the correlation for 2 observations one day apart is rho and the
Prof Brian Ripley ripley at 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)
rownames(a) = letters[1:10]
exclude = c(a,c)
Hi All,
I've a couple of questions i've been struggling with using the time
features, can anyone help? sample data
Timestampuser_id
27/05/08 22:57 763830873067 27/05/08 23:00 763830873067 27/05/08 23:01
763830873067 27/05/08 23:01 763830873067 05/06/08 11:34
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
I have a data frame. It has lots of patient information, their age, their
gender, etc etc. I need to keep all this information whilst selecting
relevant rows. So, in the example of code I provided I want to remove all
those patients who have entry N in the column with.Wcode. The dimension of
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
Dieter Menne wrote:
Maithili Shiva maithili_shiva at 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
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
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
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
On Monday 13 October 2008 09:14:23, Marie Vandresse wrote:
I have already used aidsEst() function in the micEcon package for the
estimation of elasticities.
But with the LES, I plan to estimate the minimal consumption
level...which is not possible with the AIDS model, do I?
You are right;
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 to
On Mon, 13 Oct 2008, Dieter Menne wrote:
Prof Brian Ripley ripley at 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?
It would be possible, but not I
Felix Andrews ha scritto:
thigmophobe.labels in the plotrix package tries to avoid label crashes, and
Thank you, but I choosed
There is also pointLabel() in the maptools package.
This work fine. Thank you very much.
Nicola
__
R-help@r-project.org
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
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
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
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:10,1:4]
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 library?
RCpp,
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'll use part of the
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.Width
1
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
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
__
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.
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
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!')
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 in
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,
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
(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
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
wish
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 from the
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) containing
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
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
On Oct 13, 2008, at 5:36 AM, Dieter Menne wrote:
Prof Brian Ripley ripley at 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 -
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 -
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
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.Rnw)
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
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 in
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 use
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]]
__
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
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.txt)
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.e.
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()
See ?cat for
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...
__
- Original Message -
From: Frank E Harrell Jr [EMAIL PROTECTED]
To: John Sorkin [EMAIL PROTECTED]
Cc: r-help@r-project.org; [EMAIL PROTECTED];
[EMAIL PROTECTED]
Sent: Monday, October 13, 2008 2:09 PM
Subject: Re: [R] Fw: Logistic regresion - Interpreting (SENS) and (SPEC)
John
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
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
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 intercept
x -
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
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
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
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
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
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:
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
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:
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
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
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
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))# omitting
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
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
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
1 - 100 of 125 matches
Mail list logo
