Dear List members,
I am runnning a logistic regression that includes sample weights
(expressed by a variable called WEIGHT) using glmob (from robustbase). My
model looks like this:
B12-glmrob(SERE~COMP+COM,family=binomial,weights=WEIGHT)
When I include the weights, my
On Sun, 9 Jul 2006, Daryl Manning wrote:
Apologies if this is in (one of the many) manuals somewhere... Trying
to switch to R from other stats programs.
Basically, I have a large data table I've dumped from a DB, some of
the values which are nulls '-' which I've converted to zeros. I've
Dear Justin
Ties means that you have identical values in
Year5.lm$residuals. Please remark that you can have a large
R^2, but your residuals are not normally distributed. A large
R^2 shows a strong linear relationship, but that does not say
anything about the error distribution (see example
On Mon, 2006-07-10 at 13:27 +0800, Guo Wei-Wei wrote:
Dear all,
I'm working on a data.frame named en.data, which has n cases and m columns.
I generate the correlation matrix of en.data by
cor(en.data)
I find that there is no p-value on each correlation in the correlation
matrix. I
you can use function rcor.test() from package 'ltm', e.g.,
help(rcor.test, package = ltm)
###
library(ltm)
dat - data.frame(matrix(rnorm(1000), 100, 10))
rcor.test(dat)
rcor.test(dat, method = kendall)
rcor.test(dat, method = spearman)
I hope it helps.
Best,
Dimitris
Dimitris
Dear Chi Kai,
Three years ago there was a similar thread.
At that time David Firth offered a solution for the quasipoission
problem with negative observations, see:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/16143.html
I remember that his code gave you slightly different answers than the
Hi, Gavin, your program is excellent. Thank you very much!
And I have two further questions.
1. Since it is very possible that the data contains missing value and
the program will failed against missing values, I have to delete all
the cases contained NA. Can it be done pairwisely?
2. Can the
On 10-Jul-06 Spencer Graves wrote:
[...]
Excellent overview of many string reasons for R, Spencer! Many
thanks for this well-thought-out advocacy piece.
[...]
At one time, the market dominance of both SAS and IBM was
sustained by the perception that 'nobody ever got fired'
for using them.
Hello everyone
I am very much a newcomer to R and I need to call R processing routines
from java and pass back the results to java. It seems sjava is the way
to do this, but there seems to be a lot of complexity surrounding
configuring R for this purpose. I need to do this under windows XP.
Thank you, Dimitris.
The function rcor.test() is very nice. It can pass arguments to cor()
and solve my first problem of pairwised case deletion.
Best regards,
Wei-Wei
2006/7/10, Dimitris Rizopoulos [EMAIL PROTECTED]:
you can use function rcor.test() from package 'ltm', e.g.,
Brian Lunergan wrote:
Uwe Ligges wrote:
Where did you find them? At least not in their current versions on
the official repository on CRAN master for R-2.3.x, I hope.
FracSim_0.2.zip
http://cran.stat.ucla.edu/bin/windows/contrib/2.3/
RDCOMClient_0.91-0.zip
Does anyone have any code for drawing weighted histograms (a la
Manet/Mondrian) in R?
Thanks,
Hadley
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hi
I'm having difficulty with a multiple plot. What I want is 12
plots, all of which are the same size and shape, differing
only in colour. Each one is a square, so there is an asp=1
in the plot call. I'm working in an Sweave environment so I
am free to choose the height and width of the plot.
(R)SJava is part of the Omegahat project, not the R project. In theory at
least, Omegahat has its own mailing lists.
In any case, the R posting guide asks you to contact the maintainer of a
contributed package before this list, so please do so now.
There are easier-to-setup alternatives to
[apologies for possible multiple posting]
Hi Joris
great suggestion! I never thought to use layout() in this way.
If I have
postscript(file=~/f.ps,width=5,height=8)
nf - layout(matrix(
c(00,01,00,02,00,03,
04,05,00,06,00,07,
Hello all,
I posted a question to this list last week and received no response. I am
unsure if this means no-one knows the answer or if I posed the question badly.
I'm going to assume I posed the question badly and try again. I am new to R so
it is quite likely it's a very naive question,
On Mon, 2006-07-10 at 16:22 +0800, Guo Wei-Wei wrote:
Hi, Gavin, your program is excellent. Thank you very much!
And I have two further questions.
1. Since it is very possible that the data contains missing value and
the program will failed against missing values, I have to delete all
Hi,
I am looking forward if some of you could help me with this. I have been
looking R-functions to conduct tree-ring analysis in R. I know that exists
a COFECHA version in R developed by C. Bigler. But I don't know if is there
another collection of functions to do dendroclimatological
On Mon, 2006-07-10 at 11:56 +0100, Richard Birnie wrote:
Hello all,
Hi Richard,
Sorry, I know nothing about pvclust and have never used it, but here are
a couple of general suggestions/observations.
You are asked to contact the package maintainer *not* R-Help for
questions relating to
Dear All,
Is there a better way to include NA's of a factor in the output of table()
than using as.character()?
Admittedly, I do not understand the help page for table concerning the
exclude argument applied to factors. I tried in different ways, but could
not get NA to be included in the table,
fcv - factor(c('a', NA, 'c'), exclude=NULL)
table(fcv, exclude=NULL)
---
Jacques VESLOT
CNRS UMR 8090
I.B.L (2ème étage)
1 rue du Professeur Calmette
B.P. 245
59019 Lille Cedex
Tel : 33 (0)3.20.87.10.44
Fax : 33 (0)3.20.87.10.31
If we specify exclude=NULL in factor then it need not also be
specified in table:
fcv - factor(c('a', NA, 'c'), exclude=NULL)
table(fcv)
On 7/10/06, Jacques VESLOT [EMAIL PROTECTED] wrote:
fcv - factor(c('a', NA, 'c'), exclude=NULL)
table(fcv, exclude=NULL)
Thanks for the responses,
Now I know why I got no response I'll know better in the future. Apologies all
round for any inconvenience caused that was never my intention, just
inexperience.
Richard
Dr Richard Birnie
Scientific Officer
Section of Pathology and Tumour Biology
Welcome Brenner
Hello.
Is it possible to choose the distance in the kmeans algorithm?
I have m vectors of n components and I want to cluster them using kmeans
algorithm but I want to use the Mahalanobis distance or another distance.
How can I do it in R?
If I use kmeans, I have no option to choose the
On Mon, 2006-07-10 at 12:48 +0100, Gavin Simpson wrote:
On Mon, 2006-07-10 at 16:22 +0800, Guo Wei-Wei wrote:
Hi, Gavin, your program is excellent. Thank you very much!
And I have two further questions.
1. Since it is very possible that the data contains missing value and
the
R 2.3.1
Windows XP
Question: How can I obtain the values of BIAS and STD. ERROR from a
bootstrap.
Background:
I am running a bootstrap:
result2-boot(1:400,regSEvssample,R=5000)
and obtain the following results:
Bootstrap Statistics :
originalbiasstd. error
t1* 1.876602
Look at the code of boot:::print.boot , which does the calculations.
On Mon, 10 Jul 2006, John Sorkin wrote:
R 2.3.1
Windows XP
Question: How can I obtain the values of BIAS and STD. ERROR from a
bootstrap.
Background:
I am running a bootstrap:
result2-boot(1:400,regSEvssample,R=5000)
Thank you Jacques and Gabor!
Your solution does work and it led me to try also:
table(factor(fcv, exclude=NULL)) # shows a, c, NA
Greetings,
Heinz
At 08:34 10.07.2006 -0400, Gabor Grothendieck wrote:
If we specify exclude=NULL in factor then it need not also be
specified in table:
fcv -
Hi,
I have calculated chi-square goodness of fit test,Sample coming from
Poisson distribution.
please copy this script in R run the script
The R script is as follows
## start
#
No_of_Frauds-
I have a data frame with named columns and I would like to know if it is
possible to retrieve a column name once selected:
print(colnames(df)) # assumes to print col1 col2
print.name(df$col1) # would like to print col1
print.name(df$col2) # would like to print col2
So what the print.name
Hi,
I have calculated chi-square goodness of fit test,Sample coming from
Poisson distribution.
please copy this script in R run the script The R script is as follows
## start
#
No_of_Frauds-
On Mon, 10 Jul 2006, Laurent Deniau wrote:
I have a data frame with named columns and I would like to know if it is
possible to retrieve a column name once selected:
Not really. df$col1 is a new object which does not know where it came
from.
If you wanted to do this before selection, then
On Fri, 7 Jul 2006, Valentin Dimitrov wrote:
I do not need a accelerated failure model, but a
proportional hazard model with a f0= weibull,
exponential, loglogistic or lognormal baseline
distribution. The hazard function is
lambda(t)=exp(Xi*beta)*lambda0(t),
where lambda0 is the baseline
On Sat, 8 Jul 2006, S?ren Merser wrote:
Hi
It seems that survfit() doesn't accept the argumnet 'error' as below
survfit(fit, error='greenwood')
Error in survfit.coxph(fit, error = greenwood) :
unused argument(s) (error ...)
Isn't is allowed to do that for a coxph object?
Looking
priti desai [EMAIL PROTECTED] writes:
Hi,
I have calculated chi-square goodness of fit test,Sample coming from
Poisson distribution.
please copy this script in R run the script
The R script is as follows
## start
#
Use drop = FALSE. For example using builtin data frame BOD we
can display the Time column with its heading:
BOD[, Time, drop = FALSE]
On 7/10/06, Laurent Deniau [EMAIL PROTECTED] wrote:
I have a data frame with named columns and I would like to know if it is
possible to retrieve a column
On Mon, 10 Jul 2006, Gabor Grothendieck wrote:
The problem can be reduced to this:
x - 1
x[1] - 2 # error
The following are ok:
x - 1
x[1] - 3
x - 1
x - 4
x - 1
x - 5
Does anyone know why? Is this a bug in - ?
No, it's a feature. The fact that x-5 works is arguably a bug
roger == roger koenker [EMAIL PROTECTED]
on Sun, 9 Jul 2006 12:31:16 -0500 writes:
roger On Jul 9, 2006, at 11:06 AM, Douglas Bates wrote:
Your matrix Mc should be flagged as invalid. Martin and
I should discuss whether we want to add such a test to
the validity method.
Prof Brian Ripley wrote:
On Mon, 10 Jul 2006, Laurent Deniau wrote:
I have a data frame with named columns and I would like to know if it is
possible to retrieve a column name once selected:
Not really. df$col1 is a new object which does not know where it came
from.
If you wanted
Hi,
I'm drawing a very simple plot with both axes logarithmic (default base 10).
Example:
vec=c(1,10,100,1000,1,10,100,1000)
plot(vec,vec,log=xy)
The axes on the plot now show the technical notation like 1E+3 but I would
prefer to have it the notation 10 ^3 i.e. with the
On 7/10/2006 12:43 PM, [EMAIL PROTECTED] wrote:
Hi,
I'm drawing a very simple plot with both axes logarithmic (default base 10).
Example:
vec=c(1,10,100,1000,1,10,100,1000)
plot(vec,vec,log=xy)
The axes on the plot now show the technical notation like 1E+3 but I would
You can always draw the axes by hand.
?par with axes =FALSE
?axis
?plotmath for mathematical notation in R (for exponents)
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the scientific learning
process. - George E. P.
With some help from those with expertise on this list, I managed to
produce a plot using trellis that looked like I wanted it to look.
Now, I need to take the same plot and make the lines on it color, but
I want to specify the color for the lines myself.
I've managed to make the key use the
Use col= to specify colors, e.g.
library(lattice)
x - 1:12
xyplot(x ~ x, group = gl(3,4), col = 1:3, type = l, auto.key = TRUE)
If this is not sufficiently close to your problem
1. cut your example down to a *minimal* size and
2. provide it as *self contained* and
3. *reproducible*
Carlos and Uwe:
I apologize to all, but my server has crashed and I need to order
a new computer! I'm at home now, but when I can get to my office
I'll try to send each of you a copy.
I will be submitting my packages to CRAN soon. Uwe, perhaps the
quad3d stuff should just be incorporated into
The col command worked.
Thanks!
Jamie
On 7/10/06, Gabor Grothendieck [EMAIL PROTECTED] wrote:
Use col= to specify colors, e.g.
library(lattice)
x - 1:12
xyplot(x ~ x, group = gl(3,4), col = 1:3, type = l, auto.key = TRUE)
If this is not sufficiently close to your problem
1.
Try looking at the tkwait.window function, it may be what you need. Create
your tkwindow with all its components, then call tkwait.window on the toplevel
window and the calling function will wait until that window goes away before
continuing execution.
Hope this helps,
--
Gregory (Greg)
On 7/10/06, Jamieson Cobleigh [EMAIL PROTECTED] wrote:
With some help from those with expertise on this list, I managed to
produce a plot using trellis that looked like I wanted it to look.
Now, I need to take the same plot and make the lines on it color, but
I want to specify the color for
There are a couple of different options that you can try (each
generalize a bit differently to other situations).
If you just want lines at the tick marks to go the entire width (and be
generated automatically to match the labels), then look at ?par and look
at 'tck', you can use this in the
I am a very novice R user, a social scientist (linguist) who is trying
to learn to use R after being very familiar with SPSS. Please be kind!
My concern:
I cannot figure out a way to get an accurate count of observations of
one column of data split by a factor when there are NAs in the data.
I
Wouldn't something like table(status) give you want you want? E.g.:
R status - factor(c(A, B, A, NA, A, B))
R table(status)
status
A B
3 2
Andy
From: Jenifer Larson-Hall
I am a very novice R user, a social scientist (linguist) who
is trying to learn to use R after being very familiar
Dear Dr. Graves,
Thanks for your information and kindness! I really appreciate it.
According to my findings so far, POLYCLASS doesnot seem to allow for
variables being forced into the model. Even the commercialized MARS
doesnot have this function either. However, one ad hoc approach that
they
Jenifer Larson-Hall [EMAIL PROTECTED] writes:
I am a very novice R user, a social scientist (linguist) who is trying
to learn to use R after being very familiar with SPSS. Please be kind!
My concern:
I cannot figure out a way to get an accurate count of observations of
one column of data
Hi Folks,
I'm trying to locate the source code for a (typical)
R graphics device, in order to study how it's done.
The underlying reason is that I'm thinking of trying
to create a graphics device for 'pic' (the diagram
drawing component of [g]troff).
I thought the xfig device would be a good
Dear R Users,
Is it possible to make R print the largest item in each row of a matrix X
with red font? Example:
1247
8431
...
Therefore 7 and 8 should be in red color.
I would appreciate any suggestion
Robert McFadden
That's correct. The option is not there.
You can force it in a starting model - but it can still be deleted,
At some stage Insightful implemented a commercial set of polynomial
spline functions in S-Plus that does have that option. They got it
as far as a beta-library [quoting from my own
Hello,
I am new to R and still feeling my way thru it.
I am trying to plot the values from this file below on the X-axis of a plot. I
have attached the graph to the email...the one i am trying to recreate.
Exonstart end
5'UTR 2254006022540121
1 22540122
On 7/10/2006 3:52 PM, (Ted Harding) wrote:
Hi Folks,
I'm trying to locate the source code for a (typical)
R graphics device, in order to study how it's done.
The underlying reason is that I'm thinking of trying
to create a graphics device for 'pic' (the diagram
drawing component of
Dear Gang,
I'm not clear about the form of the model that you're trying to fit,
nor what the form of the likelihood would be for such as model, but the
short answer is that sem() requires that each covariance among observed
variables be for the same observations and hence for the same number
(N)
On 10-Jul-06 Ted Harding wrote:
Hi Folks,
I'm trying to locate the source code for a (typical)
R graphics device, in order to study how it's done.
[...]
When, in R, I type xfig to see the R code, I get the
line
.Internal(XFig(file, old$paper, old$family, old$bg, old$fg,
One option is
library(R2HTML)
?HTML.cormat
The thing you're after is traffic highlighting (via CSS or HTML tags).
If HTML.cormat() doesn't do exactly what you want, modify the source
code. (By the way, I haven't used R2HTML so far so maybe there's a
more appropriate function.)
-Original
[EMAIL PROTECTED] wrote:
Hi,
I'm drawing a very simple plot with both axes logarithmic (default base 10).
Example:
vec=c(1,10,100,1000,1,10,100,1000)
plot(vec,vec,log=xy)
The axes on the plot now show the technical notation like 1E+3 but I would
prefer to have it the
dear R wizards: sorry, I am stumped. what is the parameter to just
move the location of the x-label and y-lable (not the labels on the
ticks of the x-axis and y-axis)? probably obvious, but not to me
right now... regards, /iaw
__
Does this do what you want:
x - Exonstart end
5'UTR 2254006022540121
1 2254012222540140
2 2254030322540493
3 2254155222541565
4 2254237322542519
5 2254426522544432
3'UTR 2254443322544856
y -
Hello there,
This question is relative to WindowsXP, using R 2.2.1:
I am looking for a function that allows a user to interactively choose a
directory so I can use list.files to process all the files in that
directory. I've looked at getwd, but this is not interactive. The
functions
Hi,
for testing if the slope of experimental data differs from a
given slope I'm using the function
test_regression_against_slope (see below).
I am now confronted with the problem that I have data which
requires a modelII regression (also called reduced major axes
regression (RMA) or geometric
choose.dir
it is documented on the same page as choose.files
__
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
On 7/10/2006 8:51 PM, Randy Zelick wrote:
Hello there,
This question is relative to WindowsXP, using R 2.2.1:
Time to upgrade. choose.dir (mentioned by Rich) was introduced in the
next release.
Duncan Murdoch
I am looking for a function that allows a user to interactively choose a
dear R wizards: is it possible to instruct R to save or no-save from
inside R? or does this have to be given at invokation on the
command-line?The same question applies to --no-restore-data,
although this presumably would have to be decided in a .First()
function or something like it.
on a
Thanks to all who responded. I upgraded to 2.3.1 straight away, and
choose.dir does exactly what I need.
Cheers,
=Randy=
R. Zelick email: [EMAIL PROTECTED]
Department of Biology voice: 503-725-3086
Portland State University fax:
ivo welch wrote:
dear R wizards: is it possible to instruct R to save or no-save from
inside R?
You can specify the action when you call q() or quit(), but you can't
change the default action.
or does this have to be given at invokation on the
command-line?The same question applies to
Hi there,
I tried to install gstat package to R in Linux. I follow the
instruction to use command R CMD INSTALL -1 lib pkgs. But R
shows syntax error in R CMD.
Can someone help me with the installation of gstat package?
Thank you!
Jing
__
72 matches
Mail list logo