similar to: Fitting data from a spectrophotometer.

Dear List, Please bear with a poor newbee, who might be doing everything backwards (I was brought up in pure math). I want to make a simple multi-linear regression on a set of data. I did some expreiments, and if X is a 4 by 2 array and Y is a 4 by 1 array, I can do a linear regression by lm(y~x). Now I have a tab-delimited text file with 10 rows of 300 measurements and an other file with

Dear list, I would like to call a Kerberos-authenticated web-service from within R. Curl can do it: $ curl --negotiate -u : "http://my.web.service/" so I would expect that RCurl also has the capability, but I have not been able to find the correct options to set. listCurlOptions() does not return anything with negotiate, and searching the source of RCurl, the only thing I found was

Hi I need to use rep() to get a vector out, but I have spotted something very strange. See the reproducible example below. N <- 79 seg <- 5 segN <- N / seg # = 15.8 d1 <- seg - ( segN - floor(segN) ) * seg d1 # = 1 rep(2, d1) # = numeric(0), strange - why doesn't it print one "2"? rep(2, 1) # 2, ok rep(2, d1 / 1,1) # 2, this

The Register is carrying a article written by Kevin Poulsen of Securtiy Focus, calling asterisk "..the most powerful tool for manipulating and accessing CPN data.." > http://www.theregister.co.uk/2004/07/07/hackers_gut_voip/ I hope NuFone doesn't drop asterisk-set-able callerid's after this article; i've been wanting that feature from voicepluse for a long time.

Hi All- I am running into a small problem trying to implement voicemail Directory(). I'm sure it is a simple thing, but I can't figure out where the problem lies. I can get into the directory without a problem and can look up users by their last names, however I hit a snag when asterisk says "if this is the person you are looking for press 1 now". When I hit 1, the attendant

Hi all. I am new to the list and new to asterisk. I have asterisk installed and running. I am using it as a voicemail server only. What I would like to do is send users to a general mailbox that will be addressed as <companyname>@asterisk and give them the option to wait for the tone and leave a message, or press 9 to dial by name. My questions are: 1. What is the best way to do

I ran into a problem today which surprised me, and after two hours of reading and testing I am more confused than ever, so I feel the need for a sanity check. It would be good to have a clear and authoritative article about this somewhere; I thought it would be in either "Using Samba" or in the HOWTO, but I don't see it there. I am using Samba in several installations in homes and

I want a list of the number of times some factor levels appear together, similar to the following SQL statement: SELECT A, B, COUNT(C) FROM TBL GROUP BY A, B How do I do that with a data.frame in R? Thanks, Jesper Frickmann Statistician, Quality Control Novozymes North America Inc. 77 Perry Chapel Church Road Franklinton, NC 27525 USA Tel. +1 919 494 3266 [[alternative HTML version deleted]]

Hello R-Users, I am trying to label points on an xYplot. I can label the points, but then I lose the error bars. I suspect that I should be using panel.xYplot rather than panel.xyplot, but that didn't work either... Can anybody help me out? Example code is below. First xYplot has error bars, second tries to label the points and the bars disappear. Thanks, John a=c(1:4)

Hello, I am trying to run a logistic regression with random effects on proportional data in glmmBUGS. I am a newcomer to this package, and wondered if anyone could help me specify the model correctly. I am trying to specify the response variable, /yseed/, as # of successes out of total observations... but I suspect that given the error below, that is not correct. Also, Newsect should be a

Hello R-users, I have code that gives me the important variables from an analysis. I need to input these variables into a different analysis. To do this, I need to modify them slightly... 1) remove all numbers at the end of the variables, 2) remove all periods. I tried to do it with the awkward code below. It works to remove all the numbers, but when I try to remove the period everything

Yes it is. Example from my program: workdir <- "S:\\Stat\\C03202\\R-output\\" ... sink(file=paste(workdir, assayname, "-1.txt", sep="")) ... win.metafile(filename = paste(workdir, assayname, "-2%02d.wmf", sep="")) You get the idea! Kind regards, Jesper Frickmann Statistician, Quality Control Novozymes North America Inc. Tel. +1 919 494

Hello, I want to build a function to call up a column of a data.frame by the names of the columns. I have column names that are sequentially named (col1, col2, etc.). How do I change a character expression into something that will be understood as a data.frame column. For example: example<-data.frame(cbind(col1=1:10, col2=21:30, col3=41:50)) call.fun<-function(t){

Hello, I am trying to plot a curve over points plotted with se's in xYplot (see example below). I can get Figure 1 below to plot the data with error. However, I keep getting a the error message "Error using packet 1 object "y" not found" Can anyone see what I am doing wrong? Thanks! John ------------------------------- a=0.002; b=31.7; c=0.51

O'Reilly ang "Teach yourself samba in 24 hours" I have tried to get swat up running and configure samba therefrom. I have added /etc/services swat tcp 901 and because I cannot find = inetd.conf I have instead added in xinetd.conf (apparently red hat 7.1 don't have inetd.conf anymore = but xinetd.conf) =20 I have made the appriopriete changes there too /etc/xinetd.conf added =

Hello, I am trying to make a xyplot plot with points that are different symbols. I want to call the symbol type (pch) from a column in my dataframe. Here is a simplified example. In my real example I also have groups, which I have not included here. This example doesn't change the symbols or colors. Any help you can provide would be appreciated. Thanks, John x<-c(1:12)

Firstly, I am running beta2 as AD member. Auth through Winbind. Upgrading from beta1 to beta2 using RPM on RH 9 I noticed that Winbind will not start and smb will not stop using Init scripts. Seems to me RH9 /etc/init.d/winbind need the "winbind uid" parameter ? Here is the /var/log/samba/log.winbind output. [2003/07/04 13:54:52, 1] nsswitch/winbindd.c:main(811) winbindd version

Hello, I am using glmer() from lmer(lme4) to run generalized linear mixed models. However, I am having a problem extracting the standard errors for the fixed effects. I have used: summary(model)$coef fixed.effects(model) coef(model) to get out the parameter estimates, but do not seem able to extract the se's. Anybody have a solution? Thanks, John

Hello, I have been using mgcv to run GAM hurdle models, analyzing presence/absence data with GAM logistic regressions, and then analyzing the data conditional on presence (e.g. without samples with no zeros) with GAMs with a negative binomial distribution. It occurs to me that using the negative binomial distribution on data with no zeros is not right, as the negative binomial allows zeros.

Dear list, using the packages Hmisc and lattice i produced some nice xYplots. However, since the data range of the conditioning variable is very big, i need to define more than one y-scale for the plot (in some panels you just see a flat line of data points very close to the x-axis), e.g. different y-axis scales for the different rows of the plot. Is there a way to do so? Thanks, Henning