On Mon, 2005-11-21 at 18:16 +0100, Matthieu Cornec wrote:
> Hello,
> How do you change the size of the caracters (tiny, small) using xtable?
> It works out for print.xtable when typing
> print.xtable(xtable(mydata),size="small")
> but I do not see any results when doing
> xtable(mydata,size="small
On Saturday 19 November 2005 22:09, Gabor Grothendieck wrote:
> Getting back to your original question of using apply, solving the LP
> gives us the number of components in any minimal solution and
> exhaustive search of all solutions with that many components can
> be done using combinations from
On 11/19/05, Adrian DUSA <[EMAIL PROTECTED]> wrote:
> On Saturday 19 November 2005 19:17, Patrick Burns wrote:
> > [snip...] One cheat would be to do the LP problem
> > multiple times with the rows of your matrix randomly
> > permuted. Assuming you keep track of the real rows,
> > you could th
On Saturday 19 November 2005 20:51, Ted Harding wrote:
> [..snip...]
> There is bound to be a good algorithm out there somewhere
> for finding a "minimal coveriung set" but I don't know it!
> Best wishes to all,
> Ted.
I found this presentation very explicit:
http://www.cs.ualberta.ca/~amaral/cour
Dear Ted,
On Saturday 19 November 2005 20:51, Ted Harding wrote:
> [...snip...]
> There is bound to be a good algorithm out there somewhere
> for finding a "minimal coveriung set" but I don't know it!
>
> Comments?
>
> Best wishes to all,
> Ted.
My case is probably a subset of your general algori
On 19-Nov-05 Adrian Dusa wrote:
> On Saturday 19 November 2005 17:24, Gabor Grothendieck wrote:
>> [...snip...]
>> Although the above is not wrong I should have removed the
>> rbind which is no longer needed and simplifying it further,
>> as it seems that lp will do the rep for you itself for
>> ce
On Saturday 19 November 2005 19:17, Patrick Burns wrote:
> [snip...] One cheat would be to do the LP problem
> multiple times with the rows of your matrix randomly
> permuted. Assuming you keep track of the real rows,
> you could then get a sense of how many solutions there
> might be.
Thanks
I suspect that the answer is that finding all solutions
will be hard. L1 regression is a special case of LP.
I learned how to move around the corners of the
solution space, and could easily find all of the solutions
in the special case of a two-way table. However,
sometimes there were a lot of so
On Saturday 19 November 2005 17:24, Gabor Grothendieck wrote:
> [...snip...]
> Although the above is not wrong I should have removed the rbind
> which is no longer needed and simplifying it further, as it seems
> that lp will do the rep for you itself for certain arguments, gives:
>
> lp("min", rep
On 11/19/05, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> On 11/19/05, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> > Try minizing 1'x subject to 1 >= x >= 0 and m'x >= 1 where m is your mtrx
> > and ' means transpose. It seems to give an integer solution, 1 0 1,
> > with linear programming
On 11/19/05, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> Try minizing 1'x subject to 1 >= x >= 0 and m'x >= 1 where m is your mtrx
> and ' means transpose. It seems to give an integer solution, 1 0 1,
> with linear programming even in the absence of explicit integer
> constraints:
>
> library(
Try minizing 1'x subject to 1 >= x >= 0 and m'x >= 1 where m is your mtrx
and ' means transpose. It seems to give an integer solution, 1 0 1,
with linear programming even in the absence of explicit integer
constraints:
library(lpSolve)
lp("min", rep(1,3), rbind(t(mtrx), diag(3)), rep(c(">=", "<="
On 11/19/2005 8:00 AM, Adrian DUSA wrote:
> Dear list,
>
> I have a problem with a toy example:
> mtrx <- matrix(c(1,1,0,1,1,1,0,1,1,0,0,1), nrow=3)
> rownames(ma) <- letters[1:3]
>
> I would like to determine which is the minimum combination of rows that
> "covers" all columns with at least a 1
Wen Zhang wrote:
> I have three files of data which are available at http://zhangw.com/
> R/, varied at the number of data. I tried to use R to analyze using
> shapiro.test, ks.test, and t.test. t.test ran as expected, however,
> when I run shapiro.test and ks.test commands, error message alw
On Thu, 2005-11-17 at 17:09 +0100, Matthieu Cornec wrote:
> Hello,
> I am importing the following file
> ;aa;bb;cc
> 1988;12;12;12
> 1989;78;78;12
> 1990;78;78;12
> 1991;78;78;12
> 1992;78;78;12
> 1993;78;78;12
> 1994;78;78;12
>
> data<-read.csv2(
Please do not post to both BioConductor and R.
On Thu, 2005-11-10 at 09:51 -0700, Nayeem Quayum wrote:
> Hello everybody,
> I am trying to use mas5 to normalize some array data and using mas5 and
> mas5calls. But I received these warning message. If anybody can explain the
> problem I would real
Luciana Alves asked
> Dear Sirs,
> Could you please be so kind as to send us some information on
residuals in
> multinomial logistic models?
Peter L. Flom, PhD
Assistant Director, Statistics and Data Analysis Core
Center for Drug Use and HIV Research
National Development and Research Institutes
hi, hi all,
> Dear Sirs,
> Could you please be so kind as to send us some information on residuals in
> multinomial logistic models?
here are some references to multinomial models:
Agresti A.(1996) An Introduction to Categorical Data Analysis
Agresti A. (2002) Categorical Data Analysis, 2nd Ed
1) There is no 'R 2.0'. What version did you mean?
2) We cannot reproduce your script (no data files), and JPEGs are not
allowed on R-help: see http://www.r-project.org/mail.html. (PNGs are,
though).
3) You have given us no indication of what the problems are nor in which
lines.
Please give
David Zhao wrote:
> Hi there,
>
>
> Could somebody help me disect this legacy R script I inherited at work, I
> have two questions:
> 1. I've tried to upgrade our R version from 1.6.2 (yeah, I know), to R 2.0,
> but some of the lines in this script are not compatible with R 2.0, could
> someone
I don't think AICc is available in a standard script. A couple of
years ago, I modified the stepAIC script from the MASS package to use
AICc. If and when I get time for this again, I plan to redo it using a
more complete Bayesian analysis.
However, if you want to use AICc,
On Mon, 7 Nov 2005, Nelson, Gary (FWE) wrote:
> I am wondering if anyone might know why the package clim.pact won't
> install properly. I have tried many URL sites and the same thing
> happens. I get the error message below. I also tried downloading the
> ZIP from the CRAN site and extracting th
1. Have you considered specifying the random effect as a list? This
is mentioned in the help file and discussed in Pinheiro and Bates (2000)
Mixed-Effects Models in S and S-Plus (Springer). On p. 40, they discuss
the following example that might help you:
lme(pixel~day+day^2, data=
Hi Sean,
I am sorry for that. I am using the following code for
doing MDS:
groups = list( Day1C=c(1), Day1T=c(2,3,4),
Day2C=c(5,6,7,8,9,10,11,12), Day2T=c(14,15,16,17,18))
len = length(groups)
dist = dist(data, method="euclidean")
data.mds=isoMDS(dist)
plot(data.mds$points,type = "n",main="MDS p
On 11/2/05 6:42 AM, "A Ezhil" <[EMAIL PROTECTED]> wrote:
> Hi All,
>
> I am trying to apply MDS for 4 groups in my data. The
> groups are:
>
> groups = list( Day1C=c(9), Day1T=c(7,8,10),
> Day2C=c(1,2,3,6,11,13,14,15), Day2T=c(4,5,12,16,17,18)
> )
>
> When I do the MDS plot the group1 appears
If I understand you correctly, you cannot do this with a data.frame, which
must be rectangular with equal numbers of entries (columns) in each row. See
?.data.frame and read "An Introduction to R" for these basics. You could
make the 3rd column for exhaustive = NA (or maybe an empty string, "") I
There are several bootstrap packages available at CRAN that probably
provide an elegant solution, but simply permuting the matrix is pretty easy
> data <- matrix(1:100,nrow=5) # matrix of 5 rows and 20 columns
> x <- data[sample(1:5),] # permute the rows
> y <- x[,sample(1:20)] # permute the c
> "booop" == booop booop <[EMAIL PROTECTED]>
> on Wed, 2 Nov 2005 13:22:55 -0800 (PST) writes:
booop> Could anybody help me in mixing the matrix values
booop> randomly (first rows and then columns)... for eg
booop> putting the 1 st row in the place of 3 rd and 3rd in
b
Hi
although I do not know anything about hier.part package I try few
comments
- see posting guide as it suggest to try to present a toy example
which shows your problem
- are there some error messages or the result is not as you expect?
- what is TEMP - seems to me that you need to define it
If you wish to apply labels to identified points: I am sure there is a better
way, but what I have done is use
locator(n) #n is how many points you wish to identify. Then click on the
points. The coordinates of each point are returned.
Then I use
identify(x, y, labels = My Label)
...for each
Hi,
try environment(Fn12)$y
environment(Fn12)$x
Nolwenn
**
Nolwenn Le Meur, PhD
Fred Hutchinson Cancer Research Center
Computational Biology
1100 Fairview Ave. N., M2-B876
P.O. Box 19024
Seattle, WA 98109-1024
On Tue, 1 Nov 2005, klebyn wrote:
>
>
> Hi
Just a bit of nitpicking: I believe the preferred way is to see if the
result of try() _inherits_ the try-error class. This applies to all S3
classes. I.e., the relevant line should be something like:
if (inherits(answer, "try-error")) ...
Andy
> From: David Scott
>
> On Mon, 31 Oct 2005
On Mon, 31 Oct 2005, McGehee, Robert wrote:
> It sounds like you want `try` with the argument `silent = TRUE`. This
> will allow you to keep running your program without errors. If you want
> to check if the line had an error, you can error control by seeing if
> the class of the resulting object
It sounds like you want `try` with the argument `silent = TRUE`. This
will allow you to keep running your program without errors. If you want
to check if the line had an error, you can error control by seeing if
the class of the resulting object is "try-error". For example, let's say
I wanted to ma
"Bill Shipley" <[EMAIL PROTECTED]> writes:
> Hello. I am having a problem setting up a self-starting function for
> use in nonlinear regression (and eventually in the mixed model version).
> The function is a non-rectangular hyperbola - called "NRhyperbola" -
> which is used for fitting leaf phot
Dear Simon,
The population partial correlation rho[12|3...p] is 0 when the regression
coefficient beta[2] for x[2] from the regression of x[1] on x[2] ... X[p] is
0. Thus, the usual t-test for a regression coefficient also tests that the
partial correlation is 0.
Now, the sample partial correlati
Ciao Emanuele,
you could give a look to this contribute on fitting
distributions with R, maybe it could be helpful to
you:
http://cran.r-project.org/doc/contrib/Ricci-distributions-it.pdf
Regards,
Vito
Emanuele Mazzola wrote
Hello to everybody,
I'd like to submit a problem I'm dealing with
Hi Jose,
I need to make a small correction in my code -
mod$data$B works for glm objects, but not lm or aov objects. For
those use, mod$model$B.
On Oct 17, 2005, at 6:26 AM, Martin Henry H. Stevens wrote:
> Hi Jose,
> I am just beginning to plumb the depths of lattice, but perhaps my
> recent
Hi Jose,
I am just beginning to plumb the depths of lattice, but perhaps my
recent experience can help.
I recently figured out (with encouragement from the list) how to plot
predicted values from a model into the appropriate panel. I am
certain that what I have done can be done better, but th
Hi
You are looking for functions
panel.*
especially panel.lmline, but I wondered if you can use linear for one
panel and quadratic for other panels. You could use a structure
provided in examples in xyplot help page to try to achieve what
you want. I do not have instant solution to your proble
The Matrix package is under active development and the documentation
has not caught up with the code. Examples of usage can be found in
the tests subdirectory of the source package. At present we are
concentrating on the class hierarchy and writing methods and test
cases for those methods. Becau
[mailto:[EMAIL PROTECTED]
Sent: 13 October 2005 11:33
To: rob foxall (IFR)
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] Help with Matrix package
The first thing is to ensure that you are using an optimized BLAS. On
Windows, use Goto's BLAS if you have it (is not currently availabl
The first thing is to ensure that you are using an optimized BLAS. On
Windows, use Goto's BLAS if you have it (is not currently available and
redistribution is not allowed) or one of the pre-built ATLAS-based
Rblas.dll on CRAN or (best of all) optimize your own build of ATLAS.
The Matrix packa
It's because you've asked R to create an object of size 5 * 4 * 12488 *
1000 * 8 bytes, or 1905.5 MB, and R failed to get that much memory from
the operating system.
Andy
> From: Gao Fay
>
> Hi there,
>
> When I run the following, why does it give a error like that?
>
> > res2<-array(0,c(5,4,1
On Sat, 1 Oct 2005, Jim Hurd wrote:
>
> Which provides data in DTA (STATA), XPT (SAS), and POR (SPSS) formats all
> of which I have tried to read with the foreign package but I am not able to
> load any of them. I have 2 gb of RAM, but R crashes when the memory gets
> just over 1 GB. I am using Wi
Shengzhe Wu <[EMAIL PROTECTED]> wrote:
> In the Unix environment, I open a window by x11(). May I specify the
> position of this window by specifying the position of the top left of
> the window as in Windows environment?
I use "xwit" (version 3.4), a system command which manipulates existing
X w
Marc Schwartz (via MN) wrote:
>
> I don't believe so. In general, under Unix/Linux, the Window Manager
> determines window positioning upon startup unless the application
> overrides this behavior. Some applications let you specify application
> window positioning via command line 'geometry' argu
On Mon, 2005-09-26 at 17:45 +0200, Shengzhe Wu wrote:
> Hello,
>
> In the Unix environment, I open a window by x11(). May I specify the
> position of this window by specifying the position of the top left of
> the window as in Windows environment? Or some other parameters can be
> used to do that?
Dimitris Rizopoulos med.kuleuven.be> writes:
>
> you don't need L-BFGS-B in this case since you can easily
> re-parameterize you problem, i.e., you can always write:
>
> pi_i = exp(a_i) / sum(exp(a_i)), with e.g., a_1 = 0
>
> and thus maximize with respect to a_i's.
>
> I hope it helps.
>
>
you don't need L-BFGS-B in this case since you can easily
re-parameterize you problem, i.e., you can always write:
pi_i = exp(a_i) / sum(exp(a_i)), with e.g., a_1 = 0
and thus maximize with respect to a_i's.
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistica
Hi Manoj,
AFAIK there is no way in R to optimize functions under general linear
equality constraints (for inequality constraints there is constrOptim). In
your case you could reparametrize such that you estimate weights w(1) ...
w(n-1) and simply compute w(n) from those.
Even so, it would be great
RSiteSearch("Dagum")
Francisco
>From: "Uri Iskin" <[EMAIL PROTECTED]>
>To:
>Subject: [R] help with estimating parameters with nls
>Date: Tue, 20 Sep 2005 15:44:07 -0300
>
>Dear helpeRs,
>
>I have a vector containing values of incomes and I would like to estimate
>the three parameters of a Dagum
Have you received a reply to your question? I haven't seen one.
I got from Google a reference that gave a definition(a); it doesn't
seem like it would be hard to program. However, before I did that, I
think I would review the contents of the circular and CircStats
package
Krishna <[EMAIL PROTECTED]> writes:
> Hi all
>
> i am facing a peculiar problem for data input using read.table which i
> never faced previously.
>
> i have a data file by name abnew.txt with two coloumns data as depicted below.
>
> A B
> 420 422
> 314 321
(with A left-justified and B righ
Hi
if you still have Excel file and can open it
open Excel
select what you want to read, including header
press ctrl-C
then in R
a<-read.delim("clipboard")
and to write from R
write.table(tab, "clipboard", sep = "\t", row.names = F)
open Excel
press ctrl-V
Looks like Excel left in your txt fil
Rogers, James A [PGRD Groton] wrote:
> Jose -
>
> Before implementing SNK and Duncan's, you may want to be aware of some
> criticisms of these methods:
>
>>From Hsu (1996),
>
> "Newman-Keuls multiple range test is not a confident inequalities method and
> cannot be recommended."
>
> "Duncan's
Jose -
Before implementing SNK and Duncan's, you may want to be aware of some
criticisms of these methods:
>From Hsu (1996),
"Newman-Keuls multiple range test is not a confident inequalities method and
cannot be recommended."
"Duncan's multiple range test is not a confident inequalities metho
Hi Jim,
Many thanks for the function!
--
Jose Claudio Faria
Brasil/Bahia/UESC/DCET
Estatistica Experimental/Prof. Adjunto
mails:
[EMAIL PROTECTED]
[EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: 73-3634.2779
__
R-help@stat.math.ethz.ch mailing list
http
On 9/11/05, Jose Claudio Faria <[EMAIL PROTECTED]> wrote:
> Gabor Grothendieck wrote:
>
> > On 9/10/05, Jose Claudio Faria <[EMAIL PROTECTED]> wrote:
> >
> >>Dear R-list,
> >>
> >>Could anybody tell me how to make one matrix as the below:
> >>
> >> [,1] [,2] [,3] [,4] [,5] [,6]
> >>[1,]-
Gabor Grothendieck wrote:
> On 9/10/05, Jose Claudio Faria <[EMAIL PROTECTED]> wrote:
>
>>Dear R-list,
>>
>>Could anybody tell me how to make one matrix as the below:
>>
>> [,1] [,2] [,3] [,4] [,5] [,6]
>>[1,]-23456
>>[2,]2-2345
>>[3,]32
On 9/10/05, Jose Claudio Faria <[EMAIL PROTECTED]> wrote:
> Dear R-list,
>
> Could anybody tell me how to make one matrix as the below:
>
> [,1] [,2] [,3] [,4] [,5] [,6]
> [1,]-23456
> [2,]2-2345
> [3,]32-234
> [4,]4
Is this what you want?
make.odd.matrix<-function(matsize) {
newmat<-matrix(0,nrow=matsize,ncol=matsize)
for(i in 1:matsize) {
for(j in 1:matsize)
newmat[i,j]<-ifelse(i==j,NA,abs(i-j)+1)
}
return(newmat)
}
Jim
__
R-help@stat.math.ethz.ch
howzit dom
i just saw your mail now!
i tried the following and got an error:
> u=v=seq(1,3)
> u
[1] 1 2 3
> v
[1] 1 2 3
> f=function(x,y){max(x+y-1,0)}
> z=outer(u,v,f)
Error in outer(u, v, f) : dim<- : dims [product 9] do not match the
length of object [1]
it seems s if r is not performing el
On Mon, 5 Sep 2005, [EMAIL PROTECTED] wrote:
> A trap easy to fall into!
>
So easy, in fact, that it's a FAQ.
-thomas
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http:/
On 05-Sep-05 Dominique Katshunga wrote:
> Dear helpeRs,
> I seem to be a little bit confused on the result I am getting from the
> few codes below:
>> u=v=seq(0,1,length=30)
>> u
> [1] 0. 0.03448276 0.06896552 0.10344828 0.13793103 0.17241379
> [7] 0.20689655 0.24137931 0.27586207 0.31034
Dear Bjørn-Helge,
Sorry, I wrote the wrong number of observation. It should be 1184.
I saw on the book that variance is defined by sd^2. If variation is a
different concept from variance and defined by sd^2*(n-1) ? Since I
formerly took variance and variation as the same.
Thank you,
Shengzhe
Shengzhe Wu writes:
> I have a data set with 15 variables (first one is the response) and
> 1200 observations. Now I use pls package to do the plsr as below.
[...]
> Because the trainSet has been scaled before training, I think Xtotvar
> should be equal to 14, but unexpectedly Xtotvar = 16562,
On Sat, 27 Aug 2005, Shengzhe Wu wrote:
> Hello,
>
> I use qda (package "MASS") to obtain an object. If there is any
> function to plot density plot of qda object with one dimension$B!)(B
That makes no sense. qda objects do not have a density.
Please do study the background material (as you h
Hello,
I use qda (package "MASS") to obtain an object. If there is any
function to plot density plot of qda object with one dimension？
Thank you,
Shengzhe
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do r
On Fri, 26 Aug 2005, Luwis Tapiwa Diya wrote:
> I have been trying to write my own user defined function in Rpart.I
> imitated the anova splitting rule which is given as an example.In the
> work I am doing ,I am calculating the concentration index(ci) ,which
> is in between -1 and +1.So my devianc
On Fri, 26 Aug 2005, Shengzhe Wu wrote:
> I compared "posterior" of these three prediction results, they are a
> little different.
>
> The book you mentioned should be "Modern Applied Statistics with S.
> 4th edition". But this book has been borrowed out from our univeristy
> library by someone el
I compared "posterior" of these three prediction results, they are a
little different.
The book you mentioned should be "Modern Applied Statistics with S.
4th edition". But this book has been borrowed out from our univeristy
library by someone else, and I have checked the book "Pattern
Recognition
On Fri, 26 Aug 2005, Shengzhe Wu wrote:
> Thanks for your reply. Actually I called function as below.
>
> p1 = predict(object, newdata, dimen=1)
> p2 = predict(object, newdata, dimen=1, method=debiased)
> p3 = predict(object, newdata, dimen=1, method="predictive")
So why did you say something dif
On Fri, 26 Aug 2005, Shengzhe Wu wrote:
> I use lda (package: MASS) to obtain a lda object, then want to employ
> this object to do the prediction for the new data like below:
>
> predict(object, newdata, dimen=1, method=c("plug-in", "predictive",
> "debiased"))
That is not how you call it: when
> "Don" == Don MacQueen <[EMAIL PROTECTED]>
> on Thu, 25 Aug 2005 08:11:42 -0700 writes:
Don> Also, for the three dimensional graphic,
Don> help.search("3d")
Don> will lead to a reference to the cloud() function in the lattice
package.
Don> I don't remember if the lat
Sorry to reply to my own post, but in reviewing the NAMESPACE file for
lme4, it looks like Doug is perhaps planning to add additional model
object accessor methods for the lmList class, including resid() and
summary(), which are commented out now. coef() is available presently.
So when in place,
That works also (using the example in ?lmList)
library(lme4)
?lmList
fm1 <- lmList(breaks ~ wool | tension, warpbreaks)
However, one still would need to use either sapply() or lapply() as
below to get the details that Krishna is looking for.
'fm1' above is a list of models (S4 class 'lmList'
You might find image() useful, to plot the factory map and overlay the sound
levels as colors. You could also use persp(), or have a look at rgl.surface
in the rgl package on http://cran.r-project.org.
Reid Huntsinger
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] O
Also, for the three dimensional graphic,
help.search("3d")
will lead to a reference to the cloud() function in the lattice package.
I don't remember if the lattice package is installed by default. If
not, you will have to install it.
(If you're using a Mac or Windows computer, there's a menu
What about using lmList from the lme4 package?
Randy
On 8/25/05 9:44 AM, "Marc Schwartz" <[EMAIL PROTECTED]> wrote:
> Also, looking at the last example in ?by would be helpful:
>
> attach(warpbreaks)
> tmp <- by(warpbreaks, tension, function(x) lm(breaks ~ wool, data = x))
>
> # To get coeffi
Le 25.08.2005 15:54, Philippe a écrit :
>Hello,
>My name is Philippe Favrot, I'm a french occupational doctor (working in
>Luxembourg), and I'm an "R beginner".
>I would be happy you could help me about the utilization of "R".
>Recently, I measured sound levels in a plant. Before the measuring, I
Philippe a écrit :
> ... how can I draw a three dimensional graphic of the
> plant's sound level ?
Not an answer, but I find image() very useful
to visualize 3D data simply on a 2D colored image.
hih
__
R-help@stat.math.ethz.ch mailing list
https://st
Also, looking at the last example in ?by would be helpful:
attach(warpbreaks)
tmp <- by(warpbreaks, tension, function(x) lm(breaks ~ wool, data = x))
# To get coefficients:
sapply(tmp, coef)
# To get residuals:
sapply(tmp, resid)
# To get the model matrix:
sapply(tmp, model.matrix)
To get th
Look more carefully at
?lm
at the See Also section ...
X <- rnorm(30)
Y <- rnorm(30)
lm(Y~X)
summary(lm(Y~X))
Best,
Matthias
> Hi all
>
> I used a function
> > qtrregr <- by(AB, AB$qtr, function(AB) lm(AB$X~AB$Y))
>
> objective is to run a regression on quartery subsets in the
> data set
Hello.
Sorry for the delay, (not so) long vacations...
Seems to me that it has something to deal with the 'adehabitat' package
which is very useful for radio-tracking data. ('adehabitat' is available
on CRAN)
BUT the 'getareahr' function doesn't exist. Could it be a mix between
'getverticeshr
Agnes Gault wrote:
> Hello
>
> I am working on radio tracking data, with a short programme someone gave me
> and ... which should, supposedly, work ... In this programme, there is the
> function : getareahr(kern, levels = 95). But i cannot find any 'getareahr'
> in R ...
>
> could anyone help
Agnes Gault a écrit :
> Hello
>
> I am working on radio tracking data, with a short programme someone gave me
> and ... which should, supposedly, work ... In this programme, there is the
> function : getareahr(kern, levels = 95). But i cannot find any 'getareahr'
> in R ...
>
> could anyone he
You have not told us what software you used to get the results you
present. My first question is whether you are working with prices or
log(prices)? If the former, I suggest you consider the latter; price
changes tend to be much better behaved, more nearly normal, etc., on the
log
I think a question about the S family() functions is best determined by
reading the S(-PLUS) documentation. The weights component computes the
working weights for the IWLS, unsurprisingly (and as stated by
?family.object). R has mu.eta to do that, the crucial line in R's glm.fit
being
w
On Mon, 8 Aug 2005, Haibo Huang wrote:
> I had a logit regression, but don't really know how to
> handle the "Warning message: non-integer #successes in
> a binomial glm! in: eval(expr, envir, enclos)"
> problem. I had the same logit regression without
> weights and it worked out without the warni
trellis.device()
xyplot(states ~ size, groups=type, data=tmp,
panel = "panel.superpose",
panel.groups ="panel.linejoin", auto.key=TRUE)
Jamieson Cobleigh wrote:
> I have a data frame with three columns, type (a factor with two
> values: "Monolithic" and "Composi
On 8/8/05, Krishna <[EMAIL PROTECTED]> wrote:
> Hi Everyone
>
> May I request for a small help while performing the regression analysis.
>
> I would like to know is there any possibility of conducting the
> regression for different data subsets (in the same data file),
> classified on the basis o
you could use function lmList() from the nlme package, i.e.,
dat <- data.frame(y = rnorm(120), x = runif(120, -3, 3), g = rep(1:3,
each = 40))
library(nlme)
m <- lmList(y ~ x | g, data = dat)
m
summary(m)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biost
Hi Rangesh,
Perhaps I mis-understand your question, but it could be as simple as...
p <- 1-exp(-exp(s))
In R, this is vectorized such that a new vector is calculated - one
value for each value of s, so p will have the same length as s.
In the "Introduction to R" read up on vectors and how to av
On Fri, 5 Aug 2005, Peter Dalgaard wrote:
> "Grathwohl, Dominik, LAUSANNE, NRC-BAS" <[EMAIL PROTECTED]> writes:
>
>> Dear R-helpers,
>>
>> First of all I have nothing against the French language!
>> But now my problem, yesterday I installed R 2.1.1
>> and I had to experience that my RGui is speak
"Grathwohl, Dominik, LAUSANNE, NRC-BAS" <[EMAIL PROTECTED]> writes:
> Dear R-helpers,
>
> First of all I have nothing against the French language!
> But now my problem, yesterday I installed R 2.1.1
> and I had to experience that my RGui is speaking French.
> My windows locals is French (Switzer
Grathwohl, Dominik, LAUSANNE, NRC-BAS wrote:
> Dear R-helpers,
>
> First of all I have nothing against the French language!
> But now my problem, yesterday I installed R 2.1.1
> and I had to experience that my RGui is speaking French.
> My windows locals is French (Switzerland).
> I'm used to Engl
On 8/4/2005 2:25 PM, Spencer Graves wrote:
> If you are using XEmacs, it's not that easy. I'm not completely
> clear on the protocol in that case, but I think you have to first click
> in the R window, then + + . Perhaps an Emacs
> wizard will enlighten us.
>
> spencer grave
The R for Windows FAQ has many useful tips, including this one:
http://cran.r-project.org/bin/windows/base/rw-FAQ.html
It's well worth at least reading the contents, so you know what
questions it can answer.
- Martin
Duncan Murdoch wrote:
>On 8/4/2005 11:39 AM, Shengzhe Wu wrote:
>
>
>>and
If you are using XEmacs, it's not that easy. I'm not completely
clear on the protocol in that case, but I think you have to first click
in the R window, then + + . Perhaps an Emacs
wizard will enlighten us.
spencer graves
Duncan Murdoch wrote:
> On 8/4/2005 11:39 AM, S
501 - 600 of 1084 matches
Mail list logo
