similar to: glim and gls

Hello, I'm just start learning R/S-Plus, so I think this 2 doubts is going to be too easy for you... 1) I couldn't discover what is the command for a concatenation of 2 variable strings. 2) For example, if I have three variable strings, and each one has the name of a variable in a data matrix: a<-V1 b<-V2 c<-V3 , is it possible to construct a command like this:

Hi, I'm trying to fit a marginal model via GEE but I'm getting strange results and few problems. If I set the working correlation as exchangeable I'm getting the same fitting when I set as independent. Comparing to SAS results it shouldn't happen. If I try to use another working correlation (like AR-M or stat_M_dep), R just exits without giving any error message. Another doubt

Dear useRs is there any nonparametric test for the analysis of variance in a design with two within-factors (repeated measures on both factors)? Friedman is not appropriate here, therefore I am grateful for any alternative test. thanks for any hint cheers christoph --

Dear R users, I have looked in the reference Schnabel, R. B., Koontz, J. E. and Weiss, B. E. (1985) A modular system of algorithms for unconstrained minimization. _ACM Trans. Math. Software_, *11*, 419-440. cited in the nlm help. This article says that the algorithm permits the use of step selection (line search, dogleg and optimal step), analytic or finite diference gradient

After some futile searches, I decided to ask the list to see if any of the sages out there would have an answer: I have a function I wrote a few years ago in S, which calls glim numerous times. I'd like to port it to R, but glm works differently from glim, which takes as part of its input an X design matrix. I probably could write a function to convert glim to glm, but hope this

Dear Experts, In glim() of S Plus, one of the returned values is "var", the estimated variance matrix of coefficients. However, in glm() of R (there is no glim() in R), "var" is not one of the returned values. Anyone know what could I get the varience matrix of coefficients in glm() in R? As a novice in R and S+, I'd appreciate your help Sincerely, Charlie Liu

Hi all, In my continuous transition of GLIM to R I try to make a glm with binomial errors. The data file have 3 vectors: h -> the factor that is ajusted (have 3 levels) d -> number of animais alive (the response) n -> total number of animals To test proportion of alive, make d/n. In GLIM: $yvar d$ $error binomial n$ $fit +h$ scale deviance = 25.730 (change = -9.138) at cycle 4

Dear All, I wonder if anyone can advise me as to whether there is a consensus as to how the effect size should be calculated from GLIM models in R for any specified significant main effect or interaction. In investigating the causes of variation in infection in wild animals, we have fitted 4-way GLIM models in R with negative binomial errors. These are then simplified using the STEP procedure,

Hi all, Probably a simple question, but I just can't find a simple answear in the older threads or anywhere else. I've added some new vectors as columns in a data frame using cbind(). As they're all put as the last columns inte the data frame, I would like to move them to specific positions. How do you do to change the position of a column in a data frame? I know I can use

Hi all I have another GLIM question. I have been using R as well as Genstat (version 6) in order to fit GLIM models to the data (displayed below). The same models are fitted but the answers supplied by the two packages are not the same. Why? Can anyone help? A discription of the data and the type of model/s fitted can be found below. Regards Allan The

Hi all, I guess this might be an easy question, but I've searched multiple help pages without finding any answear... so now I put my trust in you! I have a data frame (36 variables and 556 observations). One column contains three factors, and I would like to divide the data frame into three new ones, based on the value of the factors, thereby having only one value for all elements of the

Hello, I am just installing XCP 0.1.1 on a server with 24GB RAM Normally XCP is a 64 bits version but * *uname -a* Linux node012 2.6.27.42-0.1.1.xs0.1.1.737.1065xen #1 SMP Fri Jan 15 16:20:16 EST 2010 i686 i686 i386 GNU/Linux (not a x86_64 version !!!!) * *cat /proc/meminfo* MemTotal: 746496 kB MemFree: 384332 kB Buffers:

Hi! I have a runing asterisk box and i want to dial to a analog. pbx using a 4FXS Welltech. Let's say that my pbx have no. 700. If i want to dial to a person in that direction i have to dial pbx prefix (ex. 700), wait for pbx to ansear with hello message and after that to dial internal number(ex. 101). It is possible to dial directly 700101 and asterisk to dial PBX prefix, wait for PBX to

OK, I still have problems to run samba as PDC. I have followed the instructions from Samba By Examples chapter 6. Also I have configured dynamic dns + dhcp. But, I still doesn't work. I have possibility to logon on Windows machine 2 times. At the 3 i have recieve erroe message, that windows can not copy profile from server. I have check that at fist 2 times I am connecting as any user from

Rowland wrote: >Can I suggest you put your smb.conf back to what it was just after the >provision and try again. If it does work (and it should), before you >start adding lines again, can I also suggest you read 'man smb_conf', >most of the lines you have added are the defaults, others do nothing on >a DC and some are just plain wrong. > >Rowland Thank you for

Hi, First, I want to precise that I'm new to linux, so please answear with simple words ;) My problem is that I have to install and configure linux (RH9) on 36 identical computers. This would take ages. What I want to do, is to configure one, and then clone the install. But thos computers have no CD, no floppy, and I didn't manage to boot them brom an usb stick. The last option is to

Hi,guys, During a HVM''s running which take a qcow format image file as its own virtual disk, the qcow image file will be always read.In the situation that its qcow format image is based on a raw format image, if nesethe backingfile ,just that raw format image file,would be read .my purpose is to cache the data that is read from the backingfile when the hvm is running . Now what I concern

Hi, In the help files in the ?mgcv package for the gam.control() function, there is an option irls.reg. The help files describe this option as: For most models this should be 0. The iteratively re-weighted least squares method by which GAMs are fitted can fail to converge in some circumstances. For example, data with many zeroes can cause problems in a model with a log link, because a mean of

Hi all, I am also confused about the manual: a. The input arguments: wt.method are the weights case weights (giving the relative importance of case, so a weight of 2 means there are two of these) or the inverse of the variances, so a weight of two means this error is half as variable? w (optional) initial down-weighting for each case. init (optional) initial values for the

Dear users of the R mailing list, I am a ph.d. student in biology working on red deer in Norway, who would like to conduct an analysis with random factor where the response is binomially distributed. This cannot be conducted in S-plus, and I was told by others that it may be possible in R. However, I soon got into trouble which I hope you can help me to solve. My model is on this form: