similar to: is zscore() deprecated or Windows only?

Hi, I am writing a basic function to extract the z scores for some linear regression coefficients: zscore<-function( y, x) { lm<-lm( y ~ x ) z <- coef(lm)/sqrt(diag(vcov(lm))) return(z) } I would like to pass a dataframe to the function as a argument so the function call changes from zscore(df$y1,df$x1) to zscore(y1,x1,data=df) but I am not sure how to reference the data

Hi, R 1.0.1, on Linux(mandrake/RedHat). I am trying to write my first function using R and I need help figuring out the error I have indicated in the Subject line above. The purpose of the function is not to "wow" anyone, but just to serve as a stepping stone so that I might write more complicated functions in the future. First I created my function named "CIfunc" and saved

Dear R Help list, My question is regarding extracting the standard error or Z-score from a cph or coxph call. My Cox model is: - modz=cph(Surv(TSURV,STATUS)~RAGE+DAGE+REG_WTIME_M+CLD_ISCH+POLY_VS, data=kidneyT,method="breslow", x=T, y=T) I've used names(modz) but can't see anything that will let me extract the Z scores for each coefficient or the standard errors in the same

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi, I just installed R (2.11.1) 64bit under Windows 7 64 bit When trying to readin files I get the following error message: > > source(file="C:\\Users\\me\\Documents\\My PhD\\Modelling\\R\\Scripts\\Functions\\functions 20jan10.r") Error in file(file, "rt") : cannot open the connection In addition: Warning message: In

Hi all, I have played a bit with the "reshape" package and function along with "melt" and "cast", but I feel I still don't have a good handle on how to use them efficiently. Below I have included a application of "reshape" that is rather clunky and I'm hoping someone can offer advice on how to use reshape (or melt/cast) more efficiently. #For this

what does this mean in R-1.8.1 release notes? o median() no longer `works' for odd-length factor variables.

Hi- I'm wrestling with some of my data apparently not being called into a GLM or an LM. I'm looking at factors affecting fish annual catch rates (ie. CPUE) over 30 years. Two of the factors I'm using are sea surface temperature and sea surface temperature anomaly. A small sample of my data is below: CPUE Year Vessel_ID Base_Port Boat_Lgth Planing SST Anomaly 0.127

Forwarded Message: > To: Steve Arthur <sarthur at crick.protogene.com> > From: Thomas Lumley <thomas at biostat.washington.edu> > Subject: Re: [R] couldn't find function "..." > Date: Wed, 5 Jul 2000 16:16:25 -0700 (PDT) > ----- > On Wed, 5 Jul 2000, Steve Arthur wrote: > > > Hi, > > > > R 1.0.1, on Linux(mandrake/RedHat).

Hi All: I was looking for, but could not locate in the packages, or in the R archive searches if there exists an R implementation of Cuzick's test of trend. The test is described as follows: An extension of the Wilcoxon rank-sum test is developed to handle the situation in which a variable is measured for individuals in three or more (ordered) groups and a non-parametric test for trend

I have the following function for which I need to find the root of a: f <- function(R,a,c,q) sum((1 - (1-R)^a)^(1/a)) - c * q To give context for the problem, this is a psychometric issue where R is a vector denoting the percentage of students scoring correct on test item i in class j, c is the proportion correct on the test by student k, and q is the number of items on the test in total. I

Greetings- a reviewer for a paper of mine noted an anomaly in some models I ran using glmmPQL (from the MASS package). Specifically, the models are three-level hierarchical probit models estimated using PQL under R. The anomaly is that the log-likelihoods decrease (or, alternatively -2logLik increases) as variables are added to the null model. This is unusual, and I'm trying to figure out

circ.plot in CircStats under R 1.8.0 yields incorrect plots (easily compare hist() and rose.diag()). However, the code in circ.plot has not changed. Has anyone tracked this down yet?

Dear all, I'm trying to make an R code for the drought severity index (DSI) developed by Philips and McGregor (1998). You can refer to the description of the algorithm on page 19 from http://dissertations.port.ac.uk/229/01/MaloneS.pdf The code is given in Excel as the following and can be found on page 60 from the same link. C7 = IF(C6<0,IF(@SUM(A6:A1)<0,C6+A6,"0"),

Full_Name: Michael W. Grant Version: 1.7.1(2003-06-16) OS: Windows 98 Second Edition 4.10.2222 A Submission from: (NULL) (24.55.165.101) I just installed version 1.7.1 tonight. In the process I installed via CRAN from the menu several packages. One installation, exactRankTests, had a distinct anomaly. A DLL attempt to change FPU control word as seem in the last line of the following listing:

I am having consistent difficulty with write access and installing applications. The errors show up as write protect errors or not enough disk space errors. 'winefile' shows the disks as being NTFS. Is this correct or is it a winefile anomaly? Are these conditions related? What is the correct file system type? How can this be managed? 'winecfg' does not seem to allow

We're running into an issue where sendmail is only receiving one e-mail, but after a client accesses the message via IMAP and moves it to another IMAP folder, the message duplicates but each message has different X-UID headers. Has anyone ever seen this type of anomaly before? Is it a bad IMAP client, or a known bug? Running 1.0.7 Thanks!

Dear R-helpers, Please have a look at the following. f1 is the same as f2 except that it has some values replaced by NA. But it's corresponding file is slightly bigger than the file containing f2. Could someone please tell me if this is an anomaly ? > load("file1") > ls() [1] "f1" > load("file2") > ls() [1] "f1" "f2" > >

How does one go about doing a repeated measure regression? The documentation I have on it (Lorch & Myers 1990) says to use linear / (subj x linear) to get your F. However, if I put subject into glm or lm I can't get back a straight error term because it assumes (rightly) that subject is a nominal predictor of some sort. In looking at LME it seems like it just does the right thing

We have encountered an anomaly that appears to be related somehow to Samba. When users are moved to the new domain hosted by Samba4.3 then the numlock on their keyboards is always turned off when they log onto the domain regardless of its previous setting. This in itself is just a nuisance. However we have an application called SnetTerm, an ssh client, which will not recognise the numlock at

Hello, trying to get this game to run on wine 1.3.16 on a debian sid amd64 box. This is what the error output when running it : Code: $ wine /media/Data/Program\ Files/Anomaly\ -\ Warzone\ Earth/AnomalyWarzoneEarth.exe err:winedevice:ServiceMain driver L"SecDrv" failed to load wine: Unhandled page fault on read access to 0x00000000 at address 0x7e084c95 (thread 0009), starting