On Thu, 8 Sep 2005, Erich Neuwirth wrote:
>
> sapply(1:length(xxx),function(x)b[xxx[x],yyy[x]])
> does what I need and produces
> [1] 1 2 1 2 1 4 3 4 3 4
>
> Is there a function taking xxx,yyy, and b as arguments
> producing the same result?
b[cbind(xxx,yyy)]
> Essentially, I am asking for a vers
On Thu, 8 Sep 2005, Shige Song wrote:
> Dear All,
>
> The "coxph" function in the "survival" package allows multiple frailty
> terms.
Um, no, it doesn't.
> In all the examples I saw, however, the frailty terms are nested.
> What will happen if I have non-nested (that is, cross-classified) frailt
On Thu, 8 Sep 2005, John Sorkin wrote:
> I have a batch of data in each line of data contains three values,
> calcium score, age, and sex. I would like to predict calcium scores as a
> function of age and sex, i.e. calcium=f(age,sex). Unfortunately the
> calcium scorers have a very "ugly distributi
On Wed, 7 Sep 2005, Clark Allan wrote:
>
> $MLE$message
> [1] "ERROR: ABNORMAL_TERMINATION_IN_LNSRCH"
>
>
> WHAT DOES THIS ERROR MESSAGE MEAN???
>
Looking at the code in optim() a little, it looks as though this error
comes when the optimiser tries to do a line search in the steepest descent
di
of the usefulness of these frailty models and I don't know how
well their properties are known. I wouldn't use them except when I was
actually interested in the variance components, and I haven't worked on
any problems like that, so I haven't investigated the issue
On Wed, 7 Sep 2005, Chris Buddenhagen wrote:
>
> 1) I have been really pleased with R as a means of doing and learning
> statistics. I work in a Spanish speaking country- and I wanted to pass on
> the benefits of R to my Spanish speaking colleagues. There are a couple of
> introductions to R in Spa
On Tue, 6 Sep 2005, [EMAIL PROTECTED] wrote:
>
> However, a related method available for 'optim' is "L-BFGS-B"
> "which allows _box constraints_, that is each variable can be
> given a lower and/or upper bound. The initial value must satisfy
> the constraints." This can be set in a parameter for '
On Mon, 5 Sep 2005, [EMAIL PROTECTED] wrote:
> A trap easy to fall into!
>
So easy, in fact, that it's a FAQ.
-thomas
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http:/
ou should know how they were
sampled.
To drop these observations you could use
obsChina <- subset(China, !is.na(psu) & !is.na(strata) & !is.na(weight0x))
and then use obsChina rather than China in the svydesign() function.
-thomas
> -Bobby
&g
On Sun, 4 Sep 2005, A Das wrote:
> Thanks, Thomas.
>Yes, that's exactly what happened: the warnings
> came first after "data(China)", and then after
> "dchina<-svydesign..." So the design object isn't
> being produced? The dataset is very large, and the
> weights were already set in Stata befo
On Sat, 3 Sep 2005, A Das wrote:
> Hi all, I've been trying to get a large (12mb) Stata
> survey database into R. I managed that, but when I
> attach survey weights, something goes wrong. The error
> message is: object dchina not found. Here's the
> script:
If that is the *first* message then som
On Sat, 3 Sep 2005, John Sorkin wrote:
> A question for R (and perhaps S and SPlus) historians.
>
> Does anyone know the reason for the inconsistency in the way that the
> action that should be taken when data are missing is specified? There
> are several variants, na.action, na.omit, "T", TRUE,
On Fri, 2 Sep 2005, Spencer Graves wrote:
> Permit a mild protest on the word "appropriate" in this context. The
> global assignment operator "<<-" provides, for my tastes, excessive
> opportunities for problems. If I define "x" someplace else and then
> call your function, it may change
On Fri, 2 Sep 2005, Brahm, David wrote:
> In a clean environment under R-2.1.0 on Linux:
>> x <- 1:5
>> x[3] <<- 9
> Error: Object "x" not found
>
> Isn't that odd? (Note x <<- 9 works just fine.)
>
Well, yes and no.
It is the result of a bug fix a version or two ago that dealt with the
case w
On Mon, 29 Aug 2005 [EMAIL PROTECTED] wrote:
> Hello,
> I'm interested in LaTeX notations in R graphs. For example, I would like to
> know
> the R code to obtain the graph here,
> http://www.r-project.org/screenshots/power.png
Here is the code:
x<-seq(-10,10,length=400)
y1<-dnorm(x)
y2<-dnorm(x
guide! http://www.R-project.org/posting-guide.html
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-
On Thu, 25 Aug 2005, Devarajan, Karthik wrote:
>
> Hello
>
> I am fitting a Cox PH model using the function coxph(). Does anyone know how
> to obtain the estimate of the covariance matrix under the null hypothesis.
> The function coxph.detail() does not seem to be useful for this purpose.
>
You c
On Fri, 19 Aug 2005, [iso-8859-1] G?ran Brostr?m wrote:
Similar strange results appear in other date-related functions. I plan to
write functions that converts, eg, "1977-01-31" to the real number 1977.084
and back. What function in R does what 'mdy.date' claims to do?
Well, mdy.date does (
On Fri, 19 Aug 2005, Fredrik Thuring wrote:
I have a few questions concerning reading of tables from R to
other programs. My main question is if it???s even possible to read a table
created in R (with the functions data.frame and save) to Excel (or
maybe SAS) and if so how does one do
On Mon, 15 Aug 2005 [EMAIL PROTECTED] wrote:
>
> Dear R wizards:
>
> under R-2.1.0:
>
> eargs <- 3:5;
> line <- paste(c("echo A B", eargs));
> cat("executing from R: '", line, "'\n");
> system(line);
>
> Oddly, only "A" and "B" are echoed, not the eargs. I had hoped that
> line would be one strin
On Sat, 13 Aug 2005, Prof Brian Ripley wrote:
> I suspect you want
>
> AdjForBase2 <- function (data, inds)
>
> and to refer to data[inds, 1] and data[inds, 2], but since your code is
> completely devoid of spaces and indentation, I have paid it little
> attention.
http://msr.uwaterloo.ca/msr2005/
ropriate functions or a line of code!
>
> Best,
>
> Shelby
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http:
-
O__ Peter Dalgaard ?ster Farimagsgade 5, Entr.B
c/ /'_ --- Dept. of Biostatistics PO Box 2099, 1014 Cph. K
(*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918
~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907
--
---
On Thu, 11 Aug 2005, Clark Allan wrote:
> hi all
>
> is there a package that undertakes subset selection but BASED ON AIC or
> any other information criteria.
>
> i've seen the "subselect" and the "leaps" package but i have not played
> around with them yet.
The leaps package finds a best model o
On Wed, 10 Aug 2005, Gorjanc Gregor wrote:
> Dear R-users!
>
> I would like to fit exponential, Weibull and log-logistic via glm() like
> functions. Does anyone know a way to do this? Bellow is a bit longer
> description of my problem.
I think you want to use survreg(). It will still work when t
On Tue, 9 Aug 2005 [EMAIL PROTECTED] wrote:
>
> dear R wizards:
>
> plot( 1, 1, ylim=(2,10), xlim=(2,10), type="n");
> rect( -1, -1, 12, 12, col=gray(0.99) );
>
> unfortunately wipes out the border axes around the plot. how do I keep
> this?
I think you meant
plot( 1, 1, ylim=c(2,10), xlim=c(
On Mon, 8 Aug 2005, Jake Michaelson wrote:
> I use Mac OS X at home and Linux at work, so the R Aqua GUI has spoiled
> me. I have not seen its equal so far (on Windows or Linux). The most
> important thing to me is how easily accessible the help and
> documentation is. I like how when I begin t
On Thu, 4 Aug 2005, Matt Crawford wrote:
> I understand that in R, for loops are not used as often as other
> languages, and am trying to learn how to avoid them. I am wondering
> if there is a more efficient way to write a certain piece of code,
> which right now I can only envision as a for loo
On Thu, 4 Aug 2005, Jose Claudio Faria wrote:
I would like to suggest that all R functions/etc like:
codes-deprecated
grid-internal
ns-alt
ns-dblcolon
ns-hooks
ns-internals
ns-lowlev
ns-reflect.Rd
tools-internal
ts-defunct
utils
On Wed, 3 Aug 2005, array chip wrote:
> Hi, I have a matrix with both positive and negative
> numbers, I would like to use image() to draw a
> heatmap. How can I can design a palette (or is there a
> function already available) that treat negative
> numbers in a blue gradient and positive numbers
On Wed, 3 Aug 2005, Vicky Landsman wrote:
>
> Dear all,
> I have a list containing 150 simulated datasets in R. Is there a way to
> convert it to the Stata format such that I will be able to apply some Stata
> functions on each dataset in the list?
For Stata your best best is probably to stack al
p@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
>
> --
> Spencer Graves, PhD
> Senior Development Engineer
> PDF Solutions, Inc.
> 333 West San Carlos Street Su
On Mon, 1 Aug 2005, Chris Wallace wrote:
> I am struggling with migrating some stata code to R. I have a data
> frame containing, sometimes, repeat observations (rows) of the same
> family. I want to keep only one observation per family, selecting
> that observation according to some other varia
On Fri, 29 Jul 2005, Dieter Menne wrote:
> John Sorkin grecc.umaryland.edu> writes:
>
>>
>> I am trying to model data in which subjects are followed through time to
>> determine if they fall, or do not fall. Some of the subjects fall once,
>> some fall several times. Follow-up time varies from su
On Fri, 29 Jul 2005, John Sorkin wrote:
> I am trying to model data in which subjects are followed through time to
> determine if they fall, or do not fall. Some of the subjects fall once,
> some fall several times. Follow-up time varies from subject to subject.
> I know how to model time to the f
p@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
>
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
_
On Wed, 27 Jul 2005, McGehee, Robert wrote:
>
> However, the sink() does not seem to be capturing the error messages as
> I would have hoped. That is, if the R-script is {print("abc");
> stop("def")}, the print output is captured to the file, the stop error
> message is not, and instead sent to th
On Wed, 27 Jul 2005, Peter Dalgaard wrote:
> Adaikalavan Ramasamy <[EMAIL PROTECTED]> writes:
>
>> I think it is doing what is supposed to do but I never used read.spss,
>> so take this with a pinch of salt.
>>
>> In R when you use as.integer on a factor, the one with the lowest level
>> gets a va
On Tue, 26 Jul 2005, Garrett Fox wrote:
> I am trying to do logistic regression with a categorical predictor variable
> with the glm() function, family=binomial. Using glm() I would like to be
> able to calculate the confidence intervals of all three possible odds ratios
> for a factor (the facto
On Mon, 25 Jul 2005, Wladimir Eremeev wrote:
> Hello, all
>
> Please, consider the following pieces of code.
>
> 1.
> v<-0.5
> text(x=2,y=2,eval(substitute(expression(bold(S==V)),list(V=formatC(v,format="f",digits=2)
>
> This plots "S=0.5" in bold. Both "S" and "0.5" are bold.
>
> 2.
> v<-0.5
On Mon, 25 Jul 2005, Petr Mandys wrote:
> Thank you!
>
> I'm using this:
>
> require(NADA)
>
> var_sf=survfit(Surv(a, !b))
> var_result=new("cenfit", survfit=var_sf)
>
> Class cenfit has methods to get mean a variance from survfit result. Is
> this correct?
>
NADA is designed for left censoring.
On Fri, 22 Jul 2005, Gabor Grothendieck wrote:
>
> I think legend accepts a list argument directly so that could be
> simplified to just:
>
> a<-7
> plot(1)
> L <- list(bquote(alpha==.(a)),bquote(alpha^2+1==.(a^2+1)))
> legend("topleft",legend=L)
Except that it wouldn't then work: the mathemat
On Fri, 22 Jul 2005, Dan Bolser wrote:
> On Fri, 22 Jul 2005, Gabor Grothendieck wrote:
>
>> Try as.expression(bquote(...whatever...))
>
> Sob, wimper, etc.
a<-7
plot(1)
legend("topleft",legend=do.call("expression",
list(bquote(alpha==.(a)),bquote(alpha^2+1==.(a^2+1
On Thu, 21 Jul 2005 [EMAIL PROTECTED] wrote:
> Hi,
>
> I have a question about the output for variance of random effect from a gamma
> frailty model using coxph in R. Is it the vairance of frailties themselves or
> variance of log frailties? Thanks.
>
For a Gamma frailty model it is the variance
ng-guide.html
>
> --
> Spencer Graves, PhD
> Senior Development Engineer
> PDF Solutions, Inc.
> 333 West San Carlos Street Suite 700
> San Jose, CA 95110, USA
>
> [EMAIL PROTECTED]
> www.pdf.com <http://www.pdf.com>
> Tel: 408-938-4420
> Fax: 408-280-7915
&
On Tue, 19 Jul 2005, [iso-8859-1] [EMAIL PROTECTED] wrote:
>
> beta0,eta0,n are initialized as single integer, temp is a 1 dimension array
> containing 9 integer.
>
> integrate(function(beta)
> ((beta/(eta0)^beta)^n)*prod(temp^(1-beta)*exp(-sum(temp^beta)/(eta^beta)))*(1/(sqrt(2*pi))*exp(((beta-
= PHI" onto a graphic
>> i.e. " hat(phi)=2 "
>> how does one do this?
>>
>
> legend(-5, .3, substitute(hat(phi) == PHI, list(PHI = PHI)))
>
or
legend(-5, .3, bquote(hat(phi) == .(PHI)))
-thomas
Thomas Lumley Assoc.
On Mon, 18 Jul 2005, Weiwei Shi wrote:
> Hi,
> I have a dataset with 2194651x135, in which all the numbers are 0,1,2,
> and is bar-delimited.
>
> I used the following approach which can handle 100,000 lines:
> t<-scan('fv', sep='|', nlines=10)
> t1<-matrix(t, nrow=135, ncol=10)
> t2<-t(t1)
On Sat, 16 Jul 2005, Thomas Lumley wrote:
>
> Yes, but you don't need to go via the baseline. The survival curves for
> any two covariate vectors z1 and z2 are related by
>
> S(t; z1)= S(t; z2)^(z1-z2)
>
Actually
S(t; z1)=S(t;z2) ^(beta'(z1-z2))
On Sat, 16 Jul 2005, Bernardo Rangel Tura wrote:
> At 10:11 12/7/2005, [EMAIL PROTECTED] wrote:
>
>> hi all
>>
>> why does R do this:
>>
>> (-8)^(1/3)=NaN
>>
>> the answer should be : -2
>
Yes and no.
The problem is that the reciprocal of 3 is not exactly representable as a
floating point numbe
On Sat, 16 Jul 2005, Kylie-Anne Richards wrote:
> Thank you for your help.
>
>> In any case, to specify f.pom You need it to be a factor with the same set
>> of levels. You don't say what the lowest level of pom is, but if it is,
>> s
On Fri, 15 Jul 2005, David Groos wrote:
> Thank you-all very much for your help, your responses and help has been
> very encouraging. The following doesn't close the case but it tables
> it...
>
> First I copied Ken's code into my R Console and...it worked great!
> That was baffling as it looked
On Fri, 15 Jul 2005, Kylie-Anne Richards wrote:
>
> FIRST Q: The default uses the mean of 'vo' and mean of 'po', but what is it
> using for the factors?? Is it the sum of the coef of the factors divided by
> the number of factors??
It uses the mean of each factor variable. The $means component of
/dsc2005/program.html
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
On Thu, 14 Jul 2005, Duncan Murdoch wrote:
> On 7/14/2005 7:19 AM, Ginters wrote:
>> Why does memory need so much (1.6 GB) space? How can I enlarge it? Is it
>> possible to allocate a part of memory used to the hard drive? Or, is the
>> trouble only with my script?
>
> This sounds like a problem w
On Wed, 13 Jul 2005, EJ Nikelski wrote:
>
> Your suggestion that the unprintable characters represent UTF-8 encoded
> Unicode left and right double quotes also appears correct. Now, although
> the suggested work-around may well help, the foreign package does seem
> to be creating a corrupted file.
On Mon, 11 Jul 2005, Adaikalavan Ramasamy wrote:
> I am not sure if there is an easy way around this. An ugly hack is to
> make a copy the function "survival:::plot.cox.zph" and make your
> modified function. But there are others in the list who might know
> neater solutions.
If you then send a p
On Fri, 8 Jul 2005, Florian Hahne wrote:
> Dear listers,
>
> I am trying to compute the exact conditional test given strata margins
> of a 2 by 2 by K array using the mantelhaen.test function to get a
> common odds ratio estimate.
> The estimate for the test on the following data is 0, which in my
hen discovered to my surprise that the chi-square
> test for independence of the significance of the F-tests indicated
> that they were not independent. I was stumped by this but fortunately
> Thomas Lumley came to my rescue with an explanation. There is no
> reason why the results of the
On Tue, 5 Jul 2005, Chuck Cleland wrote:
>
> In SPSS, open the syntax file "warpbreaks.sps". This syntax will
> read in the data and apply variable and value labels to variables that
> were factors in R. I find you often need to edit the syntax file
> created by write.foreign just a bit in SPSS
On Tue, 5 Jul 2005, Erin Hodgess wrote:
> Dear R:
>
> In the foreign library, there is a function "write.foreign".
>
> When this is used to write to an SPSS file, there are
> actually 2 files as output: a data file and a code file.
>
> How would you coordinate these to become an SPSS sav file,
> p
On Tue, 5 Jul 2005, Sean O'Riordain wrote:
> Hi!
>
> I'm trying to use lm(y~x) amongst others in an automated way; I've
> gone through the section on indexing in R-lang and I've looked MASS4.
> How do I find out more about the structure of the returned object? In
> perl I can look at object struc
On Mon, 4 Jul 2005, E. Michael Foster wrote:
> I'm a long-time STATA user and a R newbie. I'm doing ok, but I'm addicted
> to STATA macro variables. Is there something like a macro variable in R?
>
> Specifically, I'd like to be able to do something like
>
> for (i in 1:3) {
> .
>
On Wed, 22 Jun 2005, Dan Bebber wrote:
> I would like to extract predicted failure times from a
> coxph model in library(survival). However, none of the
> prediction options ("lp", "risk", "expected", "terms")
> seem to bear any relationship to failure time.
>
> Perhaps I am asking the wrong quest
On Fri, 17 Jun 2005, Mike R wrote:
> k = c(1:9)
> if( length( which(k==3) ) ){ print("contained") }else{ print("not contained")
> }
>
> is therre a simple way to test if a vector/list contains a particular value?
value %in% vector
> more generally, is the a documentaion page that lists/descri
On Fri, 17 Jun 2005, Johanna Sundvik wrote:
> However, this "Mean1" is categorical when it should be real numbers.
>
>> Mean1
> [1] 4.4332 8.5113 35.1624 9.1693 2.974 65.1578 43.2241 3.1278 5.3364
> Levels: 2.974 3.1278 35.1624 4.4332 43.2241 5.3364 65.1578 8.5113 9.1693
>
> Why R does not un
On Thu, 16 Jun 2005, Mark Hempelmann wrote:
> Dear WizaRds,
>
> I am struggling to compute correctly a cluster sampling design. I want
> to do one stage clustering with different parametric changes:
>
> Let M be the total number of clusters in the population, and m the
> number sampled. Let
gt;> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>>
>
> --
> Brian D. Ripley, [EMAIL PROTECTED]
> Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
> University of Oxford, Tel: +44 1865 272861
On Wed, 15 Jun 2005, Spencer Graves wrote:
> 1. The function "KalmanLike" seems to change its inputs AND
> PREVIOUSLY MADE copies of the inputs.
The strange thing is that it changes its inputs. It would be expected that
this also changed previously made copies. The copies will just be
On Mon, 13 Jun 2005, Kerry Bush wrote:
>
> Suppose I fit the following model:
>
>> library(gam)
>
>> fit <- gam(y~x1+x2+s(x3),family=binomial)
>
> and then I use
>
>> fitf$coef
> x1x2 s(x3)
> 4.1947460 2.7967200 0.0788252
>
> are the coefficients for x1 and x2 the estimated
n,
>
>linkinv = linkinv, variance = variance, dev.resids = dev.resids,
>
>aic = aic, mu.eta =mu.eta, initialize = initialize,
>
>validmu = validmu), class = "family")
> }
>
> Thank you for your kind attention.
> Emmanuel
>
>
> EMMANUEL NJI ABATIH
&g
On Mon, 13 Jun 2005, BXC (Bendix Carstensen) wrote:
> The key to solving your problem is that read.spss per default
> gives you a *list* and not a *dataframe* (can anyone explain this
> choice of default?).
>
The reason for the default is speed.
It's possible that the default could be changed to
On Fri, 10 Jun 2005, Jagarlamudi, Choudary wrote:
> Hi Group,
>
> I performed the following commands to build my package in R 2.0 under Windows
> XP
> I got all my tools from Dr. Duncan Mudroch's website.
> I did a RCMD build dnal and it built a tar file for me.
> I did a RCMD INSTALL dnal and it
On Fri, 10 Jun 2005, Hanke, Alex wrote:
> Dear All,
> I'm having just a little terminology problem, relating the language used in
> the Hosmer and Lemeshow text on Applied Survival Analysis to that of the
> help that comes with the survival package.
>
> I am trying to back out the values for the b
On Fri, 10 Jun 2005, Pavel Khomski wrote:
> hello,
>
> the question concerning to the memory used and g.c. after having removed
> objects. What is wrong?
Nothing is wrong.
> bevor
> ---
>
> > gc()
> used (Mb) gc trigger (Mb) max used (Mb)
> Ncells 31
On Fri, 10 Jun 2005, Remigijus Lapinskas wrote:
> Dear All,
>
> A very simple question:
>
> > library(survival)
> > fit <- coxph(Surv(time, status) ~ x, data=aml)
> > survfit(fit)
> Call: survfit.coxph(object = fit)
>
> n events median 0.95LCL 0.95UCL
> 23 18 30 18
On Thu, 9 Jun 2005, McGehee, Robert wrote:
I'm seeing some inconsistent behavior when re-assigning values in a data
frame. The first assignment turns all of the 0s in my data frame to 2s,
the second fails to do so.
df1 <- data.frame(a = c(NA, 0, 3, 4))
df2 <- data.frame(a = c(NA, 0, 0, 4))
df1
On Thu, 9 Jun 2005 [EMAIL PROTECTED] wrote:
I only have an other question:
since h(t) = h0(t) exp(B1*X1+ B2*X2 + B3*X3 + B4*X4)
represent the hazard at time t.
In a linear prediction,
what Value = B1*(X1-mean(X1)) + B2*(X2-mean(X2)) +
represent?
coxph() parametrizes the model so th
ibull has to be parametrized as a location-scale family as described
above.
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
__
R-help@stat.math.ethz.ch mai
On Tue, 7 Jun 2005, David Whiting wrote:
When I have encountered this error message in the past seems to have
resulted from a blank/empty level in a factor or an empty character.
For example:
I've sent a fixed version to CRAN.
The C standard says that the return value of fwrite() is equal to
On Tue, 7 Jun 2005, Liaw, Andy wrote:
The easiest way, perhaps, is to assign the internal functions to frame 1:
Another easy way is to use the MC() function in section 3.3.1 of the R
FAQ, which makes function closures in a way that works in S-PLUS.
-thomas
ite.dta(x, file = "temp.dta") :
a binary write error occurred
My work-around at the time was to go through the data replacing ""
with something else that I could then deal with later.
--
David Whiting
University of Newcastle upon Tyne, UK
Thomas Lumley
On Mon, 6 Jun 2005, Jean Eid wrote:
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
On Fri, 3 Jun 2005, Sander Oom wrote:
Dear R users,
I have received a table in the following format:
id a b c1 c2 d1 d2
1 1 1 65 97 78 98
2 1 2 65 97 42 97
3 2 1 65 68 97 98
4 2 2 65 97 97 98
Factors of the design are: a, b, and e, where e has levels c a
On Thu, 2 Jun 2005, hadley wickham wrote:
An environment is a hash table, and copying occurs only on modification,
when any language would have to copy in this context.
Yes, I'm aware that copying only occurs on modification. However, it
was my understanding that
a <- list(a =1)
a$b <- 4
wo
On Thu, 2 Jun 2005, hadley wickham wrote:
cache results. Other languages usually offer a dictionary data type
which I can use as an efficient way to dynamically cache already
calculated results - what's the best way to do this in R?
It really depends on what sort of data you want to cache - i
not to do that. Now, if
you do have write permission for that file and R is just confused we have
a different problem, and one that needs more details.
Does the problem happen with other filenames or other directories, for
example?
-thomas
Thomas Lumley Assoc.
On Wed, 1 Jun 2005 [EMAIL PROTECTED] wrote:
Dear all,
I wrote the following batch script on a iMac, and ran it on a linux
mosix cluster.
tu <- read.table("cage.mm5.tags.rna_lib.CAA-CAJ.tu-reshape.table")
tu_reshaped <- t(reshape(tu[1:50,], direction="wide", timevar="tu",
idvar=c("rna","lib"))
On Wed, 1 Jun 2005, luan_sheng wrote:
hello,everyone. I have one question:
example 1
x=numeric(0)
y=5
print(x+y)
numeric(0)
example 2
x=numeric(1)
y=5
print(x+y)
[1] 5
why the print(x+y) is numeric(0) at the first example, but the result is 0
at the second example?
numeric(0) is a zer
On Tue, 31 May 2005, Stefan Pohl wrote:
Dear list,
I want o fit a shared gamma frailty model with the frailty specification
in the survival package.
I have partly left-truncated data and time-dependent covariates. Is it
possible to combine these two things in the frailty function. Or are th
On Tue, 31 May 2005, Matt McCall wrote:
*** malloc[477]: Deallocation of a pointer not malloced: 0x143cf5c0; This
could be a double free(), or free() called with the middle of an allocated
block; Try setting environment variable MallocHelp to see tools to help debug
And this is a double fre
other hand, does definitely sound like a bug.
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/li
On Tue, 31 May 2005, Duncan Murdoch wrote:
Mäkinen Jussi wrote:
Dear All,
I luckily found the following feature (or problem) when tried to apply
ifelse-function to an ordered data.
test <- c(TRUE, TRUE, TRUE, FALSE, FALSE, FALSE, FALSE)
ifelse(test, 0, 1:4)
[1] 0 0 0 4 1 2 3
As Dim
On Tue, 31 May 2005, Bliese, Paul D LTC USAMH wrote:
We have been using Least Angle Regression (lars) to help identify
predictors in models where the outcome is continuous. To do so we have
been relying on the lars package. Theoretically, it should be possible
to use the lars procedure within
On Mon, 30 May 2005, huang min wrote:>
My intention is to invert the covariance matrix to perform some
algorithm which is common in the estimating equations like GEE.
In that case there is no benefit in being able to invert very extreme
covariance matrices. The asymptotic approximations to the
On Mon, 30 May 2005, Thomas Yee wrote:
Hello,
It seems that if glm used a namespace then the conflict would be avoided?
No. glm does use a namespace, so this can't be true. Remember that R
passes arguments by value, and consider
glm(y~x, family=poisson())
The namespace ensures that f
On Fri, 27 May 2005 [EMAIL PROTECTED] wrote:
Sorry, still confused. If I dont have fpc's ready in my dataset (calculate
myself?) that means that R will use the weight of an individual for each of his
repeated observations. But is that then still correct? The "cluster" individual
is ignored and e
On Fri, 27 May 2005 [EMAIL PROTECTED] wrote:
Thank you for your reply.
Does that mean that in order to take in account the repeated measures I denote
these as another cluster in R?
Yes, but unless you have multistage finite population corrections to put
in the design object only the first
of random effect.
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
[EMAIL PROTECTED] University of Washington, Seattle
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE d
On Thu, 26 May 2005, Mark Hempelmann wrote:
Dear WizaRds,
Working through sampling theory, I tried to comprehend the concept of
stratification and apply it with Survey to a small example. My question is
more of theoretic nature, so I apologize if this does not fully fit this
board's intenti
401 - 500 of 1258 matches
Mail list logo