similar to: Supplemental groups in samba

I am trying to fit a linear mixed model of the form y_ij = X_ij \beta + delta_i + e_ij where e_ij ~N(0,s^2_ij) with s_ij known and delta_i~N(0,tau^2) I looked at the ecme routine in package:pan, but this routine does not allow for different Vi (variance covariance matrix of the e_i vector) matrices for each cluster. Is there an easy way to fit this model in R or should I bite the bullet and

Hi, I have a share in my smb.conf file that I want to share to members of a supplemental group, but I keep getting a network error. here is the share: [mp3s] comment = MPEG3 Encoded Audio path = /usr/share/mp3 public = no browsable = no valid users = +mp3users writable = yes create mask = 0644 I have a group called mp3users that inside my /etc/group file looks like this:

I've downloaded and installed R 1.9.0 under OS X 10.3.3. I then removed my previous R 1.8.0 installation by following the instructions in the readme file, e.g., 1. drag the "StartR" icon into the Trash; 2) go inside your Home directory and search for the Library folder. Inside the Library folder in your Home directory you'll find a RAqua directory: drag the RAqua

I've installed R 1.8.1 on OSX 10.3 (Panther) using the RAqua.pkg and all indications were that the installation was successful. However, after double clicking the R icon in the Applications folder, nothing appears to happen. The following message appears on the console: Mac OS X Version 10.3 (Build 7B85) 2003-12-11 08:45:31 -0500 prefs written dyld:

I've installed R 1.8.1 on OSX 10.3 (Panther) using the RAqua.pkg and all indications were that the installation was successful. However, after double clicking the R icon in the Applications folder, nothing appears to happen. The following message appears on the console: Mac OS X Version 10.3 (Build 7B85) 2003-12-11 08:45:31 -0500 prefs written dyld:

Hello, I'm trying to install a package that requires a Fortran compiler (Hmisc) using R CMD INSTALL. I downloaded the package source onto my Desktop, unzipped it, and then typed: R CMD INSTALL /Users/brianbeckage/Desktop/Hmisc * Installing *source* package 'Hmisc' ... ** libs g77 -fno-common -g -O2 -c cidxcn.f -o cidxcn.o g77 -fno-common -g -O2 -c cidxcp.f -o cidxcp.o g77

Suggestions from both Prof Brian Ripley and Gabor Grothendiec solved my problem. From Prof Ripley: >This is a problem of your OS: your example works on all of mine. > >Can't you change the origin in chron? If not, you should certainly be >able to do > >as.POSIXct(strptime(as.character(testDATES), "%m/%d/%y")) This still resulted in NA's on my system. >

How can I extract K (number of parameters) from an AIC calculation, both to report K itself and to calculate AICc? I'm aware of the conversion from AIC -> AICc, where AICc = AIC + 2K(K+1)/(n-K-1), but not sure of how K is calculated or how to extract that value from either an AIC or logLik calculation. This is probably more of a basic statistics question than an R question, but I thank

I'm having trouble loading Ryacas. I've downloaded and extracted Ryacas 0.2-9 (also tried Ryacas ) and yacas 1.0.63, have the latest version of R and have tried the following (this works for installing other packages): install.packages("Ryacas") library(Ryacas) install.packages("yacas") library(yacas) and get this: Warning message: package ‘yacas’ is not available

Is it possible to give a temporary object a name that varies with each run of a foreloop? For example, I want to fill a matrix every time I run a loop, and I want a new matrix with each run, with an appropriate new name. i.e.: for(i in 1:5){... matrix.i<-some values ...} so that in the end I would have: matrix.1 matrix.2 matrix.3 matrix.4 matrix.5 Thanks, Ben Osborne -- Botany Department

When reading a data set into R using readLines, na.strings="-99.99" is ignored. Is there an equivalent command for readLines? Alternatively, either immediately after reading into R or once the data set has been converted to a data frame, what is the appropriate command (or appropriate use of is.na) to convert my -99.99s to NAs? Thanks, Ben Osborne -- Botany Department University of

What is the difference between the two R-squareds returned for a linear regression by summary(lm...)? When might one report multiple vs. adjusted R-squared? Thank you, Ben Osborne -- Botany Department University of Vermont 109 Carrigan Drive Burlington, VT 05405 benjamin.osborne at uvm.edu phone: 802-656-0297 fax: 802-656-0440

Hi all. I have currently a server: cat /etc/redhat-release CentOS release 5.7 (Final) uname -a Linux host.domain.com 2.6.18-274.18.1.el5 #1 SMP Thu Feb 9 12:45:44 EST 2012 x86_64 x86_64 x86_64 GNU/Linux I have there a filesystem mounted: /dev/vg0/paczki /home/paczki-workdir ext4 defaults,noatime 0 0 on which df gives strange output: LANG=C df -h

I'm attempting to do model selection with AIC, using a glm and a lognormal distribution, but: fit1<-glm(BA~Year,data=pdat.sp1.65.04, family=gaussian(link="log")) ## gives the same result as either of the following: fit1<-glm(BA~Year,data=pdat.sp1.65.04, family=gaussian) fit1<-lm(BA~Year,data=pdat.sp1.65.04) fit1 #Coefficients: #(Intercept) Year2004 # -1.6341

Does anyone know of a list of functions that R already "knows?" I can't seem to find this anywhere in the help documentation. For example, I want to count the number of occurences of a certain value in a column of a data frame: What do I have to do to tell R to "Count?" Thanks, Ben Osborne -- Botany Department University of Vermont 109 Carrigan Drive Burlington, VT 05405

I want to organize outputs from several regressions into a handy table. When I try the following, each of my "fit_s" is replaces instead of read. Is there a way to read from the regression summaries that does not require writing separate lines of code for each? -Ben Osborne > fit1<-lm(dBA.spp16$sp2.dBA.ha~dBA.spp16$sp1.dBA.ha) >

Hello, thank you very much for your previous answers about the C++ code. I am interested in the application of the Gibbs Sampler in the IRT models, so in the function MCMCirt1d and MCMCirtkd. I've found the C++ source codes, as you suggested, but I cannot find anything about the Gibbs Sampler. All the files are for the Metropolis algorithm. Maybe I am not able to read them very well, by the

1) When adding additional data sets to a plot using "plot" followed by "lines", is there a way to automate the scaling of the axes to allow for all data sets to fit within the plot area? 2) I attempted to solve this by setting xlim=c(min(c(data1,data2,data3)),max(c(data1,data2,data3))) however, there are some NAs and Infs in these data sets, and min(data1) and max(data1) both

Hi all, Earlier today I posted this question on s-news, so apologies to some for the duplication. > Please put aside your snobbery about Microsoft products for a moment. > > I am fitting population models to annual survey data for trout. For those of you familiar with ecological models, I am working in the Lefkovitch matrix framework; for those unfamiliar with that shorthand, the

I have a complex function that I want to maximize (I have multiplied this function by -1 so that it becomes a minimization problem in the code below). This function has two equality constraints. I get the programs to run but the answer isn't correct because, when it does converge, at least one of the constraints is violated. Any suggestions? Code below Violated constraint (an easy check):