HDTenEightyP/NetNews huggingface dataset:
From: c.oneill@brahma.trl.oz.au (Chris O'Neill)
Newsgroups: aus.snow
Subject: RE: Snow report Hotham
Date: 16 Aug 92 05:45:04 GMT
X-To: aus-snow@unimelb.edu.au
X-Vms-To: IN%"aus-snow@unimelb.edu.au"

>
>
>In article <djm.713421735@phobos>, djm@kbs.citri.edu.au (David Mills) writes:
>> Having just spent a week at Mt. Hotham, I'd have to say it is probably the
>> worst snow cover I've seen in August, since 1982.
>
>Earlier this year I posted that I had observed a good indicator of the coming
>snow season. This was based on moisture in the ground, or in this case the lack
>of it - the drought is still widespread.
>
>Those who had intricate theories why this (and any other predictor) would never
>work must be now wondering. After all we had the best start possible - early
>snow and lots of it.

Interesting idea of "the best start possible" - two feet on the ground in
the middle of June.  Not really in the same class as 1960 - two metres on
the ground in the second week of May, dipping below two metres by the end
of May but going back above two metres by the second week of June and
staying above two metres until the end of October.  It peaked above three
metres.

Chris O'Neill
From: dougcc@brt.deakin.edu.au (Douglas Miller)
Newsgroups: aus.snow
Subject: Re: Snow reports going downhill?
Date: 18 Aug 92 19:48:56 GMT
Organization: C&CS, Deakin University

In article <-199208180018.AA19179@mineng.dmpe.csiro.au>, old hacker
<thornton@dmpe.csiro.au> writes:

> 		Firstly, does anyone else feel that this information service
> 		has gone down-hill (more or less literally) since it became
> 		a News topic ?

Well, I just happen to have some stats here:

   1989:           10 snow reports
   1990:           34 snow reports
   1991:           70 snow reports
   1992 (so far):  29 snow reports

so I don't think your theory is sustained (given that the 1992 season is
only just stated, and hence is going to be shorter).

And as to your allusion about the influx of fixed-heel resort skiers ---
yes, there are more of them.  I reckon this is a good thing, as all snow
reports of of some value to everyone (e.g., a free-heeler can extrapolate
from a Falls Creek report to gauge conditions at Mt Nelse, and a resort
skier could extrapolate from a Munyang report to gauge conditions at Blue
Cow).

> 		Can I suggest that a separate aus.snow.waffle.waffle topic be
> 		started for those who wish talk about whether they are clever
> 		enough to be able to predict snow or not (and similar Citizen
> 		Band radio behaviour) as distinct from report it !

Well, I don't think this is necessary (yet).  However, it would help if
everbody prefixed snow reports with the words "Snow report:".
Newsgroups: aus.snow
From: lloyd@cs.monash.edu.au (lloyd allison)
Subject: Re: Snow report: Lake Mountain
Organization: Computer Science, Monash University, Australia
Date: Mon, 17 Aug 1992 08:12:08 GMT

dougcc@brt.deakin.edu.au (Douglas Miller) writes:

>Sunday 16-AUG-1992
>There was about 25cm of snow off trails, and 10cm hard-packed on trails. 
>The cover on the trails was adequate, apart from a few small washouts, that
>were easy enough to avoid.  Despite this, ARC signs claimed conditions
>ranging from "excellent" (Echo Flat trail) to "fair" (Jubilee trail).  In
>fact all trails were about "good" standard.
>I bought my first "groomed trail" ticket for this day (having no real
>choice).  I couldn't see much evidence of recent grooming.  Long Heath
>trail was clearly ungroomed, and a sign at the start said so.

Have you noticed "inflation" of ARC reports this year?
Up to last year,  excellent meant very good, very good meant good, etc.
This year  excellent seems to mean good or worse,
good means poor (mud) etc.
Are they desperate for the money?
Is it the state of the Vic' economy?

Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Wed, 12 Aug 1992 10:57:58 GMT

Received: by utstat; Wed Aug 12 05:37 EDT 1992
Via: uk.ac.glasgow.geology; Wed, 12 Aug 1992 10:36:56 +0100
Received: from slioch.geology.glasgow.ac.uk by suilven.geology.glasgow.ac.uk;
          Wed, 12 Aug 92 10:34:37 BST
