Nagu wrote:
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
On 12/1/05, Arne Henningsen [EMAIL PROTECTED] wrote:
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
On Fri, 02 Dec 2005 10:38:21 +1000,
Ross Darnell (RD) wrote:
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.
You can always renew the Sinput and Soutput environments to
[EMAIL PROTECTED] wrote:
[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 +
I have 130 or so rows of questionnaire data. It is 67 columns wide
mostly of categorical nature (and often yes/no).
I have been using daisy, pam and plot to view cluster plots with subsets
of columns. With different subsets, there is great variation in the
degree of overlap of the
Jim == Jim Lemon [EMAIL PROTECTED]
on Thu, 01 Dec 2005 21:34:51 -0500 writes:
Jim 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
On Friday 02 December 2005 09:00, you wrote:
On 12/1/05, Arne Henningsen [EMAIL PROTECTED] wrote:
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.systemfi
t .
html
On Thu, 01 Dec 2005 11:34:01 -0600,
Qunfeng (Q) wrote:
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,
Hi
I totally agree with Martin because when I see boxplot I immediately
expect median in the middle and all other parts defined accordingly.
It is possible to use
bp - boxplot(..., plot=F)
and then to change the median values in bp to means and IQRs to SD
and everything to anything else but
I am using R in Windows. I see that I will have to use batch processes with
R. I will have to read and write text files, and run some R code; probably
some external code too. I have never done scripting. Is there any document
that explains simple steps with examples? I also have heard that Python
Hi
look at the following session, in which I have a dataframe,
and I want to extract the second row, without the first column.
Everything works as expected until the last line, where I set
the names of x to NULL, and get a non-desired object (I
want c(4,3).).
Three questions:
(1) why is
On 12/2/05, Arne Henningsen [EMAIL PROTECTED] wrote:
On Friday 02 December 2005 09:00, you wrote:
On 12/1/05, Arne Henningsen [EMAIL PROTECTED] wrote:
On Friday 02 December 2005 02:50, August Berg wrote:
I am confused by the following description in
RH == Robin Hankin [EMAIL PROTECTED] writes:
RH Three questions:
RH (1) why is as.vector(a[2,-1]) not a vector?
Did you read the help file of as.vector() ?
a - data.frame(male=c(T,T,F),mass=c(1,4,3),height=c(4,3,2))
a
male mass height
1 TRUE1 4
2 TRUE4 3
3
Evgeniy Kachalin wrote:
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)
On Friday 02 December 2005 09:58, you wrote:
On 12/2/05, Arne Henningsen [EMAIL PROTECTED] wrote:
On Friday 02 December 2005 09:00, you wrote:
On 12/1/05, Arne Henningsen [EMAIL PROTECTED] wrote:
On Friday 02 December 2005 02:50, August Berg wrote:
I am confused by the following
On 12/1/05, Arne Henningsen [EMAIL PROTECTED] wrote:
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
Hi
or maybe
x - as.numeric(a[2,-1])
can suit to your needs.
Petr
On 2 Dec 2005 at 17:14, Berwin A Turlach wrote:
From: Berwin A Turlach [EMAIL PROTECTED]
Date sent: Fri, 2 Dec 2005 17:14:42 +0800
To: Robin Hankin [EMAIL PROTECTED]
Dear list,
I'm currently trying to develop a model to assess clam yield potential in a
lagoon. I'm using the zeroinfl function of the pscl package to fit a
Zero-inflated negative binomial model, given the high occurrence of zero
counts.
I don't understand from the sentence in the pscl guide
On Friday 02 December 2005 10:21, August Berg wrote:
On 12/1/05, Arne Henningsen [EMAIL PROTECTED] wrote:
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.systemfi
t .
Hello,
the k-means {stats} help and the HartiganWon paper say nothing about
the way random sets works (parameter nstart). I would expect to get
the different results for each random initial set but I always obtain
only one result: how is it selected?
Charles Raux
Hi,
using mca function in library(MASS) I obtained the following result:
miacm=mca(factor.variables,abbrev=TRUE,nf=11)
miacm
Call:
mca(df = factor.variables, nf = 11, abbrev = TRUE)
Multiple correspondence analysis of 1000 cases of 3 factors
Correlations 0.605 0.599 0.586 0.577 0.577 0.577
Arne Henningsen wrote:
On Friday 02 December 2005 10:21, August Berg wrote:
On 12/1/05, Arne Henningsen [EMAIL PROTECTED] wrote:
On Friday 02 December 2005 02:50, August Berg wrote:
I am confused by the following description in
The key things here are:
1) a data frame is a list with some attributes (unlike a
matrix which is (generally) an atomic vector with attributes).
2) 'as.vector' removes attributes.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the
On 12/1/05, Arne Henningsen [EMAIL PROTECTED] wrote:
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
I have data from several sensors that recorded data at hourly intervals
during seven months. I want to separate daily variation from the trend,
and also be able to zoom in on only one month of data.
I have not been able what functions to use, I can not figure out from
the help for 'ts' how to
The source of an R package can be checked by the command
R CMD check package_name
(on a Unix/Linux machines; on Windows machines this is slightly
different).
In fact, it is exactly the same command on Windows.
Uwe Ligges
Why it did not work for me?
library( systemfit )
R CMD
On Fri, 2 Dec 2005, Charles Raux wrote:
the k-means {stats} help and the HartiganWon paper say nothing about
the way random sets works (parameter nstart). I would expect to get
the different results for each random initial set but I always obtain
only one result: how is it selected?
The code
The source of an R package can be checked by the command
R CMD check package_name
(on a Unix/Linux machines; on Windows machines this is slightly
different).
In fact, it is exactly the same command on Windows.
Uwe Ligges
Why it did not work for me?
library( systemfit )
R CMD
Hi,
I usually use lme from the nlme library. Now I have read an article about
lmer in Rnews and lmer seemed to me more comfortable to use. Unfortunately,
I didn't find out how to use covariance structures (e. g. corSymm(),
corAR1()). Is there a way to use them similarly as in lme ? Is it
August Berg wrote:
The source of an R package can be checked by the command
R CMD check package_name
(on a Unix/Linux machines; on Windows machines this is slightly
different).
In fact, it is exactly the same command on Windows.
Uwe Ligges
Why it did not work for me?
library(
R 2.1.1 on Windows 2K
I hope one and all will allow both an R question and a general stats question:
(1) Is there any function that will perform non-parametric tests such as a sign
test or a signed rank test? I know I could program both but I would prefer not
re-inventing the wheel.
(2)
Stephan:
It is my understanding that the corClass functions are not able to be
used with lmer() at this time.
Harold
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Stephan Moratti
Sent: Friday, December 02, 2005 6:58 AM
To: r-help@stat.math.ethz.ch
A lot is personal preference. I use Perl since it is very good at doing any
of the preprocessing that I need to setup data for R; e.g., reading a text
file, extracting the data and then setting it up in a format for R. I can
then invoke R from the script. Here is a sample of come Perl code that
On Fri, 02 Dec 2005 12:35:46 +0100 Henrik Andersson wrote:
I have data from several sensors that recorded data at hourly
intervals during seven months. I want to separate daily variation
from the trend, and also be able to zoom in on only one month of data.
I have not been able what
Read R-helpers,
I'm relatively new to R and trying to jump in feet first. I've been
able to learn a lot from on-line and printed documentation, but
here's one question to which I can't find an answer. Well, it's a
question with a couple parts. Thanks in advance for any direction
Hi,
Does anybody have a good tip of how to treat bimodal data to perform
statistical analyses? My data set ranges from -1 to 1 (any values are posssible
in between) and most data are either close to -1 or close to 1. They are the
results of a two choice experiment where individuals
Hi Simone,
It depends on what you want from the data ...
A bimodal distribution can typically be modelled with a mixture
distribution, where the assumption is that each of the modes corresponds to
an underlying process generating those data.
hth, ingmar
Hi,
Does anybody have a good tip of
Hi
On 2 Dec 2005 at 6:40, Jason Miller wrote:
To: r-help@stat.math.ethz.ch
From: Jason Miller [EMAIL PROTECTED]
Date sent: Fri, 2 Dec 2005 06:40:48 -0600
Subject:[R] plot and factors
Read R-helpers,
I'm relatively new to R
In addition to prior suggestion look at:
http://cran.r-project.org/doc/contrib/Ricci-refcard-ts.pdf
http://zoonek2.free.fr/UNIX/48_R/15.html
http://cran.r-project.org/src/contrib/Views/
On 12/2/05, Henrik Andersson [EMAIL PROTECTED] wrote:
I have data from several sensors that recorded data at
On 02-Dec-05 Simone Immler wrote:
Hi,
Does anybody have a good tip of how to treat bimodal data to
perform statistical analyses? My data set ranges from -1 to 1
(any values are posssible in between) and most data are either
close to -1 or close to 1. They are the results of a two choice
I do nearly all my scripting in R if its part of an R project. That keeps
things simple since everything is in one language. See ?BATCH,
?commandArgs, ?readLines, ?system, ?grep, ?regex (and commands
pointed to by it). If your script deals with XML, see the R XML package.
If R is not involved
See Wilcox.test and Wilcox_test functions.
Jarek Tuszynski
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of John Sorkin
Sent: Friday, December 02, 2005 7:13 AM
To: r-help@stat.math.ethz.ch
Subject: [R] sign and sign rank tests
R 2.1.1 on Windows 2K
I
Your problem sounds like it could be modeled with logistic regression
whereby the propensity for one result or another is linked to the factors
that control it. Logistic regressions are a special case of generalized
linear models. Look at ?glm
Charles Annis, P.E.
[EMAIL PROTECTED]
phone:
Hi!
First of all: I'm a newbie to both statistics and R, so please be
patient with me... I do however, like R because I've been programming
(pascal, IDL, perl, C etc) and designing models since -92, but never
related to statistics.
Ok, here we go:
I've got a set of 15 people, all of them
Dear All,
I am looking for function's code to implement the
Hosmer-Lemeshow gof test for survival data.
Could anyone share it?
TIA
Gioavnni
PS: I know that this test has low power, depends from
dr. Giovanni Parrinello
Section of Medical Statistics
Department of Biotecnologies
Viale
Hi guys,
Im struggling with this:
I have a dataset in which I have 9 measures and for each measure id like to
compute the kendall tau correlation using cor.test with
each of 20 covariates. Is there any systematic way to do this without writing
a for loop?
Thanks in advance,
Hello, R users,
I have a long list with over 20,000 elements (say foo), now I want to make a sub
list of 30 elements from this list using the component names such as subname-c(
a, b, d, f, g, h, o, q,...). One way to do it is to make rest of
elements in the list NULL, that is too silly. Is there
Dear R-users,
We expect to develop statistic procedures and environnement for the
computational analysis of our experimental datas. To provide a proof of
concept, we plan to implement a test for a given experiment.
Its design split data into 10 groups (including a control one) with 2
mesures
[EMAIL PROTECTED] writes:
Hello, R users,
I have a long list with over 20,000 elements (say foo), now I want to make a
sub
list of 30 elements from this list using the component names such as
subname-c(
a, b, d, f, g, h, o, q,...). One way to do it is to make rest
of
elements in the
2. Have you tried lmer in library(lme4)? The syntax will be
different, but it's a different algorithm and can handle problems that
crash lme. This may or may not apply to you. For documentation on
this, see Douglas Bates (2005) Fitting linear mixed models in R. R
News,
Using R itself has the advantage that you won't need to learn a new
language. You may even learn R better as a result.
Using Python could also work well. I have found it worth learning, at
least the basics. There a wider range of documentation (books)
available for Python than for R, and
I am looking at using R to analyze time series data containing a tidal
component. I need to remove the tidal signal to extract the time
series of the phenomena I seek to study. A browse of R-project search
engines has not been too fruitful? I've found 'hoa' and 'Rwave', but
need further
(Posted also to debian-user)
In my office network I have access to a debian powerpc server and 2 freebsd
6 servers (actually one of them is my notebook).
Experienced user of the statistical software R, I have now a go at parallel
computation via (r-)pvm and snow under R to enhance the performance
In the robust library of S+, it is possible to compute Rocke's
constrained M-estimator of scatter and location. Can't find anything
within R that does the same thing. Can someone tell me whether any package
exists for computing it?
Thanks.
Rand Wilcox
Professor
Dept of
Hi,
I am running R on my iBook (OS X 10.4) and whenever I try and plot a
model that has been run, I get a grey background. This makes the plot
difficult to read on bw printouts. I have tried changing the background to
white (or any other color) without any success. Here's a current example
On 12/2/05, Peter Rich [EMAIL PROTECTED] wrote:
Hi,
I am running R on my iBook (OS X 10.4) and whenever I try and plot a
model that has been run, I get a grey background. This makes the plot
difficult to read on bw printouts. I have tried changing the background
to
white (or any other
Dear group,
could some one help me how i can split a string into
pieces and access them.
For example:
test = this is a test
typeof(test)
[1] character
strsplit(test,split='is')
[[1]]
[1] th a test
Now I want to get test into : a vector and another
object as list
Also,
Rand == Rand Wilcox [EMAIL PROTECTED]
on Fri, 02 Dec 2005 10:42:48 -0800 (PST) writes:
Rand In the robust library of S+, it is possible to compute
Rand Rocke's constrained M-estimator of scatter and
Rand location. Can't find anything within R that does the
Rand same thing.
Dear R user,
My input data, positive and negative, is complete without missing data.
After running ksmooth( ) , I receive for $y , many NAs. What could be the
reasons and how can I receive complete output?
So many thanks for any help.
Amir Safari
Hello,
I am interfacing R and C code using the function .Call.
Is there any way to print to file form C?
I.e. is there something analogous to:
fopen(file1, w);
fprintf(file1, %f, ... );
fclose(file1);
which would not cause an error in the excution?
The Writing R Extension manual does not
As I get more familiar with R I tend to find the need for massaging data
with other scripts decreasing.
I still use python as a front end to some R tasks and it is a great
language to have in your personal arsenal (as is R).
The kind of decision making process for me goes something like:
o
On 02/12/05, paul sorenson [EMAIL PROTECTED] wrote:
I usually dump data from Python to R in CSV format. I call R scripts
from Python and about the only trick I use here is to read in an R
script template and perform string variable expansion (interpolation in
Perl) before sending it to an R
(Haven't seen an anwer to this yet; maybe I missed it.)
klebyn wrote:
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
Hello,
I am creating a 3D scatter plot with the cloud function. I ask it to
group my rows by a Treatment column (there are 4 different groups). I
need to know which points in the plot correspond to which rows (which
treatment). By default it assigns pchs and cols to each group. When I
try to
Dear Rand,
unfortunately, right now none of the R packages can compute Rocke's
constrained M-estimator of scatter and location. I have implemented this
estimator as function covMTB in rrcov - see
http://www.dst.unive.it/rsr/Todorov-treviso.pdf , slide 35 - but I have
still some open problems to
On 12/2/05, Michelle DePrenger-Levin [EMAIL PROTECTED] wrote:
Hello,
I am creating a 3D scatter plot with the cloud function. I ask it to
group my rows by a Treatment column (there are 4 different groups). I
need to know which points in the plot correspond to which rows (which
treatment). By
Date: Fri, 18 Nov 2005 13:26:42 +0200
From: [iso-8859-1] Pirjet? Antti [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Subject: [R] R-News 5/2, Bayesian Model Averaging, a detail
The article on BMA (Bayesian model averaging) presents most valuable tools for
model selection,
but I find one detail
I have not seen a reply to this, so I will offer a couple of feeble
comments, hoping that you might find them useful even though they don't
answer your question. It is difficult for me to understand your
question, because my SAS experience is so limited.
1. R is open
68 matches
Mail list logo