I do not know exactly what you are looking for but it seems that you are
writing the column names (which become row names) when transposing the
data. So to fix this try using write.table(..., sep=,, row.names=F)
Jean
Daniel Gruner wrote:
Hello:
I need to take a species-sample matrix and
You might want to read
?scan and pay attention to what= argument
Jean
mark salsburg wrote:
How do I manipulate the read.table function to read in only the 2nd
column???
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing
If you know what the dependent variable is called and is the same for
all regressions see
?formula
Or you can do glm(ABC~. , ...,data=dat)
PS. the way you are calling the formula you are taking the first element
in the first column and the first element in the second column. So that
is not
If believe that merge is trying to put first whatever cells that are
nonempty. For example if you instead did
tmp2 - data.frame(col1 = factor(c(C, D, E, F,A), levs),
col2 = 1:5)
tmp2
col1 col2
1C1
2D2
3E3
4F4
5A5
merge(tmp2, tmp1, all.y = TRUE,
Hi
I do not get why the symbols function produces warnings when axes=F is
added. The following example illustrate this
symbols(0,10, inches=T, circles=1, axes=F, xlab=, ylab=)
Warning message:
parameter axes could not be set in high-level plot() function
I augmented symbols and added the
Hi,
In my experience pdftotext did not do a very good job at this because it
screws up the formatting of tables. This of course depends on what
program the pdf document was originally constructed with. What I found
most appealing is the use of cut and paste into xemacs or emacs and use
M-x
I do not totally understand your question as well. You seem to want a
descriptive statistic about a unitary number. What is the sd of a number?
or any other descriptive statictic. Maybe you mean for the columns or rows
or it could be that these are t-stats or z-stats that you need to get
p_values
If your loop is from 1:n then you can do the following. suppose you call
the resulslts results 1:n using assign or something. so like this
for(i in 1:n){
assign(results_, i, sep=), lm(bla bla))
save(get(
for(i in 1:10){
+ temp - paste(results_, i, sep=)
+ assign(temp, rnorm(i))
+
I am sure all this work but If you want exaclty the output to be the way
you mentioned do this
temp-read.table(yourfile, as.is=T, header=T)
temp1-temp[, 1:3]
temp2-temp[, c(1,2,4)]
colnames(temp1)[3]-X
colnames(temp2)[3]-X
temp3-merge(temp1, temp2, all=T)
temp3$type-toupper(substr(temp3$X, 1,2))
You can turn the identity vector of the firms into a factor and do lm
Jean
On Mon, 5 Sep 2005, Tobias Muhlhofer wrote:
Hi, all!
Anyone know an easy way to specify the following model.
Panel dataset, with stock through time, by firm.
I want to run a model of y on a bunch of
is the
factor of all firm identifiers, I could just go
lm(y ~ x + firm)
and that will implicitly include a dummy for each level of factor firm,
thus making this a fixed effects (aka LSDV) model?
T
Jean Eid wrote:
You can turn the identity vector of the firms into a factor and do lm
want to call at the end, in a directory
# of its own.
patchCode(path=hot/onFinally/, removeAfter=TRUE)
# End of your batch code
Cheers
Henrik
Duncan Murdoch wrote:
Jean Eid wrote:
This is probably a weird question but I need to know if there is a way...
I run an R batch job without
just one more thing. Why are you defining the function over and over
again...
On Wed, 31 Aug 2005, Hathaikan Chootrakool wrote:
I was wondering why this loop doesn't work!
for (i in 1:k)
fnTr[i] - function (p) 0.5* sum ( n*log(2*pi) - log(sd(i)^2)
This is probably a weird question but I need to know if there is a way...
I run an R batch job without saving the variables at each step to the
disk. Is there a way to invoke another session of R and link it to the
same environment for read only.
The problem is that I am running optim with
why not write.table with sep=\t
On Sat, 27 Aug 2005, Duncan Golicher wrote:
Could anyone help with what should be a simple task? I have data as a
fixed format (fortran) table. I have no trouble getting it into R using
read.table. Each column is separated by a space, including the first
Although it is not as fancy as all other responses, I usually just print
the value of i first and see where it stopped. Of course this assumes you
it is stored in the main env. If it is inside a function and it is failing
I usually use the - operator to get it to the main env. and print it.
Here
for website will be greatly appreciated
Thank you
Jean Eid
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
bates
and you should be able to find a copy.
Jean Eid wrote:
Dear all,
My question is concerning the line
This is adequate for small files, but for anything more complicated we
recommend using the facilities of a language like perl to pre-process
the file.
in the import/export
Thank you Gabor,
Jean
On Tue, 16 Aug 2005, Gabor Grothendieck wrote:
On 8/16/05, Jean Eid [EMAIL PROTECTED] wrote:
Dear all,
My question is concerning the line
This is adequate for small files, but for anything more complicated we
recommend using the facilities of a language like
can you give the line you are using to import the data.
On Sun, 14 Aug 2005 [EMAIL PROTECTED] wrote:
Well, you were right. When I check with str() or summary(),
the data frame appears as having a factor for that column.
However, if I ask is.factor(C2), I get FALSE, and if I ask
sapply(nrow(matrix):1, function(x) matrix[x,])
On Mon, 8 Aug 2005, Jake wrote:
Quick question: how can I reverse the order of the rows in a matrix?
i.e. make the last row first and the first row last, etc.?
__
R-help@stat.math.ethz.ch mailing
You can do the following without resorting to a hard coded loop
sapply( paste(q35, letters[1:grep(r, letters)], sep=), function(x)
ifelse(temp[, x]%in%1,temp[, sub(5, 3, x)],NA)
as the following example shows
temp - matrix(sample(c(0,1), 360, replace=T), nrow=10)
colnames(temp) - c(paste(q33,
if labda is the elements of the vector and you know what kexp is , you can
use apply
apply(your_vector, 1, function(x) 1- exp^(kexp^(-x)))
HTH
Jean
On Wed, 3 Aug 2005, Rangesh Kunnavakkam wrote:
I have a large vector of around 12597 elements and I wish to calculate
p-value for each
Here is one way this can be done
do.call(rbind, by(dat, list(dat$fam) ,function(x) {
+ if(NROW(x)1) return(x[which.max(x$wt),])
+ else return(x)}
+ ))
and it returns
fam wt keep
1 1 1.01
2 2 1.01
3 3 0.61
4 4 0.41
hth,
On Mon, 1 Aug 2005, Chris Wallace wrote:
I
I just installed it on a Debian 2.6.8.1 and the following R
platform i386-pc-linux-gnu
arch i386
os linux-gnu
system i386, linux-gnu
status
major2
minor0.1
year 2004
month11
day 15
language R
By the way are you using apt-get install vcd. and if so why? just use
You can do the following (don't know it this is the most efficient way but
it works)
temp-read.table(your file to read the data, header=T)
temp1-table(temp)
plot(temp$x, temp$y, cex=0)
text(as.numeric(rownames(temp1)), as.numeric(colnames(temp1)), temp1)
HTH
On Wed, 13 Jul 2005, Kerry Bush
Hope everyone id doing great ..
Just need some clarification over the limit of write.dta. I have some
coauthors that use stata and I need to send them my data in .dta format.
the data.frame is 41706x229 and I get the following
Error in write.dta(Panel, file = STATADATA/Panel.dta, version = 7) :
did not clarify the system and such ...
platform i386-pc-linux-gnu
arch i386
os linux-gnu
system i386, linux-gnu
status
major2
minor0.1
year 2004
month11
day 15
language R
The foreign package version is 0.8-4
On Mon, 6 Jun 2005, Jean Eid wrote:
Hope
Thank you all for pointing out the name of the variables.. I do apologize
for not catching that,..
have a good day,
Jean
On Fri, 20 May 2005, Prof Brian Ripley wrote:
On Fri, 20 May 2005, Jean Eid wrote:
Hi all, hope you having a nice day,
I ahve this weird results with identical
Hi all, hope you having a nice day,
I ahve this weird results with identical (probably I am not understanding
correctly what it does ...)
I have these two data frames and I issue :
identical(temp, temp1)
[1] FALSE
However, these data frames are Nx2 and when I issue:
identical(temp[,2],
Why can't you just solve x1 in terms of x2 and plug it in..
ie. min g(const-x2, x2, x3)
Jean
On Fri, 20 May 2005 [EMAIL PROTECTED] wrote:
Hello,
I've got to compute a minimization equation under an equality constraint
(Min g(x1,x2,x3) with x1+x2=const). The Constroptim function does not
?table
On Fri, 20 May 2005, Wladimir Eremeev wrote:
Dear Ravi,
From a vector, I want to get the unique values and the counts of these
unique values in the vector. For example, x-c(2,1,2,1,4,2,1,4,1,1)
try
hist(x,plot=FALSE,breask=unique(x))$counts
[1] 5 3 0 0 0 2
--
Best regards
and get its expectation and match that to data. I thought of doing a
bootstrap but it will take so much time. I guess if this is the only way,
then it has to be done.
Jean
On Wed, 23 Mar 2005, Spencer Graves wrote:
Have you considered bootstrap or Monte Carlo?
spencer graves
Jean Eid
Hi AlL,
I ahve this problem that my objective function is discontinous in the
paramaters and I need to use methods such as nelder-mead to get around
this. My question is: How do i compute standard errors to a problem that
does not have a gradient?
Any literature on this is greatly appreciated.
do you mean n choose k which is a built in function see
?choose
On Sat, 19 Mar 2005, Mario Morales wrote:
En español (In Spanish)
Necesito calcular la en numeros de combinaciones de n cosas
tomando k al tiempo.
Como hago eso en R ???
Yo escribí mi propia función pero pienso que de
try and read the ess manual a bit and look at the folowing reference card
http://stat.ethz.ch/ESS/refcard.pdf
On Tue, 8 Mar 2005, Mahdi Osman wrote:
Hi all,
I have got a dataframe in coma delimted text format. My ESS and R processes
are working well and active. I can read R help files and
I have this somewhat large data set that is given to me in a fixed width
format. The file itself is already 100MB (Maybe R can actually handle
this but I am trying to gain some experience in postgres and RODBC). I
am using postgres to preprocess the file and connect to the database through
the
I have the same problem as well. I installed atalas3-sse2 and still hasd
the same problem. Here's the log of messages
* Installing *source* package 'quantreg' ...
** libs
g77 -fPIC -g -O2 -c akj.f -o akj.o
g77 -fPIC -g -O2 -c boot.f -o boot.o
g77 -fPIC -g -O2 -c boundc.f -o boundc.o
g77
that is because your first object is a data.frame but when you transposed
it you turned it into a matrix. so doing
mat1 - data.frame(matrix(rnorm(20*1532), ncol=20))
mat2 - t(mat1)
dim(sample(mat1, 5, replace=T))
[1] 15325
dim(sample(mat2, 5, replace=T))
NULL
length(sample(mat2, 5,
It seems that sample picks columns when the object is turned into a
data.frame. I do not knoe why it is doing that
Is this something that was meant and not documented or something?
Jean
On Fri, 11 Feb 2005, Adaikalavan Ramasamy wrote:
See below.
On Thu, 2005-02-10 at 14:53 -0800, T.
?order
On Fri, 11 Feb 2005, Luis Ridao Cruz wrote:
R-help,
Long ago I asked how to sort a data frame according to one of the
columns.
But Iwas wondering ho I could double sort , i.e, sort according to
one column and another one afterwards.
It is quite easy in Excel but I am unable to
I have a problem either understanding what loess is doing or that loess
has a problem itself.
As the x-axis variables become more concentrated on a particular point,the
estimated loess tends to zero. the examples below show what i am
talking about, why is that? my intution tells me
that it
list1 - list(x=c(1,2), y=c(3,4), z=c(4,5))
matrix(unlist(list1), nrow=length(list1), byrow=T)
HTH
On Thu, 3 Feb 2005, michael watson (IAH-C) wrote:
Hi
Sorry to ask such a basic question. I have a list, each element of
which is a vector of two values. What I actually want is a matrix with
Why do you use png? I think that yu need X11 to save a png file. why not
use a postscript?
This is what ?png says
R can be compiled without support for either or both of these devices:
this will be reported if you attempt to use them on a system where they
are not supported. They will not be
LS2 9JT
tel: +44 113 343 1596
fax: +44 113 343 6716
mail: [EMAIL PROTECTED]
On Wed, 2 Feb 2005, Jean Eid wrote:
Why do you use png? I think that yu need X11 to save a png file. why not
use a postscript?
This is what ?png says
R can be compiled without support for either
as.vector(coef(nlsobj))
as.vector(resid(nlsobj))
On Wed, 2 Feb 2005, Chung Chang wrote:
Hi,
I have a question about how to get the residuals and estimations
of the parameters in my program.
For example,
x-c(1,2,3,4,5,6)
y-c(2.9, 1.24, 1.71, 2.989358, 1.455979, 1.4)
nls(y ~
I have this awkward problem with trellis (lattice). I am trying to
generate some plots through loops but the .eps file is empty. When I
generate them in a list and print them outside the loop all is fine. this
is an example below:( nothing shows up in foo.eps, but all show up in
foo1.eps)
R
reshape(XX, direction=long, varying=list(c(one1, one2, one3),
c(two1, two2, two3)), v.names=c(one, two), times=c(1, 2,
3))
# I get the correct output.
I there a way not to generate the XX$two3 above
Thank you
Jean
On Wed, 26 Jan 2005, Thomas Lumley wrote:
On Wed, 26 Jan 2005, Jean
Hi,
Are there any plans to do multi line commenting? like /*...*/
Jean
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
try
sink(filename)
cat(paste(gsub('/', '', 'c:/dir1/dir2/file.ext'),\n,gsub('/',
'', 'c:/dir1/dir2/file.ext'),\n ))
sink()
the idea is to use cat,
instead of paste you can do a for loop where you cat the string into the
sinked file
On Wed, 19 Jan 2005, Joerg Klausen wrote:
Dear
when tcltk is called and a C-c C-c is issued a message of Illegal
instruction is issued and the R process is killed. To reproduce please
do the following
require(tcltk)
C-c C-c
It will say Illegal instruction
I use a debian 2.4 kernel with R 2.0.1 tcltck version 2.0.1.
I use tcltk to have gui
Can't you turn the lists into data frames issue unique and
force them back to lists. Here's the code:
L - list(c(a1,a3,a4), c(a1,a4,a5), c(a1,a5,a6))
M - list(c(a1,a3,a4), c(a2,a4,a5), c(a1,a5,a6), c(a7,
a1, a4))
LL - as.data.frame(I(L))
MM - as.data.frame(I(M))
X -
Here's tow examples that do what you want (the way I understand them).
layout(matrix(c(1,0,1,0), byrow = TRUE))
plot(10:100)
par(xpd=NA)
co - par(usr)
legend(co[1], -10, legend=c(Hello long time no see,\n , and My day
was great,\n, How was yours), text.width=95, col=red,,
text.col=c(red,blue,
You can do the following also
X - matrix(c(1,2,1,2,1,3,1,4), ncol=2)
Y-unique(X)
Y[which.max(apply(Y, 1,function(i) sum(apply(matrix(1:nrow(X)),
1,function(x) identical(X[x,], i[1:2]),]
I do not know what your strategy is when there are multiple maxima i.e two
different rows appear at
I think the polr is in MASS and indeed has a probit method.
Jean,
On Fri, 3 Dec 2004, Charlie Sprenger wrote:
Hello All,
I'm trying to run a multinomial probit on a dataset with 28 data
points and five levels (0,1,2,3,4) in the latent choice involving
response variable.
I downloaded the
wouldn't it be return[2:NROW(return)]*marketcap[1:(NROW(return)-1)]
(note that return is also a function in R so maybe you should stay away
from calling your variable return. Anyways if you have a data frame you
can add another variable to it but attach an NA to the first element
(since you are
see ?assign and ?get
i,e instead ob blank below have something like assign(paste(v, i,
sep=), a[i:N])
and if you need to loop over calling them say get(paste(v, i, sep=))
of course typing v1 will call the variable...
Jean
On Tue, 30 Nov 2004, Tobias Muhlhofer wrote:
I am trying to define a
sub(.*/, , abc/foo)
[1] foo
Jean
On Thu, 18 Nov 2004, Vadim Ogranovich wrote:
Hi,
I am looking for a way to extract an interesting part of the match to
a regular expression. For example the pattern [./](*.) matches a
substring that begins with either . or / followed by anything. I am
You can use setdiff if you only need the unique values of a that are not
in b. If you want all values you can use
a[a%in%setdiff(a,b)]
There are also intersection, union etc...
see
?setdiff
On Thu, 11 Nov 2004, Alexander Sokol wrote:
Hello,
I have the following very simple problem:
Say I
Now I understand,
R gives numbers zero to about 6 digits and Stata gives zero to about 30
digits. The intercepts are the same in both packages.
Thank you,
Jean,
On Thu, 11 Nov 2004, Thomas Lumley wrote:
On Thu, 11 Nov 2004, Jean Eid wrote:
Thank you Thomas for your answer
Dear All,
I have been struggling to understand why for the housing data in MASS
library R and stata give coef. estimates that are really different. I also
tried to come up with many many examples myself (see below, of course I
did not have the set.seed command included) and all of my
`random'
Dear Thomas,
Where you also able to replicate the second example? (the exaample
that I turned the housing data into numerical variables) That is the one
that my estimates differ.
Jean,
On Wed, 10 Nov 2004, Thomas Lumley wrote:
On Wed, 10 Nov 2004, Jean Eid wrote:
Dear All,
I have been
I am trying to automate an outide of the plot region table that has the
correlation coef. of the various variables in the plot. I am currently
using bquote to get the greek/latex representation of $\rho$. what I want
so what I have is the following
rhoa - cor(x, data1$No.Msa.Hosp)
rhon -
:
On Tuesday 19 October 2004 17:12, Jean Eid wrote:
Is there any function like par(new=T) for lattice. I want to plot a
histogram in percentages on the right hand side and also superimpose the
densityplot with its density scale on the lhs. so far I am only able to do
this
histogram( temp[,2]~ temp
I figured out both questions,.. The idea is to set them with
trellis.par.set(). However, I still would like to know if there is a
tutorial for lattice package.
Jean,
On Thu, 21 Oct 2004, Jean Eid wrote:
Deepayan,
Thank you so much,... works like a charm. However, I have two more
questions
If I understand you correctl, you just want to call the point of (0,10) to
be (0,6). if so have a look at text() function text(0,10, pos=2) or
something like that will do. The reason why you are getting the error
because your 0:10 vector is larger than 0:6 vector.
Hope this helps,
P.S. why are
Is there any function like par(new=T) for lattice. I want to plot a
histogram in percentages on the right hand side and also superimpose the
densityplot with its density scale on the lhs. so far I am only able to do
this
histogram( temp[,2]~ temp[,1],nint=100,type=desnity,
xlab =
I have been having problems with these two 'libraries' since I installed
2.0.0.
I have built a package with couple of functions so that I can load it at
startup every time R is booted. The problem is that I have the following
error every time I call the library
Loading required package: tclk
: tclk
in a DESCRIPTION file for one of your packages.
On Tue, 12 Oct 2004, Jean Eid wrote:
I have been having problems with these two 'libraries' since I installed
2.0.0.
I have built a package with couple of functions so that I can load it at
startup every time R is booted. The problem
there does not seem to be a package tcltk on CRAN for 2.0.0.
I have successfully installed the same package for 1.9.1.
In essence I require the package for a GUI interface to setwd.
All work fine with 1.9.1.
This is on a Linux Debian unstable kernel 2.4.20
The build version of R 2.0.0 (issuing
Is there a dataload utility for linux. The link in genstat is down but I
managed to find the utility at:
http://gurukul.ucc.american.edu/econ/gaussres/UTILITYS/DATALOAD.HTM
but this is a dos/windows version.
Thank you
Jean
__
[EMAIL PROTECTED]
to delete
all characters before a particular pattern.
If anyone has some help file for this, it will be greatly appreciated.
Jean Eid
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http
Thank you all for the help, specially Gabor that is exactly what I needed.
A few examples that do the same thing is very helpful in understanding the
structure of the call.
Thank you again,
Jean
On Thu, 23 Sep 2004, Gabor Grothendieck wrote:
Jean Eid jeaneid at chass.utoronto.ca writes
you might want to turn your data into a matrix. You get much much faster
for loops doing that.
Jean,
On Fri, 6 Aug 2004 [EMAIL PROTECTED] wrote:
Dear R-helpers,
I'm dealing with large datasets, say tables of 60 000 times 12 or so, and
some of the functions are (too ) slow and I'm therefore
I am trying to put legends underneath the plot (in the outer margins). Is
there an easy way to do this. I have been tinkering with split..screen but
I could not make it work.
Thank in advance
Jean
__
[EMAIL PROTECTED] mailing list
You have to have postscript(foo.eps, ..) before you issue the plot
commands and isuue a dev.off() at the end.
If you want to copy the plot graph to a postscript driver use
dev.copy2eps(file=foo.eps)
?dev.print
?postscript
Hope this helps.
On Fri, 30 Jul 2004, marzban wrote:
Hi,
I'm
, it will crash and on reboot it will check
the disk for errors.
I have also used gc() after I removed everything in the env. This did not
make a difference.
can someone guide me through a solution.
Thank you
Jean Eid
University of Toronto.
__
[EMAIL
you
Jean Eid
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
I have done that with no problems...
Jean,
On Tue, 18 May 2004 [EMAIL PROTECTED] wrote:
I use linux debian testing for which the latest debianized version of R
is 1.8.1. Therefore I installed:
r-base-core_1.8.1-0.cran.1_i386.deb
r-base-dev_1.8.1-0.cran.1_all.deb
Dear All,
I am trying to solve a Generalized Method of Moments problem which
necessitate the gradient of moments computation to get the
standard errors of estimates.
I know optim does not output the gradient, but I can use numericDeriv to
get that. My question is: is this the best function to do
something better, but this
has produced acceptable answers for me in several cases, and I did not
push it beyond that.
hope this helps.
Jean Eid wrote:
Dear All,
I am trying to solve a Generalized Method of Moments problem which
necessitate the gradient of moments computation to get
I just did this yesterday. you do not have an ODBC driver ( as it says in
the log). you need to (as root) use apt-get install unixodbc and apt-get
install unixodbc-dev.
The unixodbc package alone will not work, you need to install unixodbc-dev
(development),
Hope this helps.
I have a question
you can do this,
tmpmatrix - matrix(unlist((lapply(as.data.frame(x[, 3:ncol(x)]),
function(t)
return(rep(t[!is.na(t)], 9))) )), byrow=F, ncol=8)
matrix - cbind(matrix[, 1:2], tmpmatrix)
On Thu, 22 Apr 2004, Luis Rideau Cruz wrote:
Hi all,
I want to fill in this matrix vectors (by column)
Hi,
Is there a function that does percentile-percentile plot. I do not mean
the qqplot. I need to plot the percentiles rather than points themselves.
I am hoping for a plot that tells me that the x percentile of one data set
corresponds to the y percentile of the other. for example a point on
Thank you Duncan,
On Tue, 6 Apr 2004, Duncan Murdoch wrote:
On Tue, 6 Apr 2004 13:33:58 -0400, Jean Eid
[EMAIL PROTECTED] wrote :
Hi,
Is there a function that does percentile-percentile plot. I do not mean
the qqplot. I need to plot the percentiles rather than points themselves.
I am
When I originally posted my question under ( C-c C-c does not kill data
editor), I thought that it was a problem with ESS rather than R. However,
I have tried to do the same from a command line R (R 1.8.1 (2003-11-21).
on a debian testing i386 machine) and I got the same problem.
Specifically I
Dear all,
This might not be the best place to post this (maybe the ESS help is
better) but I will try anyways.
I am running ESS, Xemacs, R 8.1 on a debian testing dist (i386 arch). when
I envoke the edit(data) command, a data editor appears and all is fine.
However executing C-c C-c will only
I don't know what these files are. so depending on whether you want to
call them with 400 different names or just have one data set for the 400
files. but in either case you can do a for loops on the directory. i.e.
put the 400 files in a seperate directory and setwd(to that directory)
do this
sum(diag(X))
On Mon, 23 Feb 2004, li xian wrote:
How to get the sum of the diag of matrix?
Thanks!
-
[[alternative HTML version deleted]]
__
[EMAIL PROTECTED] mailing list
Dear All,
I have the following data.frame
`data.frame': 1563 obs. of 4 variables:
$ Model :Class 'AsIs' chr [1:40] Astro Astro Astro Astro
$ Make :Class 'AsIs' chr [1:40] Chevrolet Chevrolet
$ Production : num 11219 12384 1082 5409 5458 ...
$ date :`POSIXlt', format:
I want to be able to exctract a matrix from a data frame that contains repeated
measurements for individuals. i.e I want to exctract the second observation for each
individual in the study.
Is there a way of doing this?
I guess what I am trying to ask is if there as a multidimensional version
Dear all,
I am looking for a function that can scatter plot the residuals obtained
from a longitudinal model i.e. plot e_{i,j} e_{i,k} for
all j k = 1,..n ( I have 7 observations for each subject). something
similar to the pairs() function.
I can do it the long way by constructing a residual
I am working with a model that I have to estimate a nonparametric
function. The model is partial linear i.e.
Y=X$\beta$ + f(z) + $\epsilon$
I am using the ' double residual methods' Robinson (1988) Speckman (1988)
where I estimate a nonparametric function for each of the parametric
variables in
Please do forget my questions as they are really trivial and I do not kno
what I was thinking of.
Thank you thomas for clarifying my cloudy head today.
jean,
On Mon, 22 Sep 2003, Jean Eid wrote:
I am working with a model that I have to estimate a nonparametric
function. The model is partial
the xtable library to
turn output into latex tables and do not want to manually delete each and
every single design word in the coefficients name vector.
Thank you so much for any feedback,
Jean Eid
__
[EMAIL PROTECTED] mailing list
https
to change the configure file without making sure I am doing the right thing,
Any help is appreciated,
Jean Eid
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
: Martin Maechler [EMAIL PROTECTED]
To: Jean Eid [EMAIL PROTECTED]
Cc: [EMAIL PROTECTED]
Sent: Tuesday, July 08, 2003 12:04 PM
Subject: Re: [R] readline and R installation
One way is setting the (or adding to an existing one)
environment variable LD_LIBRARY_PATH to your directory
(.../local/R
I have been struggling to find some informaation on what lm exactly does.
I know it uses the QR decomp. However, I was recently faced with a
somewhat badly scaled matrix and summary(lm) said
Coefficients: ( 4 not defined because of singularities)
does anyone know how lm chooses these 4
NA NA
Jean does make reference to summing parts of a matrix, apparently in a
loop, so it might help to know what's being done besides the column sums.
Regards,
John
At 09:18 PM 6/14/2003 -0400, Jonathan Baron wrote:
On 06/14/03 20:51, Jean Eid wrote:
Dear R users,
I am looking
* 202 matrix.
thanks,
Jean Eid
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
100 matches
Mail list logo