From: Colin Farrow <colin@geology.glasgow.ac.uk>
Date: Wed, 12 Aug 92 10:38:19 BST
Message-Id: <6055.9208120938@slioch.geology.glasgow.ac.uk>
To: s-news@utstat.toronto.edu
Subject: Re: filled circles


Debby
dfs@bellcore.com asks

I'd like to plot filled circles in S.  Is there an extension
to symbols() or polygon(), or a different function entirely,
that would let me do that?

How about

filled.circle <- function(cx, cy, r, ...)
{
	# cx, cy, coordinates for centre; r is radius
        z <- (0:360 * pi)/180
        x <- sin(z) * r + cx
        y <- cos(z) * r + cy
        polygon(x, y, ...)
        return(invisible())
}
hence
> par(usr=c(0,100,0,100),pty="s")
> filled.circle(20,20,10,col=2)

Colin Farrow,Dept of Geology,Univ of Glasgow,Glasgow,G12 8QQ
Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Fri, 14 Aug 1992 20:18:11 GMT

Received: by utstat; Fri Aug 14 14:53 EDT 1992
Received: by atc.boeing.com (5.57) 
	id AA23493; Fri, 14 Aug 92 11:48:02 PDT
Received: from trifle.bcs.eca by espresso.bcs.eca (4.1/SMI-4.1)
	id AA05835; Fri, 14 Aug 92 11:53:31 PDT
Date: Fri, 14 Aug 92 11:53:31 PDT
From: tjoelker@espresso.boeing.com (R. Tjoelker 865-3197)
Message-Id: <9208141853.AA05835@espresso.bcs.eca>
To: s-news@utstat.toronto.edu
Subject: predict.lm



I'm working with Splus3.0 on the Sun.

While trying out "lm", "predict" and "pointwise" to
get a confidence bound for the predicted mean value for a
simple linear regression equation, the answers do not
seem to agree with those calculated by hand.

I think there may be an error in the code for "predict.lm"

Looking more closely at the code for "predict.lm", I see
this near the end of the function:

        if(se.fit) {
                fit.summary <- summary.lm(object)
>>>>            pred <- Build.terms(x, coefs, fit.summary$cov * fit.summary$
>>>>                    sigma, asgn, collapse = type != "terms")
                pred$residual.scale <- fit.summary$sigma
                pred$df <- object$df.resid
        }

