similar to: Algorithm did not converge

Hello World, I have a function that makes pairwise comparisons between two strings. I would like to apply this function to my data (which consists of columns with different strings) in the way that it compares the first with the second entry, and then the third with the fourth, and then the fifth with the sixth, and so on down each column... So (2x-1) and (2x) would be the different entries to be

Hi! Does anybody know a string function that would calculate how many characters two strings share? I.e. ("Hello World","Hello Peter") would be 7. Thanks. Laetitia

Dear R users, I am trying to fit a glm with a distribution free family, link = log and variance = constant*mu. I guess I have to use the quasi family but the choices of variance are restricted to constant or mu or mu^2..., I don't know the way to choose the variance that I need, i.e. constant*mu. If you have any ideas or advice, please tell me. Thanks, Laetitia Mestdagh Laetitia Mestdagh

Hello everyone, How can I specify in tcltk file browser the initial directory to "MyComputer" in Windows where Drives and Partition are accessible? And just a little question if anyone knows, is there a way to use the function choose.files under windows to select a directory? Thanks a lot for your help. Laetitia.

Hi, I've had a little problem for several weeks now. It is annoying and therefore I will ask for help now: When I write a script with several iterations, I make it write out a text file to save the data during the run. For example I write: if (i %% 25) write.table(output,"temporary_output.txt") Later on, when I read in this output and want to calculate things, R complains that

Hi R people! I am looking for some suggestions writing an if-else function. The idea is to characterize different plots containing counts of variables (here parasites). If a plot has a count equal or higher than 4 for any parasite the function should return a 1 else a 0. Later I can loop the function over all plots. Here I have a little subset of my data: VariablePAR Plot1

Hi, I have R installed under a Mandrake linux system and I don't have shell utilities any more under my R console such as completion when writing a file path, back and forth in the history, bindkeys... Moreover when I quit R by saving, no .Rhistory file is created while the .Rdata is. I don't get how this work, I thought it was based on the user unix shell but it does not seem. Any

Hello, Is there a simple and fast function that returns a vector of the number of replications for each object of a vector ? For example : I have a vector of IDs : ids <- c( "ID1", "ID2", "ID2", "ID3", "ID3","ID3", "ID5") I want the function returns the following vector where each term is the number of replicates for the

Hi, Can anybody show me how to extract all columns in my dataset that are polymorphic? Or phrased in another way I would like to delete all columns that have no more than one letter in it (that are monomorphic). Thank you. Laetitia -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: seqCol.txt URL:

In the previous versions of R (2.6.1), when a vector of NA was given to the functions 'sd' or 'var' with parameter na.rm = TRUE, it used to return NA. Now (2.7.1) it returns an ERROR : Example in 2.6.1: > sd(c(NA, NA, NA, NA), na.rm = TRUE) [1] NA Example in 2.7.1: > sd(c(NA, NA, NA, NA), na.rm = TRUE) Error in var(x, na.rm = na.rm) : paires d'éléments

I am unfamiliar with R and I’m trying to do few statistical things like GLM and GAM with it. I hope my following questions will be clear enough: My datas ( y(i,j ))are run off triangles for example : J=1 J=2 J=3 I=1 1 2 3 I=2 4 5 I=3 6 My model is : E[y(i,j)] =m(i,j) Var[y(i,j)] =constant *m(i,j) Log(m(i,j)) = eta (i,j) eta (i,j) = c + alpha(i)

Dear Wu Gong and Peter Ehlers, thank you very much for your help debugging my script. Now I have a general following up question: Is there a straightforward way to rearrange the following dataset so that all first letters of each column will be combined in one column, all the second letters in a second column, all the third ones in a third column and so on, resulting in 7 columns, i.e. for

Hi! I need some help to finish my script. I have two tables that I combine randomly to produce a third table. This I do for hundreds of iterations. In the output file I get all the simulated tables after each other. It looks like this (in this case 3 iterations): output file: [[1]] [,1] [,2] [,3] [1,] "GM030005" "WI920024" "CCCC" [2,]

Hi All, I run function "kmeans" to cluster a matrix. But when the matrix size is big enough, it keeps saying "did not converge in 10 iterations". Could you explain what it means and if the result is wrong? And the interesting thing is sometimes this warning happens when the sample size is around 51200 x 6, sometimes it happens around 30000 x 6. Does the warning related

Hi everyone, I want to write a function that modify directly variables passed as parameters (the equivalent in C language of *ptr/&ptr) so that I don't have to return a list and to reaffect all my variables. Is it possible to do so in R? Thanks a lot. Laetitia Marisa.

Hello, When I launch a script under windows it does not display sequentially my cat calls or maybe the console is not refreshed at every line of my script. For instance with that code cat("\n\n================== IMPORT DATA FROM FILE ===================\n\n") fileschosen <- choose.files(caption="Select gpr files", filters = matrix(c("genepix

I am using the Ben bolker's R package "bbmle" to estimate the parameters of a binomial mixture distribution via Maximum Likelihood Method. For some data sets, I got the following warning messages: *Warning: optimization did not converge (code 1: ) There were 50 or more warnings (use warnings() to see the first 50)* Also, warnings() results the following: *In 0:(n - x) : numerical

After running a model for a while and seeing that it did not converge yet, how can I continue to run, ie not starting anew, the model? I know if I manually/interactively use winbugs, this is possible anytime, but how can I do this in r2winbugs, so that my existing sim$sims.array and other stuff in the object that bugs() returns gets extended? Thanks Toby

Hello R users, I'm trying to do simulations for comparing cox and weibull I have come across this problem: Warning messages: 1: In survreg.fit(X, Y, weights, offset, init = init, controlvals = control, : Ran out of iterations and did not converge 2: In survreg.fit(X, Y, weights, offset, init = init, controlvals = control, : Ran out of iterations and did not converge what i did is

Someone can help me? I tried several things and always don't converge # Model library(sem) dados40.cov <- cov(dados40,method="spearman") model.dados40 <- specify.model() F1 -> Item11, lam11, NA F1 -> Item31, lam31, NA F1 -> Item36, lam36, NA F1 -> Item54, lam54, NA F1 -> Item63, lam63, NA F1 -> Item65, lam55, NA F1 -> Item67, lam67, NA F1 ->