Hard as it is for me to imagine, the ggobi windows stay open and functional
while R (in emacs) has crashed in the background after throwing the error
messages. As a perhaps naive Linux user, I thought that if a parent process
crashed any processes it spawned would crash too. I guess not in this
On May 6, 2008, at 10:30 PM, Prof Brian Ripley wrote:
What were H W? For png() they are (by default) in pixels, for
pdf() in inches.
You haven't told us your OS, but I guess Mac OS. Please update to R
2.7.0: that offers you two new png() devices for higher-quality
plots, and various
Hi,
I have a dataframe wf existing of a header with different labels and beneath
the values of those labels :
wf:
label1 label2 ...
0,450,21
0,100,45
I have a list
fl - c(label2,label3,..)
Isn't possible to use the list elements in the list in order to subtract
values
Thanks Martin,
I also posted the question on the bioconductor list but I have no reply yet.
In the meanwhile I found out that instead of saying
d=list(data.matrix2,y,censored)
I should specify the arguments:
d=list(x=data.matrix2,y=y,censoring.status=censored)
Strange, huh? Anyway, it solves
Hi,
I want to find solution of function : f(x,y) = x'Cx - a under constraints :
0 x,y p
0 x-y q
where a, p,q are given constants and x = (x, y) and C is a 2X2 matrix (given)
Can anyone suggest me any R function to do that?
On Tue, May 6, 2008 at 11:11 PM, Mark Kimpel [EMAIL PROTECTED] wrote:
Hard as it is for me to imagine, the ggobi windows stay open and
functional while R (in emacs) has crashed in the background after throwing
the error messages. As a perhaps naive Linux user, I thought that if a
parent
Hi,
I'm trying to use the i-best software. Does anyone have experience from this, I
can't get it to work with my data.
Here is the code:
T1 - read.delim(file=S://SEDIM//Yvonne//2_5//T1.txt,col.names=
c(Dye/Sample_Peak, Sample_File_Name, Size, Height, Area_in_Point,
Area_in_BP, Data_Point,
Hello Kate,
Sorry I've only seen that thread. I've not been using A2R for some time,
and never really found the time to release it properly, in CRAN for
example. I'll try to allocate some time to that.
If you have the sufficient tools installed
(http://www.murdoch-sutherland.com/Rtools/),
Spencer Graves:
Bates' condemnation of R^2 has merit, but I would not go as far as
he did in the comment cited below (dated 13 Aug 2000). A standard
definition of R^2 is as follows:
R^2 = (1 - var(prediction error) / var(obs)).
I can name several different ways of getting a negative R^2
What you want is:
wf[[fl[1]]]
'fl' is a vector, so you should only be using a single '[' for indexing.
On Wed, May 7, 2008 at 3:20 AM, Dirkheld [EMAIL PROTECTED] wrote:
Hi,
I have a dataframe wf existing of a header with different labels and beneath
the values of those labels :
wf:
Forgot to send one copy to R help. Sorry
Megh Dal [EMAIL PROTECTED] wrote: Date: Wed, 7 May 2008 02:45:09 -0700 (PDT)
From: Megh Dal [EMAIL PROTECTED]
Subject: Re: [R] Solution of function
To: Berwin A Turlach [EMAIL PROTECTED]
Hi Berwin,
Thanks for having look on my problem. However on
Thanks, I had indeed several versions of R and removed the older ones and it
worked.
I suspect you have more than one version of R installed and are mixing
them up. Those symbols have been in package stats for quite a while.
Try starting R with --vanilla, and if that works, clean out
T1 - read.delim(file=S://SEDIM//Yvonne//2_5//T1.txt,col.names=
c(Dye/Sample_Peak, Sample_File_Name, Size, Height,
Area_in_Point, Area_in_BP, Data_Point, Begin_Point,
Begin_BP, End_Point, End_BP, Width_in_Point, Width_in_BP,
User_Comments, User_Edit))
T1 - subset(T1, Size 1000 Size 50)
Jorge Ivan Velez wrote:
I think the function could be better but try this:
# Function: M is your matrix and n MUST be an integer0
mat.pow-function(M,n) {
result-M
if(n1){
for ( iter in 2:n) result-M%*%result
result
}
else {result}
result
}
Hello!!
I would like create a window that has diferent element as:
http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/checkboxes.html
http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/radiobuttons.html
I know as make it, but I don´t know as I could (align the diferent elemnts
to left,
[EMAIL PROTECTED] wrote:
T1 - read.delim(file=S://SEDIM//Yvonne//2_5//T1.txt,col.names=
c(Dye/Sample_Peak, Sample_File_Name, Size, Height,
Area_in_Point, Area_in_BP, Data_Point, Begin_Point,
Begin_BP, End_Point, End_BP, Width_in_Point, Width_in_BP,
User_Comments, User_Edit))
T1 -
Beautiful! Thanks Jim.
DaveT.
-Original Message-
From: jim holtman [mailto:[EMAIL PROTECTED]
Sent: May 6, 2008 07:33 PM
To: Thompson, David (MNR)
Subject: Re: [R] list manipulation
The reason for the NULLs is that is the output of the lapply you are
executing. If you don't want to see
Hello everyone,
Very quick question: How should I cite the use of R in a publication ?
Thanks in advance.
Sebastien
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
citation()
To cite R in publications use:
R Development Core Team (2008). R: A language and environment for
statistical computing. R Foundation for Statistical Computing,
Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org.
A BibTeX entry for LaTeX users is
@Manual{,
Sébastien wrote:
Hello everyone,
Very quick question: How should I cite the use of R in a publication ?
Thanks in advance.
Sebastien
See the FAQs:
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Citing-R
This is also referenced in the startup banner in an R console:
...
R is a
Thank you Prof Ripley for your answer.
The characteristic function is the inverse Fourier transform of the
distribution function. The characteristic function of a normaly
distributed random variable is exp(-t^2/2).
The fft is a discrete Fourier transforn, not a continuous one.
This is
Hello
I have to solve a multiple linear regression. Most programs like Excel
or Mathlab only support 5-10 dimensions. Now I have installed CRAN and
I have no clue what to do next. At the moment I am entering my data
into an excelsheet (for quick copy- paste). The Y-array will be 20
columns
On Tue, May 6, 2008 at 11:44 PM, ronggui [EMAIL PROTECTED] wrote:
library(ggplot2)
(p- qplot(mpg, wt, data=mtcars))
What I am doing is to set color of the ticks to hide them.
grid.gedit(gPath(xaxis, ticks), gp=gpar(col=white))
It should be a better way to achieve the purpose. Thanks.
The last example in ?fisher.test is not a 2x2 table, in fact it uses levels
with a natural ordering similar to the original question. Why would this not
be applicable to the situation?
From: [EMAIL PROTECTED] [EMAIL PROTECTED] On Behalf Of David
Megh Dal wrote:
Hi,
I want to find solution of function : f(x,y) = x'Cx - a under constraints :
0 x,y p
0 x-y q
where a, p,q are given constants and x = (x, y) and C is a 2X2 matrix
(given)
Can anyone suggest me any R function to do that?
Not
Multiple linear regression is handled by the function lm() in the
default installation of R. This takes inputs as lm(y~x1+x2+x3).
If you're going to be using R regularly, there are several books which
cover the basic statistical analyses available in R (and then some),
including those by
PS Perhaps it'd still be really useful to be able to change text direction
in labels.
Hadley, what do you think?
Agreed. I've added it to my ggplot2 customisation to do list.
Thanks,
Hadley
--
http://had.co.nz/
__
R-help@r-project.org
First step would be to read the manual on the R site
R Data Import/Export describes the import and export
facilities available either in R itself or via
packages which are available from CRAN
Then if that does not solve the problem you need to
explain in detail what the problems are, preferably
hello,
i have been using cor.test() for calculating the correlation coefficient and p
values for some data. however, since the data consist of two dichotomous
sequences (actually just binary data), i understand that simply using the
pearson correlation is not sufficient. however, having done
Greg Snow [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
The last example in ?fisher.test is not a 2x2 table, in fact it uses
levels with a natural ordering similar to the original question.
Why would this not be applicable to the situation?
Apologies. Clearly I misunderstood the R
i usually import data from exel, using read.table or read.csv (which implies
that i have to save exel files as .txt or .csv)
JM
El Miércoles, 7 de Mayo de 2008 11:25, John Kane escribió:
First step would be to read the manual on the R site
R Data Import/Export describes the import and export
Abhijit Dasgupta wrote:
Multiple linear regression is handled by the function lm() in the
default installation of R. This takes inputs as lm(y~x1+x2+x3).
If you're going to be using R regularly, there are several books which
cover the basic statistical analyses available in R (and then some),
Yemi Oyeyemi wrote:
Hi everyone, please I'm having problem importing data from Stata and
excel. Help me out.
Thanks
You don't provide...
a) the code that you've tried
b) the error message that relates to the problem you are having
...without these people have little information
import stata data should be straight. take a look at foreign package
On Wed, May 7, 2008 at 9:30 AM, Yemi Oyeyemi [EMAIL PROTECTED] wrote:
Hi everyone, please I'm having problem importing data from Stata and excel.
Help me out.
Thanks
-
[[elided Yahoo
What do you mean by 'the sequence option'?
The authot of lm.ridge
On Wed, 7 May 2008, Rodrigo Briceño wrote:
Dear R users. I have a doubt about the use of the sequence option on
Ridge regression. I'm trying to understand the use of this option when
variables are highly linear correlated. I'm
Greetings R users, maybe there is someone who can help
me with this problem:
I define a function optim.fun and want as output the
sum of squared errors between predicted and measured
values, as follows:
optim.fun - function (ST04, SM08b, ch2no, a, b, d, E)
{
predR -
Hi Simon and all,
I'm pretty sure that you are correct about this. I think it is a
misconception to say that the fisher exact test is only for a 2 by 2 table.
It is presented that way in textbooks because, for a 2x2 table, it is easy
to perform. For larger tables, it becomes complex quickly due
Hi,
I am using lmer to analyze habitat selection in wolverines using the
following model:
(me.fit.of -
lmer(USED~1+STEP+ALT+ALT2+relM+relM:ALT+(1|ID)+(1|ID:TRKPT2),data=vdata,
control=list(usePQL=TRUE),family=poisson,method=Laplace))
Here, the habitat selection is calaculated using a
Nice to meet you Brian. The question is about the numbers that appears
after lambda (0,0, 0.1, 0.0001). I know that seq is a set of values used
for testing which value fits best. But I'm not sure if I need to put
whatever I think or what. I tried also with a set of 5 values and I get
an error.
Hi,
The unique function is easy to understand and use. Beyond that, I want to get
also the frequency of repetition of each individual row in a data frame
Let me explain with an example :
x-data.frame(a=c(1,2,3,1,2),b=c(2,3,4,2,3),c=c(10,20,30,10,20))
xu-unique(x)
We have,
x
a b c
1 1 2 10
2 2
Hi Ravi,
Try this:
x-data.frame(a=c(1,2,3,1,2),b=c(2,3,4,2,3),c=c(10,20,30,10,20))
cbind(unique(x),lapply(x,table)$c)[,-4]
a b c Freq
1 1 2 102
2 2 3 202
3 3 4 301
HTH,
Jorge
On Wed, May 7, 2008 at 12:11 PM, ravi [EMAIL PROTECTED] wrote:
Hi,
The unique function is easy to
ravi -
This may get you started
count.reps - function(df) {
hash - do.call(paste, c(df, sep = \r))
cbind(unique(df), Freq = unclass(table(hash)))
}
test - data.frame(a = rep(1:10, 2), b = rep(1:10, 2))
count.reps(test)
Best,
Erik Iverson
ravi wrote:
Hi,
The unique function is easy to
ravi wrote:
Hi,
The unique function is easy to understand and use. Beyond that, I want to get
also the frequency of repetition of each individual row in a data frame
Let me explain with an example :
x-data.frame(a=c(1,2,3,1,2),b=c(2,3,4,2,3),c=c(10,20,30,10,20))
xu-unique(x)
We have,
x
a b
Hello,
As I mentioned in the previous message we are developing solution for
wholesale companies to analyze their sales transactions by associative
rules. I would very much appreciated if the community could give us
some hint of what is a typical PC configuration of a professional
statist
I think I should be clear exactly what I want :
take following example :
a = b = seq(1, 5, by=500)
v = matrix(0, nrow=length(a), ncol=length(a))
for (i in 1:length(a))
{
for (j in 1:length(a))
{
d = c(17989*a[i], -18109*b[j])
v[i,j] = t(d) %*%
Hi,
I have a simple R script for printing arguments
cat printargs.R EOF
args = commandArgs()
print(args)
q()
EOF
To run this script, first I set PATH to ~/src/R-2.6.2 and execute
R --no-save printargs.R
I want to run this script with different R versions by
On 07-May-08 16:31:23, Erik Iverson wrote:
ravi -
This may get you started
count.reps - function(df) {
hash - do.call(paste, c(df, sep = \r))
cbind(unique(df), Freq = unclass(table(hash)))
}
test - data.frame(a = rep(1:10, 2), b = rep(1:10, 2))
count.reps(test)
That doesn't
?fixef
gets you the coefficient vector, from which you can make your predictions.
-- Bert Gunter
Genentech
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of May, Roel
Sent: Wednesday, May 07, 2008 7:23 AM
To: r-help@r-project.org
Subject: [R] predict lmer
Sorry, my reply below may be too terse. You'll need to also construct the
appropriate design matrix to which to apply the fixef() results to.
If newDat is a data.frame containing **exactly the same named regressor and
response columns** as your original vdata dataframe, and if me.fit.of is
your
Another way to produce the data frame is
subset(as.data.frame(table(x)),Freq0)
- Phil Spector
Statistical Computing Facility
Department of Statistics
Jeff and others:
Paul Mielke and Ken Berry at Colorado State University developed very fast
Fortran routines for doing Fisher exact tests (as well as some others like
Chi-square and Zelterman's statistics). Their book (Permutation Methods:
A Distance Function Approach, 2nd ed. 2007.
Try
get(paste(wf$,fl[[1]],sep=))
See ?get
Julian
Dirkheld wrote:
Hi,
I have a dataframe wf existing of a header with different labels and beneath
the values of those labels :
wf:
label1 label2 ...
0,450,21
0,100,45
I have a list
fl - c(label2,label3,..)
Isn't
Hi,
From R 2.6, I would like to update to R2.7. I would like to have some tips on
the recommended method of installing the latest versions of an entire list of
packages in R2.7 - i.e. all the packages that I have presently installed in
R2.6.
I am hoping that there is an easier method than
Hi,
I am trying to figure how to use dlm with constant terms
(possibly time-dependent) added to both equations
y_t = c_t + F_t\theta_t + v_t
\theta_t = d_t + G_t\theta_{t-1} + w_t,
in the way that S-PLUS Finmetrics does?
Is there any straightforward way to transform the above to
Hi Endre,
I'm not sure if my configuration is typical. I'm running XP SP2, 4GB RAM and
Intel Core 2-Duo Processor of 2.4 GHz with 200GB HDD. My R version is R
2.7.0 Patched.
HTH,
Jorge
On Wed, May 7, 2008 at 12:49 PM, Endre Domiczi [EMAIL PROTECTED] wrote:
Hello,
As I mentioned in the
Thank you for the rapid response. I have no trouble in following your
example, but I was specificly interested in calculating QAIC and QAICc from
a glm fitted with the family=quasibinomial option. This provides an
estimate of dispersion, but does not contain an explicit value for AIC or
Hi R People:
Is there a package for automatic differentiation, please?
thanks in advance,
Erin
--
Erin Hodgess
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
mailto: [EMAIL PROTECTED]
__
On Tue, May 6, 2008 at 12:32 PM, Mark Kimpel [EMAIL PROTECTED] wrote:
I am running 64-bit Ubuntu 8.04 and when I invoke rggobi the interactive
graph displays but R crashes. See my sessionInfo() and a short example
below. Ggobi and rggobi installed without complaints. Mark
sessionInfo()
Uninstalling and reinstalling ggobi via Synaptic solved the problem, at
least for the demo data mtcars. Rotation works fine. No crashes on exit.
Thanks for the good advice.
Mark
On Wed, May 7, 2008 at 2:12 PM, Paul Johnson [EMAIL PROTECTED] wrote:
On Tue, May 6, 2008 at 12:32 PM, Mark Kimpel
There are a few different methods discussed in the README to the batchfiles
distribution. Home page at:
http://batchfiles.googlecode.com
On Wed, May 7, 2008 at 12:51 PM, ravi [EMAIL PROTECTED] wrote:
Hi,
From R 2.6, I would like to update to R2.7. I would like to have some tips
on the
Hi
I would like to know how to pivot a table that will sum the number of plants
(a or b) for each location (x,y,z)
I have read on the listserve similar questions but which involve summing up
numbers rather than factors. I have also read about the R package reshape
on the listserve but wanted to
Is this what you want:
x
location plant
1x a
2x b
3x b
4y b
5y b
6y a
7z a
8z a
table(x$plant,x$location)
x y z
a 1 1 2
b 2 2 0
On Wed, May 7, 2008 at 2:56 PM, juanita choo [EMAIL
assuming what you have is a dataframe called x
try:
table(Plant=x$plant,Location=x$location)
Location
Plant x y z
a 1 1 2
b 2 2 0
or simply:
table(x$plant,x$location)
x y z
a 1 1 2
b 2 2 0
thanks
juanita choo wrote:
Hi
I would like to know how to pivot a table
Hi,
I am using coxph with weights to represent sampling fraction of subjects.
Our simulation results show that the robust SE of beta systematically
under-estimate the empirical SD of beta.
Does anyone know how the robust SE are estimated in coxph using weights?
Is there any analytical formula
What is the command for 2D fft? Or 2D convolution and or auto-correlation?
Thanks,
Paul
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and
Please check out these upcoming courses offered by XLSolutions
Corporation.
**Payments due after the class
R/Splus Advanced Programming course: May 29-30, 2008 in San Francisco -
4 seats left.
R/Splus Advanced Programming course: May 29-30, 2008 in Boston - 2
seats left.
I have just installed the 64 bit version of R, using yum. The version
is: 2.6.2-1.fc7.1.x86_64.
I installed zoo without any major problem and the same with quadprog (a
few warnings). However, when I came to install tseries I get the
following:
install.packages()
Warning in install.packages() :
Hi,
Does anyone know if R has a function that is similar to lag.plot but instead
of auto-correlation, it plots cross-correlation with lags?
Thanks,
--
Tom
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Yesterday I posted the following question to the help list. Thanks to John
Fox (copied below) who pointed out the solution.
Original question:
I have come across a result that I cannot explain, and am hopingthat someone
else can provide an answer. A student fitted a mixed model usingthe
On Wed, May 7, 2008 at 11:40 AM, Mark Kimpel [EMAIL PROTECTED] wrote:
Uninstalling and reinstalling ggobi via Synaptic solved the problem, at
least for the demo data mtcars. Rotation works fine. No crashes on exit.
Thanks for the good advice.
I am pretty sure what Paul is referring to is
Hi!
Try ?ccf
HTH,
Erin
On 5/7/08, tom soyer [EMAIL PROTECTED] wrote:
Hi,
Does anyone know if R has a function that is similar to lag.plot but instead
of auto-correlation, it plots cross-correlation with lags?
Thanks,
--
Tom
[[alternative HTML version deleted]]
On Wed, 2008-05-07 at 21:09 +0100, Michael Steinbeck-Reeves wrote:
I have just installed the 64 bit version of R, using yum. The version
is: 2.6.2-1.fc7.1.x86_64.
Could it be that you don't have the development headers for R? Try:
yum install R-devel
as root in a console and then try
One question that arises is: at what level is the prediction desired?
Within a given ID:TRKPT2 level? Within a given ID level? At the
marginal level (which Bert's code appears to produce). Also, there is
the question: how confident can you be in your predictions. This
thread discusses
Hi,
I'm new to R and this mailing list, so I will attempt to state my question as
appropriately as possible.
I am running R version 2.7 with Windows XP and have recently been exploring the
use of the function regsubsets in the leaps package in order to perform
all-subsets regression.
Date: Wed, 07 May 2008 13:03:14 -0400
From: Stoyanov, Tsvetan (MSCIBARRA) [EMAIL PROTECTED]
Importance: normal
Priority: normal
Thread-topic: dlm with constant terms
Thread-index: AciwZDz3Kl800OdmR4CQkEC8PUECZg==
Hi,
I am trying to figure how to use dlm with constant terms
(possibly
There may already be tools for gravity model APPLICATION, which is
simply an exercise in normalization. Perhaps if I describe the gory
details, someone will recognize the problem by another name.
A singly constrained gravity model is simply an exercise in normalizing
the rows of a matrix
Hi,
I'm using Sweave with xtable to generate Latex tables. My problem is that some
of the tables are too wide, and I would like to use the tiny font for these
tables, but not for the others nor for the text. I can't work out how to
accomplish that.
The Latex code below does what I want it
On 07/05/2008 7:34 PM, Peter Jepsen wrote:
Hi,
I'm using Sweave with xtable to generate Latex tables. My problem is that some of the tables are too wide, and I would like to use the tiny font for these tables, but not for the others nor for the text. I can't work out how to accomplish that.
R-users
E-mail: r-help@r-project.org
but I was specificly interested in calculating QAIC and QAICc from
a glm fitted with the family=quasibinomial option.
If you use family=quasibinomial(link = logit) in glm(),
the program will be:
function ()
{
xx - c(1,2,3,4,5,6,7,8,9,10)
yy -
hello all,
Quick question, how do I get the p value out of the anova?
Thanks,
Paul
pb-aov(as.numeric(diff[5,16:33]) ~ grF)
summary(pb)
Df Sum SqMean Sq F value Pr(F)
grF 3 2.7860e+10 9.2867e+09 4.2236 0.02534 *
Residuals 14 3.0783e+10 2.1988e+09
---
Signif.
Hi,
I tried to run SVD on a 500,000* 500,000 matrix and i get a message that it
can not allocate a vector of length 270 mb
doe snayone know how to solve this problem? any ideas on other softwares
where I can do this?
I appreciate your help
thanks
--
View this message in context:
Yea the anova object seems to be odd. It's not S4 so that's why I tried
originally the attr() funtion but
summary(pb)$Pr(F)
Error: unexpected '' in summary(pb)$Pr(
summary(pb)$Pr
NULL
summary(pb)@Pr(F)
Error: unexpected '' in summary(pb)@Pr(
summary(pb)@Pr(F)
Error: no slot of name Pr
Hi Paul,
Perhaps
# Data set
set.seed(123)
x=rnorm(100)
y=x+2*rnorm(100)
# ANOVA
AOV=anova(lm (y ~ x))
p1=AOV$Pr(F)[1]
p1
or
# ANOVA 2
AOV1=aov(y ~ x)
p2=unlist(summary(AOV1))[Pr(F)1]
p2
HTH,
Jorge
On Wed, May 7, 2008 at 8:47 PM, Paul Benton [EMAIL PROTECTED] wrote:
hello all,
Quick
What works, amazingly, is
summary(pb)$Pr(F) [note the quotes]
Abhijit
On Wed, 07 May 2008 18:17:09 -0700
H. Paul Benton [EMAIL PROTECTED] wrote:
Yea the anova object seems to be odd. It's not S4 so that's why I tried
originally the attr() funtion but
summary(pb)$Pr(F)
Error:
yep That worked Jorge. Thanks!
summary(pb)$Pr(F)1
NULL
summ-summary(pb)
class(unlist(summ))
[1] numeric
unlist(summ)[Pr(F)1]
Pr(F)1
0.02533637
names(summ)
NULL
Jorge Ivan Velez wrote:
Hi Paul,
Perhaps
# Data set
set.seed(123)
x=rnorm(100)
y=x+2*rnorm(100)
# ANOVA
Sorry to bother with this topic, but I'm still not clear about the
meaning of the value set that is used on lambda values. Is there a
correct way of doing that? My doubt is how to choose those 3 values
that appear in the example.
Thanks.
On Wed, May 7, 2008 at 9:29 AM, Prof Brian Ripley [EMAIL
http://www.stat.pitt.edu/stoffer/tsa2/Examples.htm
tom soyer wrote:
Hi,
Does anyone know if R has a function that is similar to lag.plot but
instead
of auto-correlation, it plots cross-correlation with lags?
Thanks,
--
Tom
[[alternative HTML version deleted]]
David Winsemius [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
raymond chiruka [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
i am trying to carry out a categorical data analysis but my problem
is that when in i use the chi squared test some of my expected
values are less than
Erin Hodgess [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
Is there a package for automatic differentiation, please?
Have you looked at Goedman, Grothendieck, Hojsgaard, and Pinkus' Ryacas?
From the Ryacas documentation:
Analytical derivatives of functions can be evaluated with the D()
No, computer algebra systems such as yacas do symbolic differentiation.
Automatic differentiation is a numerical technique that can be used to
find derivatives of functions that can be implemented as computer
programs, through successive uses of the chain rule on the computer code
itself. It would
Hi,
I have a data matrix in which there are 1000 rows by 30 columns. The first
column of each row is a numeric indicating the number of trips taken to a
particular location with location attributes in the following column entries for
that row.
I want to repeat each row based on the number
Thanks to all!
Sincerely,
Erin
On Wed, May 7, 2008 at 10:35 PM, Simon Blomberg [EMAIL PROTECTED] wrote:
No, computer algebra systems such as yacas do symbolic differentiation.
Automatic differentiation is a numerical technique that can be used to
find derivatives of functions that can be
Hi Marion,
Try this:
set.seed(123)
mydf=data.frame(trips=rpois(10,5), matrix(rnorm(10*5),ncol=5))
mydf
sapply(mydf[,-1],rep,mydf[,1])
HTH,
Jorge
On Wed, May 7, 2008 at 11:41 PM, [EMAIL PROTECTED] wrote:
Hi,
I have a data matrix in which there are 1000 rows by 30 columns. The
first
On Wed, 7 May 2008, Rodrigo Briceño wrote:
Sorry to bother with this topic, but I'm still not clear about the
meaning of the value set that is used on lambda values. Is there a
correct way of doing that? My doubt is how to choose those 3 values
that appear in the example.
I think you have not
94 matches
Mail list logo
