similar to: Fwd: R: Studdy Missing Data, differentiate between a percent with in the valid answers and with in the different missing answers

Trying to select a subset of cases (rows of data) I encountered several problems: Firstly, because I did not read the help to read.spss() thoroughly enough, I treated the data read as a data frame. For example, dr2000 <- read.spss('myfile.sav') d <- subset(dr2000,RBINZ99 > 0) and thus received an error message (Object "RBINZ99" not found), because dr2000 is not a

I'm sorry, but my previous message, as often happens, some brackets were wrong: Here are the correct formulas: sd = 1/n * sum((x-mean(x))^2) # (1) sd = 1/(n-1) * sum((x-mean(x))^2) # (2) This also occured in the last paragraph. Dirk ************************************************* Dr. Dirk Enzmann Institute of Criminal Sciences Dept. of Criminology Edmund-Siemers-Allee 1 D-20146

Just out of curiosity (and I hope that my question does not confuse those who intuitively understood the par(new=TRUE) statement right in the beginning): I would like to know whether beginners (learning R) had the same difficulty as I had: Intuitively I thought that par(new=TRUE) would draw a new plot (and not into an already existing plot) and that par(new=FALSE) would not draw a new plot

I wonder whether the functions for skewness and kurtosis in the e1071 package are based on correct formulas. The functions in the package e1071 are: # -------------------------------------------- skewness <- function (x, na.rm = FALSE) { if (na.rm) x <- x[!is.na(x)] sum((x - mean(x))^3)/(length(x) * sd(x)^3) } # -------------------------------------------- and #

After loading the package "Matrix" (version 0.995-12), using the summary function with an lme (package nlme version 3.1-75) object results in an error message saying Fehler in dim(x) : kein Slot des Namens "Dim" f?r dieses Objekt der Klasse "correlation" (translated: 'Error in dim(x) : no slot of the name "Dim" for this object of class

Although I tried to find an answer in the manuals and archives, I cannot solve this (please excuse that my English and/or R programming skills are not good enough to state my problem more clearly): I want to write a function with an indeterminate (not pre-defined) number of arguments and think that I should use the "..." construct and the match.call() function. The goal is to write

Using ranef() (package nlme, version 3.1-75) with an 'lme' object I can obtain random effects for intercept and slope of a certain level (say: 1) - this corresponds to (say level 1) "residuals" in MLWin. Maybe I'm mistaken here, but the results are identical. However, if I try to get the standardized random effects adding the paramter "standard=T" to the

I receive R-help messages in digest form that makes it difficult to answer a post. I noticed that my answer is not added to the thread (instead, a new thread is started) although I use the same subject line (starting with "Re: ") as the original post. Is there a solution (I prefer the digest to separate messages for several reasons and don't want to change my email reader)? The

Hi R experts I'm trying to emigrate from SPSS to R, thou I have some problems whit getting R to distinguish between the different kind of missing. I want to distinguish between data that are missing because a respondent refused to answer and data that are missing because the question didn't apply to that respondent. In other words I wante to create data values where I control

I hesitate to post this question in the light of recent threads, indeed I have hesitated for several weeks, however I have come to a full stop and really need some help if I am going to progress. I am a new user of R for medical statistics. I have attempted to read all the relevant documents, but would welcome any suggestions as to what I have missed. I am trying to contruct "table 1"

The update of VR by downloading VR_7.1-6.zip and using install.packages (from local zip files) fails with the following error message: Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file `VR/DESCRIPTION' Other packages can be installed without problems, except of dse_2003.4-1 with a similar error message. Why? Operating System:

Can someone show me how to specify zero inflated poisson and zero inflated negative poisson models in R? I would like to replicate the example given in Long (1997: Regression Models for Categorical and Limited Dependent Variables) in Chapter 8.5 (pp. 242-247). TIA Dirk ************************************************* Dr. Dirk Enzmann Criminological Research Institute of Lower Saxony

Is it possible to avoid the loop in the following function (or make the function otherwise more efficient) and can someone point me to a possible solution? (It would be great if hours could be reduced to seconds :-). # --------------------------------------------- RanEigen=function(items=x,cases=y,sample=z) { X=matrix(rnorm(cases*items),nrow=cases,byrow=F) S=crossprod(X-rep(1,cases) %*%

Hi, I'm using R 2.5.1 and gstat 0.9-39 library and I'm working with the jura data set provided by gstat library. I tried to plot a graph of metal concentrations (let's say Cd) with the command spplot, but I realized that the default lags are equally distributed between Min and Max. I did: library(gstat) data(jura) jura.pred.xy <- jura.pred names(jura.pred.xy)[1:2] <-

Hi Dovecot fails to chdir() into users home directory giving the following error: Oct 9 16:43:49 jura dovecot: dovecot: Fatal: chdir(/a/home/mina) failed: Permission denied (euid=1118 0(mina) egid=11332(devel) missing +x perm: /a/home) Oct 9 16:43:49 jura dovecot: dovecot: child 14441 (imap) returned error 89 (Fatal failure) The reason seems to be that user primary group is "devel"

Enough small fixes have pilled up to make it worth a release: A Google patch adds sched_setaffinity, sched_getaffinity support. Openembedded uses kexec_load(). Gentoo folks add a Kbuild fix. ipconfig no longer wild guesses a nameserver when none is provided by the DHCP server. strndup() and unlinkat() saw fixes for various problems. codingstyle cleanup in kinit and tools. git repository:

Hello All, Is the application SetMusicOnHold always followed by the Dial application? I'm studdying the info on Google and Wikki about SetMusicOnHold, and I'm not sure I understand what I found. Thanks, Anon

Hello. We are running individual software firewalls (as provided by the respective OSs) on: MAC OSX Panther MSWindows XP Pro SuSE 8.2 LINUX Machines. In addition, the MSWindows XP Pro machine is running Norton's Software Firewall (in conjunction with Antivirus version 10.0.1.13) I am trying to configure SAMBA 2.2.8a from the SuSE LINUX system and am unable to see either the MAC OSX

samba@samba.org wrote: > SAMBA Digest 2150 > > For information on unsubscribing see http://samba.org/listproc/ > Topics covered in this issue include: > > 1) Autoconnect on restarts doesn't work > by =?iso-8859-1?Q?Str=F6mberg_Peter?= <Peter.Stromberg@scania.com> > 2) Re: How to backup from Win to Linux using smb? >

i use the CrossTable (frome gregmic package) function to get such a table as below. but the percentage of the non-NA levels(here 1,2,3,4,5) is not totally 100%. is there any way to get a table that percentage of the non-NA levelsis totally 100%,as the SPSS' valid percentage.thank you! Cell Contents |-------------------------| | Count | | Row Percent |