Eric == Eric C Jennings [EMAIL PROTECTED]
on Wed, 30 Nov 2005 21:39:12 -0800 writes:
Eric Hey there
Eric I have two vectors:
Eric y- c(0.4, 0.0, 0.2, -0.2, -0.6, 0.2, 0.0, 0.0, 0.4, 0.4, 0.2)
hmm, this is *one* vector , not two !
Eric In the vector y, I want to access
UweL == Uwe Ligges [EMAIL PROTECTED]
on Thu, 01 Dec 2005 08:37:37 +0100 writes:
UweL [EMAIL PROTECTED] wrote:
hi all,
I load a table with headers that enable me to acces it by the column
names:
tab-read.table(blob/data.dat,h=T)
attach(tab)
Martin Maechler wrote:
UweL == Uwe Ligges [EMAIL PROTECTED]
on Thu, 01 Dec 2005 08:37:37 +0100 writes:
UweL [EMAIL PROTECTED] wrote:
hi all,
I load a table with headers that enable me to acces it by the column
names:
tab-read.table(blob/data.dat,h=T)
Rick == Rick Bilonick [EMAIL PROTECTED]
on Wed, 30 Nov 2005 23:11:07 -0500 writes:
Rick I'm trying to fit a generalized mixed effects model to a data set
where
Rick each subject has paired categorical responses y (so I'm trying to use
a
Rick binomial logit link). There are
paul sorenson wrote:
I can't seem to track down R functions to calculate Internal Rate of
Return and NPV?
Thanks for the answers people. Comparing the answers with what Excel
pops out shows just how assumptions can vary. In particular whether the
first payment is at T0 or T1.
cheers
#
Let us say I have a time series I want to forecast. I have decided I am
going to use ppr, nnet and svm. I calibrate my model with these 3 algorithms
(let us assume I use all the data for the calibration; I do not distinguish
here between training and testing: all data are used for training, just
Gabor Grothendieck [EMAIL PROTECTED] writes:
What would be nice would be an R routine that automatically
implements this flowchart.
I'd recommend learning about p.adjust and the multcomp package
*instead* of following the flowchart.
Another approach to correcting for multiple testing is
Dear R-brains,
I'm rather new to state-space models and would benefit from the extra
confidence in using the excellent package sspir.
In a one-factor model, If I am trying to do a simple regression where
I assume the intercept is constant and the 'Beta' is changing, how do
I do that? How do i
On Thu, 2005-12-01 at 10:13 +, Prof Brian Ripley wrote:
On Thu, 1 Dec 2005, Martin Maechler wrote:
Rick == Rick Bilonick [EMAIL PROTECTED]
on Wed, 30 Nov 2005 23:11:07 -0500 writes:
Rick I'm trying to fit a generalized mixed effects model to a data set
where
Rick each
On 11/30/05, Scott Story [EMAIL PROTECTED] wrote:
I am trying to write a function that will solve a simple system of
nonlinear equations for the parameters that describe the beta
distribution (a,b) given the mean and variance.
mean = a/(a+b)
variance = (a*b)/(((a+b)2) * (a+b+1))
Any help as
YNoel == NOEL Yvonnick [EMAIL PROTECTED]
on Thu, 01 Dec 2005 14:42:44 +0100 writes:
YNoel On 11/30/05, Scott Story [EMAIL PROTECTED] wrote:
I am trying to write a function that will solve a simple system of
nonlinear equations for the parameters that describe the beta
Dear Sunil,
There are 7*8/2 = 28 raw moments among the 7 observed variables. Of these,
6*7/2 = 21 are used for the moments among the 6 fixed-exogenous variables,
leaving 28 - 21 = 7 df. You model has 11 free parameters. So df for the
model = 11 - 7 = -4.
Some additional comments:
If you're
Sorry perhaps I should have posted again. It wasnt a package error
instead it was my misunderstanding of what the package was asking for.
For future information the colAUC function requires a dataframe with at
least two datasets as well as an outcome set. It can be 'tricked' into
drawing an ROC
Hi,
I'm trying to get maximum likelihood estimates of \alpha, \beta_0 and
\beta_1, this can be achieved by solving the following three equations:
n / \alpha + \sum\limits_{i=1}^{n} ln(\psihat(i)) -
\sum\limits_{i=1}^{n} ( ln(x_i + \psihat(i)) ) = 0
\alpha \sum\limits_{i=1}^{n} 1/(psihat(i)) -
[EMAIL PROTECTED] wrote:
Hi,
I'm trying to get maximum likelihood estimates of \alpha, \beta_0 and
\beta_1, this can be achieved by solving the following three equations:
n / \alpha + \sum\limits_{i=1}^{n} ln(\psihat(i)) -
\sum\limits_{i=1}^{n} ( ln(x_i + \psihat(i)) ) = 0
\alpha
It is a bug, I will have to fix. I assumed that data will have column names
I can use for legend. So temporary fix is to provide column names:
library(caTools)
a-matrix(rnorm(100), 100,1)
b-rbinom(100,1,0.7)
colnames(a) = test;
colAUC(a,b,plotROC=TRUE)
Thanks for reporting it.
Jarek Tuszynski
Dear all,
Trying to produce 4 maImage plots (marray package) on the same device (2 on
the top and 2 on the bottom) with the layout() function or the split.screen()
function, we are facing the following problem:
it seems that maImage() does nt care about any of these 2 functions, and
plots only
Dear R-helpers,
I made the package depmix for fitting hidden markov models.
After updating R to version 2.2.0 on my PC I got the following error when
installing the package from a local directory. One of the C files includes
the R.h file which apparently causes problems.
I call
R CMD INSTALL
Ingmar Visser wrote:
Dear R-helpers,
I made the package depmix for fitting hidden markov models.
After updating R to version 2.2.0 on my PC I got the following error when
installing the package from a local directory. One of the C files includes
the R.h file which apparently causes problems.
hello,
is R able to suppress tick labels (not tick marks)? i know there is a way
around this with axes=F and then draw new axes, but it would be easier to
suppress them in the first place.
--
Sebastian Leuzinger
Institute of Botany, University
Hello to all users and wizards.
I am regulary using 'boxplot' function or its analogue - 'bwplot' from
the 'lattice' library. But they are, as far as I understand, totally
flawed in functionality: they miss ability to select what they would
draw 'in the middle' - median, mean. What the box
Dear Uwe,
I made the package depmix for fitting hidden markov models.
After updating R to version 2.2.0 on my PC I got the following error when
installing the package from a local directory. One of the C files includes
the R.h file which apparently causes problems.
I call
R CMD INSTALL
Boxplots were invented by John W. Tukey and I think should be
counted among the top small but smart achievements from the
20th century. Very wisely he did *not* use mean and standard deviations.
Even though it's possible to draw boxplots that are not boxplots
(and people only recently explained
I'm no wizard but looking at ?boxplot I think you should try ?bxp.
HTH,
Jean-Christophe.
2005/12/1, Evgeniy Kachalin [EMAIL PROTECTED]:
Hello to all users and wizards.
I am regulary using 'boxplot' function or its analogue - 'bwplot' from
the 'lattice' library. But they are, as far as I
All,
I've written some functions that use a list and a
list of sub-lists and I'm running into memory problems,
even after changing memory.limit. Does it make
any difference to the handling of memory if I use
simple vectors and matrices instead of the list and
list of sub-lists? I suspect
On Thu, 1 Dec 2005, Sebastian Leuzinger wrote:
is R able to suppress tick labels (not tick marks)? i
know there is a way around this with axes=F and then
draw new axes, but it would be easier to suppress them
in the first place.
Something wrong with setting them to null string?
Dear R-helpers,
When installing a source package I can suppress the compilation of .chm
files by using the --docs=normal option. Is it also possible to suppress
the creation and checking of .chm files when calling R CMD check ?
best, ingmar
__
on 12/1/2005 4:05 AM Martin Maechler said the following:
UweL == Uwe Ligges [EMAIL PROTECTED]
on Thu, 01 Dec 2005 08:37:37 +0100 writes:
UweL [EMAIL PROTECTED] wrote:
hi all,
I load a table with headers that enable me to acces it by the column
names:
Dear Tariq Khan
The initial conditions m0 and C0 can be specified according to your needs. If
you are a Bayesian (as in WestHarrison 1997), you will use m0 and C0 to
express your prior information. If you use a vague prior, you will give a high
weight to your observations in the beginning,
Martin Maechler пишет:
Boxplots were invented by John W. Tukey and I think should be
counted among the top small but smart achievements from the
20th century. Very wisely he did *not* use mean and standard deviations.
Even though it's possible to draw boxplots that are not boxplots
(and
Dear List:
I have created some code to simulate data from a complex sample where
5000 students are nested in 50 schools. My code returns a dataframe with
a variable representing student achievement at a single time point. My
actual code for creating this is below.
What I would like to do is
You asked about the behavior of KalmanSmooth(...)$var with missing
values. I also got very counterintuitive results from your excellent,
reproducible example using R 2.2.0 under Windows XP:
* With no missing values, KalmanSmooth(...)$var = 2, independent of
the data and
On 1 Dec 2005, at 17:58, Sebastian Leuzinger wrote:
hello,
is R able to suppress tick labels (not tick marks)? i know there is a
way
around this with axes=F and then draw new axes, but it would be easier
to
suppress them in the first place.
You mean the numbers below or beside each
On Thu, 2005-12-01 at 10:19 -0600, Paul Roebuck wrote:
On Thu, 1 Dec 2005, Sebastian Leuzinger wrote:
is R able to suppress tick labels (not tick marks)? i
know there is a way around this with axes=F and then
draw new axes, but it would be easier to suppress them
in the first place.
I'd like to add two comments to Martin's sensible response.
1. I've seen several intro-stats textbooks that define a
boxplot to have whiskers to the extreme data values
and then define Tukey's boxplot as a modified boxplot.
I wish authors wouldn't do that.
2. I've also seen boxplots used for
I think the question was about the __tick mark__ labels not the __axis__
labels. The standard way of dealing with this is, as Sebastian said, to draw
customized axes. However, par(lab=c(1,5,5)) reduces the number of tick mark
labels to 2 at the range of the axes, apparently, if that helps.
Sebastian Leuzinger wrote:
hello,
is R able to suppress tick labels (not tick marks)? i know there is a way
around this with axes=F and then draw new axes, but it would be easier to
suppress them in the first place.
Not really suppressing them, but you could you do the following:
On 12/1/05, Evgeniy Kachalin [EMAIL PROTECTED] wrote:
Martin Maechler пишет:
Boxplots were invented by John W. Tukey and I think should be
counted among the top small but smart achievements from the
20th century. Very wisely he did *not* use mean and standard deviations.
Even though
Morning,
I've downloaded the precompiled R 2.1.1 version and am using Windows XP
on my office workstation. As mentioned previously, I've resorted to batch
jobs to avoid the hanging that occurs when I try to plot the 3690 length
vector of data. If it's warranted, I can do a build from the source
Hi All,
I have two time series, each has length 354. I tried to calculate the
coherency^2 between them, but the value I got is always 1. On a website,
it says: Note that if the ensemble averaging were to be omitted, the
coherency (squared) would be 1, independent of the data. Does any of
1. To evalute the significance of the random variable (a random
effect?) using 'lmer', have you considered fitting models with and
without that effect, as in the example with 'example(lmer)'?
2. Regarding 'predict.lmer', I tried the following:
predict(fm1)
Error in
Hello,
I am new to the R package. After I use R to perform the hierarchical
clustering, I am only interested in retrieving the leaf nodes that share
the last common ancestors. As illustrated below, I'd like to retrieve (B,
C) as a cluster and then (D, E) as another cluster.Any chance to
On Thu, 2005-12-01 at 19:40 +0300, Evgeniy Kachalin wrote:
Martin Maechler пишет:
Boxplots were invented by John W. Tukey and I think should be
counted among the top small but smart achievements from the
20th century. Very wisely he did *not* use mean and standard deviations.
Even
On 12/1/2005 11:24 AM, Ingmar Visser wrote:
Dear R-helpers,
When installing a source package I can suppress the compilation of .chm
files by using the --docs=normal option. Is it also possible to suppress
the creation and checking of .chm files when calling R CMD check ?
I'd guess that using
All--
Would someone kindly post the reference to Tukey's formula for a boxplot
without whiskers?
I am looking at his book Exploratory Data Analysis from 1977. The
index includes box-and-whisker plot but not boxplot. On page 39-40
construction of the plot is described, including the
At office, using the internal LAN at my disposal, I'm having a go at parallel
computing - to begin with - with pvm, rpvm snow.
The two boxes are as follows
Remote machine uffbsd:
CPU: Intel(R) Pentium(R) 4 CPU 2.00GHz (1994.13-MHz 686-class CPU)
Origin = GenuineIntel Id = 0xf24 Stepping =
[EMAIL PROTECTED] wrote:
Hi,
I'm trying to get maximum likelihood estimates of \alpha, \beta_0 and
\beta_1, this can be achieved by solving the following three equations:
n / \alpha + \sum\limits_{i=1}^{n} ln(\psihat(i)) -
\sum\limits_{i=1}^{n} ( ln(x_i + \psihat(i)) ) = 0
\alpha
On Thu, 1 Dec 2005, Marc Schwartz (via MN) wrote:
On Thu, 2005-12-01 at 10:19 -0600, Paul Roebuck wrote:
On Thu, 1 Dec 2005, Sebastian Leuzinger wrote:
is R able to suppress tick labels (not tick marks)? i
know there is a way around this with axes=F and then
draw new axes, but it would be
On Thu, 1 Dec 2005, Ingmar Visser wrote:
Dear R-helpers,
When installing a source package I can suppress the compilation of .chm
files by using the --docs=normal option. Is it also possible to suppress
the creation and checking of .chm files when calling R CMD check ?
Part of the check is
On Thu, 2005-12-01 at 18:33 +, Prof Brian Ripley wrote:
On Thu, 1 Dec 2005, Marc Schwartz (via MN) wrote:
On Thu, 2005-12-01 at 10:19 -0600, Paul Roebuck wrote:
On Thu, 1 Dec 2005, Sebastian Leuzinger wrote:
is R able to suppress tick labels (not tick marks)? i
know there is a way
It's already there in EDA! Se pp 39-47.
-- 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. Box
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL
In the example at
http://www.stat.uiowa.edu/~luke/R/cluster/cluster.html clusterCall
divides the parameter R by the number of nodes -- what you've done is
calculate 999 bootstraps on each node, and compared the execution time
to 999 bootstraps on one node.
You probably want 'clusterCall' to be
Marc Schwartz (via MN) пишет:
On Thu, 2005-12-01 at 19:40 +0300, Evgeniy Kachalin wrote:
Martin Maechler пишет:
So I analize genetics data. I have some factor (gene variant, c(1,2,3))
and the quantitative variable corresponding to that factor. How do I
visualize this situation? Compare mean of
Alle 19:03, giovedì 01 dicembre 2005, Martin Morgan ha scritto:
In the example at
http://www.stat.uiowa.edu/~luke/R/cluster/cluster.html clusterCall
divides the parameter R by the number of nodes -- what you've done is
calculate 999 bootstraps on each node, and compared the execution time
to
Hi
Florence Combes wrote:
Dear all,
Trying to produce 4 maImage plots (marray package) on the same device (2 on
the top and 2 on the bottom) with the layout() function or the split.screen()
function, we are facing the following problem:
it seems that maImage() does nt care about any of
Marc Schwartz (via MN) пишет:
On Thu, 2005-12-01 at 19:40 +0300, Evgeniy Kachalin wrote:
Martin Maechler пишет:
So I analize genetics data. I have some factor (gene variant, c(1,2,3))
and the quantitative variable corresponding to that factor. How do I
visualize this situation?
Marc Schwartz (via MN) пишет:
Marc Schwartz (via MN) пишет:
So plotmeans is incapable of: boxplot(numerical~fact1+fact2). Is there
any way further?
I think that somehow we are talking past each other here.
plotmeans() does what it is designed to do, which is to simplify the
process of
interaction(A, B) will create a single factor made up of the combinations of
the two factors A and B. Perhaps that would let you use plotmeans.
Hope this helps,
Matt Wiener
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Evgeniy Kachalin
Sent:
Evgeniy Kachalin wrote:
Marc Schwartz (via MN) пишет:
Marc Schwartz (via MN) пишет:
So plotmeans is incapable of: boxplot(numerical~fact1+fact2). Is there
any way further?
I think that somehow we are talking past each other here.
plotmeans() does what it is designed to do, which is to
On Thu, 2005-12-01 at 23:27 +0300, Evgeniy Kachalin wrote:
Marc Schwartz (via MN) пишет:
Marc Schwartz (via MN) пишет:
So plotmeans is incapable of: boxplot(numerical~fact1+fact2). Is there
any way further?
I think that somehow we are talking past each other here.
plotmeans()
Dear List,
This is my first post, and I'm a relatively new R user trying to work out a
mixed effects model using lme() with random effects, and a correlation
structure, and have looked over the archives, R help on lme, corClasses,
etc extensively for clues. My programming experience is minimal
I have a data frame which has the following data.
data-read.table(table.txt,header=TRUE)
data
X14A_U133A_StatPairs X14A_U133A_Detection X14B_U133A_Signal
1 AFFX-BioB-5_at403.0 409.3
2 AFFX-BioB-M_at757.3 574.4
3
Frank E Harrell Jr пишет:
Evgeniy Kachalin wrote:
Marc Schwartz (via MN) пишет:
Marc Schwartz (via MN) пишет:
library(Hmisc)
library(lattice)
?panel.bpplot
bwplot(, panel=panel.bpplot)
By default, panel.bpplot shows the mean (dot) and median (line) plus
several quantiles.
Wiener, Matthew пишет:
interaction(A, B) will create a single factor made up of the combinations of
the two factors A and B. Perhaps that would let you use plotmeans.
Hope this helps,
Matt Wiener
So you think plotmeans(num~interaction(A,B)) will work? How? There is NO
'num' data for a.d,
Hello
I do not know very much about statistics (and English language too :-( ),
then I come in search of a clarification (explanation):
I found two distinct results on KURTOSIS and
I do not know which of them is the correct one.
Any aid will be welcome!
klebyn
CODE
Vasu,
On Dec 1, 2005, at 5:20 PM, Vasundhara Akkineni wrote:
i want to extract the first column of names as a matrix and pass it
to a String[] type in java. I am using RServe and was able to pass
the other two columns to double[] in java but was not able to do so
for data[,1] to a
P Ehlers wrote:
I'd like to add two comments to Martin's sensible response.
1. I've seen several intro-stats textbooks that define a
boxplot to have whiskers to the extreme data values
and then define Tukey's boxplot as a modified boxplot.
I wish authors wouldn't do that.
2. I've also
Ed Wang wrote:
Morning,
I've downloaded the precompiled R 2.1.1 version and am using Windows XP
on my office workstation. As mentioned previously, I've resorted to batch
jobs to avoid the hanging that occurs when I try to plot the 3690 length
vector of data. If it's warranted, I can do a
Kjetil Brinchmann Halvorsen wrote:
P Ehlers wrote:
I'd like to add two comments to Martin's sensible response.
1. I've seen several intro-stats textbooks that define a
boxplot to have whiskers to the extreme data values
and then define Tukey's boxplot as a modified boxplot.
I wish
Check your syntax on the bwplot call.
fa - data.frame(doz=sample(500:2000, size=500), fabp2=rep(1:20, 25))
bwplot(factor(fabp2) ~ doz, data=fa, panel=panel.bpplot)
fa.sum - summarize( fa$doz, list( fabp2 = fa$fabp2), smean.sd,
stat.name=doz)
Dotplot( factor(fabp2) ~ Cbind(doz, doz - SD, doz
Sweavers
As the title suggests I would appreciate any help to include S code in
ordinary paragraph mode. I can use the textsl font but is isn't the
same.
Thanks
Ross Darnell
__
R-help@stat.math.ethz.ch mailing list
If I understand correctly, you are asking how solve.QP interprets the
arguments. I clarified this by constructing a problem I could work
manually but that was sophisticated enough to seemingly answer your
question: minimize (x1^2+x2^2) subject to (x1+x2)=1. I computed
manually
I am confused by the following description in
http://www.maths.lth.se/help/R/.R/library/systemfit/html/hausman.systemfit.html
what does the Not run mean? if we do not load systemfit, how can we run
the following code?
## Not run: library( systemfit )
data( kmenta )
attach( kmenta )
...
I
1. Since your example is incomplete, I can't easily replicate the
phenomenon. This limits me to general comments and brainstorming on
things I might try.
2. Have you tried lmer in library(lme4)? The syntax will be
different, but it's a different algorithm and can handle
Hi,
Problem statement:
There are two variables Y, X. Y is a response from X. I want to find
a closed form formula to express Y in terms of X.
Comments:
I tried with simple polynomial basis but its not a good fit due to the
non-linear nature of Y.
Are there any set of functions or algorithm to
Dear R-helper,
Is there a command to get an array indexed 1 to T from T to 1?
For example:
a - c(1, 2, 3)
and by applying such a command I can get
a[1] = 3
a[2] = 2
a[3] = 1
Thanks a lot and best regards
Julio
Please use R's existing help system before posting to the list.
help.search('reverse') is an obvious first thing to try, don't you think
(and gives you the almost obvious answer immediately)?
-- Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED]
Austin, Matt пишет:
Check your syntax on the bwplot call.
fa - data.frame(doz=sample(500:2000, size=500), fabp2=rep(1:20, 25))
bwplot(factor(fabp2) ~ doz, data=fa, panel=panel.bpplot)
Yes, that's almost the same But there is a huge amount of data on
the graphic, too much for
{diverted back to R-help}
There are several R packages that provide plots of
mean +/- SD (or mean +/- 2*SD which is an approximate 95%
confidence interval for the case of normally distributed data)
or so called error bars.
E.g. function plotCI() in package 'gplots' and errbar() in
package
On Friday 02 December 2005 02:50, August Berg wrote:
I am confused by the following description in
http://www.maths.lth.se/help/R/.R/library/systemfit/html/hausman.systemfit.
html
what does the Not run mean?
This means that the command library( systemfit ) is not executed when the
package
Hi,
[...] the command/script window does not appear. [...]
Have you tried
Misc = Toolbars = Command
on the menu? Don't click until you get to Command, otherwise the menu closes
again for some reason.
Regards,
Andreas
- Original Message -
From: Henrik Parn [EMAIL PROTECTED]
Eric C. Jennings wrote:
Hey there
I have two vectors:
y- c(0.4, 0.0, 0.2, -0.2, -0.6, 0.2, 0.0, 0.0, 0.4, 0.4, 0.2)
In the vector y, I want to access (in the order given) all of the values in
between each of the specific values of given.
I understand subsetting with y[i], but how
anil kumar rohilla wrote:
HI LIST
i am a new R user,i am trying to make a model,which will give me
output in probability,which will take predictors and predictand serie as
input and and give me output in terms of probability(e.g below
normal,normal,above normal etc.).What is
Eric C. Jennings wrote:
To start with, pardon my mistake regarding two vectors
Yes I want all values (to two decimal digits if I can) between each of
the values given in the vector.
meaning I'm really try to do something like this:
y- c(0.4 : 0.0 : 0.2 : -0.2 : 0.6 : 0.2 : 0.0 : 0.0 : 0.4
84 matches
Mail list logo