> glm.24.pred<-predict(glm.24,newdata=nestday, type="response", SE.fit=T)
What is SE.fit? The help says se.fit. That _may_ be a problem.
However, I think the real problem is that the link function argument
includes a reference to vc.apfa$days that is appropriate for fitting, not
prediction. O
Note:
> configure:4041: checking whether the C compiler works
> configure:4059: error: cannot run C compiled programs.
so your compiler installation is broken, seriously enough that configure
can make no progress.
Also, gcc 3.4.2 is old and has a known bug that stops some R packages
using Fort
Hi,
I am trying to compile R-2.2.1 on Solaris 2.9 with a 64-bit build. Following
the instructions in "R Installation and Adminstration", I changed the
following settings in "config.site":
CC="gcc -m64"
F77="g77 -64"
CXX="g++ -m64"
LDFLAGS="-L/usr/local/lib/sparcv9 -L/usr/local/lib"
But I got th
Hi,
I am trying to compile R-2.2.1 on Solaris 2.9 with a 64-bit build. Following
the instructions in "R Installation and Adminstration", I changed the
following settings in "config.site":
[[alternative HTML version deleted]]
__
R-help@stat.math
Detlef Steuer wrote:
> Hi again,
>
...
> This new rpm R-base-2.3.0-beta should automatically resolve dependencies. At
> least it did so on my machine. I would be happy to get a report, if you try
> to install this one and find difficulties or success!.
> Most important would be to report missin
You need to tell us what compare() in S-PLUS does; i.e., what does it take
as input and what does it output. Presumably the information would be in
the help page for compare().
Telling us what the function looks like in S-PLUS is of no use, since most
likely the code cannot be used in R without
As suggested below,
sign(ii*1e9-f)
will give you the same result as the S-plus
compare(ii*1e-9,f)
-Christos
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Bellinger Instruments
P/L
Sent: Thursday, April 20, 2006 11:02 PM
To: 'Prof Brian Ripley'
Cc:
On 20 April 2006 at 22:00, zhongmiao wang wrote:
| Hello:
| I am generating a random number with rnorm(1). The generated number
| has 8 decimals. I don't want so many decimals. How to control the
| number of decimals in R?
You must have installed the professional version of R, just downgrade to t
The sign() function requires one argument,
and returns the sign of the value, where I require a number
-1, 0, or 1 depending on if the value is <, ==, or > the
first argument.
This is the Splus function when listed,
> compare
function(e1, e2)
.Internal(compare(e1, e2), "do_op", T, 18)
unfortunat
Hello:
I am generating a random number with rnorm(1). The generated number
has 8 decimals. I don't want so many decimals. How to control the
number of decimals in R?
Thanks!
Zhongmiao Wang
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/ma
An update for all:
Using the combined contributions from Mark and Dr. Ripley, I've been
(apparently) successfully formulating both GLM's and GLMM's (using
the MASS function glmmPQL) analyzing my nest success data. The beta
parameter estimates look reasonable and the top models resemble those
from
To: Martin Maechler <[EMAIL PROTECTED]>
> Indeed. Please do check the archives.
Yep. Post is there.
> Now back to the subject: Jose, I think your main contribution
> is based on "autoHotKeys" and that only works on Windoze, right?
> Michael explicitly mentioned he's working in Mac OS X.
>
> Ma
Is this what you are after?
floor(data[,model.list])
Or I just didn't understand what you are trying to accomplish?
cheers
Francisco
PS: try to avoid using names that are already reserved to a function like
"data" See ?data
>From: Chad Reyhan Bhatti <[EMAIL PROTECTED]>
>To: R-help@stat.math.
On 4/20/2006 7:02 PM, Chad Reyhan Bhatti wrote:
> Hello,
>
> I am not sure what to write in the subject line, but I would like to take
> a character string that is a variable in a data frame and apply a function
> that takes a numeric argument to this character string.
Remember that dataframes ar
Try this:
# test data
set.seed(1)
DF <- data.frame(x1 = rnorm(5), x2 = rnorm(5), x3 = rnorm(5))
DF
model.list <- c("x2", "x3")
# transform
for(v in model.list) DF[v] <- floor(DF[v])
On 4/20/06, Chad Reyhan Bhatti <[EMAIL PROTECTED]> wrote:
> Hello,
>
> I am not sure what to write in the subject
Hello,
After getting help to solve part of my problem and some delay on my
part, I am posting a more refined version to see if someone can help me
further. I am trying to autocalibrate a model in my subject area using
the snow and rgenoud packages. I want to use the key function "fn" that
is cal
Hello,
I am not sure what to write in the subject line, but I would like to take
a character string that is a variable in a data frame and apply a function
that takes a numeric argument to this character string.
Here is a simplified example that would solve my problem.
Imagine I have my data stor
I think you're asking how to plot arrows for variables in one dataset
onto the PCA plot for another dataset? You can do this with functions
from the vegan package. Particularly, vector and factor fitting can be
done with envfit() and surface fitting with ordisurf().
There's a good tutorial at:
Hello,
Im trying to fit the following model:
Dependent variable: MAXDEPTH (the maximum depth reached by a penguin during
a given dive)
Fixed effects: SUCCESSMN (an index of the individual quality of a bird),
STUDYDAY (the day of the study, from -5 to 20, with 0=Dec 20), and the
interaction
On 4/19/06, Doran, Harold <[EMAIL PROTECTED]> wrote:
> I am writing a small function to approximate an integral that cannot be
> evaluated in closed form. I am partially successful at this point and am
> experiencing one small, albeit important problem. Here is part of my
> function below.
> This
Hi everyone,
I'd like to project two pcas onto one device window.
I plot my first PCA:
biplot(prcomp(t(cerebdevmat)), var.axes=FALSE, cex=c(1,.1),
pc.biplot=TRUE)
Now I'd like to project the features of another PCA onto this graph.
Any suggestions?
I know this is easily done in MatLab but
This worked for my example data frame but not when I used it on my
actual data frame because some of the levels of x have only one element.
However, when I replace 'sample' with 'some' (from the car package),
it works like a charm!
Many thanks to Andy and Gabor for their help!
Kelly
Liaw, A
Thanx Marc and Gabor for your help.
Sachin
"Marc Schwartz (via MN)" <[EMAIL PROTECTED]> wrote:
On Thu, 2006-04-20 at 11:46 -0700, Sachin J wrote:
> Hi,
>
> How can I accomplish this in R. Example:
>
> R1 R2
> 3 101
> 4 102
> 3 102
> 18 102
> 11 101
>
> I want to find Sum(101) = 14 - i.e
On Thu, 2006-04-20 at 11:46 -0700, Sachin J wrote:
> Hi,
>
> How can I accomplish this in R. Example:
>
> R1 R2
> 3 101
> 4 102
> 3 102
> 18102
> 11101
>
> I want to find Sum(101) = 14 - i.e SUM(R1) where R2 = 101
> Sum(102) = 25
Using the built in data.frame iris sum each of the first 4 columns for
each value of the 5th column.
rowsum(iris[,-5], iris[,5])
On 4/20/06, Sachin J <[EMAIL PROTECTED]> wrote:
> Hi,
>
> How can I accomplish this in R. Example:
>
> R1 R2
> 3 101
> 4 102
> 3 102
> 18102
>
Hi,
How can I accomplish this in R. Example:
R1 R2
3 101
4 102
3 102
18102
11101
I want to find Sum(101) = 14 - i.e SUM(R1) where R2 = 101
Sum(102) = 25- SUM(R2) where R2 = 102
TIA
Sachin
-
Yes, you're exactly right. You can use the mcmcsamp() function to sample from
the posterior of an lmer object. This returns an object of mcmc class and you
can do all of your diagnostics using the coda package.
-Original Message-
From: Spencer Graves [mailto:[EMAIL PROTECTED]
Sent: Thur
Linda Lei wrote:
> Hi All,
>
>
>
> I have codes as follows to get the ecdf plots:
>
>
>
> >
> day.hos2<-c(6,4,6,6,4,6,5,4,7,5,6,6,8,6,17,9,8,4,6,3,5,8,7,12,5,10,6,4,6
> ,13,7,6,6,25,4,9,96,6,6,6,6,9,4,5,5,4,10,5,7,6)
>
>
>
> >
> day.hos3<-c(5,6,7,6,4,5,6,6,6,6,19,7,5,9,8,8,7,5,6,20,40
Hi, Harold:
Am I correct that the tool currently preferred for estimating
p-values for lmer is "mcmcsamp"?
Amélie: My favorite tool for exploring the archives is 'RSiteSearch'.
You can also get to it via www.r-project.org, but the last time I tried
to copy a web address to paste in
> "COMTE" == COMTE Guillaume <[EMAIL PROTECTED]> writes:
> I know it isn't the best way to represent data, but people
> are sometimes more interested by the look and feel than by
> the accuracy of the results...
If they aren't intersted in the results, why not just print some
rand
You didn't do anything wrong, lmer doesn't give them. And, for good reason.
I've been a bit indoctrinated by D. Bates, so let me share what I've learned.
With simple analysis of variance models with simple error structures, it is
known that the ratio of the variances follow and F distribution. H
Thanks for providing such a simple, complete example. I've never
used "deal" before, but a few minutes with your example led me to
something that might help you:
First, the documentation for "localprior" says, "x: an object of
class 'node' or 'network'." That information l
Hello,
Im trying to perform a REML analysis using the lmer() function (lme4
package). Well, it seems to work well, except that Im not getting any
p-value (see example below). Can someone tell me what I did wrong?
Thanks for your help,
Amélie
> library(gdata)
> dive <- read.xls("C
On Thu, Apr 20, 2006 at 12:29:58PM -0400, Gabor Grothendieck wrote:
> Just one other thought. There is a command line program
> called xlHtml that you can find via google that will convert xls files to csv
> and its used like this assuming you have placed it
> somewhere in your path:
>
> read.c
Better idea: Compare directly. ?qqplot
-- 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 PRO
Hi All,
I have codes as follows to get the ecdf plots:
>
day.hos2<-c(6,4,6,6,4,6,5,4,7,5,6,6,8,6,17,9,8,4,6,3,5,8,7,12,5,10,6,4,6
,13,7,6,6,25,4,9,96,6,6,6,6,9,4,5,5,4,10,5,7,6)
>
day.hos3<-c(5,6,7,6,4,5,6,6,6,6,19,7,5,9,8,8,7,5,6,20,40,5,8,7,7,5,6,13,
11,9,4,6,9,16,6,7,6)
> f<-e
I quoted the relevant part of the documentation for you. Have you actually
try to read what it says?
Sure, you don't get any error, but have you checked whether any
bootstrapping was actually done? Most of those functions are generics, thus
having the "..." argument that can take anything. Doe
All,
I have been able to successfully use the optim( ) function with
"L-BFGS-B" to find reasonable parameters for a one-compartment
open pharmacokinetic model. My loss function in this case was
squared error, and I made no assumptions about the distribution
of the plasma values. The model appear
Just one other thought. There is a command line program
called xlHtml that you can find via google that will convert xls files to csv
and its used like this assuming you have placed it
somewhere in your path:
read.csv(pipe("xlHtml -te -xc:1-10 -csv myfile.xls"))
Its handy since its just a sing
btw, if you change "myFun" to any R internal function, such as "cov", or
"corr", it can run successfully...
On 4/20/06, Michael <[EMAIL PROTECTED]> wrote:
>
> Andy,
>
> I've noticed there should be a "weight" or "frequency" somewhere... but my
> function does not need it.
>
> I have tried to cheat
x.unique$V1 gives the list of individual column's unique values. Thank you
again Andy.
Sachin
"Liaw, Andy" <[EMAIL PROTECTED]> wrote:
From:Sachin J
>
> Hi,
>
> This one is not working for me. It is listing all the rows
> instead of unique ones. My dataset has 30 odd rows and
> fol
"Liaw, Andy" <[EMAIL PROTECTED]> wrote: From:Sachin J
>
> Hi,
>
> This one is not working for me. It is listing all the rows
> instead of unique ones. My dataset has 30 odd rows and
> following is the resulting o/p
>
> [[308313]]
> [1] 126
> [[308314]]
> [1] 126
> [[308315]]
> [1] 126
>
But it is not giving me the list of unique elements. Count works fine.
Sachin
"Liaw, Andy" <[EMAIL PROTECTED]> wrote:
From:Sachin J
>
> Hi,
>
> This one is not working for me. It is listing all the rows
> instead of unique ones. My dataset has 30 odd rows and
> following is the res
Thanks Andy. That works.
Sachin
"Liaw, Andy" <[EMAIL PROTECTED]> wrote:
From:Sachin J
>
> Hi,
>
> This one is not working for me. It is listing all the rows
> instead of unique ones. My dataset has 30 odd rows and
> following is the resulting o/p
>
> [[308313]]
> [1] 126
> [[30831
From:Sachin J
>
> Hi,
>
> This one is not working for me. It is listing all the rows
> instead of unique ones. My dataset has 30 odd rows and
> following is the resulting o/p
>
> [[308313]]
> [1] 126
> [[308314]]
> [1] 126
> [[308315]]
> [1] 126
> [[308316]]
> [1] 126
> [[30
Hi,
This one is not working for me. It is listing all the rows instead of unique
ones. My dataset has 30 odd rows and following is the resulting o/p
[[308313]]
[1] 126
[[308314]]
[1] 126
[[308315]]
[1] 126
[[308316]]
[1] 126
[[308317]]
[1] 126
[[308318]]
[1] 126
[[308319]]
Hi all,
The package crossdes could contruct a complete sets of mutually orthogonal
latin squares.
The construction works for prime powers only.
I hope to know whether there is a way to construct a mutually orthogonal Lation
square for
10 or other numbers that could not be prime powers.
Thanks
--- Forwarded message follows ---
From: Petr Pikal <[EMAIL PROTECTED]>
To: [EMAIL PROTECTED]
Subject:Re: [R] Breakdown a number
Date sent: Thu, 20 Apr 2006 17:22:03 +0200
Hi
try this
fff<-function(x,breaks=c(5,10)) {
if
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Michael
> Sent: Thursday, April 20, 2006 3:50 AM
> To: R-help@stat.math.ethz.ch
> Subject: [R] Bootstrap error message: Error in
> statistic(data, original, ...) : unused argument(s) ( ...) [Broadcast]
This might help:
> x <- read.table("clipboard", colClasses=c("numeric", "character"))
> (x.unique <- lapply(x, unique))
$V1
[1] 155 138 126 123 103 143 111 156
$V2
[1] "A" "B" "C" "D"
> sapply(x.unique, length)
V1 V2
8 4
Andy
From: Sachin J
>
> Hi,
>
> I have a dataset which has bot
Hi,
I have a dataset which has both numeric and character values with
dupllicates. For example:
155 A
138 A
138 B
126 C
126 D
123 A
103 A
103 B
143 D
111 C
111 D
156 C
How can I count the number of unqiue entries without counting duplicate
entries. Also c
A flexible way for doing this is to define logical vectors for the types of
samples that you want to include or exclude. You can then use logical
negation to select the complementary set:
inSet1 <- HData$H < 1.3
inSet2 <- HData$H < 8 & Hdata$DBH > 20
HDataPart1 <- Hdata[!inSet1, ]
HDataPart2 <-
Just negate the result of %in%; e.g.,
! a %in% b
Andy
From: Ghislain Vieilledent
>
> Dear r-list,
>
> I've got a data base:
> > HData[1:10,]
>NumTree Site Species Date Age DBHH IdentTree
> 11 Queige Spruce 2002 184 49 33.5 Queige 1
> 22 Queige Fir 2002 NA 5
> a <- 1:10
> a
[1] 1 2 3 4 5 6 7 8 9 10
> c(1,2,3,11,12,13) %in% a
[1] TRUE TRUE TRUE FALSE FALSE FALSE
> ! c(1,2,3,11,12,13) %in% a
[1] FALSE FALSE FALSE TRUE TRUE TRUE
Is this it?
Gabor
On Thu, Apr 20, 2006 at 02:54:39PM +0200, Ghislain Vieilledent wrote:
> Dear r-list,
>
Dear r-list,
I've got a data base:
> HData[1:10,]
NumTree Site Species Date Age DBHH IdentTree
11 Queige Spruce 2002 184 49 33.5 Queige 1
22 Queige Fir 2002 NA 5 4.6 Queige 2
33 Queige Fir 2002 25 8 6.6 Queige 3
44 Queige Spruce 2002 1
Hi,
I have a simple problem writing a function that is to be called like
Myfunction( Column1 = low, Column2 = high, Column3 = all, Column4 = all,
data = mydata)
{.
contourplot(z ~ mydata$Column3* mydata$Column3)
.
}
Where Column1 and Column1 are the names of the dataframe mydata.
How do I par
Try:
c(sapply(1:2, function(i) [EMAIL PROTECTED]@my.slot))
On 4/20/06, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
> Hello,
>
> I don't manage to see if you have already focussed on this point in some
> previous messages so I post my question:
>
> I have a little problem with the S4 style of pr
Hello,
I don't manage to see if you have already focussed on this point in some
previous messages so I post my question:
I have a little problem with the S4 style of programming.
I tried to formalize my question: please consider the following example
that you can run I think:
#
Thak you very much indeed!
Eric.
On 4/20/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
>
> Try paste:
>
> paste("Axis", f)
>
> On 4/20/06, Eric Ferreira <[EMAIL PROTECTED]> wrote:
> > Dear R-colleagues,
> >
> > Is it possible to mix TEXT and VALUE of objects in y (or x) label of a
> plot?
> >
Dr. Paul Murrell will offer his "Graphics in R" course
online at statistics.com May 5 June 2.
This course teaches you how to produce publication-quality
statistical plots of data using R (a freely available open-
source statistical language and environment). It will
cover plots such as scatterp
Hi,
I haven't recieved any replies to my last email, so let me be a bit more
specific. I have a dataframe and it has the following structure:
Condition
Mapping SubjectABC
11 510 15
12
Detlef:
For your information. A few days ago I ordered SuSE 10, the version
corresponding to the old "Professional". It arrived today, and was
installed. When the R-base rpm was downloaded with Mozilla Firebird,
it was only to open it when the download was finished, click the
"install with YaST
I think this does what you require.
#Read your data in whatever way you wish:
d1<-data.frame(Date=c("2005/1/1","2005/2/1","2005/1/3","2005/1/4","2005/
1/7","2005/3/5"),
x=c(119,123,-110,114,11,200),
y=c(230,-125,300,-21,299,311))
d2<-data.frame(Date=c("2005/1/3"
I like to use the RODBC package for doing this. Here is my code sample:
xls <- odbcConnectExcel(fname)
rawdata.temp <- sqlFetch(xls, "rawdata", max=2800)
close(xls)
fname is the full path to the file and "rawdata" is the name of the excel
sheet I want to import. I trie
On Thu, 2006-04-20 at 10:28 +0200, Fredrik Karlsson wrote:
> Hi marc,
>
> I did not manage ctab to do this for me. Again, I am probably using it
> wrong, but I don't know that the problem is.
>
> You asked for a more illustrated example, so here goes:
>
> Take this table:
>
> > ftable(table(sam
Try paste:
paste("Axis", f)
On 4/20/06, Eric Ferreira <[EMAIL PROTECTED]> wrote:
> Dear R-colleagues,
>
> Is it possible to mix TEXT and VALUE of objects in y (or x) label of a plot?
>
> For instance:
>
> f<-2
> plot(..., ylab='Axis f', ...)
>
> where "f" means the VALUE of "f".
>
> Thany you in
f <- 2
plot(runif(10), ylab=paste("Axis", f, collapse=" "))
?paste
Eric Ferreira wrote:
> Dear R-colleagues,
>
> Is it possible to mix TEXT and VALUE of objects in y (or x) label of a plot?
>
> For instance:
>
> f<-2
> plot(..., ylab='Axis f', ...)
>
> where "f" means the VALUE of "f".
>
> T
Read them in as zoo objects (you can replace textConnection(Lines1)
with the filename) and then merge them using all = FALSE to retain
only common time points. Note that in my English locale I had
to modify your Apl to Apr.
Lines1 <- "Date x y
Jan-1,2005120 230
Jan-
Dear R-colleagues,
Is it possible to mix TEXT and VALUE of objects in y (or x) label of a plot?
For instance:
f<-2
plot(..., ylab='Axis f', ...)
where "f" means the VALUE of "f".
Thany you in advance,
Eric.
--
Barba
Departamento de Ciências Exatas
Universidade Federal de Lavras
Minas Gerais
Try using breaks[breaks <= x] in place of breaks. Your solution
always returns three components if breaks has two whereas this one
does not but you could extend the output it if that is an essential part.
On 4/20/06, Paul Roebuck <[EMAIL PROTECTED]> wrote:
> On Thu, 20 Apr 2006, Gabor Grothendie
On Thu, 20 Apr 2006, Manuel Gutierrez wrote:
> Is it possible to include a factor in an nls formula?
Yes. What do you intend by it? If you mean what it would mean for a lm
formula, you need A[a] and starting values for A.
There's an example on p.219 of MASS4.
> I've searched the help pages w
> On 19-Apr-06 Peter Ehlers wrote:
>> This discussion of 3-d pie charts comes at an opportune time. I have
>> just formulated a new theory of graphical information transfer which
>> is particularly simple in the case of 3-d pie charts.
>>
>> Let theta denote the angle between the normal to the pie
Hello everybody !
I try to obtain in R eigenvectors and canonical analysis on MANOVA results,
but I don't find how to process?
In particular, I would be interesting to obtain "standardized canonical
coefficients" of the canonical variates. There analysis give some
information on the correlatio
Rolf Turner wrote:
>
> People really ***should not*** be encouraged or abetted in
> wrong-headedness. Excel is terrible. Pie charts are terrible.
> Don't mess with them. Period.
>
Now I realise the opportunity I missed on April 1st, when I was going
to try and (anonymously) post the most
Detlef Steuer wrote:
> On Thu, 20 Apr 2006 10:51:57 +0200
.
> Sometimes it is very hard to think of the questions a new user has,
> if you do it yourself on a daily basis for years.
I understand that very well. Now I can tell you that I used the info
in the "Hacking OpenSuse" link you prov
The following should work:
> dfr.samp <- dfr[tapply(1:nrow(dfr), dfr$x, sample, 1),]
> dfr.samp
x y z
10 a 10 J
2 b 2 B
9 c 9 I
Andy
From: Kelly Hildner
>
> I don't use R much, and I have been unable to figure out how
> to get the
> subset of my data frame that I would like.
>
> For
Thanks Andrew. I am now trying but without much
success. I don't now how to give start values for the
factor?.
Could you give me an example solution with my toy
example?
a<-as.factor(c(rep(1,50),rep(0,50)))
independ<-1:100
respo<-rep(NA,100)
respo[a==1]<-(independ[a==1]^2.3)+2
respo[a==0]<-(indepe
Hello everybody !
I try to obtain in R eigenvectors and canonical analysis on MANOVA results,
but I don't find how to process?
In particular, I would be interesting to obtain "standardized canonical
coefficients" of the canonical variates. There analysis give some
information on the correlatio
Manuel,
I don't think that it works very easily. Instead, try gnls() in the
nlme package.
Cheers
Andrew
On Thu, Apr 20, 2006 at 11:18:02AM +0200, Manuel Gutierrez wrote:
> Is it possible to include a factor in an nls formula?
> I've searched the help pages without any luck so I
> guess it is n
On Thu, 20 Apr 2006, Gabor Grothendieck wrote:
> On 4/19/06, Paul Roebuck <[EMAIL PROTECTED]> wrote:
> > On Wed, 19 Apr 2006, Gabor Grothendieck wrote:
> >
> > > On 4/19/06, Paul Roebuck <[EMAIL PROTECTED]> wrote:
> > >
> > > > Isn't there a builtin method for doing this and, if so,
> > > > what i
Is it possible to include a factor in an nls formula?
I've searched the help pages without any luck so I
guess it is not feasible.
I've given it a few attempts without luck getting the
message:
+ not meaningful for factors in:
Ops.factor(independ^EE, a)
This is a toy example, my realworld case is
Gabor Csardi wrote:
> So you don't want to read that manual after all? :)
I did not say so, not at all. The point is that it is not very
helpful to say "Read the manual" when the the problem is very much one
of not knowing where to look in the manual. Instructions for
installing on a particul
Hi All,
How can extract AIC,BIC from a fitted Garch model?
--
SUMANTA BASAK.
[[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! http://ww
Detlef Steuer wrote:
> Hi,
>
> I`m the one to blame for the readme :-) and for providing the rpms.
I am sorry, but the note was not primarily directed at you, perhaps (I
am not sure) just as much at Linux in general, or Suse, or OpenSuse,
or ... In any case, I am sure you do a great job at provi
Hi marc,
I did not manage ctab to do this for me. Again, I am probably using it
wrong, but I don't know that the problem is.
You asked for a more illustrated example, so here goes:
Take this table:
> ftable(table(sample(paste("dim1_no",1:5,sep=""),10,replace=TRUE),
sample(paste("dim2_no",1:5,se
Dear r-users,
Suppose I have three datasets:
Dataset-1:
Date x y
Jan-1,2005120 230
Jan-2,2005123 -125
Jan-3,2005-110 300
Jan-4,2005114 -21
Jan-7,200511299
Mar-5,2005200 311
Dataset-2:
Date
Hi all,
I can understand that it isn't the right way to represent data, i knew that i
would hurt some people on this mailling list to discuss about pie.
As i've specified on my first message:
My only purpose of drawing 3D pie is for customer who don't have to understand
what is drawn , but only b
On 19-Apr-06 Peter Ehlers wrote:
> This discussion of 3-d pie charts comes at an opportune time. I have
> just formulated a new theory of graphical information transfer which
> is particularly simple in the case of 3-d pie charts.
>
> Let theta denote the angle between the normal to the pie cylind
Hi all:
As to vsn package,how can I transform my raw data format to the demo data of
vsn(named "kidney")?
My raw data:
namesignal dye
gene1 815.32 green
gene2 671.66 green
gene3 713.93 green
gene4 703.97 green
gene5 493.59 green
gene6 477.92 green
gene7 346.55 green
"zhongmiao wang" <[EMAIL PROTECTED]> writes:
> Hello,
> I have used nlme to fit a model, the R syntax is like
>
> fmla0<-as.formula(paste("~",paste(colnames(ldata[,9:13]),collapse="+"),"-1"))
> > fmla1<-as.formula(paste("~",paste(colnames(ldata[,14:18]),collapse="+"),"-1"))
> > fmla2<-as.formula(
Dear colleagues,
I've been swamped and fighting with error for a few hours but still
desperately having absolutely no clue:
What's wrong with my bootstraping code?
Thanks a lot!
Error Message:
> bootResults=boot(X, myFun, R=1);
Error in statistic(data, original, ..
Dear r-users,
Suppose I have three datasets:
Dataset-1:
Date x y
Jan-1,2005120 230
Jan-2,2005123 -125
Jan-3,2005-110 300
Jan-4,2005114 -21
Jan-7,200511299
Mar-5,2005200 311
Dataset-2:
Date
On Wed, 2006-04-19 at 21:18 +0200, Marco Girardello wrote:
> Hi,
>
> is there an R function which can compute the R-squared
> for a mars model fitted using the mda package?
>
> many thanks in advance
>
> best regards
>
> Marco Girardello
Julian Faraway's "Extending the Linear Model with R" Cha
Hi,
is there an R function which can compute the R-squared
for a mars model fitted using the mda package?
many thanks in advance
best regards
Marco Girardello
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEAS
94 matches
Mail list logo