Seth Pruitt wrote:
To All,
I am using the BRugs package. In running the meta file BRugsFit() with a
syntactically correct model .txt file, I see the message:
Error in FUN(X[[1]], ...) : .C(..): 'type' must be real for this format
Hi,
can you send me in a private message a reproducible
Dear fellow R-users,
I wish to analyse a lattice of presence-absence data which are spatially
autocorrelated.
For normally distributed errors I used gls {nlme} with the appropriate
corStruct-method.
Is there any method for other families (binomial and poisson)?
A method that look suitable to
apply transorm your data.frame into a matrix which can contains only one
type of values (character or numeric or logical).
data.frame are list, so they can contain various kind of data, and a
solution to your problem is given by lapply:
b=lapply(df, function(r) print(class(r['a'])))
[1]
doubleEm - function(p1,p2,p3,p4){return(p1 *p1,p2 * p2, p3 * p3, p4 * p4)}
suppressWarnings(doubleEm(1,2,3,4))
[[1]]
[1] 1
[[2]]
[1] 4
[[3]]
[1] 9
[[4]]
[1] 16
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Paul Roebuck
Sent: Wednesday, 20 July 2005
dear R-tists,
i want to graph information for the German Federal States (Bundeslaender) using
the maps package. unfortunately there is no maps for the German Bundeslaender.
does anyone have an idea / a source where to get map data that can be used in
the maps package that graphs structures
Hi All,
I want to print a square matrix of 7000 x 7000 into a text file. But I
got a error after few hours of computation...
write.table(MyDistMxDF, file = temp.csv, sep=,, quote=F)
*** malloc: vm_allocate(size=8421376) failed (error code=3)
*** malloc[2889]: error: Can't allocate
TaeHC == Tae-Hoon Chung [EMAIL PROTECTED]
on Tue, 19 Jul 2005 10:26:01 -0700 writes:
TaeHC Hi, All; I tried to use library mclust in 64-bit
TaeHC compiled R 2.0.1 but failed.
you mean package mclust
{a package is installed in a library full of packages}
TaeHC Installation
Hi Adrian,
give a look to these links concerning R DM:
http://www.togaware.com/datamining/survivor/R.html
http://sawww.epfl.ch/SIC/SA/publications/FI01/fi-sp-1/sp-1-page45.html
http://www.stats.ox.ac.uk/~vos/DataMining.html
Hoping I helped you.
Regards,
Vito
secretario academico FACEA
Please use a current version of R: they are much better at this!
write.table was rewritten in R 2.1.0 to use *much* less memory.
There _is_ a `R Data Import/Export Manual' that dicsusses this.
write.table was designed to write `tables' (data frames) not matrices.
write.matrix (package MASS) does
It probably is not a problem to leave hyperthreading on: we found little
performance difference on a P4 either way.
The Windows task manager is misleading, as `50%' is about as much as a
P4-class processor with hyperthreading can actually deliver.
On Tue, 19 Jul 2005, Lukasz Komsta wrote:
x - as.Date(21-07-2005,%d-%m-%y)
is base R, and gives you a object of class Date.
days() is not part of R, and I guess from your subject is part of chron.
Its help page does not say what it actually does, but my guess is that it
finds days of the month, *for a dates object not a Date object*.
Chun-Ying Lee [EMAIL PROTECTED] writes:
Hi,
We used known Vm and Km to simulate the data set (time, Cp) without
adding random error in there. Yes, the line looks like very close
to a straight line. But why can't we obtain the correct values with
fitting process? We used optim first
On Tue, 19 Jul 2005, Marc Schwartz (via MN) wrote:
For the TAB delimited columns, adjust the 'sep' argument to:
read.table(data.gct, skip = 2, header = TRUE, sep = \t)
The 'quote' argument is by default:
quote = \'
which should take care of the quoted strings and bring them in as a
My guess is that you have two data points which are very close together.
The normal way out is to included a nugget effect in your model, but it
can also be appropriate to leave out one of a pair of points.
On Tue, 19 Jul 2005, orkun wrote:
hello
I try to build DEM using Krig function of
hi all,
I need to normalize a function, did something exist in R who can do it for me,
instead of integrate the function then divide it by the result?
thks, i'm sorry i didn't found any information in the documentations and
statistic vocabulary in english is a pain for me.
Hi,
I am newbie in R. I have R installed on Solaris machine and I am connected
to Solaris machine from Windows using telnet. I want to create png image
file using R. But when I issue R command from telnet:-
png(test.png)
Following error occurs:
Error in X11(paste(png::, filename, sep = ),
What does class(r1) give you ? If it is data.frame, then try
exp( diff( log( as.matrix( df ) ) ) )
BTW, I made the assumption that both x and y are positive values only.
Regards, Adai
On Tue, 2005-07-19 at 16:30 +0100, Gilbert Wu wrote:
Hi Adai,
When I tried the optimized routine, I
Hallo
I am not sure what you want to achieve.
your factor has 3 levels but with only 2 different labels
hb
[1] 1 1 1 1 1 1 2 2 2
Levels: 1 1 2
but
str(hb)
Factor w/ 3 levels 1,1,2: 1 1 1 2 2 2 3 3 3
so you gave only one label to level 1 and 2. You can give the same
label to any level you
Dear members,
I have numerous arrays that are organised in a list.
For example, suppose I have 2 arrays in a list called
alist
alist - list(array(rpois(12,5), 6:8) ,
array(rpois(15,5), 10:12))
with array dimnames
dimnames(alist[[1]]) - list(LETTERS[1:6],
paste(namesd, 1:7, sep=),
Those wondering what gtc file stands for, might be interested in [1].
The original poster can see if the package 'ctc' [2] supports reading in
this format but I think Prof. Ripley's solution works too.
[1]http://www.broad.mit.edu/cancer/software/genepattern/tutorial/gp_tutorial_fileformats.html
You are probably not running a X server on your Windows machine, or if you
are, you are not forwarding X connections. Even then, this is likely to
be unsatisfactory unless you have a very good X server and a fast
connection.
Please read the help page for png, and note the See Also section
Petr Pikal [EMAIL PROTECTED] writes:
Hallo
I am not sure what you want to achieve.
your factor has 3 levels but with only 2 different labels
hb
[1] 1 1 1 1 1 1 2 2 2
Levels: 1 1 2
but
str(hb)
Factor w/ 3 levels 1,1,2: 1 1 1 2 2 2 3 3 3
so you gave only one label to level
Hello,
Several changes in the latest Rcmdr broken the SciViews compatibility. I
am working on it... but it takes longer than expected. Sorry for the
inconvenience.
Best,
Philippe
..°}))
) ) ) ) )
( ( ( ( (Prof. Philippe Grosjean
) )
James McDermott wrote:
Would the unique quadratic defined by the three points be the same
curve as the curve predicted by a quadratic B-spline (fit to all of
the data) through those same three points?
Yes, if you restrict attention to an interval between knots. You'll
need to re-evaluate it
Your alternatives are :
a) Use pdf if you can.
b) See help(bitmap) and especially the details section which tells you
about setting the R_GSCMD path to ghostscript
c) Try Xvfb. I do not know much about this. Try a google search.
Regards, Adai
On Wed, 2005-07-20 at 15:12 +0530, bhumir jhaveri
Hi Bhumir
Do you happen to have an X client installed on your Windows machine?
e.g. cygwin, winaXe, eXceed, ...
I've used winaXe with some success.
Andrew
On Wed, Jul 20, 2005 at 03:12:24PM +0530, bhumir jhaveri wrote:
Hi,
I am newbie in R. I have R installed on Solaris machine and
Here are some different approaches to this problem:
#1 just convert it to chron without using Date class.
days(chron(21-07-2005, format = d-m-y))
#2 after converting it to Date, convert it to numeric and then
# to chron to avoid the problem
dd - as.Date(21-07-2005, %d-%m-%y)
I believe the following is correct:
1.
first of all, as peter daalgaard already pointed out, your data Cp(t)
are following a straight line
very closely, i.e. 0.-order kinetics
2.
for your diff. eq. this means that you are permanently in the range cp
Km so that
dCp/dt = - Vm/Vd = const. =: -b and,
Dear Martyn
Thank you for your fast and helpful answer. It showed me to
my wrong thinking.
I confused the previous value of the Markov Chain with the
initial values to construct the envelope.
In the original C code by W. Gilks there are the arguments
xinit and xprev. The first is used to
Is there an R function to kill a process? I found one in package
fork but it is specific to UNIX and I want something that also
works on Windows. The XP console command, taskkill,
will do it so I can easily get the effect but it won't work
on other Windows systems, even 2000 and NT. I found a
Dear colleague,
Just an e-mail to know if they are people working in the field of
chemoinformatic that are using R in their work. If yes I was wondering if we
couldn't exchange tips and tricks about the use of R in this area ?
Best regards
Fred Ooms
[[alternative HTML version deleted]]
Hi
t.test assumes that your data within each group has a normal
distribution. This is not the case in your example.
I would recommend you a non parametric test like wilcox.test if
you want to compare the mean of two samples that are not normal
distributed.
see ?wilcox.test
Be careful. Your
Dear Philippe,
The development version 1.1-0 of the Rcmdr package (i.e., the one with
localization facilities) isn't yet on CRAN. Georges doesn't say which
version of the Rcmdr package he's using, but I assume 1.0-2 from CRAN, which
I believe should work with SciViews.
Regards,
John
Do you mean people working in areas of mass spectrometry (eg. SELDI,
LC-MS, MALDI etc) data analysis?
- Original Message -
From: Frédéric Ooms [EMAIL PROTECTED]
Date: Wednesday, July 20, 2005 1:44 pm
Subject: [R] Chemoinformatic people
Dear colleague,
Just an e-mail to know if they
I am more lloking at people working in the field of drug discovery performing
clustering analysis, MDS, LDA in order to classify chemicals based on computed
properties.
Fred
-Original Message-
From: S.O. Nyangoma [mailto:[EMAIL PROTECTED]
Sent: Wednesday, July 20, 2005 1:53 PM
To:
On 7/20/2005 7:25 AM, Gabor Grothendieck wrote:
Is there an R function to kill a process? I found one in package
fork but it is specific to UNIX and I want something that also
works on Windows. The XP console command, taskkill,
will do it so I can easily get the effect but it won't work
on
Dear All,
I am trying to build a randomization test for interaction
The problem is as follows: I have a set of stations where the ocurrence and
biomass of each species being investigated was recorded.
For each group of stations, the relative importance of a set of species was
calculated, as
You might want to check out Rajarshi Guha's work on connecting CDK to R. He
has an article in the recent issue of CDK News.
Andy
From: Frédéric Ooms
I am more lloking at people working in the field of drug
discovery performing clustering analysis, MDS, LDA in order
to classify chemicals
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Carsten Dormann
Sent: 20 July 2005 05:41
To: r-help@stat.math.ethz.ch
Subject: [R] analysing non-normal spatially autocorrelated data
Dear fellow R-users,
I wish to analyse a lattice of
Thanks for your answers,
I need to print this data frame into a .csv file to import it in WEKA.
Do you have better solution?
I quite a new user of WEKA I don't know if you can give it a binary
file. I think you can but it will be complicated...
David
On Jul 20, 2005, at 15:07, jim holtman
Simply gooling for writing ARFF file in R gave the following as first hit,
which is right on the WEKA page:
Miscellaneous code
[...]
Function for reading ARFF files into the R statistical package (kindly
provided by Dr Craig Struble).
Function for writing ARFF files from the R statistical
From: Christoph Buser
Hi
t.test assumes that your data within each group has a normal
distribution. This is not the case in your example.
Eh? What happen to the CLT?
Andy
I would recommend you a non parametric test like wilcox.test if
you want to compare the mean of two samples that
Hi,
I would like to add a prefix to colnames in a matrix but I can't get the prefix
option in colnames to work. What am I doing wrong?
X-matrix(NA,3,4)
colnames(X)-c(test,test,test,test)
colnames(X)-colnames(X,prefix=PREFIX.)
X
test test test test
[1,] NA NA NA NA
[2,] NA
Thank you for the information. I have contacted the RPM maintainer and
am awaiting a response.
It occurs to me that my problem could also be fixed by putting ATLAS
on my system. Are there advantages to doing that or any reasons not to?
Prof Brian Ripley wrote:
-lf77blas is part of ATLAS, so
On 7/20/2005 9:54 AM, [EMAIL PROTECTED] wrote:
Hi,
I would like to add a prefix to colnames in a matrix but I can't get the
prefix
option in colnames to work. What am I doing wrong?
The prefix argument is only used when there are no names and colnames is
generating some. You just want to
On Wed, 20 Jul 2005, Kevin E. Thorpe wrote:
Thank you for the information. I have contacted the RPM maintainer and am
awaiting a response.
It occurs to me that my problem could also be fixed by putting ATLAS on my
system. Are there advantages to doing that or any reasons not to?
It is a
If you set the affinity of the R process to processor 0, you can run another (R
or other) process with affinity set to processor 1 and get 100% usage.
Most applications can't take advantage of hyperthreading (or multiprocessors),
since they have to be specially written to do so.
It seems that
[EMAIL PROTECTED] writes:
Hi,
I would like to add a prefix to colnames in a matrix but I can't get the
prefix
option in colnames to work. What am I doing wrong?
X-matrix(NA,3,4)
colnames(X)-c(test,test,test,test)
colnames(X)-colnames(X,prefix=PREFIX.)
X
test test test test
On Wed, 20 Jul 2005 [EMAIL PROTECTED] wrote:
Hi,
I would like to add a prefix to colnames in a matrix but I can't get the
prefix
option in colnames to work. What am I doing wrong?
X-matrix(NA,3,4)
colnames(X)-c(test,test,test,test)
colnames(X)-colnames(X,prefix=PREFIX.)
X
test
Thank you for teh helpful answer. I do have one aother related
question; If I am not interested in teh result of the lapply, as I am
using it instaid of a for loop, example
lapply(df, function(r) sqlQuery(conn, paste(insert into t
values(,r['a'],',',r['b'],')')))
and I am not interested in teh
Hello,
My colleague and I would like to write Java code that invokes R to do a simple
TTest. I've included my sample java code below. I tried various alternatives
and am unable to pass a vector to the TTest method. In my investigation, I
tried to call other R methods that take vectors and
Is there a way to test if a variable is a date?
__
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
See:
http://www.burns-stat.com/pages/Tutor/spreadsheet_addiction.html
and
http://www.stat.uiowa.edu/~jcryer/JSMTalk2001.pdf
Greg Snow, Ph.D.
Statistical Data Center, LDS Hospital
Intermountain Health Care
[EMAIL PROTECTED]
(801) 408-8111
Wensui Liu [EMAIL PROTECTED] 07/19/05 03:22PM
I
Prof Brian Ripley wrote:
On Wed, 20 Jul 2005, Kevin E. Thorpe wrote:
Thank you for the information. I have contacted the RPM maintainer and am
awaiting a response.
It occurs to me that my problem could also be fixed by putting ATLAS on my
system. Are there advantages to doing that or any
On Wed, 20 Jul 2005, Uwe Ligges wrote:
Prof Brian Ripley wrote:
On Wed, 20 Jul 2005, Kevin E. Thorpe wrote:
Thank you for the information. I have contacted the RPM maintainer and am
awaiting a response.
It occurs to me that my problem could also be fixed by putting ATLAS on
my
If you mean of class Date, use inherits(x, Date).
On Wed, 20 Jul 2005, Omar Lakkis wrote:
Is there a way to test if a variable is a date?
--
Brian D. Ripley, [EMAIL PROTECTED]
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford,
Dear Andy
Since the power of the t-test decreases when there are
discrepancies in the data to the normal distribution and there
is only a small loss of power if the data is normal distributed,
the only reason to use the t.test is his simplicity and the
easier interpretation.
Generally I'd prefer
Your surest bet is to look into S+, not R, since the former does
'integrate' with Microsoft applications. To what extent, I don' know.
This is a desirable feature (IMO) and a reason to use S+ rather then R
(already discussed in different thread).
Try the S+ newsgroup.
-Original
There are shapefiles for Europe (and other places) at:
http://www.vdstech.com/map_data.htm (there are 16 polygons withing the one for
Germany, is that enough detail). You can read and plot the shapefiles using
the maptools package.
hope this helps,
Greg Snow, Ph.D.
Statistical Data Center,
I have a package which I use to create excel files from R.
I have not been able to produce a configure script general enough for me
to post it to cran, but I will send it to you if you like.
I use it for production jobs on our linux servers. You may have to
tinker a bit to get it to compile on
Simple question.
For a simple linear regression, I obtained the standard error of
predicted means, for both a confidence and prediction interval:
x-1:15
y-x + rnorm(n=15)
model-lm(y~x)
predict.lm(model,newdata=data.frame(x=c(10,20)),se.fit=T,interval=confidence)$se.fit
1 2
I haven't been an R lister for a bit, but I hope to enlist someone's
help here. I think this is a simple question, so I hope the answer
is not much trouble. Can you please respond directly to this email
address in addition to the list (if responding to the list is
warranted)?
I have a
This does not address the question of how to do this the way you're trying
to do it, using Omegahat. But I've had good results using Simon Urbanek's
program Rserve: http://stats.math.uni-augsburg.de/Rserve/
Hope this helps,
Matthew Wiener
Applied Computer Science Mathematics
Merck Research
John Fox wrote:
Dear Philippe,
The development version 1.1-0 of the Rcmdr package (i.e., the one with
localization facilities) isn't yet on CRAN. Georges doesn't say which
version of the Rcmdr package he's using, but I assume 1.0-2 from CRAN, which
I believe should work with SciViews.
That
Whit Armstrong whit at twinfieldscapital.com writes:
I have a package which I use to create excel files from R.
I have not been able to produce a configure script general enough for me
to post it to cran, but I will send it to you if you like.
I use it for production jobs on our linux
[EMAIL PROTECTED] writes:
Simple question.
For a simple linear regression, I obtained the standard error of
predicted means, for both a confidence and prediction interval:
x-1:15
y-x + rnorm(n=15)
model-lm(y~x)
We obtained some disturbing results from convolve() (inaccuracies and negative
probabilities). We'll try to make the context clear in as few lines as
possible...
Our function panjer() (code below) basically computes recursively the
probability mass function of a compound Poisson distribution.
Can someone please refer me to a function or method that resolves this
structuring issue:
I have two matrices with identical colnames (89), but varying number
of observations:
matrix Amatrix B
217 x 89 16063 x 89
I want to creat one
It _does_ do formatting:
write.xls(x,matrixFromR_formatted3.xls,matrix,formats=rep(#,##0.0;[
Red]#-##0.0;\\,nc))
And also writes multiple sheets to a workbook (I didn't post all the
examples).
I need those features for the reports I generate, but if the perl module
has the same features, then
Can someone please refer me to a function or method that resolves this
structuring issue:
I have two matrices with identical colnames (89), but varying number
of observations:
matrix Amatrix B
217 x 89 16063 x 89
I want to creat one
mark salsburg [EMAIL PROTECTED] writes:
Can someone please refer me to a function or method that resolves this
structuring issue:
I have two matrices with identical colnames (89), but varying number
of observations:
matrix Amatrix B
217 x 89
On 7/20/2005 12:52 PM, Vincent Goulet wrote:
We obtained some disturbing results from convolve() (inaccuracies and
negative
probabilities). We'll try to make the context clear in as few lines as
possible...
Our function panjer() (code below) basically computes recursively the
Please clarify:
Is there a column of common values that you want to use to merge the two
matrices together? Or do you just want to replicate matrix A enough times
to have the same number of rows as B? I could think this was the case, but
16063 is not a multiple of 217.
Anyway, you will have to
When you are forced to use excel (but want to really use R and just
give the result to others in excel), then there are a few options
depending on what you are trying to do. We may be able to give
better help if you can give a specific problem you are trying to
solve.
Some Ideas:
To quickly
The latter, I want to replicate A enough times to have the same rows
as B. I realize that 16063 is not a mutliple 217, that doesn't really
matter for this problem
Assuming I change the colnames of one of them, how can I create a
matrix C that is the composite of the replicated A with B.
thank you
On Wed, 20 Jul 2005, Peter Dalgaard wrote:
mark salsburg [EMAIL PROTECTED] writes:
Can someone please refer me to a function or method that resolves this
structuring issue:
I have two matrices with identical colnames (89), but varying number
of observations:
matrix A
If 16063 were an integer multiple of 217 you could do this in two steps;
first construct the left half with rep() or
L - kronecker(rep(1,16063/217),A)
and then use cbind(L,B).
However, 16063 is not a multiple of 217, so I don't know what you want to do
with the leftover 5 rows you need in L. If
Hi,
Se if this will solve your question (quick and dirty)...
A is the first matrix, B is the second...
F.CombineMat -
function (A, B)
{
ind-rep(1:nrow(A),ceiling(nrow(B)/nrow(A)))
ind-ind[1:nrow(B)]
A-A[ind,]
colnames(A)-paste('A', colnames(A), sep='')
R-cbind(A,B)
R
}
Greetings,
Pedro
At 18:24
Sundar Dorai-Raj writes:
Hi, Laura,
Would ?predict.glm be better?
plot(logarea, hempresence,
xlab = Surface area of log (m2),
ylab=Probability of hemlock seedling presence,
type=n, font.lab=2, cex.lab=1.5, axes=TRUE)
lines(logarea, predict(hemhem, logreg, response),
Thanks for the replies
1. about Rcmdr versions
Sorry, I did'nt say wich versions because I tried a lot... before asking
for help
First I tried the version of the site www.sciviews.org (0.9-14), after
I tried the Cran version (1.0-2) and at last I got the 1.1-0 version on
John Fox's web.
After
This is from a colleague who is attempting to install R 2.1.1 on a Sun
machine. I hope this is sufficient information for someone to give us
some pointers.
thanks,
Richard
Our system is a Sunfire 480 running Solaris V2.8
The LD_LIBRARY_PATH is set
Richard Valliant [EMAIL PROTECTED] writes:
This is from a colleague who is attempting to install R 2.1.1 on a Sun
machine. I hope this is sufficient information for someone to give us
some pointers.
thanks,
Richard
Our system is a Sunfire 480 running Solaris V2.8
The LD_LIBRARY_PATH
Hi,
I try to put a legend in a xyplot graphic.
xyplot(y~x|g,ylim=c(0,80),xlim=c(10,40),as.table=T,layout=c(2,3), ylab=Número
de machos capturados,xlab=expression(paste(Temperatura (,degree,C))),
key=list(corner=c(0,0),x=0, y=0, text=list(legenda),lines=list(col=cor, lwd =
espessura,
Hi,
I have a dataframe ds1.2 - 503 categorial variables
and 1 continuous response variables. I ran aregImpute
to deal with NA's and got the followig error:
fmla = terms( Response ~ . ,data=ds1.2)
ds.i = aregImpute(fmla,data=ds1.2)
Error in matrix(as.double(1), nrow = n, ncol = p,
dimnames =
Hi! I am trying to program a TURF analysis and have looked everywhere on the
web to see if I can find any clues on how to do this. I've looked at SAS
manuals as well. If anyone has ever heard of this kind of marketing analysis
and has programmed it I'd be interested in learning how. Any help
Thomas Isenbarger isen at plantpath.wisc.edu writes:
I haven't been an R lister for a bit, but I hope to enlist someone's
help here. I think this is a simple question, so I hope the answer
is not much trouble. Can you please respond directly to this email
address in addition to the
Here is an example where R is the client and Excel is the server
so that R is issuing commands to Excel. This example uses the
RDCOMClient package from www.omegahat.org:
library(RDCOMClient)
xl - COMCreate(Excel.Application) # starts up Excel
xl[[Visible]] - TRUE
See ?goodfit in package vcd.
On 7/20/05, Thomas Isenbarger [EMAIL PROTECTED] wrote:
I haven't been an R lister for a bit, but I hope to enlist someone's
help here. I think this is a simple question, so I hope the answer
is not much trouble. Can you please respond directly to this email
Check out sum.exact in the caTools package.
On 7/20/05, Vincent Goulet [EMAIL PROTECTED] wrote:
We obtained some disturbing results from convolve() (inaccuracies and negative
probabilities). We'll try to make the context clear in as few lines as
possible...
Our function panjer() (code
Thanks a lot for the input on this Professor Murdoch. I really
appreciate all the help.
Regards,
Jim
On 7/20/05, Duncan Murdoch [EMAIL PROTECTED] wrote:
James McDermott wrote:
Would the unique quadratic defined by the three points be the same
curve as the curve predicted by a quadratic
I want to do the following:
glm(y ~ x1 + x2 +...)
within a panel. Hence y, x1, and x2 all vary at the individual
level. However, there is likely correlation of these variables
within an individual, so standard errors need adjustment.
I do not want to estimate fixed effects, but do want to
Young Cho wrote:
Hi,
I have a dataframe ds1.2 - 503 categorial variables
and 1 continuous response variables. I ran aregImpute
to deal with NA's and got the followig error:
fmla = terms( Response ~ . ,data=ds1.2)
ds.i = aregImpute(fmla,data=ds1.2)
Error in matrix(as.double(1), nrow =
Dear all,
I am in the process of migrating an S programme library to R, and it
involves the family entity.
I have checked ?family but it does not give much detail of its components.
I will be very grateful if anyone can point towards sources/ways to look up
on this areas
with an aim to find
Dear All,
With the warm support of every R expert, I have built my R library
successfully.
Especially thanks: Duncan Murdoch
Gabor Grothendieck
Henrik Bengtsson
Uwe Ligges
Without your help, I will lower efficiency.
I noticed that
All R users
We observed that for the NESTED DESIGN there is a junction of the
variable effects randomic that impede the separation of the values of
these variance components.
Does this also happen in the SPLIT-PLOT method?
Is there junction between main plot error and replicate error?
There
Dear All,
I wonder whether it is still valid to use the following R code for cut. All I
have done is changed:
if (is.na(breaks) | breaks 2)
to:
if (is.na(breaks) | breaks 1)
so that it covers interval of 1?
It seems okay for my purposes but I am not sure why R specifically does
96 matches
Mail list logo
