similar to: Problem with read.spss() and as.data.frame(), or: alternative to subset()?

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 #

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

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

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

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

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

Begin forwarded message: > Subject: R: Studdy Missing Data, differentiate between a percent > with in the valid answers and with in the different missing answers > take a look at the Website from Dirk Enzmann: > http://www2.jura.uni-hamburg.de/instkrim/kriminologie/Mitarbeiter/ > Enzmann/Software/Enzmann_Software.html > He provides a function called "freq.r" - maybe

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:

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) %*%

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

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"

Hello from ealier 6.0 there is problem with synproxy in pf filter: this one 6.1-PRERELEASE #2: Wed Mar 15 02:02:37 MSK 2006 pf.conf just with single rule pass in quick on lo0 proto tcp from any to any port 22 flags S/SA synproxy state result telnet 127.0.0.1 22 Trying 127.0.0.1... Connected to 127.0.0.1. Escape character is '^]'. and it's hangs pfctl -s rules -v No ALTQ support in

Hi I finally did the merge for Xen 4.0. For this I decided to rename the source package back to xen. The packages are available via http://hermes.jura.uni-tuebingen.de/~blank/debian/xen-test. ec1482e14f2515bd6878da0d2229b725993601fa libxen-dev_4.0.0-1~experimental.1_amd64.deb 7c394f7158ca387225935f3bbe26bdc5518a1934 libxenstore3.0_4.0.0-1~experimental.1_amd64.deb

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 |

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? >

When I use the hills data in Maindonald's text (below) and apply coplot(time ~ distance | climb, data=hills, overlap=0.0, number=3) the graphs do not appear to line up with the conditional axis climb. Is this a bug or am I missing something as an R newbie? Thanks, Bob Abugov distance climb time Greenmantle 2.5 650 16.083 Carnethy 6.0 2500

=================================================================== "Change is such hard work." Billy Crystal =================================================================== Release Announcements ===================== This is is the latest stable release of Samba 3.6. Major enhancements in Samba 3.6.7 include: o Fix resolving our own "Domain Local" groups