"Build.terms" calculates standard errors as:

                var <- ((x %*% cov) * x) %*% rep(1, length(coefs))
                list( ...,  se.fit = drop(sqrt(var))

So shouldn't the code in "predict.lm" be:

>>>>            pred <- Build.terms(x, coefs, fit.summary$cov * fit.summary$
>>>>                    sigma^2, asgn, collapse = type != "terms")
                        ^^^^^^^

After making the change, the answers agree with the hand calculations.

 #  #  #  #  #  #  #  #  #  #  #  #  #  #  #  #  #  #  #  #

Rod Tjoelker
Boeing Computer Services
(tjoelker@espresso.boeing.com)


Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Sat, 15 Aug 1992 01:02:20 GMT

Received: by utstat; Fri Aug 14 20:03 EDT 1992
Received: from tajo.ucsd.edu by ucsd.edu; id AA10606
	sendmail 5.67/UCSD-2.2-sun via SMTP
	Fri, 14 Aug 92 17:03:46 -0700 for s-news@utstat.toronto.edu
Received: by tajo.UCSD.EDU (4.1/UCSDGENERIC.4)
	id AA03124 to s-news@utstat.toronto.edu; Fri, 14 Aug 92 17:10:19 PDT
Date: Fri, 14 Aug 92 17:10:19 PDT
From: cberry@tajo.UCSD.EDU (Charles C. Berry)
Message-Id: <9208150010.AA03124@tajo.UCSD.EDU>
To: s-news@utstat.toronto.edu
Subject: coercing character to numeric mode

(Using Splus 3.0)

Is there any good reason for the following?

> as.complex(real="3*nothing")	#I WOULD HOPE FOR "NA"
[1] NA
Warning messages:
  1 missing value  generated coercing to complex in: as.complex(real = "3*nothing")
> as.integer("3*nothing")	#HERE TOO, LET IT BE "NA"
[1] NA
Warning messages:
  1 missing value  generated coercing to integer in: as.integer("3*nothing")
> as.numeric("3*nothing")	#THIS ISN'T WHAT I WANTED
[1] 3

(BTW, In 1988 S, these all yield 3)

I would hope that coercing a non-numeric string which starts with a numeral to integer, complex, or numeric would yield "NA" and it does for integer and complex, but not for numeric.

Is this more than a mere bug?

 _______________________________________________________________________________
|			|							|
|Charles C. Berry	|	Department of Community and Family Medicine	|
|			|	University of California, San Diego		|
|ccberry@ucsd.edu	|	La Jolla, San Diego 92093-0622			|
|cberry@tajo.ucsd.edu	|	(619) 534-2098 (voice) (619) 534-4642 (FAX)	|
|_______________________|_______________________________________________________|
 
Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Thu, 20 Aug 1992 21:23:53 GMT

Received: by utstat; Thu Aug 20 16:03 EDT 1992
Received: by uhunix.uhcc.Hawaii.Edu (4.1/Sun690)
	id AA13038; Thu, 20 Aug 92 10:03:01 HST
Date: Thu, 20 Aug 92 10:03:01 HST
From: Jerry Wetherall  <jerryw@uhunix.uhcc.Hawaii.Edu>
Message-Id: <9208202003.AA13038@uhunix.uhcc.Hawaii.Edu>
To: s-news@utstat.toronto.edu
Subject: s-plus code for geostatistics


We're investigating spatial patterns in the abundance of fish based
on large-scale open ocean surveys.  We'd be eternally grateful for leads
on s-plus code for kriging or other geostatistical methods.

Jerry A. Wetherall
Honolulu Laboratory
National Marine Fisheries Service
2570 Dole Street
Honolulu, Hawaii 96822 U.S.A.
jerryw@uhunix.uhcc.hawaii.edu
Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Wed, 19 Aug 1992 15:22:20 GMT

Received: by utstat; Wed Aug 19 09:33 EDT 1992
Received: from Uctvax.UCT.AC.ZA by Uctvax.UCT.AC.ZA (PMDF #2607 ) id
 <01GNRR8ZTGTS8WW5ZE@Uctvax.UCT.AC.ZA>; Wed, 19 Aug 1992 15:32:41 -0200
Date: 19 Aug 1992 15:32:41 -0200
From: TREVOR@Uctvax.UCT.AC.ZA
To: S-news@utstat.toronto.edu
Message-id: <01GNRR8ZTQGY8WW5ZE@Uctvax.UCT.AC.ZA>
X-Envelope-to: S-news@utstat.toronto.edu
X-VMS-To: S_NEWS
MIME-version: 1.0
Content-transfer-encoding: 7BIT

Rod Tjoelker wrote:
>While trying out "lm", "predict" and "pointwise" to
>get a confidence bound for the predicted mean value for a
>simple linear regression equation, the answers do not
>seem to agree with those calculated by hand.
>
>I think there may be an error in the code for "predict.lm"

He then went into detail about the error. The crucial
line is:
So shouldn't the code in "predict.lm" be:

>>>>            pred <- Build.terms(x, coefs, fit.summary$cov * fit.summary$
>>>>                    sigma^2, asgn, collapse = type != "terms")
                        ^^^^^^^

He is quite correct. In fact we picked up this error and
the April 92 version of AT&T's S reflects the correction.

Bob Treder of Statsci found a bug in the way summaryi() reports
the non-linear F-statistic for "gaussian" gam models. 
The following correction should be made to summary.gam():
    aod[snames, 3] <- nl.chisq
should be changed to:
	aod[snames, 3] <- if(!chiorf) nl.chisq/nldf else nl.chisq

Trevor Hastie 
trevor@research.att.com



Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Wed, 19 Aug 1992 08:33:33 GMT

Received: by utstat; Wed Aug 19 03:15 EDT 1992
Received: from [131.174.80.1] by hermes.stat.wisc.edu; Wed, 19 Aug 92 02:15:53 -0500
Received: by wn1.sci.kun.nl (5.57/2.1) on NUNET
 	id AA16002; Wed, 19 Aug 92 09:18:11 +0200
Received: from eos.cs.kun.nl by cs.kun.nl (4.1/SMI-3.2)
	id AA05637; Wed, 19 Aug 92 09:18:09 +0200
Message-Id: <9208190718.AA05637@cs.kun.nl>
To: S-news@stat.wisc.edu
Cc: fransp@cs.kun.nl (Frans Panken)
Subject: Real time Simulation with S
Date: Wed, 19 Aug 92 09:18:15 +0200
From: fransp@eos.cs.kun.nl

Hello world,

I want to simulate some stochastic processes 
(e.g. the waiting time distribution function 
of customers, when they arrive in a queue) 
and wonder if there is anybody who has 
experiances doing this in S.

Frans Panken
The Netherlands.
E-mail : fransp@cs.kun.nl 
Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Wed, 12 Aug 1992 20:23:54 GMT

Received: by utstat; Wed Aug 12 14:33 EDT 1992
Received: from CU.NIH.GOV ([128.231.64.7]) by relay.cs.toronto.edu with SMTP id <150177>; Wed, 12 Aug 1992 14:30:27 -0400
To:	s-news@utstat.toronto.edu
From:	TCF@CU.NIH.GOV
Date:	Wed, 12 Aug 1992 14:27:49 -0400
Subject:  Receiving duplicate S-NEWS messages
Acknowledge-to:    JJ6@NIHCU.BITNET
X-Acknowledge-to:  JJ6@NIHCU.BITNET
Message-Id: <92Aug12.143027edt.150177@relay.cs.toronto.edu>

Please stop my subscription to S-NEWS under these initials
(tcf@cu.nih.gov), as I am now receiving identical copies
via JJ6.  Thank you.

Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Wed, 19 Aug 1992 15:10:41 GMT

Received: by utstat; Wed Aug 19 09:47 EDT 1992
Received: from Uctvax.UCT.AC.ZA by Uctvax.UCT.AC.ZA (PMDF #2607 ) id
 <01GNRR8ZTGTS8WW5ZE@Uctvax.UCT.AC.ZA>; Wed, 19 Aug 1992 15:32:41 -0200
Date: 19 Aug 1992 15:32:41 -0200
From: TREVOR@Uctvax.UCT.AC.ZA
To: S-news@utstat.toronto.edu
Message-id: <01GNRR8ZTQGY8WW5ZE@Uctvax.UCT.AC.ZA>
X-Envelope-to: S-news@utstat.toronto.edu
X-VMS-To: S_NEWS
MIME-version: 1.0
Content-transfer-encoding: 7BIT

Rod Tjoelker wrote:
>While trying out "lm", "predict" and "pointwise" to
>get a confidence bound for the predicted mean value for a
>simple linear regression equation, the answers do not
>seem to agree with those calculated by hand.
>
>I think there may be an error in the code for "predict.lm"

He then went into detail about the error. The crucial
line is:
So shouldn't the code in "predict.lm" be:

>>>>            pred <- Build.terms(x, coefs, fit.summary$cov * fit.summary$
>>>>                    sigma^2, asgn, collapse = type != "terms")
                        ^^^^^^^

He is quite correct. In fact we picked up this error and
the April 92 version of AT&T's S reflects the correction.

Bob Treder of Statsci found a bug in the way summaryi() reports
the non-linear F-statistic for "gaussian" gam models. 
The following correction should be made to summary.gam():
    aod[snames, 3] <- nl.chisq
should be changed to:
	aod[snames, 3] <- if(!chiorf) nl.chisq/nldf else nl.chisq

Trevor Hastie 
trevor@research.att.com



Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Wed, 19 Aug 1992 20:34:23 GMT

Received: by utstat; Wed Aug 19 14:44 EDT 1992
Received: from vms3.macc.wisc.edu by hermes.stat.wisc.edu; Wed, 19 Aug 92 13:44:22 -0500
Received: from URIACC by WISCMACC.BitNet; Wed, 19 Aug 92 13:44 CDT
Received: by URIACC (Mailer R2.07) id 0548; Wed, 19 Aug 92 14:43:35 EDT
Message-Id: <22081913441823@WISCMACC.BitNet>
Date: Wed, 19 Aug 92 14:37:12 EDT
From: Marshall Feldman <MARSH%URIACC.BitNet@vms.macc.wisc.edu>
Subject: Replacing cells in a data frame (was matrix)
To: S-News Mailing List <S-NEWS@stat.wisc.edu>

I told you it was a stupid question.  Seriously, thanks to everyone
who suggested something like:
 
m [ m == -9999] <- NA
 
where m is the matrix.
 
The reason I asked this stupid question is that I was working with a
data frame and was able to perform most matrix operations on it.  I
forgot I was working with a data frame.
 
So far the best I've been able to do with the data frame is:
 
      for (i in 1:ncol(m))
          m[[i]] [(1:nrow(m)) [ m [[i]] == -9999 ] <- NA
 
Awfully cumbersome and compute intensive.  Any ideas for a more efficient
way?  Would
 
      m <- data.matrix(m)
      m [ m == -9999 ] <- NA
      m <- data.frame(m)
 
be more efficient or less?
 
Marsh Feldman
Community Planning                      Phone: 401/792-2248
204 Rodman Hall                           FAX: 401/792-4395
University of Rhode Island           Internet: marsh@uriacc.uri.edu
Kingston, RI 02881-0815
Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Tue, 18 Aug 1992 21:22:19 GMT

Received: by utstat; Tue Aug 18 15:48 EDT 1992
Received: from iso.wwz.unibas.ch by swiba9.unibas.ch with SMTP (PP) 
          id <7181-0@swiba9.unibas.ch>; Tue, 18 Aug 1992 21:48:29 +0200
Received: by iso.wwz.unibas.ch (5.57/Ultrix3.0-C) id AA23949;
          Tue, 18 Aug 92 21:47:57 +0200
Message-Id: <9208181947.AA23949@iso.wwz.unibas.ch>
To: S News List <s-news@utstat.toronto.edu>
Subject: assigning to a list element
Date: Tue, 18 Aug 92 21:47:55 N
From: andreas@iso.wwz.unibas.ch


When usingthe assign function, it seems to me that a list element can't be
assigned a value. This is mentione din th emanual as

       assign("a$b",1)

       In the last example, note that this is distinctly not  the
       way  to  create  the  component b of object a, as a$b <- 1
       would do.

but that won't help me further.  Any way to work around this ?
Example : 

> z <- list(x=4)
> z
$x:
[1] 4

> assign("z$x",33)
> z
$x:
[1] 4

So what's the problem ? Thanks for any help.
Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Thu, 13 Aug 1992 21:00:31 GMT

Received: by utstat; Thu Aug 13 15:43 EDT 1992
Date: Thu, 13 Aug 92 15:43:21 EDT
From: zaslavsk@hustat.harvard.edu (Alan M. Zaslavsky)
Message-Id: <9208131943.AA08276@hustat.harvard.edu>
To: jmj@dopee.mlb.semi.harris.com, s-news@utstat.toronto.edu
Subject: Re: Problem with lm and get

> Can someone help me with the following:
> 
> Why does
> 
> 	lm(get("parameter1") ~ . , dataframe)
> 
> include parameter1 as a term in the analysis and thus invalidate 
> the results?  Is there anyway to model a parameter using get
> without having to list all the other parameters to be fitted?
> I tried
> 
> 	lm(get("parameter1" ~ . - get("parameter1") , dataframe) 
> 
> but it doesn't work either.  I'm presently using Splus 3.0.
> 
> Jeff Johnston
> Harris Semiconductor
> jmj@mlb.semi.harris.com
> 
> 
The lm() function is not recognizing get("par1") as par1 when it figures out
what to exclude from the . on the right hand side.  I assume that the reason
for this construct is that you have a parameter specified by a character
string in a variable that you want to be the left hand side.  How about
wrapping the following in a function:

myformula_function(whatever)
{
formula1_dummy~.
formula1[[2]]_as.name(whatever)
formula1
}

lm(myformula("parameter1"), dataframe)

 
Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Mon, 17 Aug 1992 18:06:22 GMT

Received: by utstat; Mon Aug 17 12:50 EDT 1992
Received: from ATHENA.MIT.EDU by hermes.stat.wisc.edu; Mon, 17 Aug 92 11:50:16 -0500
Received: from E51-007-3.MIT.EDU by Athena.MIT.EDU with SMTP
	id AA02662; Mon, 17 Aug 92 12:50:11 EDT
From: apk@Athena.MIT.EDU
Received: by e51-007-3 (5.57/4.7) id AA06346; Mon, 17 Aug 92 12:50:08 -0400
Message-Id: <9208171650.AA06346@e51-007-3>
To: s-news@stat.wisc.edu
Subject: please add me to the s-news list
Date: Mon, 17 Aug 92 12:50:07 EDT

Subject:please add me to the s-news list
apk@athena.mit.edu
Thanks!
Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Thu, 13 Aug 1992 15:17:39 GMT

Received: by utstat; Thu Aug 13 09:38 EDT 1992
Received: from dopee.mlb.semi.harris.com by slopoke.mlb.semi.harris.com (4.0/SMI-4.0)
	id AA13847; Thu, 13 Aug 92 09:38:03 EDT
Received: by dopee.mlb.semi.harris.com (4.1/SMI-4.1)
	id AA09989; Thu, 13 Aug 92 09:38:02 EDT
Date: Thu, 13 Aug 92 09:38:02 EDT
From: jmj@dopee.mlb.semi.harris.com (Jeff Johnston)
Message-Id: <9208131338.AA09989@dopee.mlb.semi.harris.com>
To: s-news@utstat.toronto.edu
Subject: Problem with lm and get

Can someone help me with the following:

Why does

	lm(get("parameter1") ~ . , dataframe)

include parameter1 as a term in the analysis and thus invalidate 
the results?  Is there anyway to model a parameter using get
without having to list all the other parameters to be fitted?
I tried

	lm(get("parameter1" ~ . - get("parameter1") , dataframe) 

but it doesn't work either.  I'm presently using Splus 3.0.

Jeff Johnston
Harris Semiconductor
jmj@mlb.semi.harris.com

Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Wed, 19 Aug 1992 16:41:16 GMT

Received: by utstat; Wed Aug 19 11:15 EDT 1992
Received: from vms3.macc.wisc.edu by hermes.stat.wisc.edu; Wed, 19 Aug 92 10:15:52 -0500
Received: from URIACC by WISCMACC.BitNet; Wed, 19 Aug 92 10:15 CDT
Received: by URIACC (Mailer R2.07) id 8844; Wed, 19 Aug 92 11:14:48 EDT
Message-Id: <22081910152900@WISCMACC.BitNet>
Date: Wed, 19 Aug 92 11:13:25 EDT
From: Marshall Feldman <MARSH%URIACC.BitNet@vms.macc.wisc.edu>
Subject: Replacing cells in a matrix
To: S-News Mailing List <S-NEWS@stat.wisc.edu>

This is probably a stupid question, but I can't figure out a good way
do to this.
 
I have a large matrix with some cells set to -9999.  I want to
replace this with NA's.  How do I do this on an individual cell basis?
 
Marsh Feldman
Community Planning                      Phone: 401/792-2248
204 Rodman Hall                           FAX: 401/792-4395
University of Rhode Island           Internet: marsh@uriacc.uri.edu
Kingston, RI 02881-0815
Newsgroups: aus.stats.s
From: s-sender@utstat.toronto.edu
Subject: news from S-news
X-Gateway: dmssyd.syd.dms.CSIRO.AU
Organization: CSIRO Div. Maths & Stats, CSIRO Div. Applied Physics
Date: Fri, 14 Aug 1992 10:35:10 GMT

Received: by utstat; Fri Aug 14 04:39 EDT 1992
X400-Received: by mta chx400.switch.ch in /PRMD=switch/ADMD=arcom/C=CH/;
               Relayed; Fri, 14 Aug 1992 10:38:54 +0200
X400-Received: by /PRMD=switch/ADMD=arcom/C=CH/; Relayed;
               Fri, 14 Aug 1992 10:38:28 +0200
X400-Received: by /PRMD=switch/ADMD=arcom/C=ch/; Relayed;
               Fri, 14 Aug 1992 11:37:01 +0200
Date: Fri, 14 Aug 1992 11:37:01 +0200
X400-Originator: molinari@cumuli.ethz.ch
X400-Recipients: s-news@utstat.toronto.edu
X400-MTS-Identifier: [/PRMD=switch/ADMD=arcom/C=ch/;920814103701]
X400-Content-Type: P2-1984 (2)
Content-Identifier: 1080
Conversion: Prohibited
From: Luciano Molinari <molinari@cumuli.ethz.ch>
To: s-news <s-news@utstat.toronto.edu>

I am looking for an MS-DOS utility to dearchive and decompress Unix tar.Z
files.
Can anybody help me?
Thanks,
L. Molinari
molinari@cumuli.ethz.ch
The Children's Hospital
Steinwiesstr.75
CH-8032 Zurich

Newsgroups: aus.stats.s