Dear R users,I am looking for an easy (i.e., direct) way in obtaining the F and
p values from the intercept in anovas with within-subject designs. Data are
from a psychophysics expeirment where I am using d' (d-prime) values obtained
from 3 modalities of presentation in each subject.I would
]
(801) 408-8111
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Sid Kouider
Sent: Thursday, August 30, 2007 4:33 AM
To: r-help@stat.math.ethz.ch
Subject: [R] (no subject)
Dear R users,I am looking for an easy (i.e., direct) way in
obtaining
I don't think, this has been answered:
I'm trying to run a 3-way within-subject anova in lme with 3
fixed factors (Trust, Sex, and Freq), but get stuck with handling
the random effects. As I want to include all the possible random
effects in the model, it would be something more or less
Hi Lorenz,
I really appreciate your comments.
If I understand correctly, you want to include the interactions
between the random and fixed terms?
Yes that is exactly I wanted to model.
This is done like:
model.lme - lme(Beta ~ Trust*Sex*Freq,
random =
If I understand correctly, you want to include the interactions
between the random and fixed terms?
Yes that is exactly I wanted to model.
This is done like:
model.lme - lme(Beta ~ Trust*Sex*Freq,
random = ~Trust*Sex*Freq|Subj, Model)
But this needs a lot of
I'm trying to run a 3-way within-subject anova in lme with 3 fixed
factors (Trust, Sex, and Freq), but get stuck with handling the
random effects. As I want to include all the possible random effects
in the model, it would be something more or less equivalent to using aov
fit.aov -
mara serrano teachermara at yahoo.com writes:
hi,
i'm new to R and i'm trying to port a quattro pro spreadsheet into R.
spreadsheets have optional lower and
upper limit parameters on the beta distribution function. i would like to know
how to incorporate this
with R's pbeta function.
hi,
i'm new to R and i'm trying to port a quattro pro spreadsheet into R.
spreadsheets have optional lower and upper limit parameters on the beta
distribution function. i would like to know how to incorporate this with R's
pbeta function.
thanks in advance,
mara.
Lawrence D. Brenninkmeyer [EMAIL PROTECTED] writes:
I am trying to find a way to diffuse GIS data on a European map. I have a
dataset consisting of particular locations scattered across Europe,
along with magnitude and value information.
The dataset is a spatial point pattern, where the
Dear Sir/Madam,
I am running a R program for my SNP data. There are some errors when I run
glm model in Hapassoc software, sometimes it is over the memory and
sometimes the matrix is singular. I want to ingore these errors and excute
the next statement. Otherwise, I have a big big trouble. Do
It's a FAQ (Question 7.32)
-thomas
On Mon, 23 Jul 2007, [EMAIL PROTECTED] wrote:
Dear Sir/Madam,
I am running a R program for my SNP data. There are some errors when I run
glm model in Hapassoc software, sometimes it is over the memory and
sometimes the matrix is singular. I want
Hi!
When I use lmer (lme4 package) it accept the correlation structure like lme or
glmmPQL but in the output I have no information about it. Can anyone help me
how can I get it.
Thanks
Salomé
__
R-help@stat.math.ethz.ch mailing list
Hi,
I want to know how to increase the limit size of R object. For example, i
want to store a matrix (1 rows and 8000 columns)into an object x.
However, it gives me an error message 'can't allocate a vector of size
396.7 Mb'. But there is no problem to store a matrix of 5000*3000 matrix
into
Hi All,
I want to automatically generate a number of data frames, each with an
automatically generated name and an automatically generated number of
rows. The number of rows has been calculated before and is different for
all data frames (e.g. c(4,5,2)). The number of columns is known a priori
Is this what you want to do:
auto.length - c(12,15,6)
for(i in 1:3) {
+ nam - paste(auto.data,i, sep=.)
+ assign(nam, as.data.frame(matrix(1:auto.length[i], ncol=3)))
+ }
auto.data.1
V1 V2 V3
1 1 5 9
2 2 6 10
3 3 7 11
4 4 8 12
auto.data.2
V1 V2 V3
1 1 6 11
2 2 7 12
3 3 8 13
Jonathan
the BACCO bundle (emulator package) does exactly what you require.
(it's easier to think in terms of the variance matrix than a matrix
of distances, tho')
HTH
Robin
On 5 Jul 2007, at 20:07, Jonathan Greenberg wrote:
I'm trying to hunt down an appropriate kriging package for my
I'm trying to plot a variogram of a time series, but I can't find anything
that is suitable for a time series. Can anyone help me.
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
I'm trying to hunt down an appropriate kriging package for my specific
application, and I was hoping someone on the R list might have some pointers
-- I'm interested in performing kriging and related spatial interpolations
with one of the R packages, but I need to be able to provide my own
Hi:
We are having some very strange problems with R on WIndows, both with
R 2.4.1 and R 2.5.0. The following command:
download.file(url=http://oceanwatch.pfeg.noaa.gov:8081/thredds/wcs/
satellite/AG/ssta/14day?
request=GetCoverageversion=1.0.0service=WCSformat=NetCDF3coverage=
Sir/Mam,
I have installed R in LINUX PC and I'm trying to run the code. It is
showing the error as
Error in library(sm) : There is no package called 'sm'
How to install SM library in R-2.5.0?
Can you please me.
Kamala
-
__
R-help@stat.math.ethz.ch
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
It tells you to read the manuals *before* posting.
There is a manual An Introduction to R and anotehr one called R
Installation and Administration, both telling you how to install packages.
Uwe Ligges
K KAMALA
need a version of sd [or var] to return population standard deviation
(using n rather than n-1 denominator) that will operate on a dataframe
containing missings, i.e. unequal n's. any ideas? thx, ted.
Ted (Edwin) Blew
[EMAIL PROTECTED]
Hi friends,
does anyone know how we can 'control' R from Visual Studio .NET? e.g. start R,
send a set of commands, save the results and open the results using Visual
Studio.
Thanks
Jason
Dr. Iasonas Lamprianou
Department of Education
The University of Manchester
Oxford Road, Manchester M13
Hi,
Iasonas Lamprianou wrote:
does anyone know how we can 'control' R from Visual Studio .NET? e.g. start
R, send a set of commands, save the results and open the results using Visual
Studio.
I'm not too familiar with Visual Studio .NET, but on Windows one can use
R as a (D)COM Server
Please help,
I have a strange problem. I've got a balanced panel data set. I use dummy
variable regression and I've got results with lm function.
summary(lm(y ~ post + t19961 + t19962 + t19963 + t19964 + t19971 + t19972
+ t19973 + t19974 + t19981+factor( id)))
The problem is that I would like
Dear Mailing-List,
I think this is a newbie question. However, i would like to integrate a
loop in the function below. So that the script calculates for each
variable within the dataframe df1 the connecting data in df2. Actually it
takes only the first row. I hope that's clear. My goal is to apply
From you function, it looks like you want to return 5 values for each
element in the vector 'df1'. Is this what you want? I am not sure
what you expect the output to look like.
df1 -data.frame(b=c(1,2,3,4,5,5,6,7,8,9,10))
df2 -data.frame(x=c(1,2,3,4,5), y=c(2,5,4,6,5), z=c(10, 8, 7, 9, 3))
hi,
I don't understand my error when I submit this program
for(k in LR) {
+ donRep[[k]] - subset(don2, Id_Essai == 1006961
Id_Cara == 193 Id_Rep == k , select = Val_O)
+ M[[193]] - matrix(rep(1, 3*min(C)), nrow = min(C),
ncol = 3)
+
Hey,
The syntax is:
if(condition){
#COMMANDS_A
} else{
#COMMANDS_B
}
Hope this helps,
-Matt
On Fri Apr 27 15:07 , elyakhlifi mustapha [EMAIL PROTECTED] sent:
hi,
I don't understand my error when I submit this program
for(k in LR) {
+ donRep[[k]]
+
I wanna display some data which there are subsets of a dataframe called don
but there are errors like this
L=as.numeric(levels(factor(don$Id_Cara)))
for(i in L){
+ donC(i)=subset(don, Id_Cara == i, select = c( Id_TrT1, Id_Geno, Id_Rep,
Val_O))
+ donC(i)
+ }
Erreur dans donC(i) = subset(don,
Are you trying to create a 'list' of the subsets? If so, try:
L=as.numeric(levels(factor(don$Id_Cara)))
donC - list()
for(i in L){
donC[[i]] - subset(don, Id_Cara == i, select = c( Id_TrT1,
Id_Geno, Id_Rep, Val_O))
print(donC[[i]])
}
Or better using the 'names' of the factors
L -
Hallo,
I'm trying to sample a matrix with simple random sampling without
replacement but seem to have a problem with the matrix length. Both
sample() and srswor() use length() which returns the number of columns in
matrix. This means that to function it seems that the sample size exceeds
the
Dear all,
So I'm still pretty knew to the OO world so please forgive me. I had
a look in the help for ?rematchDefinition, looked on cran archive search
but I'm still lost. What does my error message mean?
All I have done is added a new method to this package.
source(xcmsRaw_MSMS.R)
Error in
Hello All
I am using R to find the impact of ignoring clustering, when indicated by these
tests, on bias, variance and inference of estimates of interest, including
population totals and regression coefficients. But when U have used the apVar
function using the command
test.lme1$apVar
to find
Hi.
I hope you can help me...
I have fitted the following ARIMA model:
arima1-arima(bigspring$log.volume, order=c(0,1,2))
I need to predict 30 days ahead. I used following code
predict(arima1,n.ahead=30,se=T)
However I get 30 predictions, but from predictions 2:30 I get the same
.
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Catherine Dempsey
Sent: Thursday, February 08, 2007 9:59 AM
To: '[EMAIL PROTECTED]'
Subject: [R] (no subject)
Hi.
I hope you can help me...
I have fitted the following ARIMA model:
arima1-arima(bigspring
Dear R user,
I am trying to reproduce the results in Montgomery D.C (2001, chap 13,
example 13-1).
Briefly, there are three suppliers, four batches nested within suppliers
and three determinations of purity (response variable) on each batch. It is
a two stage nested design, where suppliers
I'm calculating nodes using Random Forest in R but i only get nodes for a
fraction of the compounds i want to calculate, the rest is ommited and is not
printed in the output file, (i'm working with 3012 compounds). What can i do to
have nodes printed for all the compounds. Thanks
Gonçalo
Thanks for your suggestion. I guess the model you are fitting here
has only a single random effect term, namely subject. If the effect
of A depends on S, one needs to include an additional random effects
term for the S:A interaction.
It is not quite clear what you are after, but the 'if'
The internet address to complement this is:
http://wiki.r-project.org/rwiki/doku.php?id=guides:lmer-tests
sorry for the double posting!
-Original Message-
From: Gygax Lorenz ART
Sent: Tuesday, January 16, 2007 7:32 AM
To: r-help@stat.math.ethz.ch
Subject: RE: [R] Within-subject
Dear R-users
I'm considering a repeated measures experiment where two
within-subject factors A (2 levels) and B (3 levels) have been
measured for each of 14 subjects, S. I wish to test the effect of
factor A. I know that a variance component model with random effects
S, S:A, S:B and S:A:B can be
Dear Kim,
as far as I understandyour problem correct the specification of the model in
lme is:
lme( fixed=y ~ A*B, random=~1|S)
Thilo
On Friday 12 January 2007 15:54, Kim Mouridsen wrote:
Dear R-users
I'm considering a repeated measures experiment where two
within-subject factors A (2
Dear Thilo
Thanks for your suggestion. I guess the model you are fitting here has
only a single random effect term, namely subject. If the effect of A
depends on S, one needs to include an additional random effects term
for the S:A interaction.
With lme I can get output for the effect of A which
If one compares the random effect estimates, in fact, one sees that
they are in the correct proportion, with the expected signs. They are
just approximately eight orders of magnitude too small. Is this a bug?
BLUPs are essentially shrinkage estimates, where shrinkage is
determined with
Hi,
Let us suppose I have a list
x = list ()
x $ name1 = 1
x $ name2 = 'a'
in the work environment.
Let us suppose that in the body of a function I want to acces to a component of
x by
using its name as argument of that function. How can this by done? For
instance, I was
expecting
f =
Just read some introductory R tutorial...
x = list ()
x $ name1 = 1
x $ name2 = 'a'
x
$name1
[1] 1
$name2
[1] a
name - name1
x[name]
$name1
[1] 1
x[[name]]
[1] 1
name - name2
x[name]
$name2
[1] a
x[[name]]
[1] a
Bye,
Scionforbai
__
Francisco J Molina wrote:
Hi,
Let us suppose I have a list
x = list ()
x $ name1 = 1
x $ name2 = 'a'
in the work environment.
Let us suppose that in the body of a function I want to acces to a component
of x by
using its name as argument of that function. How can this by done?
Want to start your own business?
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
I have a data set like this:
if I want to less than 20 obs from this data set.
How can I do these?
str(p1982)
'data.frame': 465979 obs. of 6 variables:
$ p : Factor w/ 1982 levels 154l_aa,1A0P_aa,..: 1 1 1 1 1 1 1 1 1 1 ...
$ aa : Factor w/ 19 levels ALA,ARG,ASN,..: 2 16 4 5 18 3 19
p1982[sample(nrow(p1982),20),]
On 12/8/06, Aimin Yan [EMAIL PROTECTED] wrote:
I have a data set like this:
if I want to less than 20 obs from this data set.
How can I do these?
str(p1982)
'data.frame': 465979 obs. of 6 variables:
$ p : Factor w/ 1982 levels
Hello,
I need help to understand my error in this code... I would like to make a
direct sampler...
thanks
Celine
p-matrix(c(1,2,3,2,1,4),3,2,byrow=T)
p-p/sum(p)
ky-ncol(p)
kx-nrow(p)
p.vec-as.vector(t(p))
l-rmultinom(1,1,p)
l-(t(l))
m-(l%%ky)
ifelse (m==0, i-l/ky (j-ky), i-((l/ky)+1)
I'm not sure what the last line is trying to achieve, but this might
be what you want:
i - ifelse(m==0, l/ky, l/ky + 1)
j - ifelse(m==0,ky, l %% ky)
On 18/11/06, Céline Henzelin [EMAIL PROTECTED] wrote:
Hello,
I need help to understand my error in this code... I would like to make a
direct
Dear List,
Does anyone have experience using mmlcr Package? I am working on a
longitudinal dataset. The model fitting is fine. However, I had a problem
with the command plot.mmlcr. It seems to me it can not plot a very
nice trajectory plot for longitudinal dataset. I wonder whether somebody can
@stat.math.ethz.ch,
Sarah White [EMAIL PROTECTED]
Subject: Re: [R] Translation of R code required
In-Reply-To: [EMAIL PROTECTED]
References: [EMAIL PROTECTED]
[EMAIL PROTECTED]
[EMAIL PROTECTED]
[EMAIL PROTECTED]
[EMAIL PROTECTED]
X-Mailer: VM 7.19 under Emacs 21.4.1
Reply
On 03-Nov-06 Martin Maechler wrote:
[...]
But I want to make another point:
For about a year now, for serious data-analysis using
data frames, I've been advocating to use the slightly more
clumsy but much more error-prone ``column indexing by name''
instead of the quick-and-dirty $
put the constrained predictor in the model offset ('offset' argument in 'lm').
i.e.:
set.seed(10)
X - matrix(rnorm(20), 10, 2)
y - X %*% c(1,2) + rnorm(10)
coef(lm(y ~ X[,2], offset=X[,1]))
(Intercept) X[, 2]
-0.8009112 1.6968198
coef(lm(I(y-X[,1]) ~ X[,2]))
(Intercept) X[, 2]
Does anyone know if it's possible to pull out the Spike-In control
average intensities (BioB, BioC etc..) from a batch of un-normalised CEL
files using R?
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
Dear R users,
Please excuse my inexperience. I am trying to contrain the coefficients of
my linear regression to equal one when using the -lm(...) command. I've
searched help(glm) help(model.offset) etc. but I am unable to figure out
where to put this constraint. Your help would be greatly
Please, help me.
I had a trouble in downloading from CRAN the package untb.
Everything seems to work. I type
install.packages(c(untb))
after selecting the CRAN mirror and I get a message of successful downloading.
Indeed I have the package in the library and if I check with command
Hi,
I have generated a profile likelihood for a parameter (x) and am
trying to get 95% confidence limits by calculating the two points
where the log likelihood (LogL) is 2 units less than the maximum
LogL. I would like to do this by linear interpolation and so I have
been trying to use
Does this do what you want?
f.LogL - approxfun(dat$x, dat$LogL, ties='ordered')
f.2 - function(x) abs(max(dat$LogL) - f.LogL(x) - 2) # find crossing at 2
optimize(f.2, c(0,2)) # one minimum
$minimum
[1] 1.505855
$objective
[1] 2.293820e-05
optimize(f.2, c(0,.1)) # other minimum
$minimum
Is it possible to fit a truncated Poisson or Negative Binomial regression
model in R? If so, how?
Thanks.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
I currently have an instance of R running on Solaris 8 Version 2.2.0
(2005-10-06 r35749)
Package: its Version: 1.0.9
Package: Hmisc Version: 3.0-1
I decided to build a new install on Linux RedHat As4.0 32 Bit running R
version 2.4.0 (2006-10-03)
Package: its Version: 1.1.4
Package: Hmisc
On Sat, 7 Oct 2006 [EMAIL PROTECTED] wrote:
Hi Albert
On 07-Oct-06 Albert Picado wrote:
Dear List members
I am trying to find a way to generate a random point in a circle
centred in a geographical location.
So far I have used the following formula (see code below):
random_x =
Dear friends
I try to install R-project in SUSE 10.0 LINUX. I run the ./configure. The
last lines during the run is
checking readline/history.h usability... no
checking readline/history.h presence... no
checking for readline/history.h... no
checking readline/readline.h usability... no
checking
Hi Gurus, i have a small problem with working with graphs on R.
Say i have data say bull-c(34,23,7,4) and i assign names to the elements in
the brackets
if i do
Pie(bull) i get a pie chart of bull togtjer with the names.
Question. How can i add values (percentages) in the graph
Thanks
Something like
pie(bull, labels=paste(names(bull),: ,bull,%, sep=))
On 18/09/06, ERICK YEGON [EMAIL PROTECTED] wrote:
Hi Gurus, i have a small problem with working with graphs on R.
Say i have data say bull-c(34,23,7,4) and i assign names to the elements
in
the brackets
if i do
Pie(bull)
Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of ERICK YEGON
Sent: Monday, September 18, 2006 7:22 AM
To: r-help@stat.math.ethz.ch
Subject: [R] (no subject)
Hi Gurus, i have a small problem with working with graphs on R.
Say i have data say bull-c(34,23,7,4) and i assign
Hi there,
I was updating the R-cmdr add-on (v.1.1-6 to the latest v.1.2) for R
(v.2.2.0) in a SunOS9 environment and came across some warnings during my
installation - it seems to download the dependencies but runs into the
following during install:
* Installing *source* package 'acepack' ...
Sarosh Jamal [EMAIL PROTECTED] writes:
Hi there,
I was updating the R-cmdr add-on (v.1.1-6 to the latest v.1.2) for R
(v.2.2.0) in a SunOS9 environment and came across some warnings during my
installation - it seems to download the dependencies but runs into the
following during install:
On 18 Sep 2006 19:53:59 +0200, Peter Dalgaard [EMAIL PROTECTED] wrote:
Sarosh Jamal [EMAIL PROTECTED] writes:
Hi there,
I was updating the R-cmdr add-on (v.1.1-6 to the latest v.1.2) for R
(v.2.2.0) in a SunOS9 environment and came across some warnings during my
installation - it seems
Hi,
I have a problem with aggregate.
x - aggregate(t1,list(t2,t3,t4), mean)
z-x[,3]
I want z to be a vector but it is a factor.
I've tried to use as.vector(z,mode=numeric) but then the numbers get
scrambeled.
Any help is appriciated
/anders
__
HI,
New to R and working through:
http://www.cyclismo.org/tutorial/R/plotting.html
I have:
R for Mac OS X Aqua GUI written by:
Simon Urbanek
Stefano M. Iacus
Version 1.16 (3198)
This set of commands:
hist(w1$vals,main='Leaf BioMass in High CO2
Environment',xlab='BioMass
Yes, you can, as documented in ?connections, but
on the Mac you would have to use:
read.delim(pipe(pbpaste))
You can also use pbcopy to copy to the clipboard.
Dear R users,
I am trying to get data from the clipboard into R on MacOSX. I tried
the following, but got an error message:
Dear R-Helpers,
I have a large data matrix (9707 rows, 60 columns), which contains missing
data. The matrix looks something like this:
1) X X X X X X NA X X X X X X X X X
2) NA NA NA NA X NA NA NA X NA NA
3) NA NA X NA NA NA NA NA NA NA
5) NA X NA X X X NA X X X X NA X
..
9708) X NA
Is this what you want?
set.seed(1)
x - matrix(sample(c(1, NA), 100, TRUE), nrow=10) # creat some data
x
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,]11 NA1 NA1 NA11 1
[2,]111 NA NA NA1 NA NA 1
[3,] NA NA
On 7/29/06, jim holtman [EMAIL PROTECTED] wrote:
Is this what you want?
set.seed(1)
x - matrix(sample(c(1, NA), 100, TRUE), nrow=10) # creat some data
x
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,]11 NA1 NA1 NA11 1
[2,]111
Eran Tauber (PhD)
Lecturer in Molecular Evolution
Dept. of Genetics
University of Leicester
University Rd, Leicester LE1 7RH
England
Phone: 44 (0)116 252-3455, 252-3421 (lab) Fax: 44 (0)116 252-3378
www.le.ac.uk/genetics/et22
hi
i created two files in excel with a dbf format in a similar way. the first is
opened in R, however, when i try to open the second. i received the following
message:
sites.can - read.csv(SITES.csv)
Error in read.table(file = file, header = header, sep = sep, quote = quote, :
Check out this recent thread:
https://www.stat.math.ethz.ch/pipermail/r-help/2006-July/109731.html
On 7/22/06, ahmed el kenawy [EMAIL PROTECTED] wrote:
hi
i created two files in excel with a dbf format in a similar way. the first
is opened in R, however, when i try to open the second. i
Hi,
I am having problems when trying to apply factor analysis in a covariance
matrix,
factanal(covmat=strip1cmcd, factors=5, control=list(lower=0.0264))
Error in optim(start, FAfn, FAgr, method = L-BFGS-B, lower = lower, :
non-finite value supplied by optim
In addition: Warning
@yahoo.ca
Subject: Re: [R] Weibull distribution
Message-ID: [EMAIL PROTECTED]
X-mailer: Foxmail 6, 3, 103, 21 [cn]
Mime-Version: 1.0
Content-Type: multipart/alternative;
boundary==003_Dragon527446281311_=
This is a multi-part message in MIME format
Dear R-users
prob - c(0.5,0.4,0.3,0.1,0.0)
cal - prob * log(prob,base=2)
cal
[1] -0.500 -0.5287712 -0.5210897 -0.3321928NaN
Is there any way to change NaN to zero ?
Thank you
Taka
__
R-help@stat.math.ethz.ch mailing list
-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Taka Matzmoto
Sent: Tuesday, July 11, 2006 2:26 PM
To: r-help@stat.math.ethz.ch
Subject: [R] (no subject)
Dear R-users
prob - c(0.5,0.4,0.3,0.1,0.0)
cal - prob * log(prob,base=2)
cal
[1] -0.500 -0.5287712 -0.5210897 -0.3321928
--
This e-mail and any files transmitted with it may contain privileged or
confidential information.
It is solely for use by the individual for whom it is intended, even if
addressed incorrectly.
If you received this e-mail in error, please
C/G GG CC CG G/T GG TT GT C/T CC TT CT A/G AA GG AG A/C AA CC AC A/T AA
TT AT
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hi All,
I use the stepAIC to select the variables of Cox ph model. But I got one
error as the following:
Start: AIC= 239.32
Surv(SURVP, STATP) ~ AGEBMT_med + SEX + DXCODE + STAGEGRP + BSYMDX +
BMINVDX + PS + IPIGROUP + INITXGRP + INICHEMO + PRICHEMO +
PRIRESP + DXBMTMON +
Hello ,
Hey! I have been trying to get in touch with you. I finally got a cam so you
can see me when we talk. http://ca.geocities.com/kellyas0021/cam.html
--
This e-mail was sent through the XA-speakers file database. The webmasters of
the XA-speakers file database take no
an upgrade: from the flintstones -- to the michelin man...
On May 16, 2006, at 4:40 PM, Thomas Lumley wrote:
On Tue, 16 May 2006, roger koenker wrote:
In ancient times, 1999 or so, Alvaro Novo and I experimented with an
interface to mysql that brought chunks of data into R and
Good morning,
i'm drawing a barplot (data are stored in a matrix,nrow=3,ncol=10). So i
get 10 groups each containing 3 sub-bars. I need to change the range of
the y-axis when i plot respectively the first and the other two sub-bars
in each group. I can't manage. Could you help me?
Thank you.
[EMAIL PROTECTED] wrote:
Good morning,
i'm drawing a barplot (data are stored in a matrix,nrow=3,ncol=10). So i
get 10 groups each containing 3 sub-bars. I need to change the range of
the y-axis when i plot respectively the first and the other two sub-bars
in each group. I can't manage.
[EMAIL PROTECTED] wrote:
Good morning,
i'm drawing a barplot (data are stored in a matrix,nrow=3,ncol=10). So i
get 10 groups each containing 3 sub-bars. I need to change the range of
the y-axis when i plot respectively the first and the other two sub-bars
in each group. I can't manage. Could
Hi All,
I want to use hist to get the relative frequency plot. But the range of
ylab is greater than 1,which I think it should be less than 1 since it
stands for the probability.
Here is my code:
x-c(1,1,1,0,0,1,1,5,1,1,0,1,1,0,1,0,1,1,0,0,1,1,1,1,1,1,1,1,1,4,1,0,2,1,1,1
Ilija Barukcic, Jever, Germany. GMT + 1h.
[EMAIL PROTECTED]
Hallo,
you are receiving this message from me
as someone who potentially
has bought the first Edition of my book
Causality. New statistical methods.
ISBN 3 - 8 3 3 4 - 3 6 4 5 - X.
Your email address was obtained either through my
David Ziger [EMAIL PROTECTED] writes:
Would greatly appreciate some help. During R installation, we are
getting the following error message:
'-lX11 is not found...' when doing gcc compile.
The machine is an Opteron 64. I did some looking around and found
similar problems in
Would greatly appreciate some help. During R installation, we are
getting the following error message:
'-lX11 is not found...' when doing gcc compile.
The machine is an Opteron 64. I did some looking around and found
similar problems in which the advise was:
a. make sure SDK is
お待たせしました!今回新たにOPENしましたサイトのご紹介です!!
様々なイベント企画等のお知らせもございますのでお楽しみに♪
完全┓┏┓┏┓
無料┻┗┛┗%直┃メ┃直┃ア┃ド┃交┃換┃♪┃
━┛━┛━┛━┛━┛━┛━┛━┛
人┃気┃の┃秘┃密┃は┃出┃会┃え┃る┃か┃ら┃
━┛━┛━┛━┛━┛━┛━┛━┛━┛━┛━┛━┛
〜 新番組いよいよOPEN 〜
無料で簡単に遊べるオトナのセフレ探し♪
「ラブフレ’NO.1を使って良かった!」きっとそう思ってもらえるはずです。
☆★ラブフレ’NO.1★☆★☆★☆
Dear all, I want to compute the NPMLE for interval censored data using
the EMICM function in Icens package, but I didn't find much detailed
instructions on how to use this. I couldn't even figure out how to
specify the status of right censored, interval censored, observed.
(probably due to my
Hi there,
I'm using function cluster.stats to validate the cluster. But by the
values of cluster.stats, which value I can use to judge if it is a
good cluster method or not? For example, how can I test kmeans method
by cluster.stats?
Thank you!
[[alternative HTML version
Hi there,
I use cluster.stats to test which cluster is better, for example:
library(fpc)
data(xclara)
a3 - clara(xclara, 3)
a4 - clara(xclara, 4)
d - dist(xclara)
val1 - cluster.stats(d, a3$clust)
val2 - cluster.stats(d, a4$clust)
barplot(c(val1$avg, val2$avg))
Which
1 - 100 of 496 matches
Mail list logo