similar to: 2 simple doubts

Hello, I would like to know if there is any package that allow us to fit Generalized Linear Models via Maximum Likelihood and Linear Models using Generalized Least Squarse in R as the functions glim and gls, respectively, from S-Plus. Also, anybody know if there is any package that fit Log-Linear Models using Generalized Least Squares? Any help will be very useful. Thanks, -- Frederico

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

Hello together i want to execute different commands that are saved in strings (eg: "d <- 4+7" and some more complicated things). How can i execute/evaluate them. I know i could save them to a temp-file a read this with source. But is there a way without a file in the middle? I know also the 'assign' command but I think it would only the simply cases. thanks and gruess joerg

Hello all: I'm sure this is quite simple, but I tried several options and I still can't get I right. I've got a data.frame and I want to delete all columns of which the sum is zero. Any ideas? Thanks, Frederico Mestre [[alternative HTML version deleted]]

Hi all This seems like a very easy question, but I have not found a simple answer to it: If you have a dataframe, how to you plot two of then in different colors depending on a third variable? All help appreciated! /Fredrik -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hello everybody, a question that connect to the question of Frederik Karlsons about 'how to stand. betas' With the stand. betas i can compare the influence of the different explaning variables. What do i with the betas of factors? I can't use the solution of JohnFox, because there is no sd of an factor. How can i compare the influence of the factor with the influence of the numeric

Hi for all, i 've installed asterisk with isdn trunk with alcatel pabx. When alcatel users dial for external numbers, a channel on asterisk is allocated for dial. When we call to an number that is an IVR the digits isn't recognized by IVR. In sip.conf i putted dtmfmode as rfc... and info, inband is only for 64k codecs, and still don't work. How can i resolve this issue ?? Thanks.

Hello: While running R doing the analysis of my data I (using packages such as BIOMOD or e1071) get the following error as a result of several of my analysis: Error: cannot allocate vector of size 998.5 Mb In addition: Warning messages: 1: In array(c(rep.int(c(1, numeric(n)), n - 1L), 1), d, dn) : Reached total allocation of 4095Mb: see help(memory.size) 2: In array(c(rep.int(c(1,

Hi guys, I have an Alcatel 4200 with a ISDN board pluged in the asterisk server with TE110P. Input calls VOIP Proider ---> Asterisk ---> Alcatel Output Calls VOIP Proider <--- Asterisk <--- Alcatel In alcatel phones, users should dial 2 for take a line tone and can dial. At this point start my problems: 1. When users dial 2 on phone (alcatel) they don't received a dial tone,

Hello: What I want to do is quite simple, but I can't find a way. I have a data frame with several points (x and y coords). I want to add another column with cluster membership. For example aggregate all the points that stand within a distance of 40 from each other. I've tried using "nncluster" from the package nnclust, but the results are not correct, for some

Hello all: Is there any way to add points to a point pattern, while keeping a given minimum distance amongst this new points and the pre-existing points? Thanks, Frederico [[alternative HTML version deleted]]

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I have a stock market trading values time series. What's the best way to extract the "last day of month" values. I looked at function window() but doesn't appear suitable for this since it expects regular dates. Thank you. lukas - -- Lukas Kubin lukas.kubin at permonik.com phone: 00420603836180 -----BEGIN PGP SIGNATURE-----

Hello, I urgently need a histogram of the vector geo1$ds1 but: > hist(geo1$ds1) Error in hist(geo1$ds1) : unused argument(s) ( ...) I can use density, summary, boxplot and stem with this vector but not hist ? Any suggestions ? Thanks in advance, Detlef Groth > density(geo1$ds1) Call: density(x = geo1$ds1) Data: geo1$ds1 (1412 obs.); Bandwidth 'bw' = 0.2805

I have a set of data frames ds1, ds2, ... each having the same columns and column names: ds1 <- data.frame(x=c(1,2,3,4), y=c(5,6,7,8)) ds1 <- data.frame(x=c(9,10,11,12), y=c(13,14,15,16)) ... and I would like to combine them into just one data frame like ds <- rbind(ds1, ds2, ...) Because there are so many of them, I will have to use a character array nms <-

Hi, I've been trying to get the logon script to with in XP. This same logon script works in Win98. I finally discovered that it will work if I change logon drive = to logon drive = z: Now I'm really confused. The Netlogon is path is /home/samba and shared as \\server\netlogon. Z: is mapped as \\server\home\. This means the path to the logon script is Z:\samba\logon.bat Not

In this simple function, how can I pass strings for index and column names to the function? I've posted this type of question before and received no response. Maybe this example will be easier to understand and troubleshoot. ds <- function(myds, vec) {myds[[vec]]*2} ds1 <- c(X=list(1:10), Y=list(11:20)) ds(get("ds1"),get("Y")) khobson at odot.org Kenneth Ray

Using R version 2.4.1 (2006-12-18) on Windows, I have a dataset which resembles this: id att1 att2 att3 1 1 1 0 2 1 0 0 3 0 1 1 4 1 1 1 ratings <- data.frame(id = c(1,2,3,4), att1 = c(1,1,0,1), att2 = c(1,0,0,1), att3 = c(0,1,1,1)) I would like to get a cross tab of counts of co-ocurrence, which might resemble this:

Hello R users, When trying to fit a nonlinear mixed model to a respiration time series, I get the following error message: Error in recalc.varFunc(object[[i]], conLin) : dims [product 30] do not match the length of object [34] In addition: Warning message: In conLin$Xy * varWeights(object) : longer object length is not a multiple of shorter object length Below is an example that generates